Kinetic advantage of intrinsically disordered proteins in coupled folding-binding process: a critical assessment of the "fly-casting" mechanism

The transition state for coupled folding and binding of a disordered DNA binding domain resembles the unbound state

The basic zippers (bZIPs) are one of two large eukaryotic families of transcription factors whose DNA binding domains are disordered in isolation but fold into stable α-helices upon target DNA binding. Here, we systematically disrupt pre-existing helical propensity within the DNA binding region of the homodimeric bZIP domain of cAMP-response element binding protein (CREB) using Ala-Gly scanning and examine the impact on target binding kinetics. We find that the secondary structure of the transition state strongly resembles that of the unbound state. The residue closest to the dimerization domain is largely folded within both unbound and transition states; dimerization apparently propagates additional helical propensity into the basic region. The results are consistent with electrostatically-enhanced DNA binding, followed by rapid folding from the folded zipper outwards. Fly-casting theory suggests that protein disorder can accelerate binding. Interestingly however, we did not observe higher association rate constants for mutants with lower levels of residual structure in the unbound state.

Correlating Disordered Activation Domain Ensembles with Gene Expression Levels

Transcription factor proteins bind to specific DNA promoter sequences and initiate gene transcription. In eukaryotes, most transcription factors contain intrinsically disordered activation domains (ADs) that regulate their transcriptional activity. Like other disordered protein regions, ADs do not have a fixed three-dimensional structure and instead exist in an ensemble of conformations. Disordered ensembles contain sequence-encoded structural preferences which are often linked to their function. We hypothesize this link exists between the structural preferences of disordered AD ensembles and their ability to induce gene expression. To test this, we used FRET microscopy to measure the ensemble dimensions of two activation domains, HIF-1α and CITED2, in live cells, and correlate this structural information with transcriptional activity. We find that point mutations that expanded the HIF-1α ensemble increased transcriptional activity, while those that compacted it reduced activity. Conversely, CITED2 showed no correlation between ensemble dimensions and activity. Our results reveal a sequence-dependent relationship between AD ensemble dimensions and their transcriptional activity.

Impact of Nonnative Interactions on the Binding Kinetics of Intrinsically Disordered p53 with MDM2: Insights from All-Atom Simulation and Markov State Model Analysis

Intrinsically disordered proteins (IDPs) lack a well-defined tertiary structure but are essential players in various biological processes. Their ability to undergo a disorder-to-order transition upon binding to their partners, known as the folding-upon-binding process, is crucial for their function. One classical example is the intrinsically disordered transactivation domain (TAD) of the tumor suppressor protein p53, which quickly forms a structured α-helix after binding to its partner MDM2, with clinical significance for cancer treatment. However, the contribution of nonnative interactions between the IDP and its partner to the rapid binding kinetics, as well as their interplay with native interactions, is not well understood at the atomic level. Here, we used molecular dynamics simulation and Markov state model (MSM) analysis to study the folding-upon-binding mechanism between p53-TAD and MDM2. Our results suggest that the system progresses from the nascent encounter complex to the well-structured encounter complex and finally reaches the native complex, following an induced-fit mechanism. We found that nonnative hydrophobic and hydrogen bond interactions, combined with native interactions, effectively stabilize the nascent and well-structured encounter complexes. Among the nonnative interactions, Leu25p53-Leu54MDM2 and Leu25p53-Phe55MDM2 are particularly noteworthy, as their interaction strength is close to the optimum. Evidently, strengthening or weakening these interactions could both adversely affect the binding kinetics. Overall, our findings suggest that nonnative interactions are evolutionarily optimized to accelerate the binding kinetics of IDPs in conjunction with native interactions.

Recent advances in de novo computational design and redesign of intrinsically disordered proteins and intrinsically disordered protein regions

In the early 2000s, the concept of "unstructured biology" has emerged to be an important field in protein science by generating various new research directions. Many novel strategies and methods have been developed that are focused on effectively identifying/predicting intrinsically disordered proteins (IDPs) and intrinsically disordered protein regions (IDPRs), identifying their potential functions, disorder based drug design etc. Due to the range of functions of IDPs/IDPRs and their involvement in various debilitating diseases they are of contemporary interest to the scientific community. Recent researches are focused on designing/redesigning specific IDPs/IDPRs de novo. These de novo design/redesigns of IDPs/IDPRs are carried out by altering compositional biases and specific sequence patterning parameters. The main focus of these researches is to influence specific molecular functions, phase behavior, cellular phenotypes etc. In this review, we first provide the differences of natively folded and natively unfolded or IDPs with respect to their potential energy landscapes. Here, we provide current understandings on the different computational design strategies and methods that have been utilized in de novo design and redesigns of IDPs and IDPRs. Finally, we conclude the review by discussing the challenges that have been faced during the computational design/design attempts of IDPs/IDPRs.

Unraveling the Behavior of Intrinsically Disordered Protein c-Myc: A Study Utilizing Gaussian-Accelerated Molecular Dynamics

The protein c-Myc is a transcription factor that remains largely intrinsically disordered and is known to be involved in various biological processes and is overexpressed in various cancers, making it an attractive drug target. However, intrinsically disordered proteins such as c-Myc do not show funnel-like basins in their free-energy landscapes; this makes their druggability a challenge. For the first time, we propose a heterodimer model of c-Myc/Max in full length in this work. We used Gaussian-accelerated molecular dynamics (GaMD) simulations to explore the behavior of c-Myc and its various regions, including the transactivation domain (TAD) and the basic helix-loop-helix-leucine-zipper (bHLH-Zipper) motif in three different conformational states: (a) monomeric c-Myc, (b) c-Myc when bound to its partner protein, Max, and (c) when Max was removed after binding. We analyzed the GaMD trajectories using root-mean-square deviation (RMSD), radius of gyration, root-mean-square fluctuation, and free-energy landscape (FEL) calculations to elaborate the behaviors of these regions. The results showed that the monomeric c-Myc structure showed a higher RMSD fluctuation as compared with the c-Myc/Max heterodimer in the bHLH-Zipper motif. This indicated that the bHLH-Zipper motif of c-Myc is more stable when it is bound to Max. The TAD region in both monomeric and Max-bound states showed similar plasticity in terms of RMSD. We also conducted residue decomposition calculations and showed that the c-Myc and Max interaction could be driven mainly by electrostatic interactions and the residues Arg299, Ile403, and Leu420 seemed to play important roles in the interaction. Our work provides insights into the behavior of c-Myc and its regions that could support the development of drugs that target c-Myc and other intrinsically disordered proteins.

Encounter Complexes Between the N-terminal of Neurotensin with the Extracellular Loop 2 of the Neurotensin Receptor 1 Steer Neurotensin to the Orthosteric Binding Pocket

Neurotensin (NT) is a linear disordered peptide that activates two different class A GPCRs, neurotensin receptor 1 (NTS1) and NTS2. Resolved structures of the complex of the C-terminal fragment of NT, NT8-13, with NTS1 shows the peptide takes a well-defined structure in the bound state. However, the mechanisms underlying NT recognition of NTS1, and the conformational transition of NT upon binding NTS1 is an open question that if answered may aid discovery of highly selective drugs and reveal potential secondary binding sites on the surface of the receptor. Herein we investigated the interactions guiding NT to the orthosteric binding pocket of NTS1 by combining NMR experiments with kinetic analysis of the binding pathway using stopped-flow fluorescence and mutagenesis on both NT and NTS1. We show the presence of transient structures in the middle part of NT that kinetically regulate the binding of NT to NTS1. Moreover, our results indicate that the binding pathway of NT onto NTS1 is mediated via electrostatic interactions between the N-terminal region of NT with the extracellular loop 2 of NTS1. These interactions induce backbone conformational changes in neurotensin similar to the bound-state neurotensin, suggesting that the N-terminal region of NT and these interactions should be considered for development of selective drugs against NTS1.

Intrinsically disordered region of talin's FERM domain functions as an initial PIP2 recognition site

Focal adhesions (FAs) mediate the interaction of the cytoskeleton with the extracellular matrix in a highly dynamic fashion. Talin is a central regulator, adaptor protein, and mechano-sensor of FA complexes. For recruitment and firm attachment at FAs, talin's N-terminal FERM domain binds to phosphatidylinositol 4,5-bisphosphate (PIP2)-enriched membranes. A newly published autoinhibitory structure of talin-1, where the known PIP2 interaction sites are covered up, lead us to hypothesize that a hitherto less examined loop insertion of the FERM domain acts as an additional and initial site of contact. We evaluated direct interactions of talin-1 with a PIP2 membrane by means of atomistic molecular dynamics simulations. We show that this unstructured, 33-residue-long loop strongly interacts with PIP2 and can facilitate further membrane contacts, including the canonical PIP2 interactions, by serving as a flexible membrane anchor. Under force as present at FAs, the extensible FERM loop ensures talin maintains membrane contacts when pulled away from the membrane by up to 7 nm. We identify key basic residues of the anchor mediating the highly dynamic talin-membrane interaction. Our results put forward an intrinsically disordered loop as a key and highly adaptable PIP2 recognition site of talin and potentially other PIP2-binding mechano-proteins.

The cofactor-dependent folding mechanism of Drosophila cryptochrome revealed by single-molecule pulling experiments

The link between cofactor binding and protein activity is well-established. However, how cofactor interactions modulate folding of large proteins remains unknown. We use optical tweezers, clustering and global fitting to dissect the folding mechanism of Drosophila cryptochrome (dCRY), a 542-residue protein that binds FAD, one of the most chemically and structurally complex cofactors in nature. We show that the first dCRY parts to fold are independent of FAD, but later steps are FAD-driven as the remaining polypeptide folds around the cofactor. FAD binds to largely unfolded intermediates, yet with association kinetics above the diffusion-limit. Interestingly, not all FAD moieties are required for folding: whereas the isoalloxazine ring linked to ribitol and one phosphate is sufficient to drive complete folding, the adenosine ring with phosphates only leads to partial folding. Lastly, we propose a dCRY folding model where regions that undergo conformational transitions during signal transduction are the last to fold.

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.

Intrinsic protein disorder uncouples affinity from binding specificity

Intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) of proteins often function by molecular recognition, in which they undergo induced folding. Based on prior generalizations, the idea prevails in the IDP field that due to the entropic penalty of induced folding, the major functional advantage associated with this binding mode is "uncoupling" specificity from binding strength. Nevertheless, both weaker binding and high specificity of IDPs/IDRs rest on limited experimental observations, making these assumptions more speculations than evidence-supported facts. The issue is also complicated by the rather vague concept of specificity that lacks an exact measure, such as the Kd for binding strength. We addressed these issues by creating and analyzing a comprehensive dataset of well-characterized ID/globular protein complexes, for which both the atomic structure of the complex and free energy (ΔG, Kd ) of interaction is known. Through this analysis, we provide evidence that the affinity distributions of IDP/globular and globular/globular complexes show different trends, whereas specificity does not connote to weaker binding strength of IDPs/IDRs. Furthermore, protein disorder extends the spectrum in the direction of very weak interactions, which may have important regulatory consequences and suggest that, in a biological sense, strict correlation of specificity and binding strength are uncoupled by structural disorder.

Conformational Heterogeneity and Frustration of the Tumor Suppressor p53 as Tuned by Punctual Mutations

The conformational heterogeneity of the p53 tumor suppressor, the wild-type (p53wt) and mutated forms, was investigated by a computational approach, including the modeling and all atoms of the molecular dynamics (MD) simulations. Four different punctual mutations (p53R175H, p53R248Q, p53R273H, and p53R282W) which are known to affect the DNA binding and belong to the most frequent hot-spot mutations in human cancers, were taken into consideration. The MD trajectories of the wild-type and mutated p53 forms were analyzed by essential dynamics to extract the relevant collective motions and by the frustration method to evaluate the degeneracy of the energy landscape. We found that p53 is characterized by wide collective motions and its energy landscape exhibits a rather high frustration level, especially in the regions involved in the binding to physiological ligands. Punctual mutations give rise to a modulation of both the collective motions and the frustration of p53, with different effects depending on the mutation. The regions of p53wt and of the mutated forms characterized by a high frustration level are also largely involved in the collective motions. Such a correlation is discussed also in connection with the intrinsic disordered character of p53 and with its central functional role.

Intrinsically disordered proteins in viral pathogenesis and infections

Disordered proteins serve a crucial part in many biological processes that go beyond the capabilities of ordered proteins. A large number of virus-encoded proteins have extremely condensed proteomes and genomes, which results in highly disordered proteins. The presence of these IDPs allows them to rapidly adapt to changes in their biological environment and play a significant role in viral replication and down-regulation of host defense mechanisms. Since viruses undergo rapid evolution and have a high rate of mutation and accumulation in their proteome, IDPs' insights into viruses are critical for understanding how viruses hijack cells and cause disease. There are many conformational changes that IDPs can adopt in order to interact with different protein partners and thus stabilize the particular fold and withstand high mutation rates. This chapter explains the molecular mechanism behind viral IDPs, as well as the significance of recent research in the field of IDPs, with the goal of gaining a deeper comprehension of the essential roles and functions played by viral proteins.

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.

Molecular Simulations of Intrinsically Disordered Proteins and Their Binding Mechanisms

Intrinsically disordered proteins (IDPs) lack well-defined secondary or tertiary structures in solution but are found to be involved in a wide range of critical cellular processes that highlight their functional importance. IDPs usually undergo folding upon binding to their targets. Such binding coupled to folding behavior has widened our perspective on the protein structure-dynamics-function paradigm in molecular biology. However, characterizing the folding upon binding mechanism of IDPs experimentally remains quite challenging. Molecular simulations emerge as a potentially powerful tool that offers information complementary to experiments. Here we present a general computational framework for the molecular simulations of IDP folding upon binding processes that combines all-atom molecular dynamics (MD) and coarse-grained simulations. The classical all-atom molecular dynamics approach using GPU acceleration allows the researcher to explore the properties of the IDP conformational ensemble, whereas coarse-grained structure-based models implemented with parameters carefully calibrated to available experimental measurements can be used to simulate the entire folding upon binding process. We also discuss a set of tools for the analysis of MD trajectories and describe the details of the computational protocol to follow so that it can be adapted by the user to study any IDP in isolation and in complex with partners.

Physics of biomolecular recognition and conformational dynamics

Biomolecular recognition usually leads to the formation of binding complexes, often accompanied by large-scale conformational changes. This process is fundamental to biological functions at the molecular and cellular levels. Uncovering the physical mechanisms of biomolecular recognition and quantifying the key biomolecular interactions are vital to understand these functions. The recently developed energy landscape theory has been successful in quantifying recognition processes and revealing the underlying mechanisms. Recent studies have shown that in addition to affinity, specificity is also crucial for biomolecular recognition. The proposed physical concept of intrinsic specificity based on the underlying energy landscape theory provides a practical way to quantify the specificity. Optimization of affinity and specificity can be adopted as a principle to guide the evolution and design of molecular recognition. This approach can also be used in practice for drug discovery using multidimensional screening to identify lead compounds. The energy landscape topography of molecular recognition is important for revealing the underlying flexible binding or binding-folding mechanisms. In this review, we first introduce the energy landscape theory for molecular recognition and then address four critical issues related to biomolecular recognition and conformational dynamics: (1) specificity quantification of molecular recognition; (2) evolution and design in molecular recognition; (3) flexible molecular recognition; (4) chromosome structural dynamics. The results described here and the discussions of the insights gained from the energy landscape topography can provide valuable guidance for further computational and experimental investigations of biomolecular recognition and conformational dynamics.

Effects of electrostatic interactions on global folding and local conformational dynamics of a multidomain Y-family DNA polymerase

The Y-family DNA polymerases specialize in translesion DNA synthesis, which is essential for replicating damaged DNA. The Y-family polymerases, which are made up of four stable domains, exhibit extensive distributions of charged residues, and are responsible for the tight formation of the protein-DNA complex. However, it is still unclear how the electrostatic interactions influence the conformational dynamics of the polymerases. Here, we focus on the case of a prototype Y-family DNA polymerase, Dpo4. Using coarse-grained models including a salt-dependent electrostatic potential, we investigate the effects of the electrostatic interactions on the folding process of Dpo4. Our simulations show that strong electrostatic interactions result in a three-state folding of Dpo4, consistent with the experimental observations. This folding process exhibits low cooperativity led by low salt concentration, where the individual domains fold one by one through one single pathway. Since the refined folding order of domains in multidomain proteins can shrink the configurational space, we suggest that the electrostatic interactions facilitate the Dpo4 folding. In addition, we study the local conformational dynamics of Dpo4 in terms of fluctuation and frustration analyses. We show that the electrostatic interactions can exaggerate the local conformational properties, which are in favor of the large-scale conformational transition of Dpo4 during the functional DNA binding. Our results underline the importance of electrostatic interactions in the conformational dynamics of Dpo4 at both the global and local scale, providing useful guidance in protein engineering at the multidomain level.

The disordered regions of the methyltransferase SETD2 govern its function by regulating its proteolysis and phase separation

SETD2 is an important methyltransferase that methylates crucial substrates such as histone H3, tubulin, and STAT1 and also physically interacts with transcription and splicing regulators such as Pol II and various hnRNPs. Of note, SETD2 has a functionally uncharacterized extended N-terminal region, the removal of which leads to its stabilization. How this region regulates SETD2 half-life is unclear. Here we show that SETD2 consists of multiple long disordered regions across its length that cumulatively destabilize the protein by facilitating its proteasomal degradation. SETD2 disordered regions can reduce the half-life of the yeast homolog Set2 in mammalian cells as well as in yeast, demonstrating the importance of intrinsic structural features in regulating protein half-life. In addition to the shortened half-life, by performing fluorescence recovery after photobleaching assay we found that SETD2 forms liquid droplets in vivo, another property associated with proteins that contain disordered regions. The phase-separation behavior of SETD2 is exacerbated upon the removal of its N-terminal segment and results in activator-independent histone H3K36 methylation. Our findings reveal that disordered region-facilitated proteolysis is an important mechanism governing SETD2 function.

Intrinsically disordered proteins and membranes: a marriage of convenience for cell signalling?

The structure-function paradigm has guided investigations into the molecules involved in cellular signalling for decades. The peripheries of this paradigm, however, start to unravel when considering the co-operation between proteins and the membrane in signalling processes. Intrinsically disordered regions hold distinct advantages over folded domains in terms of their binding promiscuity, sensitivity to their particular environment and their ease of modulation through post-translational modifications. Low sequence complexity and bias towards charged residues are also favourable for the multivalent electrostatic interactions that occur at the surfaces of lipid bilayers. This review looks at the principles behind the successful marriage between protein disorder and membranes in addition to the role of this partnership in modifying and regulating signalling in cellular processes. The HVR (hypervariable region) of small GTPases is highlighted as a well-studied example of the nuanced role a short intrinsically disordered region can play in the fine-tuning of signalling pathways.

Cancer-Associated Mutations Perturb the Disordered Ensemble and Interactions of the Intrinsically Disordered p53 Transactivation Domain

Intrinsically disordered proteins (IDPs) are key components of regulatory networks that control crucial aspects of cell decision making. The intrinsically disordered transactivation domain (TAD) of tumor suppressor p53 mediates its interactions with multiple regulatory pathways to control the p53 homeostasis during the cellular response to genotoxic stress. Many cancer-associated mutations have been discovered in p53-TAD, but their structural and functional consequences are poorly understood. Here, by combining atomistic simulations, NMR spectroscopy, and binding assays, we demonstrate that cancer-associated mutations can significantly perturb the balance of p53 interactions with key activation and degradation regulators. Importantly, the four mutations studied in this work do not all directly disrupt the known interaction interfaces. Instead, at least three of these mutations likely modulate the disordered state of p53-TAD to perturb its interactions with regulators. Specifically, NMR and simulation analysis together suggest that these mutations can modulate the level of conformational expansion as well as rigidity of the disordered state. Our work suggests that the disordered conformational ensemble of p53-TAD can serve as a central conduit in regulating the response to various cellular stimuli at the protein-protein interaction level. Understanding how the disordered state of IDPs may be modulated by regulatory signals and/or disease associated perturbations will be essential in the studies on the role of IDPs in biology and diseases.

Introducing intrinsic disorder reduces electrostatic steering in protein-protein interactions

Protein-protein interactions underlie many critical biology functions, such as cellular signaling and gene expression, in which electrostatic interactions can play a critical role in mediating the specificity and stability of protein complexes. A substantial portion of proteins are intrinsically disordered, and the influences of structural disorder on binding kinetics and thermodynamics have been widely investigated. However, whether the effect of electrostatic steering depends on structural disorder remains unexplored. In this work, we addressed the consequence of introducing intrinsic disorder in the electrostatic steering of the E3/Im3 complex using molecular dynamics simulation. Our results recapitulated the experimental observations that the responses of stability and kinetics to salt concentration for the ordered E3/Im3 complex were larger than those for the disordered E3/Im3 complex. Mechanistic analysis revealed that the native contact interactions involved in the encounter state and the transition state were essentially identical for both ordered and disordered E3. Therefore, the observed difference in electrostatic steering between ordered E3 and disordered E3 may result from their difference in conformation rather than their difference in binding mechanism. Because charged residues are frequently involved in protein-protein interactions, our results suggest that increasing structural disorder is expected to generally modulate the effect of electrostatic steering.

Solution structure of the anticancer p28 peptide in biomimetic medium

The p28 peptide derived from Pseudomonas aeruginosa azurin shows an anticancer activity after binding to p53 protein and is currently in Phase I of clinical trials. We have studied its structure in water and in a biomimetic media and show that the peptide is unstructured in water but when studied in a biomimetic medium assumes a structure very similar to the one observed in azurin, suggesting a high propensity of this peptide to maintain this secondary structure. Analysis of p28 sequences from different bacterial species indicates conservation of the secondary structure despite amino acid replacement in different positions, suggesting that others, similar peptides could be tested for binding to p53.

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

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

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.

Structural and dynamic characterization of the C-terminal tail of ErbB2: Disordered but not random

ErbB2 (or HER2) is a receptor tyrosine kinase overexpressed in some breast cancers and associated with poor prognosis. Treatments targeting the receptor extracellular and kinase domains have greatly improved disease outcome in the last 20 years. In parallel, the structures of these domains have been described, enabling better mechanistic understanding of the receptor function and targeted inhibition. However, the ErbB2 disordered C-terminal cytoplasmic tail (CtErbB2) remains very poorly characterized in terms of structure, dynamics, and detailed functional mechanism. Yet, it is where signal transduction is triggered via phosphorylation of tyrosine residues and carried out via interaction with adaptor proteins. Here, we report the first description, to our knowledge, of the ErbB2 disordered tail at atomic resolution using NMR, complemented by small-angle x-ray scattering. We show that although no part of CtErbB2 has any fully populated secondary or tertiary structure, it contains several transient α-helices and numerous transient polyproline II helices, populated up to 20 and 40%, respectively, and low but significant compaction. The presence of some structural elements suggests, along the lines of the results obtained for EGFR (ErbB1), that they may have a functional role in ErbB2's autoregulation processes. In addition, the transient formation of polyproline II helices is compliant with previously suggested interactions with SH3 domains. All in all, our in-depth structural study opens perspectives in the mechanistic understanding of ErbB2.

Influence of sequence length and charged residues on Swc5 binding with histone H2A-H2B

SWR is a member of chromatin remodeler family and participates the replacement of histone H2A with H2A.Z. One of the SWR subunits, Swc5, has an intrinsically disordered region and binds to H2A-H2B dimer. Though the binding structure of Swc5 and H2A-H2B has been resolved recently, it is still challenging to investigate the binding mechanism as well as the role of the charge interactions between Swc5 and H2A-H2B. Here we developed a coarse-grained structure-based model and performed molecular dynamics simulations to investigate the binding processes of two Swc5 regions with different lengths (swc5-a and swc5-b) to H2A-H2B. The simulation results suggest a different role of electrostatic interactions between swc5-a/swc5-b and H2A-H2B on binding. The electrostatic interactions between swc5-a/swc5-b and H2A-H2B can not only accelerate the initial capture step of binding, but can also trap the swc5-a/swc5-b at the wrong binding site on H2A. Besides, the conserved DEF/Y-2 motif of Swc5 is important for the binding affinity and the recognition with H2A-H2B at the initial step. Both swc5-a and swc5-b undergo a structural shift before reaching the final bound state. This theoretical study provides important details and the underlying physical mechanisms of the binding processes of swc5-a/swc5-b and H2A-H2B.

Phosphorylation-Dependent Conformations of the Disordered Carboxyl-Terminus Domain in the Epidermal Growth Factor Receptor

The epidermal growth factor receptor (EGFR), a receptor tyrosine kinase, regulates basic cellular functions and is a major target for anticancer therapeutics. The carboxyl-terminus domain is a disordered region of EGFR that contains the tyrosine residues, which undergo autophosphorylation followed by docking of signaling proteins. Local phosphorylation-dependent secondary structure has been identified and is thought to be associated with the signaling cascade. Deciphering and distinguishing the overall conformations, however, have been challenging because of the disordered nature of the carboxyl-terminus domain and resultant lack of well-defined three-dimensional structure for most of the domain. We investigated the overall conformational states of the isolated EGFR carboxyl-terminus domain using single-molecule Förster resonance energy transfer and coarse-grained simulations. Our results suggest that electrostatic interactions between charged residues emerge within the disordered domain upon phosphorylation, producing a looplike conformation. This conformation may enable binding of downstream signaling proteins and potentially reflect a general mechanism in which electrostatics transiently generate functional architectures in disordered regions of a well-folded protein.

Functions of short lifetime biological structures at large: the case of intrinsically disordered proteins

Although for more than a century a protein function was intimately associated with the presence of unique structure in a protein molecule, recent years witnessed a skyrocket rise of the appreciation of protein intrinsic disorder concept that emphasizes the importance of the biologically active proteins without ordered structures. In different proteins, the depth and breadth of disorder penetrance are different, generating an amusing spatiotemporal heterogeneity of intrinsically disordered proteins (IDPs) and intrinsically disordered protein region regions (IDPRs), which are typically described as highly dynamic ensembles of rapidly interconverting conformations (or a multitude of short lifetime structures). IDPs/IDPRs constitute a substantial part of protein kingdom and have unique functions complementary to functional repertoires of ordered proteins. They are recognized as interaction specialists and global controllers that play crucial roles in regulation of functions of their binding partners and in controlling large biological networks. IDPs/IDPRs are characterized by immense binding promiscuity and are able to use a broad spectrum of binding modes, often resulting in the formation of short lifetime complexes. In their turn, functions of IDPs and IDPRs are controlled by various means, such as numerous posttranslational modifications and alternative splicing. Some of the functions of IDPs/IDPRs are briefly considered in this review to shed some light on the biological roles of short-lived structures at large.

A Switch between Two Intrinsically Disordered Conformational Ensembles Modulates the Active Site of a Basic-Helix-Loop-Helix Transcription Factor

We report a conformational switch between two distinct intrinsically disordered subensembles within the active site of a transcription factor. This switch highlights an evolutionary benefit conferred by the high plasticity of intrinsically disordered domains, namely, their potential to dynamically sample a heterogeneous conformational space housing multiple states with tailored properties. We focus on proto-oncogenic basic-helix-loop-helix (bHLH)-type transcription factors, as these play key roles in cell regulation and function. Despite intense research efforts, the understanding of structure-function relations of these transcription factors remains incomplete as they feature intrinsically disordered DNA-interaction domains that are difficult to characterize, theoretically as well as experimentally. Here we characterize the structural dynamics of the intrinsically disordered region DNA-binding site of the vital MYC-associated transcription factor X (MAX). Integrating nuclear magnetic resonance (NMR) measurements, molecular dynamics (MD) simulations, and electron paramagnetic resonance (EPR) measurements, we show that, in the absence of DNA, the binding site of the free MAX₂ homodimer samples two intrinsically disordered conformational subensembles. These feature distinct structural properties: one subensemble consists of a set of highly flexible and spatially extended conformers, while the second features a set of "hinged" conformations. In this latter ensemble, the disordered N-terminal tails of MAX₂ fold back along the dimer, forming transient long-range contacts with the HLH-region and thereby exposing the DNA binding site to the solvent. The features of these divergent substates suggest two mechanisms by which protein conformational dynamics in MAX₂ might modulate DNA-complex formation: by enhanced initial recruitment of free DNA ligands, as a result of the wider conformational space sampled by the extended ensemble, and by direct exposure of the binding site and the corresponding strong electrostatic attractions presented while in the hinged conformations.

Time-Resolved Fluorescence and Essential Dynamics Study on the Structural Heterogeneity of p53DBD Bound to the Anticancer p28 Peptide

Time-resolved fluorescence emission was combined with molecular dynamics (MD) simulations to investigate the DNA-binding domain (DBD) of the tumor suppressor p53 alone and its complex with the anticancer peptide p28 (DBD/p28). The fluorescence emission decay of the lone Trp residue, from both DBD and DBD/p28, was well-described by a stretched exponential function. Such a behavior was ascribed to heterogeneity in the Trp relaxation behavior, likely due to the coexistence of different conformational states. The increase of the stretching parameter, on passing from DBD to DBD/p28, indicates a reduced heterogeneity in the Trp146 environment for DBD/p28. Moreover, the effects of p28 on the global dynamics of DBD were analyzed by the essential dynamics method on 30 ns long MD trajectories of both DBD and DBD/p28. We found the establishment of wide-amplitude anharmonic modes throughout the DBD molecule, with a particularly high amplitude being detected in the DNA-binding region. These modes are significantly reduced when DBD is bound to p28, consistently with a structure stabilization. In summary, the results indicate that p28 binding has a strong effect on both the local and global heterogeneity of DBD, thus providing some hints to the understanding of its anticancer activity.

Investigating the trade-off between folding and function in a multidomain Y-family DNA polymerase

The way in which multidomain proteins fold has been a puzzling question for decades. Until now, the mechanisms and functions of domain interactions involved in multidomain protein folding have been obscure. Here, we develop structure-based models to investigate the folding and DNA-binding processes of the multidomain Y-family DNA polymerase IV (DPO4). We uncover shifts in the folding mechanism among ordered domain-wise folding, backtracking folding, and cooperative folding, modulated by interdomain interactions. These lead to ‘U-shaped’ DPO4 folding kinetics. We characterize the effects of interdomain flexibility on the promotion of DPO4–DNA (un)binding, which probably contributes to the ability of DPO4 to bypass DNA lesions, which is a known biological role of Y-family polymerases. We suggest that the native topology of DPO4 leads to a trade-off between fast, stable folding and tight functional DNA binding. Our approach provides an effective way to quantitatively correlate the roles of protein interactions in conformational dynamics at the multidomain level.

Presence of intrinsically disordered proteins can inhibit the nucleation phase of amyloid fibril formation of Aβ(1-42) in amino acid sequence independent manner

The molecular shield effect was studied for intrinsically disordered proteins (IDPs) that do not adopt compact and stable protein folds. IDPs are found among many stress-responsive gene products and cryoprotective- and drought-protective proteins. We recently reported that some fragments of human genome-derived IDPs are cryoprotective for cellular enzymes, despite a lack of relevant amino acid sequence motifs. This sequence-independent IDP function may reflect their molecular shield effect. This study examined the inhibitory activity of IDPs against fibril formation in an amyloid beta peptide (Aβ(1–42)) model system. Four of five human genome-derived IDPs (size range 20 to 44 amino acids) showed concentration-dependent inhibition of amyloid formation (IC₅₀ range between 60 and 130 μM against 20 μM Aβ(1–42)). The IC₅₀ value was two orders of magnitude lower than that of polyethylene-glycol and dextran, used as neutral hydrophilic polymer controls. Nuclear magnetic resonance with ¹⁵ N-labeled Aβ(1–42) revealed no relevant molecular interactions between Aβ(1–42) and IDPs. The inhibitory activities were abolished by adding external amyloid-formation seeds. Therefore, IDPs seemed to act only at the amyloid nucleation phase but not at the elongation phase. These results suggest that IDPs (0.1 mM or less) have a molecular shield effect that prevents aggregation of susceptible molecules.

Charge Interactions Modulate the Encounter Complex Ensemble of Two Differently Charged Disordered Protein Partners of KIX

Disordered proteins play important roles in cell signaling and are frequently involved in protein-protein interactions. They also have a larger proportion of charged and polar residues than their folded counterparts. Here, we developed a structure-based model and applied molecular dynamics simulations to examine the presence and importance of electrostatic interactions in the binding processes of two differently charged intrinsically disordered ligands of the KIX domain of CBP. We observed non-native opposite-charged contacts in the encounter complexes for both ligands with KIX, and this may be a general feature of coupled folding and binding reactions. The ensemble of successful encounter complexes is a diverse set of structures, and in the case of the highly charged ligand, this ensemble was found to be malleable with respect to ionic strength. There are only minor differences between encounter complex ensembles for successful and unsuccessful collisions with no key interactions that appear to make the process far more productive. The energy landscape at this early stage in the process does not appear highly funneled. Strikingly we observed many native interactions that appear to reduce chances of an encounter complex being productive. Instead it appears that collectively non-native electrostatic interactions in the encounter complex increase the likelihood of productivity by holding the proteins together long enough for folding to take place. This mechanism is more effective for the more highly charged ligand.

Nonclassical nuclear localization signals mediate nuclear import of CIRBP

The specific interaction of importins with nuclear localization signals (NLSs) of cargo proteins not only mediates nuclear import but also, prevents their aberrant phase separation and stress granule recruitment in the cytoplasm. The importin Transportin-1 (TNPO1) plays a key role in the (patho-)physiology of both processes. Here, we report that both TNPO1 and Transportin-3 (TNPO3) recognize two nonclassical NLSs within the cold-inducible RNA-binding protein (CIRBP). Our biophysical investigations show that TNPO1 recognizes an arginine-glycine(-glycine) (RG/RGG)-rich region, whereas TNPO3 recognizes a region rich in arginine-serine-tyrosine (RSY) residues. These interactions regulate nuclear localization, phase separation, and stress granule recruitment of CIRBP in cells. The presence of both RG/RGG and RSY regions in numerous other RNA-binding proteins suggests that the interaction of TNPO1 and TNPO3 with these nonclassical NLSs may regulate the formation of membraneless organelles and subcellular localization of numerous proteins.

Two Differential Binding Mechanisms of FG-Nucleoporins and Nuclear Transport Receptors

Phenylalanine-glycine-rich nucleoporins (FG-Nups) are intrinsically disordered proteins, constituting the selective barrier of the nuclear pore complex (NPC). Previous studies showed that nuclear transport receptors (NTRs) were found to interact with FG-Nups by forming an “archetypal-fuzzy” complex through the rapid formation and breakage of interactions with many individual FG motifs. Here, we use single-molecule studies combined with atomistic simulations to show that, in sharp contrast, FG-Nup214 undergoes a coupled reconfiguration-binding mechanism when interacting with the export receptor CRM1. Association and dissociation rate constants are more than an order of magnitude lower than in the archetypal-fuzzy complex between FG-Nup153 and NTRs. Unexpectedly, this behavior appears not to be encoded selectively into CRM1 but rather into the FG-Nup214 sequence. The same distinct binding mechanisms are unperturbed in O-linked β-N-acetylglucosamine-modified FG-Nups. Our results have implications for differential roles of distinctly spatially distributed FG-Nup⋅NTR interactions in the cell.

•

Identification of two differential binding mechanisms in the nuclear transport pathway

•

FG-Nup214 does not bind CRM1 via an archetypal-fuzzy complex

•

Glycosylated FG-Nups maintain their NTR-binding mechanisms

•

Linker regions of FG-Nups may have functional relevance to the binding mechanism

Features of molecular recognition of intrinsically disordered proteins via coupled folding and binding

The sequence-structure-function paradigm of proteins has been revolutionized by the discovery of intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs). In contrast to traditional ordered proteins, IDPs/IDRs are unstructured under physiological conditions. The absence of well-defined three-dimensional structures in the free state of IDPs/IDRs is fundamental to their function. Folding upon binding is an important mode of molecular recognition for IDPs/IDRs. While great efforts have been devoted to investigating the complex structures and binding kinetics and affinities, our knowledge on the binding mechanisms of IDPs/IDRs remains very limited. Here, we review recent advances on the binding mechanisms of IDPs/IDRs. The structures and kinetic parameters of IDPs/IDRs can vary greatly, and the binding mechanisms can be highly dependent on the structural properties of IDPs/IDRs. IDPs/IDRs can employ various combinations of conformational selection and induced fit in a binding process, which can be templated by the target and/or encoded by the IDP/IDR. Further studies should provide deeper insights into the molecular recognition of IDPs/IDRs and enable the rational design of IDP/IDR binding mechanisms in the future.

Intrinsically Disordered Transactivation Domains Bind to TAZ1 Domain of CBP via Diverse Mechanisms

CREB-binding protein is a multidomain transcriptional coactivator whose transcriptional adaptor zinc-binding 1 (TAZ1) domain mediates interactions with a number of intrinsically disordered transactivation domains (TADs), including the CREB-binding protein/p300-interacting transactivator with ED-rich tail, the hypoxia inducible factor 1α, p53, the signal transducer and activator of transcription 2, and the NF-κB p65 subunit. These five disordered TADs undergo partial disorder-to-order transitions upon binding TAZ1, forming fuzzy complexes with helical segments. Interestingly, they wrap around TAZ1 with different orientations and occupy the binding sites with various orders. To elucidate the microscopic molecular details of the binding processes of TADs with TAZ1, in this work, we carried out extensive molecular dynamics simulations using a coarse-grained topology-based model. After careful calibration of the models to reproduce the residual helical contents and binding affinities, our simulations were able to recapitulate the experimentally observed flexibility profiles. Although great differences exist in the complex structures, we found similarities between hypoxia inducible factor 1α and signal transducer and activator of transcription 2 as well as between CREB-binding protein/p300-interacting transactivator with ED-rich tail and NF-κB p65 subunit in the binding kinetics and binding thermodynamics. Although the origins of similarities and differences in the binding mechanisms remain unclear, our results provide some clues that indicate that binding of TADs to TAZ1 could be templated by the target as well as encoded by the TADs.

Binding induced folding: Lessons from the kinetics of interaction between NTAIL and XD

Intrinsically Disordered Proteins (IDPs) are a class of protein that exert their function despite lacking a well-defined three-dimensional structure, which is sometimes achieved only upon binding to their natural ligands. This feature implies the folding of IDPs to be generally coupled with a binding event, representing an interesting challenge for kinetic studies. In this review, we recapitulate some of the most important findings of IDPs binding-induced folding mechanisms obtained by analyzing their binding kinetics. Furthermore, by focusing on the interaction between the Measles virus N_TAIL protein, a prototypical IDP, and its physiological partner, the X domain, we recapitulate the major theoretical and experimental approaches that were used to describe binding induced folding.

Disordered linkers in multidomain allosteric proteins: Entropic effect to favor the open state or enhanced local concentration to favor the closed state?

There are many multidomain allosteric proteins where an allosteric signal at the allosteric domain modifies the activity of the functional domain. Intrinsically disordered regions (linkers) are widely involved in this kind of regulation process, but the essential role they play therein is not well understood. Here, we investigated the effect of linkers in stabilizing the open or the closed states of multidomain proteins using combined thermodynamic deduction and coarse-grained molecular dynamics simulations. We revealed that the influence of linker can be fully characterized by an effective local concentration [B]₀ . When K_d is smaller than [B]₀ , the closed state would be favored; while the open state would be preferred when K_d is larger than [B]₀ . We used four protein systems with markedly different domain-domain binding affinity and structural order/disorder as model systems to understand the relationship between [B]₀ and the linker length as well as its flexibility. The linker length is the main practical determinant of [B]₀ . [B]₀ of a flexible linker with 40-60 residues was determined to be in a narrow range of 0.2-0.6 mM, while a too short or too long length would dramatically decrease [B]₀ . With the revealed [B]₀ range, the introduction of a flexible linker makes the regulation of weakly interacting partners possible.

Designed Mutations Alter the Binding Pathways of an Intrinsically Disordered Protein

Many cellular functions, including signaling and regulation, are carried out by intrinsically disordered proteins (IDPs) binding to their targets. Experimental and computational studies have now significantly advanced our understanding of these binding processes. In particular, IDPs that become structured upon binding typically follow a dock-and-coalesce mechanism, whereby the docking of one IDP segment initiates the process, followed by on-target coalescence of remaining IDP segments. Multiple dock-and-coalesce pathways may exist, but one may dominate, by relying on electrostatic attraction and molecular flexibility for fast docking and fast coalescing, respectively. Here we critically test this mechanistic understanding by designing mutations that alter the dominant pathway. This achievement marks an important step toward precisely manipulating IDP functions.

Specific and Fuzzy Interactions Cooperate in Modulating Protein Half-Life

Protein degradation is critical for maintaining cellular homeostasis. The 20S proteasome is selective for unfolded, extended polypeptide chains without ubiquitin tags. Sequestration of such segments by protein partners, however, may provide a regulatory mechanism. Here we used the AP-1 complex to study how c-Fos turnover is controlled by interactions with c-Jun. We show that heterodimerization with c-Jun increases c-Fos half-life. Mutations affecting specific contact sites (L165V, L172V) or charge separation (E175D, E189D, K190R) with c-Jun both modulate c-Fos turnover, proportionally to their impact on binding affinity. The fuzzy tail beyond the structured b-HLH/ZIP domain (~165 residues) also contributes to the stabilization of the AP-1 complex, removal of which decreases c-Fos half-life. Thus, protein turnover by 20S proteasome is fine-tuned by both specific and fuzzy interactions, consistently with the previously proposed "nanny" model.

Position-, disorder-, and salt-dependent diffusion in binding-coupled-folding of intrinsically disordered proteins

Successful extensions of protein-folding energy landscape theory to intrinsically disordered proteins' (IDPs') binding-coupled-folding transition can enormously simplify this biomolecular process into diffusion along a limited number of reaction coordinates, and the dynamics subsequently is described by Kramers' rate theory. As the critical pre-factor, the diffusion coefficient D has direct implications on the binding kinetics. Here, we employ a structure-based model (SBM) to calculate D in the binding-folding of an IDP prototype. We identify a strong position-dependent D during binding by applying a reaction coordinate that directly measures the fluctuations in a Cartesian configuration space. Using the malleability of the SBM, we modulate the degree of conformational disorder in an isolated IDP and determine complex effects of intrinsic disorder on D varying for different binding stages. Here, D tends to increase with disorder during initial binding but shows a non-monotonic relationship with disorder in terms of a decrease-followed-by-increase in D during the late binding stage. The salt concentration, which correlates with electrostatic interactions via Debye-Hückel theory in our SBM, also modulates D in a stepwise way. The speeding up of diffusion by electrostatic interactions is observed during the formation of the encounter complex at the beginning of binding, while the last diffusive binding dynamics is hindered by non-native salt bridges. Because D describes the diffusive speed locally, which implicitly reflects the roughness of the energy landscape, we are eventually able to portray the binding energy landscape, including that from IDPs' binding, then to binding with partial folding, and finally to rigid docking, as well as that under different environmental salt concentrations. Our theoretical results provide key mechanistic insights into IDPs' binding-folding, which is internally conformation- and externally salt-controlled with respect to diffusion.

The key role of electrostatic interactions in the induced folding in RNA recognition by DCL1-A

The intrinsically disordered protein domain DCL1-A is the first report of a complete double stranded RNA binding domain folding upon binding. DCL1-A recognizes the dsRNA by acquiring a well-folded structure after engagement with its interaction partner. Despite the structural characterization of the interaction complex underlying the recognition of dsRNA has been established, the dynamics of disorder-to-order transitions in the binding process remains elusive. Here we have developed a coarse-grained structure-based model with consideration of electrostatic interactions to explore the mechanism of the coupled folding and binding. Our approach led to remarkable agreements with both experimental and theoretical results. We quantified the global binding-folding landscape, which indicates a synergistic binding induced folding mechanism. We further investigated the effect of electrostatic interactions in this coupled folding and binding process. It reveals that non-native electrostatic interactions dominate the initial stage of the recognition. Our results help improve our understanding of the induced folding of the IDP DCL1-A upon binding to dsRNA. Such methods developed here can be applied for further explorations of the dynamics of coupled folding and binding systems.

Single-Molecule Förster Resonance Energy Transfer Measurement Reveals the Dynamic Partially Ordered Structure of the Epidermal Growth Factor Receptor C-Tail Domain

Intrinsically disordered proteins (IDPs) or regions (IDRs) are thought to exhibit unique functionalities without forming ordered structures. However, these molecular mechanisms are not easily elucidated, partly because of the difficultly of measuring structural information. In this study, we applied the alternative laser excitation (ALEX) method and circular dichroism (CD) spectroscopy to investigate the structure of the C-terminal tail (CTT) domain of the human epidermal growth factor receptor (EGFR). The single-molecule distributions of Förster resonance energy transfer (FRET) obtained by ALEX under solution conditions modified by the addition of potassium chloride (KCl), urea, or guanidinium chloride (GdmCl) allowed us to separately examine the influences of charge interactions and secondary structure formation. The CD spectrum analyses indicated the types of included secondary structure. The results suggested that the structure of the CTT is influenced by secondary structure formation, which is a principally antiparallel β-sheet, rather than by charge interactions and that phosphorylation of the major Grb2-binding sites partially denatures that secondary structure. Our findings suggest that the EGFR CTT might regulate ligand binding kinetics by local β-sheet formation or by the disruption associated with phosphorylation states.

Residual Structure Accelerates Binding of Intrinsically Disordered ACTR by Promoting Efficient Folding upon Encounter

Intrinsically disordered proteins (IDPs) often fold into stable structures upon specific binding. The roles of residual structure of unbound IDPs in coupling binding and folding have been under much debate. While many studies emphasize the importance of conformational flexibility for IDP recognition, it was recently demonstrated that stabilization the N-terminal helix of intrinsically disordered ACTR accelerated its binding to another IDP, NCBD of the CREB-binding protein. To understand how enhancing ACTR helicity accelerates binding, we derived a series of topology-based coarse-grained models that mimicked various ACTR mutants with increasing helical contents and reproduced their NCBD binding affinities. Molecular dynamics simulations were then performed to sample hundreds of reversible coupled binding and folding transitions. The results show that increasing ACTR helicity does not alter the baseline mechanism of synergistic folding, which continues to follow "extended conformational selection" with multiple stages of selection and induced folding. Importantly, these coarse-grained models, while only calibrated based on binding thermodynamics, recapitulate the observed kinetic acceleration with increasing ACTR helicity. However, the residual helices do not enhance the association kinetics via more efficient seeding of productive collisions. Instead, they allow the nonspecific collision complexes to evolve more efficiently into the final bound and folded state, which is the primary source of accelerated association kinetics. Meanwhile, reduced dissociation kinetics with increasing ACTR helicity can be directly attributed to smaller entropic cost of forming the bound state. Altogether, this study provides important mechanistic insights into how residual structure may modulate thermodynamics and kinetics of IDP interactions.

The influence of intrinsic folding mechanism of an unfolded protein on the coupled folding-binding process during target recognition

Intrinsically disordered proteins (IDPs) are extensively involved in dynamic signaling processes which require a high association rate and a high dissociation rate for rapid binding/unbinding events and at the same time a sufficient high affinity for specific recognition. Although the coupled folding-binding processes of IDPs have been extensively studied, it is still impossible to predict whether an unfolded protein is suitable for molecular signaling via coupled folding-binding. In this work, we studied the interplay between intrinsic folding mechanisms and coupled folding-binding process for unfolded proteins through molecular dynamics simulations. We first studied the folding process of three representative IDPs with different folded structures, that is, c-Myb, AF9, and E3 rRNase. We found the folding free energy landscapes of IDPs are downhill or show low barriers. To further study the influence of intrinsic folding mechanism on the binding process, we modulated the folding mechanism of barnase via circular permutation and simulated the coupled folding-binding process between unfolded barnase permutant and folded barstar. Although folding of barnase was coupled to target binding, the binding kinetics was significantly affected by the intrinsic folding free energy barrier, where reducing the folding free energy barrier enhances binding rate up to two orders of magnitude. This accelerating effect is different from previous results which reflect the effect of structure flexibility on binding kinetics. Our results suggest that coupling the folding of an unfolded protein with no/low folding free energy barrier with its target binding may provide a way to achieve high specificity and rapid binding/unbinding kinetics simultaneously.