Kinetics of Folding and Binding of an Intrinsically Disordered Protein: The Inhibitor of Yeast Aspartic Proteinase YPrA

Charge Distribution Patterns of IA3 Impact Conformational Expansion and Hydration Diffusivity of the Disordered Ensemble

IA3 is a 68 amino acid natural peptide/protein inhibitor of yeast aspartic proteinase A (YPRA) that is intrinsically disordered in solution with induced N-terminal helicity when in the protein complex with YPRA. Based on the intrinsically disordered protein (IDP) parameters of fractional net charge (FNC), net charge density per residue (NCPR), and charge patterning (κ), the two domains of IA3 are defined to occupy different domains within conformationally based subclasses of IDPs, thus making IA3 a bimodal domain IDP. Site-directed spin labeling (SDSL) electron paramagnetic resonance (EPR) spectroscopy and low-field Overhauser dynamic nuclear polarization (ODNP) spectroscopy results show that these two domains possess different degrees of compaction and hydration diffusivity behavior. This work suggests that SDSL EPR line shapes, analyzed in terms of their local tumbling volume (VL), provide insights into the compaction of the unstructured IDP ensemble in solution and that protein sequence and net charge distribution patterns within a conformational subclass can impact bound water hydration dynamics, thus possibly offering an alternative thermodynamic property that can encode conformational binding and behavior of IDPs and liquid-liquid phase separations.

Hydrogen Bonding Compensation on the Convex Solvent-Exposed Helical Face of IA3, an Intrinsically Disordered Protein

Saccharomyces cerevisiae IA₃ is a 68 amino acid peptide inhibitor of yeast proteinase A (YPRA) characterized as a random coil when in solution, folding into an N-terminal amphipathic alpha helix for residues 2-32 when bound to YPRA, with residues 33-68 unresolved in the crystal complex. Circular dichroism (CD) spectroscopy results show that amino acid substitutions that remove hydrogen-bonding interactions observed within the hydrophilic face of the N-terminal domain (NTD) of IA₃-YPRA crystal complex reduce the 2,2,2-trifluoroethanol (TFE)-induced helical transition in solution. Although nearly all substitutions decreased TFE-induced helicity compared to wild-type (WT), each construct did retain helical character in the presence of 30% (v/v) TFE and retained disorder in the absence of TFE. The NTDs of 8 different Saccharomyces species have nearly identical amino acid sequences, indicating that the NTD of IA₃ may be highly evolved to adopt a helical fold when bound to YPRA and in the presence of TFE but remain unstructured in solution. Only one natural amino acid substitution explored within the solvent-exposed face of the NTD of IA₃ induced TFE-helicity greater than the WT sequence. However, chemical modification of a cysteine by a nitroxide spin label that contains an acetamide side chain did enhance TFE-induced helicity. This finding suggests that non-natural amino acids that can increase hydrogen bonding or alter hydration through side-chain interactions may be important to consider when rationally designing intrinsically disordered proteins (IDPs) with varied biotechnological applications.

Unveiling induced folding of intrinsically disordered proteins - Protein engineering, frustration and emerging themes

Intrinsically disordered proteins (IDPs) can be generally described as a class of proteins that lack a well-defined ordered structure in isolation at physiological conditions. Upon binding to their physiological ligands, IDPs typically undergo a disorder-to-order transition, which may or may not lead to the complete folding of the IDP. In this short review, we focus on some of the key findings pertaining to the mechanisms of such induced folding. In particular, first we describe the general features of the reaction; then, we discuss some of the most remarkable findings obtained from applying protein engineering in synergy with kinetic studies to induced folding; and finally, we offer a critical view on some of the emerging themes when considering the structural heterogeneity of IDPs vis-à-vis to their inherent frustration.

Intrinsic Dynamics of Protein-Peptide Unbinding

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

Understanding Binding-Induced Folding by Temperature Jump

Temperature jump is a powerful technique for the characterization of fast kinetics and can be readily employed to understand both binding and folding reactions. Here we summarize briefly a temperature-jump prototypical experiment between an intrinsically disordered protein and its physiological partner. The model used is the N_TAIL domain from Measles virus Nucleoprotein and its natural ligand, the globular P_XD domain from Measles virus Phosphoprotein. We recapitulate how to set up the experiment and how to analyze data in order to extract the kinetic parameters of the reaction.

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.

Templated folding of intrinsically disordered proteins

Much of our current knowledge of biological chemistry is founded in the structure-function relationship, whereby sequence determines structure that determines function. Thus, the discovery that a large fraction of the proteome is intrinsically disordered, while being functional, has revolutionized our understanding of proteins and raised new and interesting questions. Many intrinsically disordered proteins (IDPs) have been determined to undergo a disorder-to-order transition when recognizing their physiological partners, suggesting that their mechanisms of folding are intrinsically different from those observed in globular proteins. However, IDPs also follow some of the classic paradigms established for globular proteins, pointing to important similarities in their behavior. In this review, we compare and contrast the folding mechanisms of globular proteins with the emerging features of binding-induced folding of intrinsically disordered proteins. Specifically, whereas disorder-to-order transitions of intrinsically disordered proteins appear to follow rules of globular protein folding, such as the cooperative nature of the reaction, their folding pathways are remarkably more malleable, due to the heterogeneous nature of their folding nuclei, as probed by analysis of linear free-energy relationship plots. These insights have led to a new model for the disorder-to-order transition in IDPs termed "templated folding," whereby the binding partner dictates distinct structural transitions en route to product, while ensuring a cooperative folding.

Transition path dynamics in the binding of intrinsically disordered proteins: A simulation study

Association of proteins and other biopolymers is a ubiquitous process in living systems. Recent single-molecule measurements probe the dynamics of association in unprecedented detail by measuring the properties of association transition paths, i.e., short segments of molecular trajectories between the time the proteins are close enough to interact and the formation of the final complex. Interpretation of such measurements requires adequate models for describing the dynamics of experimental observables. In an effort to develop such models, here we report a simulation study of the association dynamics of two oppositely charged, disordered polymers. We mimic experimental measurements by monitoring intermonomer distances, which we treat as "experimental reaction coordinates." While the dynamics of the distance between the centers of mass of the molecules is found to be memoryless and diffusive, the dynamics of the experimental reaction coordinates displays significant memory and can be described by a generalized Langevin equation with a memory kernel. We compute the most commonly measured property of transition paths, the distribution of the transition path time, and show that, despite the non-Markovianity of the underlying dynamics, it is well approximated as one-dimensional diffusion in the potential of mean force provided that an apparent value of the diffusion coefficient is used. This apparent value is intermediate between the slow (low frequency) and fast (high frequency) limits of the memory kernel. We have further studied how the mean transition path time depends on the ionic strength and found only weak dependence despite strong electrostatic attraction between the polymers.

Binding-induced folding under unfolding conditions: Switching between induced fit and conformational selection mechanisms

The chemistry of protein-ligand binding is the basis of virtually every biological process. Ligand binding can be essential for a protein to function in the cell by stabilizing or altering the conformation of a protein, particularly for partially or completely unstructured proteins. However, the mechanisms by which ligand binding impacts disordered proteins or influences the role of disorder in protein folding is not clear. To gain insight into this question, the mechanism of folding induced by the binding of a Pro-rich peptide ligand to the SH3 domain of phosphatidylinositol 3-kinase unfolded in the presence of urea has been studied using kinetic methods. Under strongly denaturing conditions, folding was found to follow a conformational selection (CS) mechanism. However, under mildly denaturing conditions, a ligand concentration-dependent switch in the mechanism was observed. The folding mechanism switched from being predominantly a CS mechanism at low ligand concentrations to being predominantly an induced fit (IF) mechanism at high ligand concentrations. The switch in the mechanism manifests itself as an increase in the reaction flux along the IF pathway at high ligand concentrations. The results indicate that, in the case of intrinsically disordered proteins too, the folding mechanism is determined by the concentration of the ligand that induces structure formation.

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.

How Robust Is the Mechanism of Folding-Upon-Binding for an Intrinsically Disordered Protein?

The mechanism of interaction of an intrinsically disordered protein (IDP) with its physiological partner is characterized by a disorder-to-order transition in which a recognition and a binding step take place. Even if the mechanism is quite complex, IDPs tend to bind their partner in a cooperative manner such that it is generally possible to detect experimentally only the disordered unbound state and the structured complex. The interaction between the disordered C-terminal domain of the measles virus nucleoprotein (N_TAIL) and the X domain (XD) of the viral phosphoprotein allows us to detect and quantify the two distinct steps of the overall reaction. Here, we analyze the robustness of the folding of N_TAIL upon binding to XD by measuring the effect on both the folding and binding steps of N_TAIL when the structure of XD is modified. Because it has been shown that wild-type XD is structurally heterogeneous, populating an on-pathway intermediate under native conditions, we investigated the binding to 11 different site-directed variants of N_TAIL of one particular variant of XD (I504A XD) that populates only the native state. Data reveal that the recognition and the folding steps are both affected by the structure of XD, indicating a highly malleable pathway. The experimental results are briefly discussed in the light of previous experiments on other IDPs.

Membrane interactions of intrinsically disordered proteins: The example of alpha-synuclein

Peripheral membrane proteins associate reversibly with biological membranes that, compared to protein binding partners, are structurally labile and devoid of specific binding pockets. Membranes in different subcellular compartments vary primarily in their chemical composition and physical properties, and recognition of these features is therefore critical for allowing such proteins to engage their proper membrane targets. Intrinsically disordered proteins (IDPs) are well-suited to accomplish this task using highly specific and low- to moderate-affinity interactions governed by recognition principles that are both similar to and different from those that mediate the membrane interactions of rigid proteins. IDPs have also evolved multiple mechanisms to regulate membrane (and other) interactions and achieve their impressive functional diversity. Moreover, IDP-membrane interactions may have a kinetic advantage in fast processes requiring rapid control of such interactions, such as synaptic transmission or signaling. Herein we review the biophysics, regulation and functional implications of IDP-membrane interactions and include a brief overview of some of the methods that can be used to study such interactions. At each step, we use the example of alpha-synuclein, a protein involved in the pathogenesis of Parkinson's disease and one of the best characterized membrane-binding IDP, to illustrate some of the principles discussed.

Transition path times of coupled folding and binding reveal the formation of an encounter complex

The association of biomolecules is the elementary event of communication in biology. Most mechanistic information of how the interactions between binding partners form or break is, however, hidden in the transition paths, the very short parts of the molecular trajectories from the encounter of the two molecules to the formation of a stable complex. Here we use single-molecule spectroscopy to measure the transition path times for the association of two intrinsically disordered proteins that form a folded dimer upon binding. The results reveal the formation of a metastable encounter complex that is electrostatically favored and transits to the final bound state within tens of microseconds. Such measurements thus open a new window into the microscopic events governing biomolecular interactions.

Structural metamorphism and polymorphism in proteins on the brink of thermodynamic stability

The classical view of the structure-function paradigm advanced by Anfinsen in the 1960s is that a protein's function is inextricably linked to its three-dimensional structure and is encrypted in its amino acid sequence. However, it is now known that a significant fraction of the proteome consists of intrinsically disordered proteins (IDPs). These proteins populate a polymorphic ensemble of conformations rather than a unique structure but are still capable of performing biological functions. At the boundary, between well-ordered and inherently disordered states are proteins that are on the brink of stability, either weakly stable ordered systems or disordered but on the verge of being stable. In such marginal states, even relatively minor changes can significantly alter the energy landscape, leading to large-scale conformational remodeling. Some proteins on the edge of stability are metamorphic, with the capacity to switch from one fold topology to another in response to an environmental trigger (e.g., pH, temperature/salt, redox). Many IDPs, on the other hand, are marginally unstable such that small perturbations (e.g., phosphorylation, ligands) tip the balance over to a range of ordered, partially ordered, or even more disordered states. In general, the structural transitions described by metamorphic fold switches and polymorphic IDPs possess a number of common features including low or diminished stability, large-scale conformational changes, critical disordered regions, latent or attenuated binding sites, and expansion of function. We suggest that these transitions are, therefore, conceptually and mechanistically analogous, representing adjacent regions in the continuum of order/disorder transitions.

Unified understanding of folding and binding mechanisms of globular and intrinsically disordered proteins

Extensive experimental and theoretical studies have advanced our understanding of the mechanisms of folding and binding of globular proteins, and coupled folding and binding of intrinsically disordered proteins (IDPs). The forces responsible for conformational changes and binding are common in both proteins; however, these mechanisms have been separately discussed. Here, we attempt to integrate the mechanisms of coupled folding and binding of IDPs, folding of small and multi-subdomain proteins, folding of multimeric proteins, and ligand binding of globular proteins in terms of conformational selection and induced-fit mechanisms as well as the nucleation–condensation mechanism that is intermediate between them. Accumulating evidence has shown that both the rate of conformational change and apparent rate of binding between interacting elements can determine reaction mechanisms. Coupled folding and binding of IDPs occurs mainly by induced-fit because of the slow folding in the free form, while ligand binding of globular proteins occurs mainly by conformational selection because of rapid conformational change. Protein folding can be regarded as the binding of intramolecular segments accompanied by secondary structure formation. Multi-subdomain proteins fold mainly by the induced-fit (hydrophobic collapse) mechanism, as the connection of interacting segments enhances the binding (compaction) rate. Fewer hydrophobic residues in small proteins reduce the intramolecular binding rate, resulting in the nucleation–condensation mechanism. Thus, the folding and binding of globular proteins and IDPs obey the same general principle, suggesting that the coarse-grained, statistical mechanical model of protein folding is promising for a unified theoretical description of all mechanisms.

Toward increased concentration sensitivity for continuous wave EPR investigations of spin-labeled biological macromolecules at high fields

High-field, high-frequency electron paramagnetic resonance (EPR) spectroscopy at W-(∼94 GHz) and D-band (∼140 GHz) is important for investigating the conformational dynamics of flexible biological macromolecules because this frequency range has increased spectral sensitivity to nitroxide motion over the 100 ps to 2 ns regime. However, low concentration sensitivity remains a roadblock for studying aqueous samples at high magnetic fields. Here, we examine the sensitivity of a non-resonant thin-layer cylindrical sample holder, coupled to a quasi-optical induction-mode W-band EPR spectrometer (HiPER), for continuous wave (CW) EPR analyses of: (i) the aqueous nitroxide standard, TEMPO; (ii) the unstructured to α-helical transition of a model IDP protein; and (iii) the base-stacking transition in a kink-turn motif of a large 232 nt RNA. For sample volumes of ∼50 μL, concentration sensitivities of 2-20 μM were achieved, representing a ∼10-fold enhancement compared to a cylindrical TE011 resonator on a commercial Bruker W-band spectrometer. These results therefore highlight the sensitivity of the thin-layer sample holders employed in HiPER for spin-labeling studies of biological macromolecules at high fields, where applications can extend to other systems that are facilitated by the modest sample volumes and ease of sample loading and geometry.

Induced folding in RNA recognition by Arabidopsis thaliana DCL1

DCL1 is the ribonuclease that carries out miRNA biogenesis in plants. The enzyme has two tandem double stranded RNA binding domains (dsRBDs) in its C-terminus. Here we show that the first of these domains binds precursor RNA fragments when isolated and cooperates with the second domain in the recognition of substrate RNA. Remarkably, despite showing RNA binding activity, this domain is intrinsically disordered. We found that it acquires a folded conformation when bound to its substrate, being the first report of a complete dsRBD folding upon binding. The free unfolded form shows tendency to adopt folded conformations, and goes through an unfolded bound state prior to the folding event. The significance of these results is discussed by comparison with the behavior of other dsRBDs.

Deciphering the mechanisms of binding induced folding at nearly atomic resolution: The Φ value analysis applied to IDPs

The Φ value analysis is a method to analyze the structure of metastable states in reaction pathways. Such a methodology is based on the quantitative analysis of the effect of point mutations on the kinetics and thermodynamics of the probed reaction. The Φ value analysis is routinely used in protein folding studies and is potentially an extremely powerful tool to analyze the mechanism of binding induced folding of intrinsically disordered proteins. In this review we recapitulate the key equations and experimental advices to perform the Φ value analysis in the perspective of the possible caveats arising in intrinsically disordered systems. Finally, we briefly discuss some few examples already available in the literature.

Interplay between partner and ligand facilitates the folding and binding of an intrinsically disordered protein

Protein-protein interactions are at the heart of regulatory and signaling processes in the cell. In many interactions, one or both proteins are disordered before association. However, this disorder in the unbound state does not prevent many of these proteins folding to a well-defined, ordered structure in the bound state. Here we examine a typical system, where a small disordered protein (PUMA, p53 upregulated modulator of apoptosis) folds to an α-helix when bound to a groove on the surface of a folded protein (MCL-1, induced myeloid leukemia cell differentiation protein). We follow the association of these proteins using rapid-mixing stopped flow, and examine how the kinetic behavior is perturbed by denaturant and carefully chosen mutations. We demonstrate the utility of methods developed for the study of monomeric protein folding, including β-Tanford values, Leffler α, Φ-value analysis, and coarse-grained simulations, and propose a self-consistent mechanism for binding. Folding of the disordered protein before binding does not appear to be required and few, if any, specific interactions are required to commit to association. The majority of PUMA folding occurs after the transition state, in the presence of MCL-1. We also examine the role of the side chains of folded MCL-1 that make up the binding groove and find that many favor equilibrium binding but, surprisingly, inhibit the association process.

Why are the truncated cyclin Es more effective CDK2 activators than the full-length isoforms?

Cell cycle regulating enzymes, CDKs, become activated upon association with their regulatory proteins, cyclins. The G1 cyclin, cyclin E, is overexpressed and present in low molecular weight (LMW) isoforms in breast cancer cells and tumor tissues. In vivo and in vitro studies have shown that these LMW isoforms of cyclin E hyperactivate CDK2 and accelerate the G1-S phase of cell division. The molecular basis of CDK2 hyperactivation due to LMW cyclin E isoforms in cancer cells is, however, unknown. Here, we employ a computational approach, combining homology modeling, bioinformatics analyses, molecular dynamics (MD) simulations, and principal component analyses to unravel the key structural features of CDK2-bound full-length and LMW isoforms of cyclin E1 and correlate those features to their differential activity. Results suggest that the missing N- and C-terminal regions of the cyclin E LMW isoforms constitute the Nuclear Localization Sequence (NLS) and PEST domains and are intrinsically disordered. These regions, when present in the full-length cyclin E/CDK2 complex, weaken the cyclin-CDK interface packing due to the loss of a large number of key interface interactions. Such weakening is manifested in the decreased contact area and increased solvent accessibility at the interface and also by the absence of concerted motions between the two partner proteins in the full-length complex. More effective packing and interactions between CDK2 and LMW cyclin E isoforms, however, produce more efficient protein-protein complexes that accelerate the cell division processes in cancer cells, where these cyclin E isoforms are overexpressed.

Continuous wave W- and D-band EPR spectroscopy offer "sweet-spots" for characterizing conformational changes and dynamics in intrinsically disordered proteins

Site-directed spin labeling (SDSL) electron paramagnetic resonance (EPR) spectroscopy is a powerful tool for characterizing conformational sampling and dynamics in biological macromolecules. Here we demonstrate that nitroxide spectra collected at frequencies higher than X-band (∼9.5 GHz) have sensitivity to the timescale of motion sampled by highly dynamic intrinsically disordered proteins (IDPs). The 68 amino acid protein IA3, was spin-labeled at two distinct sites and a comparison of X-band, Q-band (35 GHz) and W-band (95 GHz) spectra are shown for this protein as it undergoes the helical transition chemically induced by tri-fluoroethanol. Experimental spectra at W-band showed pronounced line shape dispersion corresponding to a change in correlation time from ∼0.3 ns (unstructured) to ∼0.6 ns (α-helical) as indicated by comparison with simulations. Experimental and simulated spectra at X- and Q-bands showed minimal dispersion over this range, illustrating the utility of SDSL EPR at higher frequencies for characterizing structural transitions and dynamics in IDPs.

The mechanism of binding of the KIX domain to the mixed lineage leukemia protein and its allosteric role in the recognition of c-Myb

The KIX domain is a mediator of the interaction between different transcription factors. This complex function is carried out via two distinct binding sites located on opposite sides of the protein; namely, the 'c-Myb site' and the 'MLL site', named after their characteristic ligands-the transactivation domain of c-Myb and the mixed lineage leukemia protein (MLL). Both these ligands are unstructured in isolation and fold only upon binding, posing the KIX domain as an ideal candidate to explore the binding induced folding reaction of intrinsically unstructured proteins. Here, we complement the recent kinetic description on the interaction between KIX and c-Myb, by characterizing the binding kinetics between KIX and MLL, at different pH and ionic strength conditions. Furthermore, we analyze quantitatively the mechanism of allosteric communication between the topologically distinct c-Myb and MLL sites. The implications of our results are discussed in the light of previous work on other intrinsically unstructured systems.

The effects of threonine phosphorylation on the stability and dynamics of the central molecular switch region of 18.5-kDa myelin basic protein

The classic isoforms of myelin basic protein (MBP) are essential for the formation and maintenance of myelin in the central nervous system of higher vertebrates. The protein is involved in all facets of the development, compaction, and stabilization of the multilamellar myelin sheath, and also interacts with cytoskeletal and signaling proteins. The predominant 18.5-kDa isoform of MBP is an intrinsically-disordered protein that is a candidate auto-antigen in the human demyelinating disease multiple sclerosis. A highly-conserved central segment within classic MBP consists of a proline-rich region (murine 18.5-kDa sequence -T92-P93-R94-T95-P96-P97-P98-S99-) containing a putative SH3-ligand, adjacent to a region that forms an amphipathic α-helix (P82-I90) upon interaction with membranes, or under membrane-mimetic conditions. The T92 and T95 residues within the proline-rich region can be post-translationally modified through phosphorylation by mitogen-activated protein (MAP) kinases. Here, we have investigated the structure of the α-helical and proline-rich regions in dilute aqueous buffer, and have evaluated the effects of phosphorylation at T92 and T95 on the stability and dynamics of the α-helical region, by utilizing four 36-residue peptides (S72-S107) with differing phosphorylation status. Nuclear magnetic resonance spectroscopy reveals that both the α-helical as well as the proline-rich regions are disordered in aqueous buffer, whereas they are both structured in a lipid environment (cf., Ahmed et al., Biochemistry 51, 7475-9487, 2012). Thermodynamic analysis of trifluoroethanol-titration curves monitored by circular dichroism spectroscopy reveals that phosphorylation, especially at residue T92, impedes formation of the amphipathic α-helix. This conclusion is supported by molecular dynamics simulations, which further illustrate that phosphorylation reduces the folding reversibility of the α-helix upon temperature perturbation and affect the global structure of the peptides through altered electrostatic interactions. The results support the hypothesis that the central conserved segment of MBP constitutes a molecular switch in which the conformation and/or intermolecular interactions are mediated by phosphorylation/dephosphorylation at T92 and T95.

A folding-after-binding mechanism describes the recognition between the transactivation domain of c-Myb and the KIX domain of the CREB-binding protein

A large body of evidence suggests that a considerable fraction of the human proteome may be at least in part intrinsically unstructured. While disordered, intrinsically unstructured proteins are nevertheless functional and mediate many interactions. Despite their significant role in regulation, however, little is known about the molecular mechanism whereby intrinsically unstructured proteins exert their function. This basic problem is critical to establish the role, if any, of disorder in cellular systems. Here we present kinetic experiments supporting a mechanism of binding-induced-folding when the KIX domain of the CREB-binding protein binds the transactivation domain of c-Myb, an intrinsically unstructured domain. The high-resolution structure of this physiologically important complex was previously determined by NMR spectroscopy. Our data reveal that c-Myb recognizes KIX by first forming a weak encounter complex in a disordered conformation, which is subsequently locked-in by a folding step, i.e. binding precedes folding. On the basis of the pH dependence of the observed combination and dissociation rate constants we propose a plausible mechanism for complex formation. The implications of our results in the light of previous work on intrinsically unstructured systems are discussed.

Thermodynamic and kinetic analysis of peptides derived from CapZ, NDR, p53, HDM2, and HDM4 binding to human S100B

S100B is a member of the S100 subfamily of EF-hand proteins that has been implicated in malignant melanoma and neurodegenerative conditions such as Alzheimer's disease and Parkinson's disease. Calcium-induced conformational changes expose a hydrophobic binding cleft, facilitating interactions with a wide variety of nuclear, cytoplasmic, and extracellular target proteins. Previously, peptides derived from CapZ, p53, NDR, HDM2, and HDM4 have been shown to interact with S100B in a calcium-dependent manner. However, the thermodynamic and kinetic basis of these interactions remains largely unknown. To gain further insight, we screened these peptides against the S100B protein using isothermal titration calorimetry and nuclear magnetic resonance. All peptides were found to have binding affinities in the low micromolar to nanomolar range. Binding-induced changes in the line shapes of S100B backbone (1)H and (15)N resonances were monitored to obtain the dissociation constants and the kinetic binding parameters. The large microscopic K(on) rate constants observed in this study (≥1 × 10(7) M(-1) s(-1)) suggest that S100B utilizes a "fly casting mechanism" in the recognition of these peptide targets.

Laser temperature-jump spectroscopy of intrinsically disordered proteins

Laser temperature-jump methods allow an experimenter to study the kinetics and dynamics of very rapid solution-phase processes, including conformational dynamics of biomolecules on time scales of nanoseconds and microseconds. The combination of laser temperature-jump (T-jump) excitation and appropriate optical detection techniques such as fluorescence energy transfer allows the study of intramolecular and intermolecular conformational changes and interactions that occur during protein folding and binding. This article describes the application of the laser temperature-jump method to UV-visible fluorescence studies of the coupled folding and binding of intrinsically disordered proteins. We emphasize the practical aspects of instrument alignment and optimization, sample preparation, and data collection using fluorescently labeled peptides with UV laser excitation.

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.

Histone chaperones facilitate the assembly and disassembly of nucleosome by interacting with the corresponding histones or histone variants. As the biomolecules in nucleosome are highly charged, electrostatic interactions are particularly important in these processes. The experiments have explored that the histone chaperon Chz1 as an intrinsically disordered protein (IDP) can fold by binding to its histone variants H2A.Z-H2B. Here, we developed a molecular simulation program that treated electrostatic interactions with Debye-Hückel model to study the mechanism of the association. We found that the inter-chain electrostatic interactions facilitate the coupled folding and binding transitions, consistent with the kinetic experiments and microscopic structural perspectives. Furthermore, we show that the intra-chain electrostatic interactions collapse Chz1 and slow the binding rate. The collapsed structure in IDPs caused by intra-chain electrostatic interactions has been widely found in experiments and the effect in binding is well studied in our work. Our theoretical approach shed new light on the role of electrostatics on inter-chain and intra-chain interactions for IDPs' binding and is applicable to the binding/folding of other IDPs to their targets.

Rate constants and mechanisms of intrinsically disordered proteins binding to structured targets

The binding of intrinsically disordered proteins (IDPs) to structured targets is gaining increasing attention. Here we review experimental and computational studies on the binding kinetics of IDPs. Experiments have yielded both the binding rate constants and the binding mechanisms, the latter via mutation and deletion studies and NMR techniques. Most computational studies have aimed at qualitative understanding of the binding rate constants or at mapping the free energy surfaces after the IDPs are engaged with their targets. The experiments and computation show that IDPs generally gain structures after they are engaged with their targets; that is, interactions with the targets facilitate the IDPs' folding. It also seems clear that the initial contact of an IDP with the target is formed by just a segment, not the entire IDP. The docking of one segment to its sub-site followed by coalescing of other segments around the corresponding sub-sites emerges as a recurring feature in the binding of IDPs. Such a dock-and-coalesce model forms the basis for quantitative calculation of binding rate constants. For both disordered and ordered proteins, strong electrostatic attraction with their targets can enhance the binding rate constants by several orders of magnitude. There are now tremendous opportunities in narrowing the gap in our understanding of IDPs relative to ordered proteins with regard to binding kinetics.

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.

Towards the physical basis of how intrinsic disorder mediates protein function

Intrinsically disordered proteins (IDPs) are an important class of functional proteins that is highly prevalent in biology and has broad association with human diseases. In contrast to structured proteins, free IDPs exist as heterogeneous and dynamical conformational ensembles under physiological conditions. Many concepts have been discussed on how such intrinsic disorder may provide crucial functional advantages, particularly in cellular signaling and regulation. Establishing the physical basis of these proposed phenomena requires not only detailed characterization of the disordered conformational ensembles, but also mechanistic understanding of the roles of various ensemble properties in IDP interaction and regulation. Here, we review the experimental and computational approaches that may be integrated to address many important challenges of establishing a "structural" basis of IDP function, and discuss some of the key emerging ideas on how the conformational ensembles of IDPs may mediate function, especially in coupled binding and folding interactions.

The transition state of coupled folding and binding for a flexible β-finger

Flexible and fully disordered protein regions that fold upon binding mediate numerous protein-protein interactions. However, little is known about their mechanism of interaction. One such coupled folding and binding occurs when a flexible region of neuronal nitric oxide synthase adopts a β-finger structure upon binding to its protein ligand, a PDZ [PSD-95 (postsynaptic density protein-95)/Discs large/ZO-1] domain from PSD-95. We have analyzed this binding reaction by protein engineering combined with kinetic experiments. Mutational destabilization of the β-finger changed mainly the dissociation rate constant of the proteins and, to a lesser extent, the association rate constant. Thus, mutation affected late events in the coupled folding and binding reaction. Our results therefore suggest that the native binding interactions of the β-finger are not present in the rate-limiting transition state for binding but form on the downhill side in a cooperative manner. However, by mutation, we could destabilize the β-finger further and change the rate-limiting step such that an initial conformational change becomes rate limiting. This switch in rate-limiting step shows that multistep binding mechanisms are likely to be found among flexible and intrinsically disordered regions of proteins.

Residual structures, conformational fluctuations, and electrostatic interactions in the synergistic folding of two intrinsically disordered proteins

To understand the interplay of residual structures and conformational fluctuations in the interaction of intrinsically disordered proteins (IDPs), we first combined implicit solvent and replica exchange sampling to calculate atomistic disordered ensembles of the nuclear co-activator binding domain (NCBD) of transcription coactivator CBP and the activation domain of the p160 steroid receptor coactivator ACTR. The calculated ensembles are in quantitative agreement with NMR-derived residue helicity and recapitulate the experimental observation that, while free ACTR largely lacks residual secondary structures, free NCBD is a molten globule with a helical content similar to that in the folded complex. Detailed conformational analysis reveals that free NCBD has an inherent ability to substantially sample all the helix configurations that have been previously observed either unbound or in complexes. Intriguingly, further high-temperature unbinding and unfolding simulations in implicit and explicit solvents emphasize the importance of conformational fluctuations in synergistic folding of NCBD with ACTR. A balance between preformed elements and conformational fluctuations appears necessary to allow NCBD to interact with different targets and fold into alternative conformations. Together with previous topology-based modeling and existing experimental data, the current simulations strongly support an “extended conformational selection” synergistic folding mechanism that involves a key intermediate state stabilized by interaction between the C-terminal helices of NCBD and ACTR. In addition, the atomistic simulations reveal the role of long-range as well as short-range electrostatic interactions in cooperating with readily fluctuating residual structures, which might enhance the encounter rate and promote efficient folding upon encounter for facile binding and folding interactions of IDPs. Thus, the current study not only provides a consistent mechanistic understanding of the NCBD/ACTR interaction, but also helps establish a multi-scale molecular modeling framework for understanding the structure, interaction, and regulation of IDPs in general.

Intrinsically disordered proteins (IDPs) are now widely recognized to play fundamental roles in biology and to be frequently associated with human diseases. Although the potential advantages of intrinsic disorder in cellular signaling and regulation have been widely discussed, the physical basis for these proposed phenomena remains sketchy at best. An integration of multi-scale molecular modeling and experimental characterization is necessary to uncover the molecular principles that govern the structure, interaction, and regulation of IDPs. In this work, we characterize the conformational properties of two IDPs involved in transcription regulation at the atomistic level and further examine the roles of these properties in their coupled binding and folding interactions. Our simulations suggest interplay among residual structures, conformational fluctuations, and electrostatic interactions that allows efficient synergistic folding of these two IDPs. In particular, we propose that electrostatic interactions might play an important role in facilitating rapid folding and binding recognition of IDPs, by enhancing the encounter rate and promoting efficient folding upon encounter.

Side-chain interactions form late and cooperatively in the binding reaction between disordered peptides and PDZ domains

Intrinsically disordered proteins are very common and mediate numerous protein-protein and protein-DNA interactions. While it is clear that these interactions are instrumental for the life of the mammalian cell, there is a paucity of data regarding their molecular binding mechanisms. Here we have used short peptides as a model system for intrinsically disordered proteins. Linear free energy relationships based on rate and equilibrium constants for the binding of these peptides to ordered target proteins, PDZ domains, demonstrate that native side-chain interactions form mainly after the rate-limiting barrier for binding and in a cooperative fashion. This finding suggests that these disordered peptides first form a weak encounter complex with non-native interactions. The data do not support the recent notion that the affinities of intrinsically disordered proteins toward their targets are generally governed by their association rate constants. Instead, we observed the opposite for peptide-PDZ interactions, namely, that changes in K(d) correlate with changes in k(off).

Dynamics and mechanisms of coupled protein folding and binding reactions

Protein folding coupled to binding of a specific ligand is frequently observed in biological processes. In recent years numerous studies have addressed the structural properties of the unfolded proteins in the absence of their ligands. Surprisingly few time-resolved investigations on coupled folding and binding reactions have been published up to date and the dynamics and kinetic mechanisms of these processes are still only poorly understood. Especially, it is still unsolved for most systems which conformation of the protein is recognized by the ligand (conformational selection vs. folding-after-binding) and whether the ligand influences the folding kinetics. Here we review experimental methods, kinetic models and time-resolved experimental studies of coupled folding and binding reactions.

Modulation of an IDP binding mechanism and rates by helix propensity and non-native interactions: association of HIF1α with CBP

Intrinsically disordered proteins that acquire their three dimensional structures only upon binding to their targets are very important in cellular signal regulation. While experimental studies have been made on the structures of both bound (structured) and unbound (disordered) states, less is known about the actual folding-binding transition. Coarse grained simulations using native-centric (i.e. Gō) potentials have been particularly useful in addressing this problem, given the large search space for IDP binding, but have well-known deficiencies in reproducing the unfolded state structure and dynamics. Here, we investigate the interaction of HIF1α with CBP using a hierarchy of coarse-grained models, in each case matching the binding affinity at 300 K to the experimental value. Starting from a pure Gō-like model based on the native structure of the complex we go on to consider a more realistic model of helix propensity in the HIF1α, and finally the effect of non-native interactions between binding partners. We find structural disorder (i.e."fuzziness") in the bound state of HIF1α in all models which is supported by the results of atomistic simulations. Correcting the over-stabilized helices in the unbound state gives rise to a more cooperative folding-binding transition (destabilizing partially bound intermediates). Adding non-native contacts lowers the free energy barrier for binding to an almost barrierless scenario, leading to higher binding/unbinding rates relative to the other models, in better agreement with the near diffusion-limited binding rates measured experimentally. Transition state structures for the three models are highly disordered, supporting a fly-casting mechanism for binding.

Synergistic folding of two intrinsically disordered proteins: searching for conformational selection

Intrinsically disordered proteins (IDPs) lack stable structures under physiological conditions but often fold into stable structures upon specific binding. These coupled binding and folding processes underlie the organization of cellular regulatory networks, and a mechanistic understanding is thus of fundamental importance. Here, we investigated the synergistic folding of two IDPs, namely, the NCBD domain of transcription coactivator CBP and the p160 steroid receptor coactivator ACTR, using a topology-based model that was carefully calibrated to balance intrinsic folding propensities and intermolecular interactions. As one of the most structured IDPs, NCBD is a plausible candidate that interacts through conformational selection-like mechanisms, where binding is mainly initiated by pre-existing folded-like conformations. Indeed, the simulations demonstrate that, even though binding and folding of both NCBD and ACTR is highly cooperative on the baseline level, the tertiary folding of NCBD is best described by the "extended conformational selection" model that involves multiple stages of selection and induced folding. The simulations further predict that the NCBD/ACTR recognition is mainly initiated by forming a mini folded core that includes the second and third helices of NCBD and ACTR. These predictions are fully consistent with independent physics-based atomistic simulations as well as a recent experimental mapping of the H/D exchange protection factors. The current work thus adds to the limited number of existing mechanistic studies of coupled binding and folding of IDPs, and provides a first direct demonstration of how conformational selection might contribute to efficient recognition of IDPs. Interestingly, even for highly structured IDPs like NCBD, the recognition is initiated by the more disordered C-terminal segment and with substantial contribution from induced folding. Together with existing studies of IDP interaction mechanisms, this argues that induced folding is likely prevalent in IDP-protein interaction, and emphasizes the importance of understanding how IDPs manage to fold efficiently upon (nonspecific) binding. Success of the current study also further supports the notion that, with careful calibration, topology-based models can be effective tools for mechanistic study of IDP interaction and regulation, especially when combined with physics-based atomistic simulations and experiments.

A free-energy landscape for coupled folding and binding of an intrinsically disordered protein in explicit solvent from detailed all-atom computations

The N-terminal repressor domain of neural restrictive silencer factor (NRSF) is an intrinsically disordered protein (IDP) that binds to the paired amphipathic helix (PAH) domain of mSin3. An NMR experiment revealed that the minimal binding unit of NRSF is a 15-residue segment that adopts a helical structure upon binding to a cleft of mSin3. We computed a free-energy landscape of this system by an enhanced conformational sampling method, all-atom multicanonical molecular dynamics. The simulation started from a configuration where the NRSF segment was fully disordered and distant from mSin3 in explicit solvent. In the absence of mSin3, the disordered NRSF segment thermally fluctuated between hairpins, helices, and bent structures. In the presence of mSin3, the segment bound to mSin3 by adopting the structures involved in the isolated state, and non-native and native complexes were formed. The free-energy landscape comprised three superclusters, and free-energy barriers separated the superclusters. The native complex was located at the center of the lowest free-energy cluster. When NRSF landed in the largest supercluster, the generated non-native complex moved on the landscape to fold into the native complex, by increasing the interfacial hydrophobic contacts and the helix content. When NRSF landed in other superclusters, the non-native complex overcame the free-energy barriers between the various segment orientations in the binding cleft of mSin3. Both population-shift and induced-fit (or induced-folding) mechanisms work cooperatively in the coupled folding and binding. The diverse structural adaptability of NRSF may be related to the hub properties of the IDP.

Characterization of the disordered-to-α-helical transition of IA₃ by SDSL-EPR spectroscopy

Electron paramagnetic resonance (EPR) spectroscopy coupled with site-directed spin labeling (SDSL) is a valuable tool for characterizing the mobility and conformational changes of proteins but has seldom been applied to intrinsically disordered proteins (IDPs). Here, IA₃ is used as a model system demonstrating SDSL-EPR characterization of conformational changes in small IDP systems. IA₃ has 68 amino acids, is unstructured in solution, and becomes α-helical upon addition of the secondary structural stabilizer 2,2,2-trifluoroethanol (TFE). Two single cysteine substitutions, one in the N-terminus (S14C) and one in the C-terminus (N58C), were generated and labeled with three different nitroxide spin labels. The resultant EPR line shapes of each of the labels were compared and each reported changes in mobility upon addition of TFE. Specifically, the spectral line shape parameters h((+1))/h(₀), the local tumbling volume (V(L)), and the percent change of the h(₋₁) intensity were utilized to quantitatively monitor TFE-induced conformational changes. The values of h((+1)/)h(₀) as a function of TFE titration varied in a sigmoidal manner and were fit to a two-state Boltzmann model that provided values for the midpoint of the transition, thus, reporting on the global conformational change of IA₃. The other parameters provide site-specific information and show that S14C-SL undergoes a conformational change resulting in more restricted motion than N58C-SL, which is consistent with previously published results obtained by studies using NMR and circular dichroism spectroscopy indicating a higher degree of α-helical propensity of the N-terminal segment of IA₃. Overall, the results provide a framework for data analyzes that can be used to study induced unstructured-to-helical conformations in IDPs by SDSL.

Multi-scaled explorations of binding-induced folding of intrinsically disordered protein inhibitor IA3 to its target enzyme

Biomolecular function is realized by recognition, and increasing evidence shows that recognition is determined not only by structure but also by flexibility and dynamics. We explored a biomolecular recognition process that involves a major conformational change – protein folding. In particular, we explore the binding-induced folding of IA3, an intrinsically disordered protein that blocks the active site cleft of the yeast aspartic proteinase saccharopepsin (YPrA) by folding its own N-terminal residues into an amphipathic alpha helix. We developed a multi-scaled approach that explores the underlying mechanism by combining structure-based molecular dynamics simulations at the residue level with a stochastic path method at the atomic level. Both the free energy profile and the associated kinetic paths reveal a common scheme whereby IA3 binds to its target enzyme prior to folding itself into a helix. This theoretical result is consistent with recent time-resolved experiments. Furthermore, exploration of the detailed trajectories reveals the important roles of non-native interactions in the initial binding that occurs prior to IA3 folding. In contrast to the common view that non-native interactions contribute only to the roughness of landscapes and impede binding, the non-native interactions here facilitate binding by reducing significantly the entropic search space in the landscape. The information gained from multi-scaled simulations of the folding of this intrinsically disordered protein in the presence of its binding target may prove useful in the design of novel inhibitors of aspartic proteinases.

The intrinsically disordered peptide IA3 is the endogenous inhibitor for the enzyme named yeast aspartic proteinase saccharopepsin (YPrA). In the presence of YPrA, IA3 folds itself into an amphipathic helix that blocks the active site cleft of the enzyme. We developed a multi-scaled approach to explore the underlying mechanism of this binding-induced ordering transition. Our approach combines a structure-based molecular dynamics model at the residue level with a stochastic path method at the atomic level. Our simulations suggest that IA3 inhibits YPrA through an induced-fit mechanism where the enzyme (YPrA) induces conformational change of its inhibitor (IA3). This expands the definition of an induced-fit model from its original meaning that the binding of substrate (IA3) drives conformational change in the protein (YPrA). Our result is consistent with recent kinetic experiments and provides a microscopic explanation for the underlying mechanism. We also discuss the important roles of non-native interactions and backtracking. These results enrich our understanding of the enzyme-inhibition mechanism and may have value in the design of drugs.

Intrinsically disordered proteins in a physics-based world

Intrinsically disordered proteins (IDPs) are a newly recognized class of functional proteins that rely on a lack of stable structure for function. They are highly prevalent in biology, play fundamental roles, and are extensively involved in human diseases. For signaling and regulation, IDPs often fold into stable structures upon binding to specific targets. The mechanisms of these coupled binding and folding processes are of significant importance because they underlie the organization of regulatory networks that dictate various aspects of cellular decision-making. This review first discusses the challenge in detailed experimental characterization of these heterogeneous and dynamics proteins and the unique and exciting opportunity for physics-based modeling to make crucial contributions, and then summarizes key lessons from recent de novo simulations of the structure and interactions of several regulatory IDPs.

The basic helix-loop-helix region of human neurogenin 1 is a monomeric natively unfolded protein which forms a "fuzzy" complex upon DNA binding

Neuronal specification is regulated by the activity of transcription factors containing the basic helix-loop-helix motif (bHLH); these regulating proteins include, among others, the neurogenin (Ngn) family, related to the atonal family of genes. Neurogenin 1 (NGN1) is a 237-residue protein that contains a bHLH domain and is involved in neuronal differentiation. In this work, we synthesized the bHLH region of NGN1 (bHLHN) comprising residues 90-150 of the full-length NGN1. The domain is a monomeric natively unfolded protein with a pH-dependent premolten globule conformation, as shown by several spectroscopic techniques (namely, NMR, fluorescence, FTIR, and circular dichroism). The unfolded character of the domain also explains, first, the impossibility of its overexpression in several Escherichia coli strains and, second, its insolubility in aqueous buffers. To the best of our knowledge, this is the first extensive study of the conformational preferences of a bHLH domain under different solution conditions. Upon binding to two DNA E-boxes, the protein forms "fuzzy" complexes (that is, the complexes were not fully folded). The affinities of bHLHN for both DNA boxes were smaller than those of other bHLH domains, which might explain why the protein-DNA complexes were not fully folded.

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

Intrinsically disordered proteins (IDPs) are recognized to play important roles in many biological functions such as transcription and translation regulation, cellular signal transduction, protein phosphorylation, and molecular assemblies. The coupling of folding with binding through a "fly-casting" mechanism has been proposed to account for the fast binding kinetics of IDPs. In this article, experimental data from the literature were collated to verify the kinetic advantages of IDPs, while molecular simulations were performed to clarify the origin of the kinetic advantages. The phosphorylated KID-kinase-inducible domain interacting domain (KIX) complex was used as an example in the simulations. By modifying a coarse-grained model with a native-centric Gō-like potential, we were able to continuously tune the degree of disorder of the phosphorylated KID domain and thus investigate the intrinsic role of chain flexibility in binding kinetics. The simulations show that the "fly-casting" effect is not only due to the greater capture radii of IDPs. The coupling of folding with binding of IDPs leads to a significant reduction in binding free-energy barrier. Such a reduction accelerates the binding process. Although the greater capture radius has been regarded as the main factor in promoting the binding rate of IDPs, we found that this parameter will also lead to the slower translational diffusion of IDPs when compared with ordered proteins. As a result, the capture rate of IDPs was found to be slower than that of ordered proteins. The main origin of the faster binding for IDPs are the fewer encounter times required before the formation of the final binding complex. The roles of the interchain native contacts fraction (Q(b)) and the mass-center distance (DeltaR) as reaction coordinates are also discussed.

Conformational selection or induced fit: a flux description of reaction mechanism

The mechanism of ligand binding coupled to conformational changes in macromolecules has recently attracted considerable interest. The 2 limiting cases are the "induced fit" mechanism (binding first) or "conformational selection" (conformational change first). Described here are the criteria by which the sequence of events can be determined quantitatively. The relative importance of the 2 pathways is determined not by comparing rate constants (a common misconception) but instead by comparing the flux through each pathway. The simple rules for calculating flux in multistep mechanisms are described and then applied to 2 examples from the literature, neither of which has previously been analyzed using the concept of flux. The first example is the mechanism of conformational change in the binding of NADPH to dihydrofolate reductase. The second example is the mechanism of flavodoxin folding coupled to binding of its cofactor, flavin mononucleotide. In both cases, the mechanism switches from being dominated by the conformational selection pathway at low ligand concentration to induced fit at high ligand concentration. Over a wide range of conditions, a significant fraction of the flux occurs through both pathways. Such a mixed mechanism likely will be discovered for many cases of coupled conformational change and ligand binding when kinetic data are analyzed by using a flux-based approach.

Solvent friction changes the folding pathway of the tryptophan zipper TZ2

Because the rate of a diffusional process such as protein folding is controlled by friction encountered along the reaction pathway, the speed of folding is readily tunable through adjustment of solvent viscosity. The precise relationship between solvent viscosity and the rate of diffusion is complex and even conformation-dependent, however, because both solvent friction and protein internal friction contribute to the total reaction friction. The heterogeneity of the reaction friction along the folding pathway may have subtle consequences. For proteins that fold on a multidimensional free-energy surface, an increase in solvent friction may drive a qualitative change in folding trajectory. Our time-resolved experiments on the rapidly and heterogeneously folding beta-hairpin TZ2 show a shift in the folding pathway as viscosity increases, even though the energetics of folding is unaltered. We also observe a nonlinear or saturating behavior of the folding relaxation time with rising solvent viscosity, potentially an experimental signature of the shifting pathway for unfolding. Our results show that manipulations of solvent viscosity in folding experiments and simulations may have subtle and unexpected consequences on the folding dynamics being studied.