Smoothing molecular interactions: the "kinetic buffer" effect of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) widely participate in molecular recognition and signaling processes in cells by interacting with other molecules. Compared with ordered proteins, IDPs usually possess stronger intermolecular interactions in binding. As a result, the interface structure of IDPs in complexes is distinct from that of ordered-protein complexes, and this difference may have essential effect on the response to various perturbations in a cell. In this study, we examined the perturbations of intermolecular interactions and temperature on the coupled folding and binding processes of pKID to KIX domains by performing molecular dynamics simulations. By comparing a series of virtual pKID systems with various degree of disorder, we found that the complex stability and the binding kinetics of the disordered systems were less sensitive to the perturbations than the ordered systems. The origin of the lower response sensitivity of IDPs was attributed to their higher flexibility in the complex interface, which was further supported by an analysis on protein complex structures. On the basis of our simulations and results from the literature, we speculate IDPs may not only interact with their biological partners with high specificity and low affinity but also may be resistant to the perturbations in the environment and transmit signals fast and smooth. We proposed to name it the "kinetic buffer" effect.     

Probing metal ion binding and conformational properties of the colicin E9 endonuclease by electrospray ionization time-of-flight mass spectrometry

Nano-electrospray ionization time-of-flight mass spectrometry (ESI-MS) was used to study the conformational consequences of metal ion binding to the colicin E9 endonuclease (E9 DNase) by taking advantage of the unique capability of ESI-MS to allow simultaneous assessment of conformational heterogeneity and metal ion binding. Alterations of charge state distributions on metal ion binding/release were correlated with spectral changes observed in far- and near-UV circular dichroism (CD) and intrinsic tryptophan fluorescence. In addition, hydrogen/deuterium (H/D) exchange experiments were used to probe structural integrity. The present study shows that ESI-MS is sensitive to changes of the thermodynamic stability of E9 DNase as a result of metal ion binding/release in a manner consistent with that deduced from proteolysis and calorimetric experiments. Interestingly, acid-induced release of the metal ion from the E9 DNase causes dramatic conformational instability associated with a loss of fixed tertiary structure, but secondary structure is retained. Furthermore, ESI-MS enabled the direct observation of the noncovalent protein complex of E9 DNase bound to its cognate immunity protein Im9 in the presence and absence of Zn(2+). Gas-phase dissociation experiments of the deuterium-labeled binary and ternary complexes revealed that metal ion binding, not Im9, results in a dramatic exchange protection of E9 DNase in the complex. In addition, our metal ion binding studies and gas-phase dissociation experiments of the ternary E9 DNase-Zn(2+)-Im9 complex have provided further evidence that electrostatic interactions govern the gas phase ion stability.     

Nonnative interactions in coupled folding and binding processes of intrinsically disordered proteins

Proteins function by interacting with other molecules, where both native and nonnative interactions play important roles. Native interactions contribute to the stability and specificity of a complex, whereas nonnative interactions mainly perturb the binding kinetics. For intrinsically disordered proteins (IDPs), which do not adopt rigid structures when being free in solution, the role of nonnative interactions may be more prominent in binding processes due to their high flexibilities. In this work, we investigated the effect of nonnative hydrophobic interactions on the coupled folding and binding processes of IDPs and its interplay with chain flexibility by conducting molecular dynamics simulations. Our results showed that the free-energy profiles became rugged, and intermediate states occurred when nonnative hydrophobic interactions were introduced. The binding rate was initially accelerated and subsequently dramatically decreased as the strength of the nonnative hydrophobic interactions increased. Both thermodynamic and kinetic analysis showed that disordered systems were more readily affected by nonnative interactions than ordered systems. Furthermore, it was demonstrated that the kinetic advantage of IDPs ("fly-casting" mechanism) was enhanced by nonnative hydrophobic interactions. The relationship between chain flexibility and protein aggregation is also discussed.     

Ordered Disorder of the Astrocytic Dystrophin-Associated Protein Complex in the Norm and Pathology

The abundance and potential functional roles of intrinsically disordered regions in aquaporin-4, Kir4.1, a dystrophin isoforms Dp71, α-1 syntrophin, and α-dystrobrevin; i.e., proteins constituting the functional core of the astrocytic dystrophin-associated protein complex (DAPC), are analyzed by a wealth of computational tools. The correlation between protein intrinsic disorder, single nucleotide polymorphisms (SNPs) and protein function is also studied together with the peculiarities of structural and functional conservation of these proteins. Our study revealed that the DAPC members are typical hybrid proteins that contain both ordered and intrinsically disordered regions. Both ordered and disordered regions are important for the stabilization of this complex. Many disordered binding regions of these five proteins are highly conserved among vertebrates. Conserved eukaryotic linear motifs and molecular recognition features found in the disordered regions of five protein constituting DAPC likely enhance protein-protein interactions that are required for the cellular functions of this complex. Curiously, the disorder-based binding regions are rarely affected by SNPs suggesting that these regions are crucial for the biological functions of their corresponding proteins.     

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.     

A Large Intrinsically Disordered Region in SKIP and Its Disorder-Order Transition Induced by PPIL1 Binding Revealed by NMR

Intrinsically disordered proteins or protein regions play an important role in fundamental biological processes. During spliceosome activation, a large structural rearrangement occurs. The Prp19 complex and related factors are involved in the catalytic activation of the spliceosome. Recent mass spectrometric analyses have shown that Ski interaction protein (SKIP) and peptidylprolyl isomerase-like protein 1 (PPIL1) are Prp19-related factors that constitute the spliceosome B, B*, and C complexes. Here, we report that a highly flexible region of SKIP (SKIPN, residues 59–129) is intrinsically disordered. Upon binding to PPIL1, SKIPN undergoes a disorder-order transition. A highly conserved fragment of SKIP (residues 59–79) called the PPIL1-binding fragment (PBF) was sufficient to bind PPIL1. The structure of PBF·PPIL1 complex, solved by NMR, shows that PBF exhibits an ordered structure and interacts with PPIL1 through electrostatic and hydrophobic interactions. Three subfragments in the PBF (residues 59–67, 68–73, and 74–79) show hook-like backbone structure, and interactions between these subfragments are necessary for PBF·PPIL1 complex formation. PPIL1 is a cyclophilin family protein. It is recruited by SKIP into the spliceosome by a region other than the peptidylprolyl isomerase active site. This enables the active site of PPIL1 to remain open in the complex and still function as a peptidylprolyl cis/trans-isomerase or molecular chaperon to facilitate the folding of other proteins in the spliceosomes. The large disordered region in SKIP provides an interaction platform. Its disorder-order transition, induced by PPIL1 binding, may adapt the requirement for a large structural rearrangement occurred in the activation of spliceosome.     

Importance of electrostatic interactions in the association of intrinsically disordered histone chaperone Chz1 and histone H2A.Z-H2B

Histone chaperones facilitate assembly and disassembly of nucleosomes. Understanding the process of how histone chaperones associate and dissociate from the histones can help clarify their roles in chromosome metabolism. Some histone chaperones are intrinsically disordered proteins (IDPs). Recent studies of IDPs revealed that the recognition of the biomolecules is realized by the flexibility and dynamics, challenging the century-old structure-function paradigm. Here we investigate the binding between intrinsically disordered chaperone Chz1 and histone variant H2A.Z-H2B by developing a structure-based coarse-grained model, in which Debye-Hückel model is implemented for describing electrostatic interactions due to highly charged characteristic of Chz1 and H2A.Z-H2B. We find that major structural changes of Chz1 only occur after the rate-limiting electrostatic dominant transition state and Chz1 undergoes folding coupled binding through two parallel pathways. Interestingly, although the electrostatic interactions stabilize bound complex and facilitate the recognition at first stage, the rate for formation of the complex is not always accelerated due to slow escape of conformations with non-native electrostatic interactions at low salt concentrations. Our studies provide an ionic-strength-controlled binding/folding mechanism, leading to a cooperative mechanism of "local collapse or trapping" and "fly-casting" together and a new understanding of the roles of electrostatic interactions in IDPs' binding.     

Molecular Recognition Features in Zika Virus Proteome

Viruses have compact genomes that encode limited number of proteins in comparison to other biological entities. Interestingly, viral proteins have shown natural abundance of either completely disordered proteins that are recognized as intrinsically disorder proteins (IDPs) or partially disordered segments known as intrinsically disordered protein regions (IDPRs). IDPRs are involved in interactions with multiple binding partners to accomplish signaling, regulation, and control functions in cells. Tuning of IDPs and IDPRs are mediated through post-translational modification and alternative splicing. Often, the interactions of IDPRs with their binding protein partner(s) lead to transition from the state of disorder to ordered form. Such interaction-prone protein IDPRs are identified as molecular recognition features (MoRFs). Molecular recognition is an important initial step for the biomolecular interactions and their functional proceedings. Although previous studies have established occurrence of the IDPRs in Zika virus proteome, which provide the functional diversity and structural plasticity to viral proteins, the MoRF analysis has not been performed as of yet. Many computational methods have been developed for the identification of the MoRFs in protein sequences including ANCHOR, MoRFpred, DISOPRED3, and MoRFchibi_web server. In the current study, we have investigated the presence of MoRF regions in structural and non-structural proteins of Zika virus using an aforementioned set of computational techniques. Furthermore, we have experimentally validated the intrinsic disorderness of NS2B cofactor region of NS2B–NS3 protease. NS2B has one of the longest MoRF regions in Zika virus proteome. In future, this study may provide valuable information while investigating the virus host protein interaction networks.     

Molecular Mechanisms of Tight Binding through Fuzzy Interactions

Many intrinsically disordered proteins (IDPs) form fuzzy complexes upon binding to their targets. Although many IDPs are weakly bound in fuzzy complexes, some IDPs form high-affinity complexes. One example is the nonstructural protein 1 (NS1) of the 1918 Spanish influenza A virus, which hijacks cellular CRKII through the strong binding affinity (K_d ～10 nM) of its proline-rich motif (PRM^NS1) to the N-terminal Src-homology 3 domain of CRKII. However, its molecular mechanism remains elusive. Here, we examine the interplay between structural disorder of a bound PRM^NS1 and its long-range electrostatic interactions. Using x-ray crystallography and NMR spectroscopy, we found that PRM^NS1 retains substantial conformational flexibility in the bound state. Moreover, molecular dynamics simulations showed that structural disorder of the bound PRM^NS1 increases the number of electrostatic interactions and decreases the mean distances between the positively charged residues in PRM^NS1 and the acidic residues in the N-terminal Src-homology 3 domain. These results are analyzed using a polyelectrostatic model. Our results provide an insight into the molecular recognition mechanism for a high-affinity fuzzy complex.     

Disordered nucleiome: Abundance of intrinsic disorder in the DNA‐ and RNA‐binding proteins in 1121 species from Eukaryota,Bacteria and Archaea

Intrinsically disordered proteins (IDPs) are abundant in various proteomes, where they play numerous important roles and complement biological activities of ordered proteins. Among functions assigned to IDPs are interactions with nucleic acids. However, often, such assignments are made based on the guilty‐by‐association principle. The validity of the extension of these correlations to all nucleic acid binding proteins has never been analyzed on a large scale across all domains of life. To fill this gap, we perform a comprehensive computational analysis of the abundance of intrinsic disorder and intrinsically disordered domains in nucleiomes (～548 000 nucleic acid binding proteins) of 1121 species from Archaea, Bacteria and Eukaryota. Nucleiome is a whole complement of proteins involved in interactions with nucleic acids. We show that relative to other proteins in the corresponding proteomes, the DNA‐binding proteins have significantly increased disorder content and are significantly enriched in disordered domains in Eukaryotes but not in Archaea and Bacteria. The RNA‐binding proteins are significantly enriched in the disordered domains in Bacteria, Archaea and Eukaryota, while the overall abundance of disorder in these proteins is significantly increased in Bacteria, Archaea, animals and fungi. The high abundance of disorder in nucleiomes supports the notion that the nucleic acid binding proteins often require intrinsic disorder for their functions and regulation.     

Mapping the transition state for a binding reaction between ancient intrinsically disordered proteins

Intrinsically disordered protein domains often have multiple binding partners. It is plausible that the strength of pairing with specific partners evolves from an initial low affinity to a higher affinity. However, little is known about the molecular changes in the binding mechanism that would facilitate such a transition. We previously showed that the interaction between two intrinsically disordered domains, NCBD and CID, likely emerged in an ancestral deuterostome organism as a low-affinity interaction that subsequently evolved into a higher-affinity interaction before the radiation of modern vertebrate groups. Here we map native contacts in the transition states of the low-affinity ancestral and high-affinity human NCBD/CID interactions. We show that the coupled binding and folding mechanism is overall similar but with a higher degree of native hydrophobic contact formation in the transition state of the ancestral complex and more heterogeneous transient interactions, including electrostatic pairings, and an increased disorder for the human complex. Adaptation to new binding partners may be facilitated by this ability to exploit multiple alternative transient interactions while retaining the overall binding and folding pathway.     

Electrostatically accelerated coupled binding and folding of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are now recognized to be prevalent in biology, and many potential functional benefits have been discussed. However, the frequent requirement of peptide folding in specific interactions of IDPs could impose a kinetic bottleneck, which could be overcome only by efficient folding upon encounter. Intriguingly, existing kinetic data suggest that specific binding of IDPs is generally no slower than that of globular proteins. Here, we exploited the cell cycle regulator p27(Kip1) (p27) as a model system to understand how IDPs might achieve efficient folding upon encounter for facile recognition. Combining experiments and coarse-grained modeling, we demonstrate that long-range electrostatic interactions between enriched charges on p27 and near its binding site on cyclin A not only enhance the encounter rate (i.e., electrostatic steering) but also promote folding-competent topologies in the encounter complexes, allowing rapid subsequent formation of short-range native interactions en route to the specific complex. In contrast, nonspecific hydrophobic interactions, while hardly affecting the encounter rate, can significantly reduce the efficiency of folding upon encounter and lead to slower binding kinetics. Further analysis of charge distributions in a set of known IDP complexes reveals that, although IDP binding sites tend to be more hydrophobic compared to the rest of the target surface, their vicinities are frequently enriched with charges to complement those on IDPs. This observation suggests that electrostatically accelerated encounter and induced folding might represent a prevalent mechanism for promoting facile IDP recognition.     

Calorimetric dissection of colicin DNase--immunity protein complex specificity

We explore the thermodynamic strategies used to achieve specific, high-affinity binding within a family of conserved protein-protein complexes. Protein-protein interactions are often stabilized by a conserved interfacial hotspot that serves as the anchor for the complex, with neighboring variable residues providing specificity. A key question for such complexes is the thermodynamic basis for specificity given the dominance of the hotspot. We address this question using, as our model, colicin endonuclease (DNase)-immunity (Im) protein complexes. In this system, cognate and noncognate complexes alike share the same mechanism of association and binding hotspot, but cognate complexes (K(d) approximately 10(-)(14) M) are orders of magnitude more stable than noncognate complexes (10(6)-10(10)-fold discrimination), largely because of a much slower rate of dissociation. Using isothermal titration calorimetry (ITC), we investigated the changes in enthalpy (DeltaH), entropy (-TDeltaS), and heat capacity (DeltaC(p)) accompanying binding of each Im protein (Im2, Im7, Im8, and Im9) to the DNase domains of colicins E2, E7, E8, and E9, in the context of both cognate and noncognate complexes. The data show that specific binding to the E2, E7, and E8 DNases is enthalpically driven but entropically driven for the E9 DNase. Analysis of DeltaC(p), a measure of the change in structural fluctuation upon complexation, indicates that E2, E7, and E8 DNase specificity is coupled to structural changes within cognate complexes that are consistent with a reduction in the conformational dynamics of these complexes. In contrast, E9 DNase specificity appears coupled to the exclusion of water molecules, consistent with the nonpolar nature of the interface of this complex. The work highlights that although protein-protein interactions may be centered on conserved structural epitopes the thermodynamic mechanism underpinning binding specificity can vary considerably.     

Investigating the disorder–order transition of calmodulin binding domain upon binding calmodulin using molecular dynamics simulation

We have studied the conformational transition of the calmodulin binding domains (CBD) in several calmodulin‐binding kinases, in which CBD changes from the disordered state to the ordered state when binding with calmodulin (CaM). Targeted molecular dynamics simulation was used to investigate the binding process of CaM and CBD of CaM‐dependent kinase I (CaMKI–CBD). The results show that CaMKI–CBD began to form an α‐helix and the interaction free energy between CaM and CaMKI–CBD increased once CaM fully encompassed CaMKI–CBD. Two series of CaM/CBD complex systems, including the complexes of CaM with the initially disordered and the final ordered CBD, were constructed to study the interaction using molecular dynamics simulations. Our analyses suggest that the VDW interaction plays a dominant role in CaM/CBD binding and is a key factor in the disorder–order transition of CBD. Additionally, the entropy effect is not in favor of the formation of the CaM/CBD complex, which is consistent with the experimental evidence. Based on the results, it appears that the CBD conformational change from a non‐compact extended structure to compact α‐helix is critical in gaining a favorable VDW interaction and interaction free energy. Copyright © 2009 John Wiley & Sons, Ltd.     

Multisite phosphorylation and binding alter conformational dynamics of the 4E-BP2 protein

Intrinsically disordered proteins (IDPs) play critical roles in regulatory protein interactions, but detailed structural/dynamic characterization of their ensembles remain challenging, both in isolation and when they form dynamic “fuzzy” complexes. Such is the case for mRNA cap-dependent translation initiation, which is regulated by the interaction of the predominantly folded eukaryotic initiation factor 4E (eIF4E) with the intrinsically disordered eIF4E binding proteins (4E-BPs) in a phosphorylation-dependent manner. Single-molecule Förster resonance energy transfer showed that the conformational changes of 4E-BP2 induced by binding to eIF4E are non-uniform along the sequence; while a central region containing both motifs that bind to eIF4E expands and becomes stiffer, the C-terminal region is less affected. Fluorescence anisotropy decay revealed a non-uniform segmental flexibility around six different labeling sites along the chain. Dynamic quenching of these fluorescent probes by intrinsic aromatic residues measured via fluorescence correlation spectroscopy report on transient intra- and inter-molecular contacts on nanosecond-to-microsecond timescales. Upon hyperphosphorylation, which induces folding of ～40 residues in 4E-BP2, the quenching rates decreased at most labeling sites. The chain dynamics around sites in the C-terminal region far away from the two binding motifs significantly increased upon binding to eIF4E, suggesting that this region is also involved in the highly dynamic 4E-BP2:eIF4E complex. Our time-resolved fluorescence data paint a sequence-level rigidity map of three states of 4E-BP2 differing in phosphorylation or binding status and distinguish regions that form contacts with eIF4E. This study adds complementary structural and dynamics information to recent studies of 4E-BP2, and it constitutes an important step toward a mechanistic understanding of this important IDP via integrative modeling.     

Electrostatic control of calcineurin's intrinsically-disordered regulatory domain binding to calmodulin

Calcineurin (CaN) is a serine/threonine phosphatase that regulates a variety of physiological and pathophysiological processes in mammalian tissue. The calcineurin (CaN) regulatory domain (RD) is responsible for regulating the enzyme's phosphatase activity, and is believed to be highly-disordered when inhibiting CaN, but undergoes a disorder-to-order transition upon diffusion-limited binding with the regulatory protein calmodulin (CaM). The prevalence of polar and charged amino acids in the regulatory domain (RD) suggests electrostatic interactions are involved in mediating calmodulin (CaM) binding, yet the lack of atomistic-resolution data for the bound complex has stymied efforts to probe how the RD sequence controls its conformational ensemble and long-range attractions contribute to target protein binding. In the present study, we investigated via computational modeling the extent to which electrostatics and structural disorder facilitate CaM/CaN association kinetics. Specifically, we examined several RD constructs that contain the CaM binding region (CAMBR) to characterize the roles of electrostatics versus conformational diversity in controlling diffusion-limited association rates, via microsecond-scale molecular dynamics (MD) and Brownian dynamic (BD) simulations. Our results indicate that the RD amino acid composition and sequence length influence both the dynamic availability of conformations amenable to CaM binding, as well as long-range electrostatic interactions to steer association. These findings provide intriguing insight into the interplay between conformational diversity and electrostatically-driven protein-protein association involving CaN, which are likely to extend to wide-ranging diffusion-limited processes regulated by intrinsically-disordered proteins.     

Interplay between binding affinity and kinetics in protein–protein interactions

To clarify the interplay between the binding affinity and kinetics of protein–protein interactions, and the possible role of intrinsically disordered proteins in such interactions, molecular simulations were carried out on 20 protein complexes. With bias potential and reweighting techniques, the free energy profiles were obtained under physiological affinities, which showed that the bound‐state valley is deep with a barrier height of 12 ? 33 RT. From the dependence of the affinity on interface interactions, the entropic contribution to the binding affinity is approximated to be proportional to the interface area. The extracted dissociation rates based on the Arrhenius law correlate reasonably well with the experimental values (Pearson correlation coefficient R = 0.79). For each protein complex, a linear free energy relationship between binding affinity and the dissociation rate was confirmed, but the distribution of the slopes for intrinsically disordered proteins showed no essential difference with that observed for ordered proteins. A comparison with protein folding was also performed. Proteins 2016; 84:920–933. © 2016 Wiley Periodicals, Inc.