Fly-casting Mechanism Research Articles

Exploring Energy Landscapes of Intrinsically Disordered Proteins: Insights into Functional Mechanisms.

Intrinsically disordered proteins (IDPs) lack a rigid three-dimensional structure and populate a polymorphic ensemble of conformations. Because of the lack of a reference conformation, their energy landscape representation in terms of reaction coordinates presents a daunting challenge. Here, our newly developed energy landscape visualization method (ELViM), a reaction coordinate-free approach, shows its prime application to explore frustrated energy landscapes of an intrinsically disordered protein, prostate-associated gene 4 (PAGE4). PAGE4 is a transcriptional coactivator that potentiates the oncogene c-Jun. Two kinases, namely, HIPK1 and CLK2, phosphorylate PAGE4, generating variants phosphorylated at different serine/threonine residues (HIPK1-PAGE4 and CLK2-PAGE4, respectively) with opposing functions. While HIPK1-PAGE4 predominantly phosphorylates Thr51 and potentiates c-Jun, CLK2-PAGE4 hyperphosphorylates PAGE4 and attenuates transactivation. To understand the underlying mechanisms of conformational diversity among different phosphoforms, we have analyzed their atomistic trajectories simulated using AWSEM forcefield, and the energy landscapes were elucidated using ELViM. This method allows us to identify and compare the population distributions of different conformational ensembles of PAGE4 phosphoforms using the same effective phase space. The results reveal a predominant conformational ensemble with an extended C-terminal segment of WT PAGE4, which exposes a functional residue Thr51, implying its potential of undertaking a fly-casting mechanism while binding to its cognate partner. In contrast, for HIPK1-PAGE4, a compact conformational ensemble enhances its population sequestering phosphorylated-Thr51. This clearly explains the experimentally observed weaker affinity of HIPK1-PAGE4 for c-Jun. ELViM appears as a powerful tool, especially to analyze the highly frustrated energy landscape representation of IDPs where appropriate reaction coordinates are hard to apprehend.

Uncovering a membrane-distal conformation of KRAS available to recruit RAF to the plasma membrane

The small GTPase KRAS is localized at the plasma membrane where it functions as a molecular switch, coupling extracellular growth factor stimulation to intracellular signaling networks. In this process, KRAS recruits effectors, such as RAF kinase, to the plasma membrane where they are activated by a series of complex molecular steps. Defining the membrane-bound state of KRAS is fundamental to understanding the activation of RAF kinase and in evaluating novel therapeutic opportunities for the inhibition of oncogenic KRAS-mediated signaling. We combined multiple biophysical measurements and computational methodologies to generate a consensus model for authentically processed, membrane-anchored KRAS. In contrast to the two membrane-proximal conformations previously reported, we identify a third significantly populated state using a combination of neutron reflectivity, fast photochemical oxidation of proteins (FPOP), and NMR. In this highly populated state, which we refer to as "membrane-distal" and estimate to comprise ∼90% of the ensemble, the G-domain does not directly contact the membrane but is tethered via its C-terminal hypervariable region and carboxymethylated farnesyl moiety, as shown by FPOP. Subsequent interaction of the RAF1 RAS binding domain with KRAS does not significantly change G-domain configurations on the membrane but affects their relative populations. Overall, our results are consistent with a directional fly-casting mechanism for KRAS, in which the membrane-distal state of the G-domain can effectively recruit RAF kinase from the cytoplasm for activation at the membrane.

The BADC and BCCP subunits of chloroplast acetyl-CoA carboxylase sense the pH changes of the light–dark cycle

Acetyl-CoA carboxylase (ACCase) catalyzes the first committed step in the de novo synthesis of fatty acids. The multisubunit ACCase in the chloroplast is activated by a shift to pH 8 upon light adaptation and is inhibited by a shift to pH 7 upon dark adaptation. Here, titrations with the purified ACCase biotin attachment domain-containing (BADC) and biotin carboxyl carrier protein (BCCP) subunits from Arabidopsis indicated that they can competently and independently bind biotin carboxylase (BC) but differ in responses to pH changes representing those in the plastid stroma during light or dark conditions. At pH 7 in phosphate buffer, BADC1 and BADC2 gain an advantage over BCCP1 and BCCP2 in affinity for BC. At pH 8 in KCl solution, however, BCCP1 and BCCP2 had more than 10-fold higher affinity for BC than did BADC1. The pH-modulated shifts in BC preferences for BCCP and BADC partners suggest they contribute to light-dependent regulation of heteromeric ACCase. Using NMR spectroscopy, we found evidence for increased intrinsic disorder of the BADC and BCCP subunits at pH 7. We propose that this intrinsic disorder potentially promotes fast association with BC through a "fly-casting mechanism." We hypothesize that the pH effects on the BADC and BCCP subunits attenuate ACCase activity by night and enhance it by day. Consistent with this hypothesis, Arabidopsis badc1 badc3 mutant lines grown in a light-dark cycle synthesized more fatty acids in their seeds. In summary, our findings provide evidence that the BADC and BCCP subunits function as pH sensors required for light-dependent switching of heteromeric ACCase activity.

Free-energy landscape of molecular interactions between endothelin 1 and human endothelin type B receptor: fly-casting mechanism.

The free-energy landscape of interaction between a medium-sized peptide, endothelin 1 (ET1), and its receptor, human endothelin type B receptor (hETB), was computed using multidimensional virtual-system coupled molecular dynamics, which controls the system’s motions by introducing multiple reaction coordinates. The hETB embedded in lipid bilayer was immersed in explicit solvent. All molecules were expressed as all-atom models. The resultant free-energy landscape had five ranges with decreasing ET1–hETB distance: completely dissociative, outside-gate, gate, binding pocket, and genuine-bound ranges. In the completely dissociative range, no ET1–hETB interaction appeared. In the outside-gate range, an ET1–hETB attractive interaction was the fly-casting mechanism. In the gate range, the ET1 orientational variety decreased rapidly. In the binding pocket range, ET1 was in a narrow pathway with a steep free-energy slope. In the genuine-bound range, ET1 was in a stable free-energy basin. A G-protein-coupled receptor (GPCR) might capture its ligand from a distant place.

The Membrane Tether of the RAS Signaling Protein Drives Nanoclustering by Fly-Casting for Anionic Lipids

Cancer is a devastating disease that results from pathological modifications to cellular decision making processes. In healthy cells, the protein-mediated signaling networks that regulate growth and movement are tightly regulated. However, mutations that disrupt or over-activate signaling proteins can drive uncontrolled cell growth. Ras, a peripheral membrane signaling protein, harbors activating mutations in 20% of all cancers, especially those of the pancreas, colon, thyroid, and lung. In theory, the signaling output from overactive Ras mutants can be decreased by small molecules that bind and stabilize inactive states. However, crystallographic structures of Ras reveal a relatively smooth surface lacking substantial hydrophobic crevices for drug binding. Therefore, we use molecular dynamics simulations to explore the topography of Ras proteins and the dynamic accessibility of their effector-binding surfaces. Our long-term goal is the identification of potentially druggable cavities at intermolecular interfaces, including those between Ras and the cell membrane. Our molecular dynamics simulations indicate that, in Ras isoform K-Ras4B, the prenylated 19-residue C-terminal highly variable region (HVR) membrane anchor remains disordered as it pulls Ras' globular G domain closer to membranes enriched in negatively charged (anionic) lipids. Furthermore, we show that although Ras' globular G domain interacts promiscuously with its own HVR and membrane lipids, there are preferential and potentially targetable interactions. Finally, simulations of multiple HVR peptides indicate that intermolecular contacts predominantly involve the G domain-proximal (N-terminal) region of the HVR, consistent with a fly-casting mechanism that augments the self-association of Ras proteins in the context of disordered membrane patches enriched in anionic lipids. This targeted adsorption of the K-Ras4B HVR to anionic membrane surfaces may play a role in regulating the multimerization ability of downstream effector proteins such as Raf kinase.

On the role of residue phosphorylation in 14-3-3 partners: AANAT as a case study

Twenty years ago, a novel concept in protein structural biology was discovered: the intrinsically disordered regions (IDRs). These regions remain largely unstructured under native conditions and the more are studied, more properties are attributed to them. Possibly, one of the most important is their ability to conform a new type of protein-protein interaction. Besides the classical domain-to-domain interactions, IDRs follow a ‘fly-casting’ model including ‘induced folding’. Unfortunately, it is only possible to experimentally explore initial and final states. However, the complete movie of conformational changes of protein regions and their characterization can be addressed by in silico experiments. Here, we simulate the binding of two proteins to describe how the phosphorylation of a single residue modulates the entire process. 14-3-3 protein family is considered a master regulator of phosphorylated proteins and from a modern point-of-view, protein phosphorylation is a three component system, with writers (kinases), erasers (phosphatases) and readers. This later biological role is attributed to the 14-3-3 protein family. Our molecular dynamics results show that phosphorylation of the key residue Thr31 in a partner of 14-3-3, the aralkylamine N-acetyltransferase, releases the fly-casting mechanism during binding. On the other hand, the non-phosphorylation of the same residue traps the proteins, systematically and repeatedly driving the simulations into wrong protein-protein conformations.

Exploring the Role of Flexibility in Binding Kinetics and Affinity of pKID-Kix through Coarse Grained Simulations

Physical properties that distinguish Intrinsically Disordered Proteins (IDPs) from structured proteins include increased association rates and the formation of bound complexes with high specificity and low affinity. In this work, we investigate the influence of flexibility on binding kinetics and thermodynamics of IDPs through coarse grained simulations. Using pKID-KIX complex as an example, we illustrate how coupled folding and binding allows the IDP to have both fast binding and unbinding rates together with low binding affinity. In particular, we show that while flexibility moderately increases binding rate (consistent with the fly-casting mechanism), the unbinding rate is enhanced to a greater extent, resulting in a low affinity complex. Analysis of equilibrium trajectories show that a flexible protein forms contacts incrementally throughout binding and that orientational constraints are satisfied gradually; both of which encourage a lower free energy barrier controlling binding and unbinding kinetics. Our study highlights the importance that flexibility has not only on an IDP's binding rate, but also on its ability to form low affinity complexes with sufficiently fast unbinding kinetics essential to their function.

Intrinsic disorder accelerates dissociation rather than association

The intrinsically disordered protein (IDP) has distinct properties both physically and biologically: it often becomes folded when binding to the target and is frequently involved in signal transduction. The physical property seems to be compatible with the biological property where fast association and dissociation between IDP and the target are required. While fast association has been well studied, fueled by the fly-casting mechanism, the dissociation kinetics has received less attention. We here study how the intrinsic disorder affects the dissociation kinetics, as well as the association kinetics, paying attention to the interaction strength at the binding site (i.e., the quality of the "fly lure"). Coarse-grained molecular dynamics simulation of the pKID-KIX system, a well-studied IDP system, shows that the association rate becomes larger as the disorder-inducing flexibility that was imparted to the model is increased, but the acceleration is marginal and turns into deceleration as the quality of the fly lure is worsened. In contrast, the dissociation rate is greatly enhanced as the disorder is increased, indicating that intrinsic disorder serves for rapid signal switching more effectively through dissociation than association. Proteins 2016; 84:1124-1133. © 2016 Wiley Periodicals, Inc.

Different Roles of Backbone-Originated and Native-Contact-Originated Flexibilities of an Intrinsically Disordered Protein in Fast Binding and Unbinding

The intrinsically disordered protein (IDP) in the unbound state is so flexible that it does not adopt a unique structure. IDP may exploit this flexibility to function as a signaling molecule, switching on and off the signal by binding to and unbinding from the target molecule, considering that the binding/unbinding kinetics is largely affected by the flexibility. For example, it has been suggested that the flexibility enhances the binding rate by the fly-casting mechanism. However, some counter evidences have been reported, so the role of the flexibility on the fly-casting mechanism is still controversial. To address this problem, we examined the effect of two different types of flexibilities that were computationally imparted to IDP to which the one-bead (Cα) Gō-like model was applied; one is by softening the backbone rigidity and the other is by weakening the native contact. We used pKID (kinase inducible domain of CREB) as a model of IDP. pKID folds upon binding to the partner protein, KIX. We calculated the potential of mean force as a function of distance between pKID and KIX. We found that the backbone-originated flexibility and the native-contact-originated one exhibit different results. The native-contact-originated flexibility yielded larger capturing radius and faster binding, consistent with the fly-casting, but the backbone-originated one did not. The strength of the interaction between IDP and the target was also shown to play a critical role. As for unbinding, the backbone-originated flexibility destabilized the bound state more effectively, leading to faster unbinding. Therefore, IDP can exploit these two different types of flexibilities to optimize their signaling function.

Hierarchical Folding of the rRNA in the Early Assembly of the E. Coli Ribosomal Small Subunit

In order to investigate the co-transcriptional folding and assembly of the E. coli ribosomal small subunit (SSU), we applied the structure-based Go model to fold 16S rRNA fragments with increasing length. The first fragment examined is the binding domain of S20, which includes the first few local and non-local helices to be transcribed (h6, h7-h10) in vivo. The simulations show that the main barriers to the formation of the S20 binding domain are the docking of the 3' and 5' strands of the non-local h7 followed by re-orientation of the helices in the four-way junction to allow tertiary contacts between h6 and h8 to be established. These observations suggest the role of the primary binding protein S20 in stabilizing key tertiary contacts and speeding up the protein:RNA recognition during the assembly process.Transcription of rRNA then proceeds to the local five-way junction (5WJ), which nucleates the assembly of the 5' domain upon interaction with the initiator and primary binding protein S4. This complex contains the largest contiguous ribosomal structure signature, helix h16, which is a conserved feature in all bacterial 16S rRNAs. We use all-atom molecular dynamics simulations to demonstrate that the co-evolving signature in the N-terminus of S4 is intrinsically disordered and capable of accelerating the 5WJ:S4 binding process through a fly-casting mechanism. The scenario is further confirmed using the structure-based Go model for protein:RNA, with which hundreds of simultaneous folding and binding events between the 5WJ and S4 are captured. Based on analysis of both explicit solvent MD and structure-based Go trajectories, we identify multiple metastable states of the 5WJ, and various protein-guided folding pathways that compare directly with the FRET experiments.

Complex Energy Landscape of a Giant Repeat Protein

Here, we reveal a remarkable complexity in the unfolding of giant HEAT-repeat protein PR65/A, a molecular adaptor for the heterotrimeric PP2A phosphatases. The repeat array ruptures at multiple sites, leading to intermediate states with noncontiguous folded subdomains. There is a dominant sequence of unfolding, which reflects a nonuniform stability distribution across the repeat array and can be rationalized by theoretical models accounting for heterogeneous contact density in the folded structure. Unfolding of certain intermediates is, however, competitive, leading to parallel unfolding pathways. The low-stability, central repeats sample unfolded conformations under physiological conditions, suggesting how folding directs function: certain regions present rigid motifs for molecular recognition, whereas others have the flexibility with which to broaden the search area, as in the fly-casting mechanism. Partial unfolding of PR65/A also impacts catalysis by altering the proximity of bound catalytic subunit and substrate. Thus, the repeat array orchestrates the assembly and activity of PP2A.

Replica exchange molecular dynamics simulations of an α/β-type small acid soluble protein (SASP)

Small acid soluble proteins (SASPs) of α/β-type play a major role in the resistance of spore DNAs to external assaults. It has been found that α/β-type SASP exhibits intrinsic disorder on isolation, but it acquires a defined native state upon binding to DNA. This disorder to order transition is not yet understood. Other questions related to the role of the thermodynamics and structure of the individual protein in the complex formation remain elusive. Characterization of the unbound state of α/β-type SASP in experiments could be a challenging problem because of the heterogeneous nature of the ensemble. Here, computer simulations can help gain more insights into the unbound state of α/β-type SASP. In the present work, by using replica exchange molecular dynamics (REMD), we simulated an α/β-type SASP on isolation with an implicit solvent. We found that α/β-type SASP undergoes a continuous phase transition with a small free energy barrier, a common feature of intrinsically disordered proteins (IDPs). Additionally, we detected the presence of residual α-helical structures at local level and a high degree of plasticity in the chain which can contribute to the fast disorder to order transition by reducing the fly-casting mechanism.

Atomistic Description of the Folding of a Dimeric Protein

Equilibrium molecular dynamics simulations are increasingly being used to describe the folding of individual proteins and protein domains at an atomic level of detail. Isolated protein domains, however, are rarely observed in vivo, where multidomain proteins and multimeric assemblies are far more common. It is clear that the folding of such proteins is often inextricably coupled with the process of dimerization; indeed, many protein monomers and protein domains are not stable in isolation, and fold to their native structures only when stabilized by interactions with other members of a protein complex. Here, we use equilibrium molecular dynamics simulations with an aggregate simulation length of 4 ms to elucidate key aspects of the folding mechanism, and of the associated free-energy surface, of the Top7-CFr dimer, a 114-amino-acid protein homodimer with a mixed α/β structure. In these simulations, we observed a number of folding and unfolding events. Each folding event was characterized by the assembly of two unfolded Top7-CFr monomers to form a stable, folded dimer. We found that the isolated monomer is unstable but that, early in the folding pathway, nascent native structure is stabilized by contacts between the two monomer subunits. These contacts are in some cases native, as in an induced-folding model, and in other cases non-native, as in a fly-casting mechanism. Occasionally, folding by conformational selection, in which both subunits form independently before dimerization, was also observed. Folding then proceeds through the sequential addition of strands to the protein β sheet. Although the long-time-scale relaxation of the folding process can be well described by a single exponential, these simulations reveal the presence of a number of kinetic traps, characterized by structures in which individual strands are added in an incorrect order.

Binding and Folding of the Small Bacterial Chaperone HdeA

The small pH stress-sensing chaperone HdeA helps pathogenic enteric E. coli survive passage through the severely acidic environment of the mammalian stomach. Under stress conditions, HdeA transitions from an inactive folded dimer to a chaperone-active unfolded monomer to prevent the acid-induced aggregation of periplasmic proteins. Here we use a topology-based Gō-like model to delineate the relationship between dimer interface formation and monomer folding and to better understand the structural details of the chaperone activation mechanism. Free energy surfaces show that dimer interface formation and monomer folding proceed concurrently through an on-pathway dimeric intermediate in which one monomer is partially unfolded. The absence of a preexisting fully folded monomer in the proposed binding mechanism is in agreement with HdeA's rapid chaperone response. Binding between unfolded monomers exhibits an enhancement of molecular recognition reminiscent of the fly-casting mechanism. Overall, our simulations further highlight the efficient nature of HdeA's chaperone response and we anticipate that knowledge of a dimeric intermediate will facilitate the interpretation of experimental studies.

Capture and Quality Control Mechanisms for Adenosine-5′-triphosphate Binding

The catalytic events in members of the nucleotidylyl transferase superfamily are initiated by a millisecond binding of ATP in the active site. Through metadynamics simulations on a class I aminoacyl-tRNA synthetase (aaRSs), the largest group in the superfamily, we calculate the free energy landscape of ATP selection and binding. Mutagenesis studies and fluorescence spectroscopy validated the identification of the most populated intermediate states. The rapid first binding step involves formation of encounter complexes captured through a fly casting mechanism that acts upon the triphosphate moiety of ATP. In the slower nucleoside binding step, a conserved histidine in the HxxH motif orients the incoming ATP through base-stacking interactions resulting in a deep minimum in the free energy surface. Mutation of this histidine significantly decreases the binding affinity measured experimentally and computationally. The metadynamics simulations further reveal an intermediate quality control state that the synthetases and most likely other members of the superfamily use to select ATP over other nucleoside triphosphates.

Structural and kinetic effects of hypertrophic cardiomyopathy related mutations R146G/Q and R163W on the regulatory switching activity of rat cardiac troponin I

Mutations in cardiac troponin I (cTnI) that cause hypertrophic cardiomyopathy (HCM) have been reported to change the contractility of cardiac myofilaments, but the underlying molecular mechanism remains elusive. In this study, Förster resonance energy transfer (FRET) was used to investigate the specific structural and kinetic effects that HCM related rat cTnI mutations R146G/Q and R163W exert on Ca2+ and myosin S1 dependent conformational transitions in rat cTn structure. Ca2+-induced changes in interactions between cTnC and cTnI were individually monitored in reconstituted thin filaments using steady state and time resolved FRET, and kinetics were determined using stopped flow. R146G/Q and R163W all changed the FRET distances between cTnC and cTnI in unique and various ways. However, kinetic rates of conformational transitions induced by Ca2+-dissociation were universally slowed when R146G/Q and R163W were present. Interestingly, the kinetic rates of changes in the inhibitory region of cTnI were always slower than that of the regulatory region, suggesting that the fly casting mechanism that normally underlies deactivation is preserved in spite of mutation. In situ rat myocardial fiber studies also revealed that FRET distance changes indicating mutation specific disruption of the cTnIIR−actin interaction were consistent with increased passive tension.

Predictive energy landscapes for protein–protein association

We investigate protein-protein association using the associative-memory, water-mediated, structure, and energy model (AWSEM), a coarse-grained protein folding model that has been optimized using energy-landscape theory. The potential was originally parameterized by enforcing a funneled nature for a database of dimeric interfaces but was later further optimized to create funneled folding landscapes for individual monomeric proteins. The ability of the model to predict interfaces was not tested previously. The present results show that simulated annealing of the model indeed is able to predict successfully the native interfaces of eight homodimers and four heterodimers, thus amounting to a flexible docking algorithm. We go on to address the relative importance of monomer geometry, flexibility, and nonnative intermonomeric contacts in the association process for the homodimers. Monomer surface geometry is found to be important in determining the binding interface, but it is insufficient. Using a uniform binding potential rather than the water-mediated potential results in sampling of misbound structures that are geometrically preferred but are nonetheless energetically disfavored by AWSEM, as well as in nature. Depending on the stability of the unbound monomers, nonnative contacts play different roles in the association process. For unstable monomers, thermodynamic states stabilized by nonnative interactions correspond to productive, on-pathway intermediates and can, therefore, catalyze binding through a fly-casting mechanism. For stable monomers, in contrast, states stabilized by nonnative interactions generally correspond to traps that impede binding.

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.

Assembly of the Five-Way Junction in the Ribosomal Small Subunit Using Hybrid MD-Go̅ Simulations

Assembly of the bacterial ribosomal small subunit (SSU) begins with the folding of the five-way junction upon interaction with the primary binding protein S4. This complex contains the largest contiguous molecular signature, which is a conserved feature in all bacterial 16S rRNAs. In a previous study, we used all-atom molecular dynamics simulations to demonstrate that the co-evolving signature in the N-terminus of S4 is intrinsically disordered and capable of accelerating the binding process through a fly casting mechanism. In this paper, comparisons between the all-atom MD simulations and FRET experiments identify multiple metastable conformations of the naked five-way junction without the presence of S4. Furthermore, we capture the simultaneous folding and binding of the five-way junction and r-protein S4 by use of a structure-based Gō potential implemented within the framework of the all-atom molecular dynamics CHARMM force field. Different folding pathways are observed for the refolding of the five-way junction upon partial binding of S4. Our simulations illustrate the complex nature of RNA folding in the presence of a protein binding partner and provide insight into the role of population shift and the induced fit mechanisms in the protein:RNA folding and binding process.

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.