Intrinsic disorder as a mechanism to optimize allosteric coupling in proteins

Molecular mechanism of SSR128129E, an extracellularly acting, small-molecule, allosteric inhibitor of FGF receptor signaling

The fibroblast growth factor (FGF)/fibroblast growth factor receptor (FGFR) signaling network plays an important role in cell growth, survival, differentiation, and angiogenesis. Deregulation of FGFR signaling can lead to cancer development. Here, we report an FGFR inhibitor, SSR128129E (SSR), that binds to the extracellular part of the receptor. SSR does not compete with FGF for binding to FGFR but inhibits FGF-induced signaling linked to FGFR internalization in an allosteric manner, as shown by crystallography studies, nuclear magnetic resonance, Fourier transform infrared spectroscopy, molecular dynamics simulations, free energy calculations, structure-activity relationship analysis, and FGFR mutagenesis. Overall, SSR is a small molecule allosteric inhibitor of FGF/FGFR signaling, acting via binding to the extracellular part of the FGFR.

Ligand-induced dynamic changes in extended PDZ domains from NHERF1

The multi-domain scaffolding protein NHERF1 modulates the assembly and intracellular trafficking of various transmembrane receptors and ion-transport proteins. The two PDZ (postsynaptic density 95/disk large/zonula occluden 1) domains of NHERF1 possess very different ligand-binding capabilities: PDZ1 recognizes a variety of membrane proteins with high affinity, while PDZ2 only binds limited number of target proteins. Here using NMR, we have determined the structural and dynamic mechanisms that differentiate the binding affinities of the two PDZ domains, for the type 1 PDZ-binding motif (QDTRL) in the carboxyl terminus of cystic fibrosis transmembrane regulator. Similar to PDZ2, we have identified a helix-loop-helix subdomain coupled to the canonical PDZ1 domain. The extended PDZ1 domain is highly flexible with correlated backbone motions on fast and slow timescales, while the extended PDZ2 domain is relatively rigid. The malleability of the extended PDZ1 structure facilitates the transmission of conformational changes at the ligand-binding site to the remote helix-loop-helix extension. By contrast, ligand binding has only modest effects on the conformation and dynamics of the extended PDZ2 domain. The study shows that ligand-induced structural and dynamic changes coupled with sequence variation at the putative PDZ binding site dictate ligand selectivity and binding affinity of the two PDZ domains of NHERF1.

Structural and energetic basis of allostery

Allostery is a biological phenomenon of fundamental importance in regulation and signaling, and efforts to understand this process have led to the development of numerous models. In spite of individual successes in understanding the structural determinants of allostery in well-documented systems, much less success has been achieved in identifying a set of quantitative and transferable ground rules that provide an understanding of how allostery works. Are there organizing principles that allow us to relate structurally different proteins, or are the determinants of allostery unique to each system? Using an ensemble-based model, we show that allosteric phenomena can be formulated in terms of conformational free energies of the cooperative elements in a protein and the coupling interactions between them. Interestingly, the resulting allosteric ground rules provide a framework to reconcile observations that challenge purely structural models of site-to-site coupling, including (a) allostery in the absence of pathways of structural distortions, (b) allostery in the absence of any structural change, and (c) the ability of allosteric ligands to act as agonists under some circumstances and antagonists under others. The ensemble view of allostery that emerges provides insights into the energetic prerequisites of site-to-site coupling and thus into how allostery works.

The disordered C-terminal domain of human DNA glycosylase NEIL1 contributes to its stability via intramolecular interactions

NEIL1 [Nei (endonuclease VIII)-like protein 1], one of the five mammalian DNA glycosylases that excise oxidized DNA base lesions in the human genome to initiate base excision repair, contains an intrinsically disordered C-terminal domain (CTD; ~100 residues), not conserved in its Escherichia coli prototype Nei. Although dispensable for NEIL1's lesion excision and AP lyase activities, this segment is required for efficient in vivo enzymatic activity and may provide an interaction interface for many of NEIL1's interactions with other base excision repair proteins. Here, we show that the CTD interacts with the folded domain in native NEIL1 containing 389 residues. The CTD is poised for local folding in an ordered structure that is induced in the purified fragment by osmolytes. Furthermore, deletion of the disordered tail lacking both Tyr and Trp residues causes a red shift in NEIL1's intrinsic Trp-specific fluorescence, indicating a more solvent-exposed environment for the Trp residues in the truncated protein, which also exhibits reduced stability compared to the native enzyme. These observations are consistent with stabilization of the native NEIL1 structure via intramolecular, mostly electrostatic, interactions that were disrupted by mutating a positively charged (Lys-rich) cluster of residues (amino acids 355-360) near the C-terminus. Small-angle X-ray scattering (SAXS) analysis confirms the flexibility and dynamic nature of NEIL1's CTD, a feature that may be critical to providing specificity for NEIL1's multiple, functional interactions.

Importance of protein dynamics during enzymatic C-H bond cleavage catalysis

Quantum tunneling and protein dynamics have emerged as important components of enzyme function. This review focuses on soybean lipoxygenase-1, to illustrate how the properties of enzymatic C-H bond activation link protein motions to the fundamental bond making-breaking processes.

Interplay between allostery and intrinsic disorder in an ensemble

Allostery is a biological phenomenon of critical importance in metabolic regulation and cell signalling. The fundamental premise of classical models that describe allostery is that structure mediates 'action at a distance'. Recently, this paradigm has been challenged by the enrichment of IDPs (intrinsically disordered proteins) or ID (intrinsically disordered) segments in transcription factors and signalling pathways of higher organisms, where an allosteric response from external signals is requisite for regulated function. This observation strongly suggests that IDPs elicit the capacity for finely tunable allosteric regulation. Is there a set of transferable ground rules that reconcile these disparate allosteric phenomena? We focus on findings from the human GR (glucocorticoid receptor) which is a nuclear transcription factor in the SHR (steroid hormone receptor) family. GR contains an intrinsically disordered NTD (N-terminal domain) that is obligatory for transcription activity. Different GR translational isoforms have various lengths of NTD and by studying these isoforms we found that the full-length ID NTD consists of two thermodynamically distinct coupled regions. The data are interpreted in the context of an EAM (ensemble allosteric model) that considers only the intrinsic and measurable energetics of allosteric systems. Expansion of the EAM is able to reconcile the paradox that ligands for SHRs can be agonists and antagonists in a cell-context-dependent manner. These findings suggest a mechanism by which SHRs in particular, and IDPs in general, may have evolved to couple thermodynamically distinct ID segments. The ensemble view of allostery that is illuminated provides organizing principles to unify the description of all allosteric systems and insight into 'how' allostery works.

Temperature and urea have opposing impacts on polyproline II conformational bias

The native states of globular proteins have been accessed in atomic detail by X-ray crystallography and nuclear magnetic resonance spectroscopy, yet characterization of denatured proteins beyond global metrics has proven to be elusive. Denatured proteins have been observed to exhibit global geometric properties of a random coil polymer. However, this does not preclude the existence of nonrandom, local conformational bias that may be significant for protein folding and function. Indeed, circular dichroism (CD) spectroscopy and other methods have suggested that the denatured state contains considerable local bias to the polyproline II (PII) conformation. Here, we develop predictive models to determine the extent that temperature and the chemical denaturant urea modulate PII propensity. In agreement with our predictive model, PII propensity is observed experimentally to decrease with an increase in temperature. Conversely, urea appears to promote the PII conformation as determined by CD and isothermal titration calorimetry. Importantly, the calorimetric data are in quantitative agreement with a model that predicts the stability of the PII helix relative to other denatured state conformations based upon solvent accessible surface area and experimentally measured Gibbs transfer free energies. The ability of urea to promote the PII conformation can be attributed to the favorable interaction of urea with the peptide backbone. Thus, perturbing denatured states by temperature or cosolutes has subtle, yet opposing, impacts on local PII conformational biases. These results have implications for protein folding as well as for the function of signaling proteins that bind proline-rich targets in globular or intrinsically disordered proteins.

Potential conformational heterogeneity of p53 bound to S100B(ββ)

The negative regulatory domain (NRD) of the p53 tumor suppressor is intrinsically disordered. It contains several posttranslational modification (PTM) sites that are important for regulation of p53 activity. Calcium-dependent binding of dimeric S100B(ββ) to p53-NRD blocks access to these PTM sites and disrupts the p53 tetramer to inhibit p53 activation. Previous nuclear magnetic resonance (NMR) structural studies have suggested that p53-NRD folds into a stable helix upon binding to S100B(ββ). Intriguingly, despite the well-converged and stably folded nature of the NMR structure ensemble, experimentally resolved intermolecular nuclear Overhauser enhancements (NOEs) are extremely weak; most have 5- to 6-Å upper bounds, and mainly involve the C-terminal segment of p53-NRD. Such a systematic lack of strong intermolecular NOEs could suggest that the p53/S100B(ββ) interface is more dynamic than currently believed. Indeed, extensive atomistic simulations in explicit solvent (with 1.0μs total effective sampling) revealed large heterogeneity in the S100B(ββ)-bound conformation of p53-NRD. Helix unwinding at the C-terminus allows key hydrophobic residues (Leu383 and Phe385) to make more extensive intermolecular contacts, whereas the highly helical N-terminus displays substantial flexibility in packing with S100B(ββ). Importantly, the predicted heterogeneous ensemble as a whole is highly consistent with experimental intermolecular NOEs, although many conformational sub-states coexist and individual sub-states satisfy only subsets of the NOE restraints. Furthermore, the simulated ensemble provides similar shielding of key PTM sites to support p53 inhibition. This study not only provides new insights into the structural basis of the p53/S100B(ββ) recognition but also highlights the importance of recognizing dynamic complexes in structural studies of intrinsically disordered protein interactions.

Toward an understanding of the sequence and structural basis of allosteric proteins

Allostery is the most efficient means of regulating protein functions, ranging from the control of metabolic mechanisms to signal transduction pathways. Although allosteric regulation has been recognized for half a century, our knowledge is limited to the characteristics of allosteric proteins and the structural coupling of allosteric sites and modulators. In this paper, we present a comprehensive analysis of allosteric proteins that provides insight into the foundation of allosteric interactions by revealing a series of common features in the allosteric proteins. Allosteric proteins mainly appear in transferases, and phosphorylation is the most common type of modification found in allosteric proteins. Disorders related to allosteric proteins primarily comprise metabolic diseases and cancers. In general, allosteric proteins prefer to exist as monomers or even-numbered multimers. Greater stability and hydrophobicity are observed in allosteric proteins than in general proteins. Further analysis of the allosteric sites reveals a series of buried and compact pockets composed of significantly greater hydrophobic surface area than the corresponding orthosteric sites. The hydrophobicity of the allosteric sites plays a dominant role in the binding of allosteric modulators as observed in the analysis of 106 diverse allosteric protein-modulator pairs. These results may be of great significance in predicting which proteins are allosteric and in designing novel triggers to inhibit or activate proteins of interest.

Disruption of disulfide-restriction at integrin knees induces activation and ligand-independent signaling of α4β7

Control of integrin activation and signaling plays critical roles in cell adhesion, spreading, and migration. Here, we report that selective breakage of two conserved disulfide bonds located at the knees of integrin, α4C589–C594 and β7C494–C526, induced α4β7 activation. This activated α4β7 had a unique structure different from the typical extended conformation of active integrin. In addition, these activated α4β7 integrins spontaneously clustered on the cell membrane and triggered integrin downstream signaling independent of ligand binding. Although these disulfide bonds were not broken during α4β7 activation by inside-out signaling or Mn2+, they could be specifically reduced by 0.1 mM dithiothreitol, a reducing strength that could be produced in vivo under certain conditions. Our findings reveal a novel mechanism of integrin activation under specific reducing conditions by which integrin can signal and promote cell spreading in the absence of ligand.

Conformational heterogeneity within the LID domain mediates substrate binding to Escherichia coli adenylate kinase: function follows fluctuations

Proteins exist as dynamic ensembles of molecules, implying that protein amino acid sequences evolved to code for both the ground state structure as well as the entire energy landscape of excited states. Accumulating theoretical and experimental evidence suggests that enzymes use such conformational fluctuations to facilitate allosteric processes important for substrate binding and possibly catalysis. This phenomenon can be clearly demonstrated in Escherichia coli adenylate kinase, where experimentally observed local unfolding of the LID subdomain, as opposed to a more commonly postulated rigid-body opening motion, is related to substrate binding. Because "entropy promoting" glycine mutations designed to increase specifically the local unfolding of the LID domain also affect substrate binding, changes in the excited energy landscape effectively tune the function of this enzyme without changing the ground state structure or the catalytic site. Thus, additional thermodynamic information, above and beyond the single folded structure of an enzyme-substrate complex, is likely required for a full and quantitative understanding of how enzymes work.

The dark energy of proteins comes to light: conformational entropy and its role in protein function revealed by NMR relaxation

Historically it has been virtually impossible to experimentally determine the contribution of residual protein entropy to fundamental protein activities such as the binding of ligands. Recent progress has illuminated the possibility of employing NMR relaxation methods to quantitatively determine the role of changes in conformational entropy in molecular recognition by proteins. The method rests on using fast internal protein dynamics as a proxy. Initial results reveal a large and variable role for conformational entropy in the binding of ligands by proteins. Such a role for conformational entropy in molecular recognition has significant implications for enzymology, signal transduction, allosteric regulation and the development of protein-directed pharmaceuticals.

Binding of the N-terminal region of coactivator TIF2 to the intrinsically disordered AF1 domain of the glucocorticoid receptor is accompanied by conformational reorganizations

Control of gene transcription by glucocorticoid receptors (GRs) is important for many physiological processes. Like other steroid hormone receptors, the regulation of target genes by GR is mediated by two transactivation domains: activation function 1 (AF1) in the N-terminal domain and AF2 in the C-terminal ligand-binding domain (LBD). Full receptor activity requires both AF1 and -2 plus assorted coregulatory proteins. Crystal structures of the ligand-bound LBD have provided insight regarding how AF2 interacts with specific coactivators. However, despite its being the major activation domain of GRs, knowledge of AF1 structure/function has languished. This is mainly because of the highly disorganized structure of the GR N-terminal domain. This lack of AF1 structure is shared by all members of the steroid/nuclear receptor superfamily for which it has been examined and AF1 is thought to allow productive interactions with assorted cofactors via protein-induced changes in secondary/tertiary structures. To date, there are no reports of a classical coactivator altering the secondary/tertiary structure of the GR AF1 domain. Earlier, we reported an N-terminal fragment of the p160 coactivator TIF2, called TIF2.0, that binds the GR N-terminal domain and alters GR transcriptional activity. We therefore proposed that TIF2.0 binding to AF1 changes both its conformation and transcriptional activity. We now report that TIF2.0 interacts with the GR AF1 domain to increase the amount of α-helical structure in the complex. Furthermore, TIF2 coactivator activity is observed in the absence of the GR LBD in a manner that requires the AF1 domain. This contrasts with previous models where TIF2 receptor interaction domains binding to GR LBD somehow alter AF1 conformation. Our results establish for the first time that coactivators can modify the structure of the AF1 domain directly via the binding of a second region of the coactivator and suggest a molecular explanation for how coactivators increase the transcriptional activity of GR-agonist complexes.

Tissue-specific splicing of disordered segments that embed binding motifs rewires protein interaction networks

Alternative inclusion of exons increases the functional diversity of proteins. Among alternatively spliced exons, tissue-specific exons play a critical role in maintaining tissue identity. This raises the question of how tissue-specific protein-coding exons influence protein function. Here we investigate the structural, functional, interaction, and evolutionary properties of constitutive, tissue-specific, and other alternative exons in human. We find that tissue-specific protein segments often contain disordered regions, are enriched in posttranslational modification sites, and frequently embed conserved binding motifs. Furthermore, genes containing tissue-specific exons tend to occupy central positions in interaction networks and display distinct interaction partners in the respective tissues, and are enriched in signaling, development, and disease genes. Based on these findings, we propose that tissue-specific inclusion of disordered segments that contain binding motifs rewires interaction networks and signaling pathways. In this way, tissue-specific splicing may contribute to functional versatility of proteins and increases the diversity of interaction networks across tissues.

The spatial architecture of protein function and adaptation

Statistical analysis of protein evolution suggests a design for natural proteins in which sparse networks of coevolving amino acids (termed sectors) comprise the essence of three-dimensional structure and function. However, proteins are also subject to pressures deriving from the dynamics of the evolutionary process itself--the ability to tolerate mutation and to be adaptive to changing selection pressures. To understand the relationship of the sector architecture to these properties, we developed a high-throughput quantitative method for a comprehensive single-mutation study in which every position is substituted individually to every other amino acid. Using a PDZ domain (PSD95(pdz3)) model system, we show that sector positions are functionally sensitive to mutation, whereas non-sector positions are more tolerant to substitution. In addition, we find that adaptation to a new binding specificity initiates exclusively through variation within sector residues. A combination of just two sector mutations located near and away from the ligand-binding site suffices to switch the binding specificity of PSD95(pdz3) quantitatively towards a class-switching ligand. The localization of functional constraint and adaptive variation within the sector has important implications for understanding and engineering proteins.

Negative coupling as a mechanism for signal propagation between C2 domains of synaptotagmin I

Synaptotagmin I (Syt I) is a vesicle-localized protein implicated in sensing the calcium influx that triggers fast synchronous release of neurotransmitter. How Syt I utilizes its two C2 domains to integrate signals and mediate neurotransmission has continued to be a controversial area of research, though prevalent hypotheses favor independent function. Using differential scanning calorimetry and fluorescence lifetime spectroscopy in a thermodynamic denaturation approach, we tested an alternative hypothesis in which both domains interact to cooperatively disseminate binding information. The free energy of stability was determined for C2A, C2B, and C2AB constructs by globally fitting both methods to a two-state model of unfolding. By comparing the additive free energies of C2A and C2B with C2AB, we identified a negative coupling interaction between the C2 domains of Syt I. This interaction not only provides a mechanistic means for propagating signals, but also a possible means for coordinating the molecular events of neurotransmission.

Thermal unfolding of the N-terminal region of p53 monitored by circular dichroism spectroscopy

It has been estimated that 30% of eukaryotic protein and 70% of transcription factors are intrinsically disordered (ID). The biochemical significance of proteins that lack stable tertiary structure, however, is not clearly understood, largely owing to an inability to assign well-defined structures to specific biological tasks. In an attempt to investigate the structural character of ID protein, we have measured the circular dichroism spectrum of the N-terminal region of p53 over a range of temperatures and solution conditions. p53 is a well-studied transcription factor that has a proline-rich N-terminal ID region containing two activation domains. High proline content is a property commonly associated with ID, and thus p53 may be a good model system for investigating the biochemical importance of ID. The spectra presented here suggest that the N-terminal region of p53 may adopt an ordered structure under physiological conditions and that this structure can be thermally unfolded in an apparent two-state manner. The midpoint temperature for this thermal unfolding of the N-terminal region of p53 was at the near-physiological temperature of 39°C, suggesting the possibility of a physiological role for the observed structural equilibrium.

Conditional disorder in chaperone action

Protein disorder remains an intrinsically fuzzy concept. Its role in protein function is difficult to conceptualize and its experimental study is challenging. Although a wide variety of roles for protein disorder have been proposed, establishing that disorder is functionally important, particularly in vivo, is not a trivial task. Several molecular chaperones have now been identified as conditionally disordered proteins; fully folded and chaperone-inactive under non-stress conditions, they adopt a partially disordered conformation upon exposure to distinct stress conditions. This disorder appears to be vital for their ability to bind multiple aggregation-sensitive client proteins and to protect cells against the stressors. The study of these conditionally disordered chaperones should prove useful in understanding the functional role for protein disorder in molecular recognition.

Intrinsically disordered proteins: a 10-year recap

The suggestion that the native state of many proteins is intrinsically disordered (or, as originally termed, unstructured) is now integral to our general view of protein structure and function. A little more than 10 years ago, however, such challenge to the almost dogmatic 'structure-function paradigm' was pure heresy due to the overwhelming evidence that structure determines function. A decade of steady progress turned skepticism around: this 10-year recap review outlines the situation a decade ago and the major directions of the breathtaking advance achieved by experimental and computational approaches. I show that the evidence for the generality and importance of this phenomenon is now so insurmountable that it demands the inclusion of 'unstructural' biology into mainstream biology and biochemistry textbooks.

Binding of two intrinsically disordered peptides to a multi-specific protein: a combined Monte Carlo and molecular dynamics study

The unique ability of intrinsically disordered proteins (IDPs) to fold upon binding to partner molecules makes them functionally well-suited for cellular communication networks. For example, the folding-binding of different IDP sequences onto the same surface of an ordered protein provides a mechanism for signaling in a many-to-one manner. Here, we study the molecular details of this signaling mechanism by applying both Molecular Dynamics and Monte Carlo methods to S100B, a calcium-modulated homodimeric protein, and two of its IDP targets, p53 and TRTK-12. Despite adopting somewhat different conformations in complex with S100B and showing no apparent sequence similarity, the two IDP targets associate in virtually the same manner. As free chains, both target sequences remain flexible and sample their respective bound, natively -helical states to a small extent. Association occurs through an intermediate state in the periphery of the S100B binding pocket, stabilized by nonnative interactions which are either hydrophobic or electrostatic in nature. Our results highlight the importance of overall physical properties of IDP segments, such as net charge or presence of strongly hydrophobic amino acids, for molecular recognition via coupled folding-binding.

A substantial fraction of our proteins are believed to be partly or completely disordered, meaning that they contain regions that lack a stable folded structure under typical physiological conditions. This is a feature which plays a key role in their functions. For example, it allows them to have many structurally different binding partners which in turn permits the construction of the intricate signaling and regulatory networks necessary to sustain complex biological organisms such as ourselves. Whereas measuring the binding strengths of associations involving disordered proteins is routine, the binding process itself is today still not fully understood. We use two different computational models to study the interactions of a folded protein, S100B, which can bind various disordered peptides. In particular, we compare two peptides whose structures are known when in complex with S100B. Our results suggest that, although the peptides assume different structures in the bound state, there are similarities in how they associate with S100B. The possibility to computationally model the interplay between proteins is an important complement to experiments, by identifying crucial steps in the binding process. This is essential to understand, e.g., how single mutations sometimes lead to serious diseases.

A preformed binding interface in the unbound ensemble of an intrinsically disordered protein: evidence from molecular simulations

Intrinsically disordered proteins play an important role in cellular signalling, mediated by their interactions with other biomolecules. A key question concerns the nature of their binding mechanism, and whether the bound structure is induced only by proximity to the binding partner. This is difficult to answer through experiment alone because of the very heterogeneous nature of the unbound ensemble, and the probable rapid interconversion of the various unbound structures. Here we report the most extensive set of simulations on NCBD to date: we use large-scale replica exchange molecular dynamics to explore the unbound state. An important feature of the study is the use of an atomistic force field that has been parametrised against experimental data for weakly structured peptides, together with an accurate explicit water model. Neither the force field nor the starting conformations are biased towards a particular structure. The regions of NCBD that have high helical propensity in the simulations correspond closely to helices in the 'core' unbound conformation determined by NMR, although no single member of the simulated unbound ensemble closely resembles the core conformation, or either of the two known bound conformations. We have validated the results against NMR spectroscopy and SAXS measurements, obtaining reasonable agreement. The two helices which most stabilise the binding of NCBD with ACTR are formed readily; the third helix, which is less important for binding but is involved in most of the intraprotein contacts of NCBD in the bound conformation, is formed more rarely, and tends not to coexist with the other helices. These results support a mechanism by which NCBD gains the advantages of disorder, while forming binding-competent structures in the unbound state. We obtain support for this mechanism from coarse-grained simulations of NCBD with, and without, its binding partner.

Local motifs in proteins combine to generate global functional moves

Literature on the topological properties of folded proteins that has emerged as a field in its own right in the past decade is reviewed. Physics-based construction of coarse-grained models of proteins from knowledge of all-atom coordinates of the average structure is discussed. Once network is thus obtained with the node and link information, local motifs provide plethora of information on protein function. The hierarchical structure of the proteins manifested in the interrelations of local motifs is emphasized. Motifs are also related to modularity of the structure, and they quantify shifts in the landscapes upon conformational changes induced by, e.g. ligand binding. Redundancy emerges as a balance between local and global network descriptors and is related to the collectivity of the protein motions. Introducing weight on links followed by sequential removal of least cohesive contacts allows interactions in proteins to be represented as the superposition of essential and redundant sets. Lack of the former makes the network non-functional, while the latter ensures robust functioning under a wide range of perturbation scenarios.

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.

Thermodynamic dissection of the intrinsically disordered N-terminal domain of human glucocorticoid receptor

Intrinsically disordered (ID) sequence segments are abundant in cell signaling proteins and transcription factors. Because ID regions commonly fold as part of their intracellular function, it is crucial to understand the folded states as well as the transitions between the unfolded and folded states. Specifically, it is important to determine 1) whether large ID segments contain different thermodynamically and/or functionally distinct regions, 2) whether any ID regions fold upon activation, 3) the degree of coupling between the different ID regions, and 4) whether the stability of ID domains is a determinant of function. In this study, we thermodynamically characterized the full-length ID N-terminal domain (NTD) of human glucocorticoid receptor (GR) and two of its naturally occurring translational isoforms. The protective osmolyte trimethylamine N-oxide (TMAO) was used to induce folding transitions. Each of the three NTD isoforms was found to undergo a cooperative folding transition that is thermodynamically indistinguishable (based on m-values) from that of a globular protein of similar size. The extrapolated stabilities for the NTD isoforms showed clear correlation with the known activities of their corresponding GR translational isoforms. The data reveal that the full-length NTD can be viewed as having at least two thermodynamically coupled regions, a functional region, which is indispensable for GR transcriptional activity, and a regulatory region, the length of which serves to regulate the stability of NTD and thus the activity of GR. These results suggest a new functional paradigm whereby steroid hormone receptors in particular and ID proteins in general can have multiple functionally distinct ID regions that interact and modulate the stability of important functional sites.

Hybrid approaches to molecular simulation

Molecular dynamics (MD) simulation is an established method for studying the conformational changes that are important for protein function. Recent advances in hardware and software have allowed MD simulations over the same timescales as experiment, improving the agreement between theory and experiment to a large extent. However, running such simulations are costly, in terms of resources, storage, and trajectory analysis. There is still a place for techniques that involve short MD simulations. In order to overcome the sampling paucity of short time-scales, hybrid methods that include some form of MD simulation can exploit certain features of the system of interest, often combining experimental information in surprising ways. Here, we review some recent hybrid approaches to the simulation of proteins.

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.

Allosteric modulators of steroid hormone receptors: structural dynamics and gene regulation

Steroid hormones are synthesized from cholesterol primarily in the adrenal gland and the gonads and play vital roles in normal physiology, the control of development, differentiation, metabolic homeostasis, and reproduction. The actions of these small lipophilic molecules are mediated by intracellular receptor proteins. It is just over 25 yr since the first cDNA for steroid receptors were cloned, a development that led to the birth of a superfamily of ligand-activated transcription factors: the nuclear receptors. The receptor proteins share structurally and functionally related ligand binding and DNA-binding domains but possess distinct N-terminal domains and hinge regions that are intrinsically disordered. Since the original cloning experiments, considerable progress has been made in our understanding of the structure, mechanisms of action, and biology of this important class of ligand-activated transcription factors. In recent years, there has been interest in the structural plasticity and function of the N-terminal domain of steroid hormone receptors and in the allosteric regulation of protein folding and function in response to hormone, DNA response element architecture, and coregulatory protein binding partners. The N-terminal domain can exist as an ensemble of conformers, having more or less structure, which prime this region of the receptor to rapidly respond to changes in the intracellular environment through hormone binding and posttranslation modifications. In this review, we address the question of receptor structure and function dynamics with particular emphasis on the structurally flexible N-terminal domain, intra- and interdomain communications, and the allosteric regulation of receptor action.

Understanding and overcoming aminoglycoside resistance caused by N-6'-acetyltransferase

Aminoglycosides occupy a special niche amongst antibiotics in part because of their broad spectrum of action. Bacterial resistance is however menacing to render these drugs obsolete. A significant amount of work has been devoted to understand and overcome aminoglycoside resistance. This mini-review will discuss aminoglycoside-modifying enzymes (AMEs), with a special emphasis on the efforts to comprehend and block resistance caused by aminoglycoside 6'-N-acetyltransferase (AAC(6')).

Agonism/antagonism switching in allosteric ensembles

Ligands for several transcription factors can act as agonists under some conditions and antagonists under others. The structural and molecular bases of such effects are unknown. Previously, we demonstrated how the folding of intrinsically disordered (ID) protein sequences, in particular, and population shifts, in general, could be used to mediate allosteric coupling between different functional domains, a model that has subsequently been validated in several systems. Here it is shown that population redistribution within allosteric systems can be used as a mechanism to tune protein ensembles such that a given ligand can act as both an agonist and an antagonist. Importantly, this mechanism can be robustly encoded in the ensemble, and does not require that the interactions between the ligand and the protein differ when it is acting either as an agonist or an antagonist. Instead, the effect is due to the relative probabilities of states prior to the addition of the ligand. The ensemble view of allostery that is illuminated by these studies suggests that rather than being seen as switches with fixed responses to allosteric activation, ensembles can evolve to be "functionally pluripotent," with the capacity to up or down regulate activity in response to a stimulus. This result not only helps to explain the prevalence of intrinsic disorder in transcription factors and other cell signaling proteins, it provides important insights about the energetic ground rules governing site-to-site communication in all allosteric systems.

Network-based models as tools hinting at nonevident protein functionality

Network-based models of proteins are popular tools employed to determine dynamic features related to the folded structure. They encompass all topological and geometric computational approaches idealizing proteins as directly interacting nodes. Topology makes use of neighborhood information of residues, and geometry includes relative placement of neighbors. Coarse-grained approaches efficiently predict alternative conformations because of inherent collectivity in the protein structure. Such collectivity is moderated by topological characteristics that also tune neighborhood structure: That rich residues have richer neighbors secures robustness toward random loss of interactions/nodes due to environmental fluctuations/mutations. Geometry conveys the additional information of force balance to network models, establishing the local shape of the energy landscape. Here, residue and/or bond perturbations are critically evaluated to suggest new experiments, as network-based computational techniques prove useful in capturing domain movements and conformational shifts resulting from environmental alterations. Evolutionarily conserved residues are optimally connected, defining a subnetwork that may be utilized for further coarsening.

Coupled folding-binding in a hydrophobic/polar protein model: impact of synergistic folding and disordered flanks

Coupled folding-binding is central to the function of many intrinsically disordered proteins, yet not fully understood. With a continuous three-letter protein model, we explore the free-energy landscape of pairs of interacting sequences and how it is impacted by 1), variations in the binding mechanism; and 2), the addition of disordered flanks to the binding region. In particular, we focus on two sequences, one with 16 and one with 35 amino acids, which make a stable dimeric three-helix bundle at low temperatures. Three distinct binding mechanisms are realized by altering the stabilities of the individual monomers: docking, coupled folding-binding of a single α-helix, and synergistic folding and binding. Compared to docking, the free-energy barrier for binding is reduced when the single α-helix is allowed to fold upon binding, but only marginally. A greater reduction is found for synergistic folding, which in addition results in a binding transition state characterized by very few interchain contacts. Disordered flanking chain segments attached to the α-helix sequence can, despite a negligible impact on the dimer stability, lead to a downhill free-energy surface in which the barrier for binding is eliminated.

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.

The road less traveled: new views of steroid receptor action from the path of dose-response curves

Conventional studies of steroid hormone action proceed via quantitation of the maximal activity for gene induction at saturating concentrations of agonist steroid (i.e., A(max)). Less frequently analyzed parameters of receptor-mediated gene expression are EC(50) and PAA. The EC(50) is the concentration of steroid required for half-maximal agonist activity and is readily determined from the dose-response curve. The PAA is the partial agonist activity of an antagonist steroid, expressed as percent of A(max) under the same conditions. Recent results demonstrate that new and otherwise inaccessible mechanistic information is obtained when the EC(50) and/or PAA are examined in addition to the A(max). Specifically, A(max), EC(50), and PAA can be independently regulated, which suggests that novel pathways and factors may preferentially modify the EC(50) and/or PAA with little effect on A(max). Other approaches indicate that the activity of receptor-bound factors can be altered without changing the binding of factors to receptor. Finally, a new theoretical model of steroid hormone action not only permits a mechanistically based definition of factor activity but also allows the positioning of when a factor acts, as opposed to binds, relative to a kinetically defined step. These advances illustrate some of the benefits of expanding the mechanistic studies of steroid hormone action to routinely include EC(50) and PAA.

Intrinsic disorder in cell signaling and gene transcription

Structural disorder, which enables unique modes of action often associated with molecular recognition and folding induced by a partner, is widespread in eukaryotic proteomes. Due to the ensuing advantages, such as specificity without strong binding, adaptability to multiple partners and subtle regulation by post-translational modification, structural disorder is prevalent in proteins of signaling and regulatory functions, such as membrane receptors, scaffold proteins, cytoskeletal proteins, transcription factors and nuclear hormone receptors. In this review we survey the most important aspects of structural disorder, with major focus on features and advantages pertinent to signal transduction. Our major goal is to elucidate how the functional requirements of these protein classes concur with specific functional modes disorder enables.

Structural and functional analysis of domains of the progesterone receptor

Steroid hormone receptors are multi-domain proteins composed of conserved well-structured regions, such as ligand (LBD) and DNA binding domains (DBD), plus other naturally unstructured regions including the amino-terminal domain (NTD) and the hinge region between the LBD and DBD. The hinge is more than just a flexible region between the DBD and LBD and is capable of binding co-regulatory proteins and the minor groove of DNA flanking hormone response elements. Because the hinge can directly participate in DNA binding it has also been termed the carboxyl terminal extension (CTE) of the DNA binding domain. The CTE and NTD are dynamic regions of the receptor that can adopt multiple conformations depending on the environment of interacting proteins and DNA. Both regions have important regulatory roles for multiple receptor functions that are related to the ability of the CTE and NTD to form multiple active conformations. This review focuses on studies of the CTE and NTD of progesterone receptor (PR), as well as related work with other steroid/nuclear receptors.

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.

Van't Hoff global analyses of variable temperature isothermal titration calorimetry data

Isothermal titration calorimetry (ITC) can provide detailed information on the thermodynamics of biomolecular interactions in the form of equilibrium constants, KA , and enthalpy changes, ΔHA . A powerful application of this technique involves analyzing the temperature dependences of ITC-derived KA and ΔHA values to gain insight into thermodynamic linkage between binding and additional equilibria, such as protein folding. We recently developed a general method for global analysis of variable temperature ITC data that significantly improves the accuracy of extracted thermodynamic parameters and requires no prior knowledge of the coupled equilibria. Here we report detailed validation of this method using Monte Carlo simulations and an application to study coupled folding and binding in an aminoglycoside acetyltransferase enzyme.

Roles for intrinsic disorder and fuzziness in generating context-specific function in Ultrabithorax, a Hox transcription factor

Surprisingly few transcription factors drive animal development relative to the number and diversity of final tissues and body structures. Therefore, most transcription factors must function in more than one tissue. In a famous example, members of the Hox transcription factor family are expressed in contiguous stripes along the anterior/posterior axis during animal development. Individual Hox transcription factors specify all tissues within their expression domain and thus must respond to cellular cues to instigate the correct tissue-specific gene regulatory cascade. We describe how, in the Drosophila Hox protein Ultrabithorax, intrinsically disordered regions implement, regulate and co-ordinate multiple functions, potentially enabling context-specific gene regulation. The large N-terminal disordered domain encodes most of the transcription activation domain and directly impacts DNA binding affinity by the Ubx homeodomain. Similarly, the C-terminal disordered domain alters DNA binding affinity and specificity, interaction with a Hox binding protein and strongly influences both transcription activation and repression. Phosphorylation of the N-terminal disordered domain and alternative splicing of the C-terminal disordered domain could allow the cell to both regulate and co-ordinate DNA binding, protein interactions and transcription regulation. For regulatory mechanisms relying on disorder to continue to be available when Ubx is bound to other proteins or DNA, fuzziness would need to be preserved in these macromolecular complexes. The intrinsically disordered domains in Hox proteins are predicted to be on the very dynamic end of the disorder spectrum, potentially allowing disorder to persist when Ubx is bound to proteins or DNA to regulate the function of these "fuzzy" complexes. Because both intrinsically disordered regions within Ubx have multiple roles, each region may implement several different regulatory mechanisms identified in fuzzy complexes. These intrinsic disorder-based regulatory mechanisms are likely to be critical for allowing Ubx to sense tissue identity and respond by implementing a context-specific gene regulatory cascade.

N-terminal segments modulate the α-helical propensities of the intrinsically disordered basic regions of bZIP proteins

Basic region leucine zippers (bZIPs) are modular transcription factors that play key roles in eukaryotic gene regulation. The basic regions of bZIPs (bZIP-bRs) are necessary and sufficient for DNA binding and specificity. Bioinformatic predictions and spectroscopic studies suggest that unbound monomeric bZIP-bRs are uniformly disordered as isolated domains. Here, we test this assumption through a comparative characterization of conformational ensembles for 15 different bZIP-bRs using a combination of atomistic simulations and circular dichroism measurements. We find that bZIP-bRs have quantifiable preferences for α-helical conformations in their unbound monomeric forms. This helicity varies from one bZIP-bR to another despite a significant sequence similarity of the DNA binding motifs (DBMs). Our analysis reveals that intramolecular interactions between DBMs and eight-residue segments directly N-terminal to DBMs are the primary modulators of bZIP-bR helicities. We test the accuracy of this inference by designing chimeras of bZIP-bRs to have either increased or decreased overall helicities. Our results yield quantitative insights regarding the relationship between sequence and the degree of intrinsic disorder within bZIP-bRs, and might have general implications for other intrinsically disordered proteins. Understanding how natural sequence variations lead to modulation of disorder is likely to be important for understanding the evolution of specificity in molecular recognition through intrinsically disordered regions (IDRs).

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).

Structural dynamics, intrinsic disorder, and allostery in nuclear receptors as transcription factors

Steroid hormone receptors (SHRs) and nuclear receptors (NRs) in general are flexible, allosterically regulated transcription factors. The classic model is inadequate to explain all their behavior. Keys to function are their regions of intrinsic disorder (ID). Data show the dynamic structure and allosteric interactions of the three classic SHR domains: ligand-binding (LBD), DNA-binding (DBD), and N-terminal (NTD). Each responds to its ligands by stabilizing its structure. The LBD responds to classic steroidal and nonsteroidal small ligands; both may selectively modify SHR activity. The DBD responds differentially to the DNA sequences of its response elements. The NTD, with its high ID content and AF1, interacts allosterically with the LBD and DBD. Each domain binds heterologous proteins, potential allosteric ligands. An ensemble framework improves the classic model, shows how ID regions poise the SHR/NR family for optimal allosteric response, and provides a basis for quantitative evaluation of SHR/NR actions.

Binding-folding induced regulation of AF1 transactivation domain of the glucocorticoid receptor by a cofactor that binds to its DNA binding domain

Intrinsically disordered (ID) regions of proteins commonly exist within transcription factors, including the N-terminal domain (NTD) of steroid hormone receptors (SHRs) that possesses a powerful activation function, AF1 region. The mechanisms by which SHRs pass signals from a steroid hormone to control gene expression remain a central unresolved problem. The role of N-terminal activation function AF1, which exists in an intrinsically disordered (ID) conformation, in this process is of immense importance. It is hypothesized that under physiological conditions, ID AF1 undergoes disorder/order transition via inter- and intra-molecular communications, which allows AF1 surfaces to interact with specific co-regulatory proteins, critical for the final outcome of target gene expression regulated by SHRs. However, the means by which AF1 acquires functionally folded conformations is not well understood. In this study, we tested whether binding of jun dimerization protein 2 (JDP2) within the DNA binding domain (DBD) of the glucocorticoid receptor (GR) leads to acquisition of functionally active structure in its AF1/NTD. Our results show that signals mediated from GR DBD:JDP2 interactions in a two domain GR fragment, consisting of the entire NTD and little beyond DBD, significantly increased secondary/tertiary structure formation in the NTD/AF1. This increased structure formation facilitated AF1's interaction with specific co-regulatory proteins and subsequent glucocorticoid response element-mediated AF1 promoter:reporter activity. These results support the hypothesis that inter- and intra-molecular signals give a functionally active structure(s) to the GR AF1, which is important for its transcriptional activity.

The C-terminal αI domain linker as a critical structural element in the conformational activation of αI integrins

The activation of α/β heterodimeric integrins is the result of highly coordinated rearrangements within both subunits. The molecular interactions between the two subunits, however, remain to be characterized. In this study, we use the integrin α(L)β(2) to investigate the functional role of the C-linker polypeptide that connects the C-terminal end of the inserted (I) domain with the β-propeller domain on the α subunit and is located at the interface with the βI domain of the β chain. We demonstrate that shortening of the C-linker by eight or more amino acids results in constitutively active α(L)β(2) in which the αI domain is no longer responsive to the regulation by the βI domain. Despite this intersubunit uncoupling, both I domains remain individually sensitive to intrasubunit conformational changes induced by allosteric modulators. Interestingly, the length and not the sequence of the C-linker appears to be critical for its functionality in α/β intersubunit communication. Using two monoclonal antibodies (R7.1 and CBR LFA-1/1) we further demonstrate that shortening of the C-linker results in the gradual loss of combinational epitopes that require both the αI and β-propeller domains for full reactivity. Taken together, our findings highlight the role of the C-linker as a spring-like element that allows relaxation of the αI domain in the resting state and controlled tension of the αI domain during activation, exerted by the β chain.

Intrinsic disorder in the androgen receptor: identification, characterisation and drugability

The androgen receptor (AR) regulates networks of genes in response to the steroid hormones testosterone and dihydrotestosterone. The receptor protein is made up of both stably folded globular domains, involved in hormone and DNA binding, and regions of intrinsic disorder, including the N-terminal domain (NTD). The AR-NTD has a modular activation function (termed AF1) and is important for gene regulation, participating in multiple protein-protein interactions. Biophysical studies have revealed that AR-NTD/AF1 has limited stable secondary structure and conforms to a 'collapsed disordered' conformation. The AR-NTD/AF1 has the propensity to adopt an α-helical conformation in response to a natural osmolyte or a co-regulatory binding partner. The AR is a key drug target in the management of advanced prostate cancer and recently a small molecule inhibitor was identified that interacts with the NTD/AF1 and impairs protein-protein interactions and recruitment of the receptor to target genes. In this review the role of intrinsic disorder in AR function is discussed along with the potential to develop new drugs that will target the structurally plastic NTD.

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.

The native ensemble and folding of a protein molten-globule: functional consequence of downhill folding

The continually emerging functional significance of intrinsic disorder and conformational flexibility in proteins has challenged the long-standing dogma of a well-defined structure contributing to a specific function. Molten-globular states, a class of proteins with significant secondary-structure but a fluid hydrophobic core, is one such example. They have however been difficult to characterize due to the complexity of experimental data and lack of computational avenues. Here, we dissect the folding mechanism of the α-helical molten-globular protein NCBD from three fundamentally different approaches: statistical-mechanical variable barrier model, C(α)-based Gō-model and explicit water all-atom molecular dynamics (MD) simulations. We find that NCBD displays the characteristics of a one-state globally downhill folder but is significantly destabilized. Using simulation techniques, we generate a highly constrained but a heterogeneous native ensemble of the molten-globule for the first time that is consistent with experimental data including small angle X-ray scattering (SAXS), circular dichroism (CD), and nuclear magnetic resonance (NMR). The resulting native ensemble populates conformations reported in other bound-forms providing direct evidence to the mechanism of conformational selection for binding multiple partners in this domain. Importantly, our simulations reveal a connection between downhill folding and large conformational flexibility in this domain that has been evolutionarily selected and functionally exploited resulting in large binding promiscuity. Finally, the multimodel approach we employ here serves as a powerful methodology to study mechanisms and suggests that the thermodynamic features of molten-globules fall within the array of folding mechanisms available to small single-domain proteins.