The dynamic properties of a nuclear coactivator binding domain are evolutionarily conserved.

Evolution of proteins is constrained by their structure and function. While there is a consensus that the plasticity of intrinsically disordered proteins relaxes the structural constraints on evolution there is a paucity of data on the molecular details of these processes. The Nuclear Coactivator Binding Domain (NCBD) from CREB-binding protein is a protein interaction domain, which contains a hydrophobic core but is not behaving as a typical globular domain, and has been described as 'molten-globule like'. The highly dynamic properties of NCBD makes it an interesting model system for evolutionary structure-function investigation of intrinsically disordered proteins. We have here compared the structure and biophysical properties of an ancient version of NCBD present in a bilaterian animal ancestor living around 600 million years ago with extant human NCBD. Using a combination of NMR spectroscopy, circular dichroism and kinetics we show that although NCBD has increased its thermodynamic stability, it has retained its dynamic biophysical properties in the ligand-free state in the evolutionary lineage leading from the last common bilaterian ancestor to humans. Our findings suggest that the dynamic properties of NCBD have been maintained by purifying selection and thus are important for its function, which includes mediating several distinct protein-protein interactions.

Characterization of the dynamics and the conformational entropy in the binding between TAZ1 and CTAD-HIF-1α.

The interaction between the C-terminal transactivation domain of HIF-1α (CTAD-HIF-1α) and the transcriptional adapter zinc binding 1 (TAZ1) domain of CREB binding protein participate in the initiation of gene transcription during hypoxia. Unbound CTAD-HIF-1α is disordered but undergoes a disorder-to-order transition upon binding to TAZ1. We have here performed NMR side chain and backbone relaxation studies on TAZ1 and side chain relaxation measurements on CTAD-HIF-1α in order to investigate the role of picosecond to nanosecond dynamics. We find that the internal motions are significantly affected upon binding, both on the side chain and the backbone level. The dynamic response corresponds to a conformational entropy change that contributes substantially to the binding thermodynamics for both binding partners. Furthermore, the conformational entropy change for the well-folded TAZ1 varies upon binding to different IDP targets. We further identify a cluster consisting of side chains in bound TAZ1 and CTAD-HIF-1α that experience extensive dynamics and are part of the binding region that involves the N-terminal end of the LPQL motif in CTAD-HIF-1α; a feature that might have an important role in the termination of the hypoxic response.

Role of Conformational Entropy in Molecular Recognition by TAZ1 of CBP.

The globular transcriptional adapter zinc binding 1 (TAZ1) domain of CREB binding protein participates in protein-protein interactions that are involved in transcriptional regulation. TAZ1 binds numerous targets, of which many are intrinsically disordered proteins that undergo a disorder-to-order transition to various degrees. One such target is the disordered transactivation domain of transcription factor RelA (TAD-RelA), which with its interaction with TAZ1 is involved in transcriptional regulation of genes in NF-κB signaling. We have here performed nuclear magnetic resonance backbone and side-chain relaxation studies to investigate the influence of RelA-TA2 (residues 425-490 in TAD-RelA) binding on the subnanosecond internal motions of TAZ1. We find a considerable dynamic response on both the backbone and side-chain levels, which corresponds to a conformational entropy change that contributes significantly to the binding energetics. We further show that the microscopic origins of the dynamic response of TAZ1 vary depending on the target. This study demonstrates that folded protein domains that are able to interact with various targets are not dynamically passive but can have a significant role in the motional response upon target association.

A structurally heterogeneous transition state underlies coupled binding and folding of disordered proteins.

Many intrinsically disordered proteins (IDPs) attain a well-defined structure in a coupled folding and binding reaction with another protein. Such reactions may involve early to late formation of different native structural regions along the reaction pathway. To obtain insights into the transition state for a coupled binding and folding reaction, we performed restrained molecular dynamics simulations using previously determined experimental binding Φb values of the interaction between two IDP domains: the activation domain from the p160 transcriptional co-activator for thyroid hormone and retinoid receptors (ACTR) and the nuclear co-activator binding domain (NCBD) of CREB-binding protein, each forming three well-defined α-helices upon binding. These simulations revealed that both proteins are largely disordered in the transition state for complex formation, except for two helices, one from each domain, that display a native-like structure. The overall transition state structure was extended and largely dynamic with many weakly populated contacts. To test the transition state model, we combined site-directed mutagenesis with kinetic experiments, yielding results consistent with overall diffuse interactions and formation of native intramolecular interactions in the third NCBD helix during the binding reaction. Our findings support the view that the transition state and, by inference, any encounter complex in coupled binding and folding reactions are structurally heterogeneous and largely independent of specific interactions. Furthermore, experimental Φb values and Brønsted plots suggested that the transition state is globally robust with respect to most mutations but can display more native-like features for some highly destabilizing mutations, possibly because of Hammond behavior or ground-state effects.

Structure and dynamics conspire in the evolution of affinity between intrinsically disordered proteins.

In every established species, protein-protein interactions have evolved such that they are fit for purpose. However, the molecular details of the evolution of new protein-protein interactions are poorly understood. We have used nuclear magnetic resonance spectroscopy to investigate the changes in structure and dynamics during the evolution of a protein-protein interaction involving the intrinsically disordered CREBBP (CREB-binding protein) interaction domain (CID) and nuclear coactivator binding domain (NCBD) from the transcriptional coregulators NCOA (nuclear receptor coactivator) and CREBBP/p300, respectively. The most ancient low-affinity "Cambrian-like" [540 to 600 million years (Ma) ago] CID/NCBD complex contained less secondary structure and was more dynamic than the complexes from an evolutionarily younger "Ordovician-Silurian" fish ancestor (ca. 440 Ma ago) and extant human. The most ancient Cambrian-like CID/NCBD complex lacked one helix and several interdomain interactions, resulting in a larger solvent-accessible surface area. Furthermore, the most ancient complex had a high degree of millisecond-to-microsecond dynamics distributed along the entire sequences of both CID and NCBD. These motions were reduced in the Ordovician-Silurian CID/NCBD complex and further redistributed in the extant human CID/NCBD complex. Isothermal calorimetry experiments show that complex formation is enthalpically favorable and that affinity is modulated by a largely unfavorable entropic contribution to binding. Our data demonstrate how changes in structure and motion conspire to shape affinity during the evolution of a protein-protein complex and provide direct evidence for the role of structural, dynamic, and frustrational plasticity in the evolution of interactions between intrinsically disordered proteins.

The transition state structure for binding between TAZ1 of CBP and the disordered Hif-1α CAD.

Intrinsically disordered proteins (IDPs) are common in eukaryotes. However, relatively few experimental studies have addressed the nature of the rate-limiting transition state for the coupled binding and folding reactions involving IDPs. By using site-directed mutagenesis in combination with kinetics measurements we have here characterized the transition state for binding between the globular TAZ1 domain of CREB binding protein and the intrinsically disordered C-terminal activation domain of Hif-1α (Hif-1α CAD). A total of 17 Hif-1α CAD point-mutations were generated and a Φ-value binding analysis was carried out. We found that native hydrophobic binding interactions are not formed at the transition state. We also investigated the effect the biologically important Hif-1α CAD Asn-803 hydroxylation has on the binding kinetics, and found that the whole destabilization effect due the hydroxylation is within the dissociation rate constant. Thus, the rate-limiting transition state is "disordered-like", with native hydrophobic binding contacts being formed cooperatively after the rate-limiting barrier, which is clearly shown by linear free energy relationships. The same behavior was observed in a previously characterized TAZ1/IDP interaction, which may suggest common features for the rate-limiting transition state for TAZ1/IDP interactions.

Dynamics, Conformational Entropy, and Frustration in Protein-Protein Interactions Involving an Intrinsically Disordered Protein Domain.

Intrinsically disordered proteins (IDPs) are abundant in the eukaryotic proteome. However, little is known about the role of subnanosecond dynamics and the conformational entropy that it represents in protein-protein interactions involving IDPs. Using nuclear magnetic resonance side chain and backbone relaxation, stopped-flow kinetics, isothermal titration calorimetry, and computational studies, we have characterized the interaction between the globular TAZ1 domain of the CREB binding protein and the intrinsically disordered transactivation domain of STAT2 (TAD-STAT2). We show that the TAZ1/TAD-STAT2 complex retains considerable subnanosecond motions, with TAD-STAT2 undergoing only a partial disorder-to-order transition. We report here the first experimental determination of the conformational entropy change for both binding partners in an IDP binding interaction and find that the total change even exceeds in magnitude the binding enthalpy and is comparable to the contribution from the hydrophobic effect, demonstrating its importance in the binding energetics. Furthermore, we show that the conformational entropy change for TAZ1 is also instrumental in maintaining a biologically meaningful binding affinity. Strikingly, a spatial clustering of very high amplitude motions and a cluster of more rigid sites in the complex exist, which through computational studies we found to overlap with regions that experience energetic frustration and are less frustrated, respectively. Thus, the residual dynamics in the bound state could be necessary for faster dissociation, which is important for proteins that interact with multiple binding partners.

Native Hydrophobic Binding Interactions at the Transition State for Association between the TAZ1 Domain of CBP and the Disordered TAD-STAT2 Are Not a Requirement.

A significant fraction of the eukaryotic proteome consists of proteins that are either partially or completely disordered under native-like conditions. Intrinsically disordered proteins (IDPs) are common in protein-protein interactions and are involved in numerous cellular processes. Although many proteins have been identified as disordered, much less is known about the binding mechanisms of the coupled binding and folding reactions involving IDPs. Here we have analyzed the rate-limiting transition state for binding between the TAZ1 domain of CREB binding protein and the intrinsically disordered transactivation domain of STAT2 (TAD-STAT2) by site-directed mutagenesis and kinetic experiments (Φ-value analysis) and found that the native protein-protein binding interface is not formed at the transition state for binding. Instead, native hydrophobic binding interactions form late, after the rate-limiting barrier has been crossed. The association rate constant in the absence of electrostatic enhancement was determined to be rather high. This is consistent with the Φ-value analysis, which showed that there are few or no obligatory native contacts. Also, linear free energy relationships clearly demonstrate that native interactions are cooperatively formed, a scenario that has usually been observed for proteins that fold according to the so-called nucleation-condensation mechanism. Thus, native hydrophobic binding interactions at the rate-limiting transition state for association between TAD-STAT2 and TAZ1 are not a requirement, which is generally in agreement with previous findings on other IDP systems and might be a common mechanism for IDPs.

Emergence and evolution of an interaction between intrinsically disordered proteins.

Protein-protein interactions involving intrinsically disordered proteins are important for cellular function and common in all organisms. However, it is not clear how such interactions emerge and evolve on a molecular level. We performed phylogenetic reconstruction, resurrection and biophysical characterization of two interacting disordered protein domains, CID and NCBD. CID appeared after the divergence of protostomes and deuterostomes 450-600 million years ago, while NCBD was present in the protostome/deuterostome ancestor. The most ancient CID/NCBD formed a relatively weak complex (Kd∼5 µM). At the time of the first vertebrate-specific whole genome duplication, the affinity had increased (Kd∼200 nM) and was maintained in further speciation. Experiments together with molecular modeling using NMR chemical shifts suggest that new interactions involving intrinsically disordered proteins may evolve via a low-affinity complex which is optimized by modulating direct interactions as well as dynamics, while tolerating several potentially disruptive mutations.

Activation Barrier-Limited Folding and Conformational Sampling of a Dynamic Protein Domain.

Folding reaction mechanisms of globular protein domains have been extensively studied by both experiment and simulation and found to be highly concerted chemical reactions in which numerous noncovalent bonds form in an apparent two-state fashion. However, less is known regarding intrinsically disordered proteins because their folding can usually be studied only in conjunction with binding to a ligand. We have investigated by kinetics the folding mechanism of such a disordered protein domain, the nuclear coactivator-binding domain (NCBD) from CREB-binding protein. While a previous computational study suggested that NCBD folds without an activation free energy barrier, our experimental data demonstrate that NCBD, despite its highly dynamic structure, displays relatively slow folding (∼10 ms at 277 K) consistent with a barrier-limited process. Furthermore, the folding kinetics corroborate previous nuclear magnetic resonance data showing that NCBD exists in two folded conformations and one more denatured conformation at equilibrium and, thus, that the folding mechanism is a three-state mechanism. The refolding kinetics is limited by unfolding of the less populated folded conformation, suggesting that the major route for interconversion between the two folded states is via the denatured state. Because the two folded conformations have been suggested to bind distinct ligands, our results have mechanistic implications for conformational sampling in protein-protein interactions.

Coupled binding and folding of intrinsically disordered proteins: what can we learn from kinetics?

Protein or protein regions that are not forming well-defined structures in their free states under native-like conditions are called intrinsically disordered proteins. Such proteins are very common in protein-protein interactions, where their disorder apparently gives several advantages including optimal binding properties. To fully appreciate why protein disorder is advantageous for protein-protein interactions we need to understand the mechanism(s) of interaction. However, elucidating mechanisms in protein-protein interactions is usually very challenging. Here we discuss how kinetics in combination with protein engineering and structural information can be used to depict details of protein-protein interactions involving intrinsically disordered proteins.

Binding Rate Constants Reveal Distinct Features of Disordered Protein Domains.

Intrinsically disordered proteins (IDPs) are abundant in the proteome and involved in key cellular functions. However, experimental data about the binding kinetics of IDPs as a function of different environmental conditions are scarce. We have performed an extensive characterization of the ionic strength dependence of the interaction between the molten globular nuclear co-activator binding domain (NCBD) of CREB binding protein and five different protein ligands, including the intrinsically disordered activation domain of p160 transcriptional co-activators (SRC1, TIF2, ACTR), the p53 transactivation domain, and the folded pointed domain (PNT) of transcription factor ETS-2. Direct comparisons of the binding rate constants under identical conditions show that the association rate constant, kon, for interactions between NCBD and disordered protein domains is high at low salt concentrations (90-350 × 10(6) M(-1) s(-1) at 4 °C) but is reduced significantly (10-30-fold) with an increasing ionic strength and reaches a plateau around physiological ionic strength. In contrast, the kon for the interaction between NCBD and the folded PNT domain is only 7 × 10(6) M(-1) s(-1) (4 °C and low salt) and displays weak ionic strength dependence, which could reflect a distinctly different association that relies less on electrostatic interactions. Furthermore, the basal rate constant (in the absence of electrostatic interactions) is high for the NCBD interactions, exceeding those typically observed for folded proteins. One likely interpretation is that disordered proteins have a large number of possible collisions leading to a productive on-pathway encounter complex, while folded proteins are more restricted in terms of orientation. Our results highlight the importance of electrostatic interactions in binding involving IDPs and emphasize the significance of including ionic strength as a factor in studies that compare the binding properties of IDPs to those of ordered proteins.

Rigidified clicked dimeric ligands for studying the dynamics of the PDZ1-2 supramodule of PSD-95.

PSD-95 is a scaffolding protein of the MAGUK protein family, and engages in several vital protein-protein interactions in the brain with its PDZ domains. It has been suggested that PSD-95 is composed of two supramodules, one of which is the PDZ1-2 tandem domain. Here we have developed rigidified high-affinity dimeric ligands that target the PDZ1-2 supramodule, and established the biophysical parameters of the dynamic PDZ1-2/ligand interactions. By employing ITC, protein NMR, and stopped-flow kinetics this study provides a detailed insight into the overall conformational energetics of the interaction between dimeric ligands and tandem PDZ domains. Our findings expand our understanding of the dynamics of PSD-95 with potential relevance to its biological role in interacting with multivalent receptor complexes and development of novel drugs.

Dynamic and thermodynamic response of the Ras protein Cdc42Hs upon association with the effector domain of PAK3.

Human cell division cycle protein 42 (Cdc42Hs) is a small, Rho-type guanosine triphosphatase involved in multiple cellular processes through its interactions with downstream effectors. The binding domain of one such effector, the actin cytoskeleton-regulating p21-activated kinase 3, is known as PBD46. Nitrogen-15 backbone and carbon-13 methyl NMR relaxation was measured to investigate the dynamical changes in activated GMPPCP·Cdc42Hs upon PBD46 binding. Changes in internal motion of the Cdc42Hs, as revealed by methyl axis order parameters, were observed not only near the Cdc42Hs-PBD46 interface but also in remote sites on the Cdc42Hs molecule. The binding-induced changes in side-chain dynamics propagate along the long axis of Cdc42Hs away from the site of PBD46 binding with sharp distance dependence. Overall, the binding of the PBD46 effector domain on the dynamics of methyl-bearing side chains of Cdc42Hs results in a modest rigidification, which is estimated to correspond to an unfavorable change in conformational entropy of approximately -10kcalmol(-1) at 298K. A cluster of methyl probes closest to the nucleotide-binding pocket of Cdc42Hs becomes more rigid upon binding of PBD46 and is proposed to slow the catalytic hydrolysis of the γ phosphate moiety. An additional cluster of methyl probes surrounding the guanine ring becomes more flexible on binding of PBD46, presumably facilitating nucleotide exchange mediated by a guanosine exchange factor. In addition, the Rho insert helix, which is located at a site remote from the PBD46 binding interface, shows a significant dynamic response to PBD46 binding.

Distinguishing induced fit from conformational selection.

The interactions between proteins and ligands often involve a conformational change in the protein. This conformational change can occur before (conformational selection) or after (induced fit) the association with ligand. It is often very difficult to distinguish induced fit from conformational selection when hyperbolic binding kinetics are observed. In light of a recent paper in this journal (Vogt et al., Biophys. Chem., 186, 2014, 13-21) and the current interest in binding mechanisms emerging from observed sampling of distinct conformations in protein domains, as well as from the field of intrinsically disordered proteins, we here describe a kinetic method that, at least in some cases, unequivocally distinguishes induced fit from conformational selection. The method relies on measuring the observed rate constant λ for binding and varying both the protein and the ligand in separate experiments. Whereas induced fit always yields a hyperbolic dependence of increasing λ values, the conformational selection mechanism gives rise to distinct kinetics when the ligand and protein (displaying the conformational change) concentration is varied in separate experiments. We provide examples from the literature and discuss the limitations of the approach.

A frustrated binding interface for intrinsically disordered proteins.

Intrinsically disordered proteins are very common in the eukaryotic proteome, and many of them are associated with diseases. Disordered proteins usually undergo a coupled binding and folding reaction and often interact with many different binding partners. Using double mutant cycles, we mapped the energy landscape of the binding interface for two interacting disordered domains and found it to be largely suboptimal in terms of interaction free energies, despite relatively high affinity. These data depict a frustrated energy landscape for interactions involving intrinsically disordered proteins, which is likely a result of their functional promiscuity.

Helical propensity in an intrinsically disordered protein accelerates ligand binding.

Many intrinsically disordered proteins fold upon binding to other macromolecules. The secondary structure present in the well-ordered complex is often formed transiently in the unbound state. The consequence of such transient structure for the binding process is, however, not clear. The activation domain of the activator for thyroid hormone and retinoid receptors (ACTR) is intrinsically disordered and folds upon binding to the nuclear coactivator binding domain (NCBD) of the CREB binding protein. A number of mutants was designed that selectively perturbs the amount of secondary structure in unbound ACTR without interfering with the intermolecular interactions between ACTR and NCBD. Using NMR spectroscopy and fluorescence-monitored stopped-flow kinetic measurements we show that the secondary structure content in helix 1 of ACTR indeed influences the binding kinetics. The results thus support the notion of preformed secondary structure as an important determinant for molecular recognition in intrinsically disordered proteins.

The binding mechanisms of intrinsically disordered proteins.

Intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) of proteins are very common and instrumental for cellular signaling. Recently, a number of studies have investigated the kinetic binding mechanisms of IDPs and IDRs. These results allow us to draw conclusions about the energy landscape for the coupled binding and folding of disordered proteins. The association rate constants of IDPs cover a wide range (10(5)-10(9) M(-1) s(-1)) and are largely governed by long-range charge-charge interactions, similarly to interactions between well-folded proteins. Off-rate constants also differ significantly among IDPs (with half-lives of up to several minutes) but are usually around 0.1-1000 s(-1), allowing for rapid dissociation of complexes. Likewise, affinities span from pM to µM suggesting that the low-affinity high-specificity concept for IDPs is not straightforward. Overall, it appears that binding precedes global folding although secondary structure elements such as helices may form before the protein-protein interaction. Short IDPs bind in apparent two-state reactions whereas larger IDPs often display complex multi-step binding reactions. While the two extreme cases of two-step binding (conformational selection and induced fit) or their combination into a square mechanism is an attractive model in theory, it is too simplistic in practice. Experiment and simulation suggest a more complex energy landscape in which IDPs bind targets through a combination of conformational selection before binding (e.g., secondary structure formation) and induced fit after binding (global folding and formation of short-range intermolecular interactions).

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.