Remarkably fast coupled folding and binding of the intrinsically disordered transactivation domain of cMyb to CBP KIX

The transition state for coupled folding and binding of a disordered DNA binding domain resembles the unbound state

The basic zippers (bZIPs) are one of two large eukaryotic families of transcription factors whose DNA binding domains are disordered in isolation but fold into stable α-helices upon target DNA binding. Here, we systematically disrupt pre-existing helical propensity within the DNA binding region of the homodimeric bZIP domain of cAMP-response element binding protein (CREB) using Ala-Gly scanning and examine the impact on target binding kinetics. We find that the secondary structure of the transition state strongly resembles that of the unbound state. The residue closest to the dimerization domain is largely folded within both unbound and transition states; dimerization apparently propagates additional helical propensity into the basic region. The results are consistent with electrostatically-enhanced DNA binding, followed by rapid folding from the folded zipper outwards. Fly-casting theory suggests that protein disorder can accelerate binding. Interestingly however, we did not observe higher association rate constants for mutants with lower levels of residual structure in the unbound state.

High-throughput affinity measurements of direct interactions between activation domains and co-activators

Sequence-specific activation by transcription factors is essential for gene regulation1,2. Key to this are activation domains, which often fall within disordered regions of transcription factors3,4 and recruit co-activators to initiate transcription5. These interactions are difficult to characterize via most experimental techniques because they are typically weak and transient6,7. Consequently, we know very little about whether these interactions are promiscuous or specific, the mechanisms of binding, and how these interactions tune the strength of gene activation. To address these questions, we developed a microfluidic platform for expression and purification of hundreds of activation domains in parallel followed by direct measurement of co-activator binding affinities (STAMMPPING, for Simultaneous Trapping of Affinity Measurements via a Microfluidic Protein-Protein INteraction Generator). By applying STAMMPPING to quantify direct interactions between eight co-activators and 204 human activation domains (>1,500 K ds), we provide the first quantitative map of these interactions and reveal 334 novel binding pairs. We find that the metazoan-specific co-activator P300 directly binds >100 activation domains, potentially explaining its widespread recruitment across the genome to influence transcriptional activation. Despite sharing similar molecular properties (e.g. enrichment of negative and hydrophobic residues), activation domains utilize distinct biophysical properties to recruit certain co-activator domains. Co-activator domain affinity and occupancy are well-predicted by analytical models that account for multivalency, and in vitro affinities quantitatively predict activation in cells with an ultrasensitive response. Not only do our results demonstrate the ability to measure affinities between even weak protein-protein interactions in high throughput, but they also provide a necessary resource of over 1,500 activation domain/co-activator affinities which lays the foundation for understanding the molecular basis of transcriptional activation.

Commonly asked questions about transcriptional activation domains

Eukaryotic transcription factors activate gene expression with their DNA-binding domains and activation domains. DNA-binding domains bind the genome by recognizing structurally related DNA sequences; they are structured, conserved, and predictable from protein sequences. Activation domains recruit chromatin modifiers, coactivator complexes, or basal transcriptional machinery via structurally diverse protein-protein interactions. Activation domains and DNA-binding domains have been called independent, modular units, but there are many departures from modularity, including interactions between these regions and overlap in function. Compared to DNA-binding domains, activation domains are poorly understood because they are poorly conserved, intrinsically disordered, and difficult to predict from protein sequences. This review, organized around commonly asked questions, describes recent progress that the field has made in understanding the sequence features that control activation domains and predicting them from sequence.

Binding-and-Folding Recognition of an Intrinsically Disordered Protein Using Online Learning Molecular Dynamics

Intrinsically disordered proteins participate in many biological processes by folding upon binding to other proteins. However, coupled folding and binding processes are not well understood from an atomistic point of view. One of the main questions is whether folding occurs prior to or after binding. Here we use a novel, unbiased, high-throughput adaptive sampling approach to reconstruct the binding and folding between the disordered transactivation domain of c-Myb and the KIX domain of the CREB-binding protein. The reconstructed long-term dynamical process highlights the binding of a short stretch of amino acids on c-Myb as a folded α-helix. Leucine residues, especially Leu298-Leu302, establish initial native contacts that prime the binding and folding of the rest of the peptide, with a mixture of conformational selection on the N-terminal region with an induced fit of the C-terminal.

Evolutionary fine-tuning of residual helix structure in disordered proteins manifests in complex structure and lifetime

Transcription depends on complex networks, where folded hub proteins interact with intrinsically disordered transcription factors undergoing coupled folding and binding. For this, local residual structure, a prototypical feature of intrinsic disorder, is key. Here, we dissect the unexplored functional potential of residual structure by comparing structure, kinetics, and thermodynamics within the model system constituted of the DREB2A transcription factor interacting with the αα-hub RCD1-RST. To maintain biological relevance, we developed an orthogonal evolutionary approach for the design of variants with varying amounts of structure. Biophysical analysis revealed a correlation between the amount of residual helical structure and binding affinity, manifested in altered complex lifetime due to changed dissociation rate constants. It also showed a correlation between helical structure in free and bound DREB2A variants. Overall, this study demonstrated how evolution can balance and fine-tune residual structure to regulate complexes in coupled folding and binding, potentially affecting transcription factor competition.

The Mediator complex as a master regulator of transcription by RNA polymerase II

The Mediator complex, which in humans is 1.4 MDa in size and includes 26 subunits, controls many aspects of RNA polymerase II (Pol II) function. Apart from its size, a defining feature of Mediator is its intrinsic disorder and conformational flexibility, which contributes to its ability to undergo phase separation and to interact with a myriad of regulatory factors. In this Review, we discuss Mediator structure and function, with emphasis on recent cryogenic electron microscopy data of the 4.0-MDa transcription preinitiation complex. We further discuss how Mediator and sequence-specific DNA-binding transcription factors enable enhancer-dependent regulation of Pol II function at distal gene promoters, through the formation of molecular condensates (or transcription hubs) and chromatin loops. Mediator regulation of Pol II reinitiation is also discussed, in the context of transcription bursting. We propose a working model for Mediator function that combines experimental results and theoretical considerations related to enhancer-promoter interactions, which reconciles contradictory data regarding whether enhancer-promoter communication is direct or indirect. We conclude with a discussion of Mediator's potential as a therapeutic target and of future research directions.

Physics of biomolecular recognition and conformational dynamics

Biomolecular recognition usually leads to the formation of binding complexes, often accompanied by large-scale conformational changes. This process is fundamental to biological functions at the molecular and cellular levels. Uncovering the physical mechanisms of biomolecular recognition and quantifying the key biomolecular interactions are vital to understand these functions. The recently developed energy landscape theory has been successful in quantifying recognition processes and revealing the underlying mechanisms. Recent studies have shown that in addition to affinity, specificity is also crucial for biomolecular recognition. The proposed physical concept of intrinsic specificity based on the underlying energy landscape theory provides a practical way to quantify the specificity. Optimization of affinity and specificity can be adopted as a principle to guide the evolution and design of molecular recognition. This approach can also be used in practice for drug discovery using multidimensional screening to identify lead compounds. The energy landscape topography of molecular recognition is important for revealing the underlying flexible binding or binding-folding mechanisms. In this review, we first introduce the energy landscape theory for molecular recognition and then address four critical issues related to biomolecular recognition and conformational dynamics: (1) specificity quantification of molecular recognition; (2) evolution and design in molecular recognition; (3) flexible molecular recognition; (4) chromosome structural dynamics. The results described here and the discussions of the insights gained from the energy landscape topography can provide valuable guidance for further computational and experimental investigations of biomolecular recognition and conformational dynamics.

Microsecond Dynamics During the Binding-induced Folding of an Intrinsically Disordered Protein

Tau is an intrinsically disordered protein implicated in many neurodegenerative diseases. The repeat domain fragment of tau, tau-K18, is known to undergo a disorder to order transition in the presence of lipid micelles and vesicles, in which helices form in each of the repeat domains. Here, the mechanism of helical structure formation, induced by a phospholipid mimetic, sodium dodecyl sulfate (SDS) at sub-micellar concentrations, has been studied using multiple biophysical probes. A study of the conformational dynamics of the disordered state, using photoinduced electron transfer coupled to fluorescence correlation spectroscopy (PET-FCS) has indicated the presence of an intermediate state, I, in equilibrium with the unfolded state, U. The cooperative binding of the ligand (L), SDS, to I has been shown to induce the formation of a compact, helical intermediate (IL₅) within the dead time (∼37 µs) of a continuous flow mixer. Quantitative analysis of the PET-FCS data and the ensemble microsecond kinetic data, suggests that the mechanism of induction of helical structure can be described by a U ↔ I ↔ IL₅ ↔ FL₅ mechanism, in which the final helical state, FL₅, forms from IL₅ with a time constant of 50-200 µs. Finally, it has been shown that the helical conformation is an aggregation-competent state that can directly form amyloid fibrils.

Kinetic Methods of Deducing Binding Mechanisms Involving Intrinsically Disordered Proteins

There are multiple examples of protein-protein interactions involving one intrinsically disordered protein region binding to an ordered protein domain in a coupled binding and folding reaction. Similarly to protein folding studies, much effort has been devoted to understanding the mechanisms of such coupled binding and folding reactions. In this chapter, we describe how kinetics can be used to assess binding mechanisms with focus on fluorescence-monitored stopped-flow experiments. The approach can be applied more generally to any protein interaction with or without a coupled conformational change and to other kinetic techniques. Determining binding mechanisms is a great challenge and while "proving" a mechanism may be futile, it is possible to deduce the simplest scenarios, which are consistent with experimental data.

Shedding Light on the Inhibitory Mechanisms of SARS-CoV-1/CoV-2 Spike Proteins by ACE2-Designed Peptides

Angiotensin-converting enzyme 2 (ACE2) is the host cellular receptor that locks onto the surface spike protein of the 2002 SARS coronavirus (SARS-CoV-1) and of the novel, highly transmissible and deadly 2019 SARS-CoV-2, responsible for the COVID-19 pandemic. One strategy to avoid the virus infection is to design peptides by extracting the human ACE2 peptidase domain α1-helix, which would bind to the coronavirus surface protein, preventing the virus entry into the host cells. The natural α1-helix peptide has a stronger affinity to SARS-CoV-2 than to SARS-CoV-1. Another peptide was designed by joining α1 with the second portion of ACE2 that is far in the peptidase sequence yet grafted in the spike protein interface with ACE2. Previous studies have shown that, among several α1-based peptides, the hybrid peptidic scaffold is the one with the highest/strongest affinity for SARS-CoV-1, which is comparable to the full-length ACE2 affinity. In this work, binding and folding dynamics of the natural and designed ACE2-based peptides were simulated by the well-known coarse-grained structure-based model, with the computed thermodynamic quantities correlating with the experimental binding affinity data. Furthermore, theoretical kinetic analysis of native contact formation revealed the distinction between these processes in the presence of the different binding partners SARS-CoV-1 and SARS-CoV-2 spike domains. Additionally, our results indicate the existence of a two-state folding mechanism for the designed peptide en route to bind to the spike proteins, in contrast to a downhill mechanism for the natural α1-helix peptides. The presented low-cost simulation protocol demonstrated its efficiency in evaluating binding affinities and identifying the mechanisms involved in the neutralization of spike-ACE2 interaction by designed peptides. Finally, the protocol can be used as a computer-based screening of more potent designed peptides by experimentalists searching for new therapeutics against COVID-19.

Measuring and Analyzing Binding Kinetics of Coupled Folding and Binding Reactions Under Pseudo-First-Order Conditions

Many intrinsically disordered proteins (IDPs) adopt a well-defined structure upon binding to their interaction partners. Kinetic characterization is a requirement for the investigation of the dynamics and mechanisms of these folding-upon-binding reactions. Here a protocol is described for the investigation of binding kinetics of bimolecular binding and folding reactions of IDPs to their ligand partner under pseudo-first-order conditions using stopped-flow mixing and fluorescence detection.

Conformational Preferences of an Intrinsically Disordered Protein Domain: A Case Study for Modern Force Fields

Molecular simulations of intrinsically disordered proteins (IDPs) are challenging because they require sampling a very large number of relevant conformations, corresponding to a multitude of shallow minima in a flat free energy landscape. However, in the presence of a binding partner, the free energy landscape of an IDP can be dominated by few deep minima. This characteristic imposes high demands on the accuracy of the force field used to describe the molecular interactions. Here, as a model system for an IDP that is unstructured in solution but folds upon binding to a structured interaction partner, the transactivation domain of c-Myb was studied both in the unbound (free) form and when bound to the KIX domain. Six modern biomolecular force fields were systematically tested and compared in terms of their ability to describe the structural ensemble of the IDP. The protein force field/water model combinations included in this study are AMBER ff99SB-disp with its corresponding water model that was derived from TIP4P-D, CHARMM36m with TIP3P, ff15ipq with SPC/E_b, ff99SB*-ILDNP with TIP3P and TIP4P-D, and FB15 with TIP3P-FB water. Comparing the results from REST2-enhanced sampling simulations with experimental CD spectra and secondary chemical shifts reveals that the ff99SB-disp force field can realistically capture the broad and mildly helical structural ensemble of free c-Myb. The structural ensembles yielded by CHARMM36m, ff99SB*-ILDNP together with TIP4P-D water, and FB15 are also mildly helical; however, each of these force fields can be assigned a specific subset of c-Myb residues for which the simulations could not reproduce the experimental secondary chemical shifts. In addition, microsecond-timescale MD simulations of the KIX/c-Myb complex show that most force fields used preserve a stable helix fold of c-Myb in the complex. Still, all force fields predict a KIX/c-Myb complex interface that differs slightly from the structures provided by NMR because several NOE-derived distances between KIX and c-Myb were exceeded in the simulations. Taken together, the ff99SB-disp force field in the first place but also CHARMM36m, ff99SB*-ILDNP together with TIP4P-D water, and FB15 can be suitable choices for future simulation studies of the coupled folding and binding mechanism of the KIX/c-Myb complex and potentially also other IDPs.

The interaction strength of an intrinsically disordered protein domain with its binding partner is little affected by very different cosolutes

A common feature of intrinsically disordered proteins (IDPs) is a disorder-to-order transition upon binding to other proteins, which has been tied to multiple benefits, including accelerated association rates or binding with low affinity, yet high specificity. Given the balanced equilibrium concentrations of folded and unfolded state of an IDP we asked the question if changes in the chemical environment, such as the presence of osmolytes or crowding agents, have a strong influence on the interaction of an IDP. Here, we demonstrate the impact of cosolutes on the interaction of the intrinsically disordered transcription factor c-Myb and its binding partner, the kinase-inducible interaction domain (KIX) of the CREB-binding protein. Temperature jump relaxation kinetics and microscale thermophoresis were employed in order to quantify the rate constants and the binding affinity of the c-Myb/KIX complex, respectively, in the presence of various cosolutes. We find the binding free energy of the c-Myb/KIX complex only marginally modulated by cosolutes, whereas the enthalpy and entropy of the interaction are very sensitive to the respective solvent conditions. For different cosolutes we observe substantial changes in enthalpy, both favorable and unfavorable, which are going with entropy changes largely compensating the enthalpy effects in each case. These characteristics might reflect a potential mechanism by which c-Myb offsets changes in the physico-chemical environment to maintain a roughly unaltered binding affinity.

Polyelectrolyte interactions enable rapid association and dissociation in high-affinity disordered protein complexes

Highly charged intrinsically disordered proteins can form complexes with very high affinity in which both binding partners fully retain their disorder and dynamics, exemplified by the positively charged linker histone H1.0 and its chaperone, the negatively charged prothymosin α. Their interaction exhibits another surprising feature: The association/dissociation kinetics switch from slow two-state-like exchange at low protein concentrations to fast exchange at higher, physiologically relevant concentrations. Here we show that this change in mechanism can be explained by the formation of transient ternary complexes favored at high protein concentrations that accelerate the exchange between bound and unbound populations by orders of magnitude. Molecular simulations show how the extreme disorder in such polyelectrolyte complexes facilitates (i) diffusion-limited binding, (ii) transient ternary complex formation, and (iii) fast exchange of monomers by competitive substitution, which together enable rapid kinetics. Biological polyelectrolytes thus have the potential to keep regulatory networks highly responsive even for interactions with extremely high affinities.

A disordered encounter complex is central to the yeast Abp1p SH3 domain binding pathway

Protein-protein interactions are involved in a wide range of cellular processes. These interactions often involve intrinsically disordered proteins (IDPs) and protein binding domains. However, the details of IDP binding pathways are hard to characterize using experimental approaches, which can rarely capture intermediate states present at low populations. SH3 domains are common protein interaction domains that typically bind proline-rich disordered segments and are involved in cell signaling, regulation, and assembly. We hypothesized, given the flexibility of SH3 binding peptides, that their binding pathways include multiple steps important for function. Molecular dynamics simulations were used to characterize the steps of binding between the yeast Abp1p SH3 domain (AbpSH3) and a proline-rich IDP, ArkA. Before binding, the N-terminal segment 1 of ArkA is pre-structured and adopts a polyproline II helix, while segment 2 of ArkA (C-terminal) adopts a 310 helix, but is far less structured than segment 1. As segment 2 interacts with AbpSH3, it becomes more structured, but retains flexibility even in the fully engaged state. Binding simulations reveal that ArkA enters a flexible encounter complex before forming the fully engaged bound complex. In the encounter complex, transient nonspecific hydrophobic and long-range electrostatic contacts form between ArkA and the binding surface of SH3. The encounter complex ensemble includes conformations with segment 1 in both the forward and reverse orientation, suggesting that segment 2 may play a role in stabilizing the correct binding orientation. While the encounter complex forms quickly, the slow step of binding is the transition from the disordered encounter ensemble to the fully engaged state. In this transition, ArkA makes specific contacts with AbpSH3 and buries more hydrophobic surface. Simulating the binding between ApbSH3 and ArkA provides insight into the role of encounter complex intermediates and nonnative hydrophobic interactions for other SH3 domains and IDPs in general.

Influence of electrostatic forces on the association kinetics and conformational ensemble of an intrinsically disordered protein

Recent work has revealed that the association of a disordered region of a protein with a folded binding partner can occur as rapidly as association between two folded proteins. This is the case for the phosphatase calcineurin (CaN) and its association with its activator calmodulin. Calmodulin binds to the intrinsically disordered regulatory domain of CaN. Previous studies have shown that electrostatic steering can accelerate the binding of folded proteins with disordered ligands. Given that electrostatic forces are strong determinants of disordered protein ensembles, the relationship between electrostatics, conformational ensembles, and quaternary interactions is unclear. Here, we employ experimental approaches to explore the impact of electrostatic interactions on the association of calmodulin with the disordered regulatory region of CaN. We find that estimated association rate constants of calmodulin with our chosen calmodulin-substrates are within the diffusion-limited regime. The association rates are dependent on the ionic strength, indicating that favorable electrostatic forces increase the rate of association. Further, we show that charged amino acids outside the calmodulin-binding site modulate the binding rate. Conformational ensembles obtained from computer simulations suggest that electrostatic interactions within the regulatory domain might bias the conformational ensemble such that the calmodulin binding region is readily accessible. Given the prevalence of charged residues in disordered protein chains, our findings are likely relevant to many protein-protein interactions.

Charge Interactions Modulate the Encounter Complex Ensemble of Two Differently Charged Disordered Protein Partners of KIX

Disordered proteins play important roles in cell signaling and are frequently involved in protein-protein interactions. They also have a larger proportion of charged and polar residues than their folded counterparts. Here, we developed a structure-based model and applied molecular dynamics simulations to examine the presence and importance of electrostatic interactions in the binding processes of two differently charged intrinsically disordered ligands of the KIX domain of CBP. We observed non-native opposite-charged contacts in the encounter complexes for both ligands with KIX, and this may be a general feature of coupled folding and binding reactions. The ensemble of successful encounter complexes is a diverse set of structures, and in the case of the highly charged ligand, this ensemble was found to be malleable with respect to ionic strength. There are only minor differences between encounter complex ensembles for successful and unsuccessful collisions with no key interactions that appear to make the process far more productive. The energy landscape at this early stage in the process does not appear highly funneled. Strikingly we observed many native interactions that appear to reduce chances of an encounter complex being productive. Instead it appears that collectively non-native electrostatic interactions in the encounter complex increase the likelihood of productivity by holding the proteins together long enough for folding to take place. This mechanism is more effective for the more highly charged ligand.

An Early Association between the α-Helix of the TEAD Binding Domain of YAP and TEAD Drives the Formation of the YAP:TEAD Complex

The Hippo pathway is an evolutionarily conserved signaling pathway that is involved in the control of organ size and development. The TEAD transcription factors are the most downstream elements of the Hippo pathway, and their transcriptional activity is regulated via the interaction with different co-regulators such as YAP. The structure of the YAP:TEAD complex shows that YAP binds to TEAD via two distinct secondary structure elements, an α-helix and an Ω-loop, and site-directed mutagenesis experiments revealed that the Ω-loop is the "hot spot" of this interaction. While much is known about how YAP and TEAD interact with each other, little is known about the mechanism leading to the formation of a complex between these two proteins. Here we combine site-directed mutagenesis with pre-steady-state kinetic measurements to show that the association between these proteins follows an apparent one-step binding mechanism. Furthermore, linear free energy relationships and a Φ analysis suggest that binding-induced folding of the YAP α-helix to TEAD occurs independently of and before formation of the Ω-loop interface. Thus, the binding-induced folding of YAP appears not to conform to the concomitant formation of tertiary structure (nucleation-condensation) usually observed for coupled binding and folding reactions. Our findings demonstrate how a mechanism reminiscent of the classical framework (diffusion-collision) mechanism of protein folding may operate in disorder-to-order transitions involving intrinsically disordered proteins.

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

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

Features of molecular recognition of intrinsically disordered proteins via coupled folding and binding

The sequence-structure-function paradigm of proteins has been revolutionized by the discovery of intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs). In contrast to traditional ordered proteins, IDPs/IDRs are unstructured under physiological conditions. The absence of well-defined three-dimensional structures in the free state of IDPs/IDRs is fundamental to their function. Folding upon binding is an important mode of molecular recognition for IDPs/IDRs. While great efforts have been devoted to investigating the complex structures and binding kinetics and affinities, our knowledge on the binding mechanisms of IDPs/IDRs remains very limited. Here, we review recent advances on the binding mechanisms of IDPs/IDRs. The structures and kinetic parameters of IDPs/IDRs can vary greatly, and the binding mechanisms can be highly dependent on the structural properties of IDPs/IDRs. IDPs/IDRs can employ various combinations of conformational selection and induced fit in a binding process, which can be templated by the target and/or encoded by the IDP/IDR. Further studies should provide deeper insights into the molecular recognition of IDPs/IDRs and enable the rational design of IDP/IDR binding mechanisms in the future.

Coupled Binding and Helix Formation Monitored by Synchrotron-Radiation Circular Dichroism

Intrinsically disordered proteins organize interaction networks in the cell in many regulation and signaling processes. These proteins often gain structure upon binding to their target proteins in multistep reactions involving the formation of both secondary and tertiary structure. To understand the interactions of disordered proteins, we need to understand the mechanisms of these coupled folding and binding reactions. We studied helix formation in the binding of the molten globule-like nuclear coactivator binding domain and the disordered interaction domain from activator of thyroid hormone and retinoid receptors. We demonstrate that helix formation in a rapid binding reaction can be followed by stopped-flow synchrotron-radiation circular dichroism (CD) spectroscopy and describe the design of such a beamline. Fluorescence-monitored binding experiments of activator of thyroid hormone and retinoid receptors and nuclear coactivator binding domain display several kinetic phases, including one concentration-independent phase, which is consistent with an intermediate stabilized at high ionic strength. Time-resolved CD experiments show that almost all helicity is formed upon initial association of the proteins or separated from the encounter complex by only a small energy barrier. Through simulation of mechanistic models, we show that the intermediate observed at high ionic strength likely involves a structural rearrangement with minor overall changes in helicity. Our experiments provide a benchmark for simulations of coupled binding reactions and demonstrate the feasibility of using synchrotron-radiation CD for mechanistic studies of protein-protein interactions.

Tunable order-disorder continuum in protein-DNA interactions

DNA-binding protein domains (DBDs) sample diverse conformations in equilibrium facilitating the search and recognition of specific sites on DNA over millions of energetically degenerate competing sites. We hypothesize that DBDs have co-evolved to sense and exploit the strong electric potential from the array of negatively charged phosphate groups on DNA. We test our hypothesis by employing the intrinsically disordered DBD of cytidine repressor (CytR) as a model system. CytR displays a graded increase in structure, stability and folding rate on increasing the osmolarity of the solution that mimics the non-specific screening by DNA phosphates. Electrostatic calculations and an Ising-like statistical mechanical model predict that CytR exhibits features of an electric potential sensor modulating its dimensions and landscape in a unique distance-dependent manner, while DNA plays the role of a non-specific macromolecular chaperone. Accordingly, CytR binds its natural half-site faster than the diffusion-controlled limit and even random DNA conforming to an electrostatic-steering binding mechanism. Our work unravels for the first time the synergistic features of a natural electrostatic potential sensor, a novel binding mechanism driven by electrostatic frustration and disorder, and the role of DNA in promoting distance-dependent protein structural transitions critical for switching between specific and non-specific DNA-binding modes.

Consideration of Binding Kinetics in the Design of Stapled Peptide Mimics of the Disordered Proteins Eukaryotic Translation Initiation Factor 4E-Binding Protein 1 and Eukaryotic Translation Initiation Factor 4G

Protein disorder plays a crucial role in signal transduction and is key for many cellular processes including transcription, translation, and cell cycle. Within the intrinsically disordered protein interactome, the α-helix is commonly used for binding, which is induced via a disorder-to-order transition. Because the targeting of protein-protein interactions (PPIs) remains an important challenge in medicinal chemistry, efforts have been made to mimic this secondary structure for rational inhibitor design through the use of stapled peptides. Cap-dependent mRNA translation is regulated by two disordered proteins, 4E-BP1 and eIF4G, that inhibit or stimulate the activity of the m7G cap-binding translation initiation factor, eIF4E, respectively. Both use an α-helical motif for eIF4E binding, warranting the investigation of stapled peptide mimics for manipulating eIF4E PPIs. Herein, we describe our efforts toward this goal, resulting in the synthesis of a cell-active stapled peptide for further development in manipulating aberrant cap-dependent translation in human diseases.

Position-, disorder-, and salt-dependent diffusion in binding-coupled-folding of intrinsically disordered proteins

Successful extensions of protein-folding energy landscape theory to intrinsically disordered proteins' (IDPs') binding-coupled-folding transition can enormously simplify this biomolecular process into diffusion along a limited number of reaction coordinates, and the dynamics subsequently is described by Kramers' rate theory. As the critical pre-factor, the diffusion coefficient D has direct implications on the binding kinetics. Here, we employ a structure-based model (SBM) to calculate D in the binding-folding of an IDP prototype. We identify a strong position-dependent D during binding by applying a reaction coordinate that directly measures the fluctuations in a Cartesian configuration space. Using the malleability of the SBM, we modulate the degree of conformational disorder in an isolated IDP and determine complex effects of intrinsic disorder on D varying for different binding stages. Here, D tends to increase with disorder during initial binding but shows a non-monotonic relationship with disorder in terms of a decrease-followed-by-increase in D during the late binding stage. The salt concentration, which correlates with electrostatic interactions via Debye-Hückel theory in our SBM, also modulates D in a stepwise way. The speeding up of diffusion by electrostatic interactions is observed during the formation of the encounter complex at the beginning of binding, while the last diffusive binding dynamics is hindered by non-native salt bridges. Because D describes the diffusive speed locally, which implicitly reflects the roughness of the energy landscape, we are eventually able to portray the binding energy landscape, including that from IDPs' binding, then to binding with partial folding, and finally to rigid docking, as well as that under different environmental salt concentrations. Our theoretical results provide key mechanistic insights into IDPs' binding-folding, which is internally conformation- and externally salt-controlled with respect to diffusion.

Collision-Induced Unfolding Reveals Unique Fingerprints for Remote Protein Interaction Sites in the KIX Regulation Domain

The kinase-inducible domain (KIX) of the transcriptional coactivator CBP binds multiple transcriptional regulators through two allosterically connected sites. Establishing a method for observing activator-specific KIX conformations would facilitate the discovery of drug-like molecules that capture specific conformations and further elucidate how distinct activator-KIX complexes produce differential transcriptional effects. However, the transient and low to moderate affinity interactions between activators and KIX are difficult to capture using traditional biophysical assays. Here, we describe a collision-induced unfolding-based approach that produces unique fingerprints for peptides bound to each of the two available sites within KIX, as well as a third fingerprint for ternary KIX complexes. Furthermore, we evaluate the analytical utility of unfolding fingerprints for KIX complexes using CIUSuite, and conclude by speculating as to the structural origins of the conformational families created from KIX:peptide complexes following collisional activation. Graphical Abstract ᅟ.

Does Ionic Screening Lower Activation Barriers for Conformational Transitions in Proteins?

In this work, we investigated the kinetics of binding of hen egg-white lysozyme with tri- N-acetylglucosamine in aqueous solutions, at two values of pH, 3.2 and 11, as a function of ionic strength, by a stopped-flow method with tryptophyl fluorescence observation of the transients. We analyzed registered reaction progress curves by employing numerical integration of appropriate chemical master equations. We discriminated between several binding models and established that the process observed in experiments follows a two-step mechanism, composed of four elementary stages: diffusional formation of an encounter complex, dissociation of the encounter complex, conformational transition of the encounter complex to the final complex, and the reverse transformation, i.e., from the final complex to the encounter complex. We evaluated rate constants of these elementary stages and determined their dependencies on solution ionic strength. Regardless of solution pH, rate constants of both forward and reverse conformational transitions increase with an increasing ionic strength. This suggests that ionic screening of intramolecular electrostatic interactions may act to lower the activation barrier for conformational transition in proteins.

Stopped-Flow Kinetic Techniques for Studying Binding Reactions of Intrinsically Disordered Proteins

Intrinsically disordered proteins are abundant in signaling processes such as transcription. Suitable binding and unbinding rates of proteins with their partners are critical for allowing them to perform their biological roles. Understanding how these are achieved, and indeed designing strategies for intervening or modulating related biological processes, therefore requires kinetic studies. In this chapter, we describe stopped-flow-based methods for determining association and dissociation rate constants for pairs of macromolecular binding partners. We describe how to select the simplest appropriate model to describe the interaction, and highlight cases where it is possible to distinguish between induced fit and conformational selection binding mechanisms. Finally, we go on to describe methods for examining the role of electrostatic forces in binding processes, and for describing the transition state for binding processes that have folding associated with them.

Using NMR Chemical Shifts to Determine Residue-Specific Secondary Structure Populations for Intrinsically Disordered Proteins

Protein disorder is a pervasive phenomenon in biology and a natural consequence of polymer evolution that facilitates cell signaling by organizing sites for posttranslational modifications and protein-protein interactions into arrays of short linear motifs that can be rearranged by RNA splicing. Disordered proteins are missing the long-range nonpolar interactions that form tertiary structures, but they often contain regions with residual secondary structure that are stabilized by protein binding. NMR spectroscopy is uniquely suited to detect residual secondary structure in a disordered protein and it can provide atomic resolution data on the structure and dynamics of disordered protein interaction sites. Here we describe how backbone chemical shifts are used for assigning residual secondary structure in disordered proteins and discuss some of the tools available for estimating secondary structure populations with a focus on disordered proteins containing different levels of alpha helical secondary structure which are stabilized by protein binding.

Conservation of coactivator engagement mechanism enables small-molecule allosteric modulators

Transcriptional coactivators are a molecular recognition marvel because a single domain within these proteins, the activator binding domain or ABD, interacts with multiple compositionally diverse transcriptional activators. Also remarkable is the structural diversity among ABDs, which range from conformationally dynamic helical motifs to those with a stable core such as a β-barrel. A significant objective is to define conserved properties of ABDs that allow them to interact with disparate activator sequences. The ABD of the coactivator Med25 (activator interaction domain or AcID) is unique in that it contains secondary structural elements that are on both ends of the spectrum: helices and loops that display significant conformational mobility and a seven-stranded β-barrel core that is structurally rigid. Using biophysical approaches, we build a mechanistic model of how AcID forms binary and ternary complexes with three distinct activators; despite its static core, Med25 forms short-lived, conformationally mobile, and structurally distinct complexes with each of the cognate partners. Further, ternary complex formation is facilitated by allosteric communication between binding surfaces on opposing faces of the β-barrel. The model emerging suggests that the conformational shifts and cooperative binding is mediated by a flexible substructure comprised of two dynamic helices and flanking loops, indicating a conserved mechanistic model of activator engagement across ABDs. Targeting a region of this substructure with a small-molecule covalent cochaperone modulates ternary complex formation. Our data support a general strategy for the identification of allosteric small-molecule modulators of ABDs, which are key targets for mechanistic studies as well as therapeutic applications.

Uncoupling the Folding and Binding of an Intrinsically Disordered Protein

The relationship between helical stability and binding affinity was examined for the intrinsically disordered transactivation domain of the myeloblastosis oncoprotein, c-Myb, and its ordered binding partner, KIX. A series of c-Myb mutants was designed to either increase or decrease helical stability without changing the binding interface with KIX. This included a complimentary series of A, G, P, and V mutants at three non-interacting sites. We were able to use the glycine mutants as a reference state and show a strong correlation between binding affinity and helical stability. The intrinsic helicity of c-Myb is 21%, and helicity values of the mutants ranged from 8% to 28%. The c-Myb helix is divided into two conformationally distinct segments. The N-terminal segment, from K291-L301, has an average helicity greater than 60% and the C-terminal segment, from S304-L315, has an average helicity less than 10%. We observed different effects on binding when these two segments were mutated. Mutants in the N-terminal segment that increased helicity had no effect on the binding affinity to KIX, while helix destabilizing glycine and proline mutants reduced binding affinity by more than 1 kcal/mol. Mutants that either increased or decreased helical stability in the C-terminal segment had almost no effect on binding. However, several of the mutants reveal the presence of multiple conformations accessible in the bound state based on changes in enthalpy and linkage analysis of binding free energies. These results may explain the high level of sequence identity (>90%), even at non-interacting sites, for c-Myb homologues.

Disparate binding kinetics by an intrinsically disordered domain enables temporal regulation of transcriptional complex formation

Intrinsically disordered regions are highly represented among mammalian transcription factors, where they often contribute to the formation of multiprotein complexes that regulate gene expression. An example of this occurs with LIM-homeodomain (LIM-HD) proteins in the developing spinal cord. The LIM-HD protein LHX3 and the LIM-HD cofactor LDB1 form a binary complex that gives rise to interneurons, whereas in adjacent cell populations, LHX3 and LDB1 form a rearranged ternary complex with the LIM-HD protein ISL1, resulting in motor neurons. The protein-protein interactions within these complexes are mediated by ordered LIM domains in the LIM-HD proteins and intrinsically disordered LIM interaction domains (LIDs) in LDB1 and ISL1; however, little is known about how the strength or rates of binding contribute to complex assemblies. We have measured the interactions of LIM:LID complexes using FRET-based protein-protein interaction studies and EMSAs and used these data to model population distributions of complexes. The protein-protein interactions within the ternary complexes are much weaker than those in the binary complex, yet surprisingly slow LDB1:ISL1 dissociation kinetics and a substantial increase in DNA binding affinity promote formation of the ternary complex over the binary complex in motor neurons. We have used mutational and protein engineering approaches to show that allostery and modular binding by tandem LIM domains contribute to the LDB1_LID binding kinetics. The data indicate that a single intrinsically disordered region can achieve highly disparate binding kinetics, which may provide a mechanism to regulate the timing of transcriptional complex assembly.

pKID Binds to KIX via an Unstructured Transition State with Nonnative Interactions

Understanding the detailed mechanism of interaction of intrinsically disordered proteins with their partners is crucial to comprehend their functions in signaling and transcription. Through its interaction with KIX, the disordered pKID region of CREB protein is central in the transcription of cAMP responsive genes, including those involved in long-term memory. Numerous simulation studies have investigated these interactions. Combined with experimental results, these can provide valuable and comprehensive understanding of the mechanisms involved. Here, we probe the transition state of this interaction experimentally through analyzing the kinetic effect of mutating both interface and solvent exposed residues in pKID. We show that very few specific interactions between pKID and KIX are required in the initial binding process. Only a small number of weak interactions are formed at the transition state, including nonnative interactions, and most of the folding occurs after the initial binding event. These properties are consistent with computational results and also the majority of experimental studies of intrinsically disordered protein coupled folding and binding in other protein systems, suggesting that these may be common features.

Phosphorylation of the IDP KID Modulates Affinity for KIX by Increasing the Lifetime of the Complex

Intrinsically disordered proteins (IDPs) are known to undergo a range of posttranslational modifications, but by what mechanism do such modifications affect the binding of an IDP to its partner protein? We investigate this question using one such IDP, the kinase inducible domain (KID) of the transcription factor CREB, which interacts with the KIX domain of CREB-binding protein upon phosphorylation. As with many other IDPs, KID undergoes coupled folding and binding to form α-helical structure upon interacting with KIX. This single site phosphorylation plays an important role in the control of transcriptional activation in vivo. Here we show that, contrary to expectation, phosphorylation has no effect on association rates—unphosphorylated KID binds just as rapidly as pKID, the phosphorylated form—but rather, acts by increasing the lifetime of the complex. We propose that by controlling the lifetime of the bound complex of pKID:KIX via altering the dissociation rate, phosphorylation can facilitate effective control of transcription regulation.

Kinetic mechanism of controlled Fab-arm exchange for the formation of bispecific immunoglobulin G1 antibodies

Bispecific antibodies (bsAbs) combine the antigen specificities of two distinct Abs and demonstrate therapeutic promise based on novel mechanisms of action. Among the many platforms for creating bsAbs, controlled Fab-arm exchange (cFAE) has proven useful based on minimal changes to native Ab structure and the simplicity with which bsAbs can be formed from two parental Abs. Despite a published protocol for cFAE and its widespread use in the pharmaceutical industry, the reaction mechanism has not been determined. Knowledge of the mechanism could lead to improved yields of bsAb at faster rates as well as foster adoption of process control. In this work, a combination of Förster resonance energy transfer (FRET), nonreducing SDS-PAGE, and strategic mutation of the Ab hinge region was employed to identify and characterize the individual steps of cFAE. Fluorescence correlation spectroscopy (FCS) was used to determine the affinity of parental (homodimer) and bispecific (heterodimer) interactions within the C_H3 domain, further clarifying the thermodynamic basis for bsAb formation. The result is a clear sequence of events with rate constants that vary with experimental conditions, where dissociation of the K409R parental Ab into half-Ab controls the rate of the reaction.

Affinity of IDPs to their targets is modulated by ion-specific changes in kinetics and residual structure

Intrinsically disordered proteins (IDPs) are characterized by a lack of defined structure. Instead, they populate ensembles of rapidly interconverting conformations with marginal structural stabilities. Changes in solution conditions such as temperature and crowding agents consequently affect IDPs more than their folded counterparts. Here we reveal that the residual structure content of IDPs is modulated both by ionic strength and by the type of ions present in solution. We show that these ion-specific structural changes result in binding affinity shifts of up to sixfold, which happen through alteration of both association and dissociation rates. These effects follow the Hofmeister series, but unlike the well-established effects on the stability of folded proteins, they already occur at low, hypotonic concentrations of salt. We attribute this sensitivity to the marginal stability of IDPs, which could have physiological implications given the role of IDPs in signaling, the asymmetric ion profiles of different cellular compartments, and the role of ions in biology.

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

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

The contribution of intrinsically disordered regions to protein function, cellular complexity, and human disease

In the 1960s, Christian Anfinsen postulated that the unique three-dimensional structure of a protein is determined by its amino acid sequence. This work laid the foundation for the sequence–structure–function paradigm, which states that the sequence of a protein determines its structure, and structure determines function. However, a class of polypeptide segments called intrinsically disordered regions does not conform to this postulate. In this review, I will first describe established and emerging ideas about how disordered regions contribute to protein function. I will then discuss molecular principles by which regulatory mechanisms, such as alternative splicing and asymmetric localization of transcripts that encode disordered regions, can increase the functional versatility of proteins. Finally, I will discuss how disordered regions contribute to human disease and the emergence of cellular complexity during organismal evolution.

Eukaryotic transcription factors: paradigms of protein intrinsic disorder

Gene-specific transcription factors (TFs) are key regulatory components of signaling pathways, controlling, for example, cell growth, development, and stress responses. Their biological functions are determined by their molecular structures, as exemplified by their structured DNA-binding domains targeting specific cis-acting elements in genes, and by the significant lack of fixed tertiary structure in their extensive intrinsically disordered regions. Recent research in protein intrinsic disorder (ID) has changed our understanding of transcriptional activation domains from ‘negative noodles’ to ID regions with function-related, short sequence motifs and molecular recognition features with structural propensities. This review focuses on molecular aspects of TFs, which represent paradigms of ID-related features. Through specific examples, we review how the ID-associated flexibility of TFs enables them to participate in large interactomes, how they use only a few hydrophobic residues, short sequence motifs, prestructured motifs, and coupled folding and binding for their interactions with co-activators, and how their accessibility to post-translational modification affects their interactions. It is furthermore emphasized how classic biochemical concepts like allostery, conformational selection, induced fit, and feedback regulation are undergoing a revival with the appreciation of ID. The review also describes the most recent advances based on computational simulations of ID-based interaction mechanisms and structural analysis of ID in the context of full-length TFs and suggests future directions for research in TF ID.

Floppy but not sloppy: Interaction mechanism of FG-nucleoporins and nuclear transport receptors

The nuclear pore complex (NPC) forms a permeability barrier between the nucleus and the cytoplasm. Molecules that are able to cross this permeability barrier encounter different disordered phenylalanine glycine rich nucleoporins (FG-Nups) that act as a molecular filter and regulate the selective NPC crossing of biomolecules. In this review, we provide a current overview regarding the interaction mechanism between FG-Nups and the carrier molecules that recognize and enable the transport of cargoes through the NPC aiming to understand the general molecular mechanisms that facilitate the nucleocytoplasmic transport.

Impact of the amino acid sequence on the conformation of side chain lactam-bridged octapeptides

Synthetic helical peptides are valuable scaffolds for the development of modulators of protein-protein interactions involving helical motifs. Backbone-to-side chain or side chain-to-side chain constraints have been and still are intensively exploited to stabilize short α-helices. Very often, these constraints have been combined with backbone modifications induced by Cα-tetrasubstituted, β-, or γ-amino acids, which facilitate the α-peptide or α/β/γ-peptide adopting an α-helical conformation. In this work, we investigated the helical character of octapeptides that were cyclized by a Lys-Asp-(i,i + 4)-lactam bridge. We started with two sequences extracted from the helix-loop-helix region of the Id proteins, which are inhibitors of cell differentiation during development and in cancer. Nineteen analogs containing the lactam bridge at different positions and displaying different amino acid core triads (i + 1,2,3) as well as outer residues were prepared by solid-phase methodology. Their conformation in water and water/2,2,2-trifluoroethanol mixtures was investigated by circular dichroism (CD) spectroscopy. The cyclopeptides could be grouped in helix-prone and non-helix-prone structures. Both the amino acid core triad (i + 1,2,3) and the pendant residues positively or negatively affected the formation of a helical structure. Computational studies based on the NMR-derived helical structure of a cyclopeptide containing Aib at position (i + 2) of the triad were generally in agreement with the secondary structure propensity of the cyclopeptides observed by CD spectroscopy. In conclusion, the Lys-Asp-(i,i + 4)-lactam bridge may succeed or fail in the stabilization of short helices, depending on the primary structure. Moreover, computational methods may be valuable tools to discriminate helix-prone from non-helix-prone peptide-based macrolactams. Copyright © 2017 European Peptide Society and John Wiley & Sons, Ltd.

Role of non-native electrostatic interactions in the coupled folding and binding of PUMA with Mcl-1

PUMA, which belongs to the BH3-only protein family, is an intrinsically disordered protein (IDP). It binds to its cellular partner Mcl-1 through its BH3 motif, which folds upon binding into an α helix. We have applied a structure-based coarse-grained model, with an explicit Debye-Hückel charge model, to probe the importance of electrostatic interactions both in the early and the later stages of this model coupled folding and binding process. This model was carefully calibrated with the experimental data on helical content and affinity, and shown to be consistent with previously published experimental data on binding rate changes with respect to ionic strength. We find that intramolecular electrostatic interactions influence the unbound states of PUMA only marginally. Our results further suggest that intermolecular electrostatic interactions, and in particular non-native electrostatic interactions, are involved in formation of the initial encounter complex. We are able to reveal the binding mechanism in more detail than is possible using experimental data alone however, and in particular we uncover the role of non-native electrostatic interactions. We highlight the potential importance of such electrostatic interactions for describing the binding reactions of IDPs. Such approaches could be used to provide predictions for the results of mutational studies.

Mechanistic roles of protein disorder within transcription

Understanding the interactions of proteins involved in transcriptional regulation is critical to describing biological systems because they control the expression profile of the cell. Yet sadly they belong to a less well biophysically characterized subset of proteins; they frequently contain long disordered regions that are highly dynamic. A key question therefore is, why? What functional roles does protein disorder play in transcriptional regulation? Experimental data exemplifying these roles are starting to emerge, with common themes being enabling complexity within networks and quick responses. Most recently a role for disorder in mediating phase transitions of membrane-less organelles has been proposed.

Binding Mechanisms of Intrinsically Disordered Proteins: Theory, Simulation, and Experiment

In recent years, protein science has been revolutionized by the discovery of intrinsically disordered proteins (IDPs). In contrast to the classical paradigm that a given protein sequence corresponds to a defined structure and an associated function, we now know that proteins can be functional in the absence of a stable three-dimensional structure. In many cases, disordered proteins or protein regions become structured, at least locally, upon interacting with their physiological partners. Many, sometimes conflicting, hypotheses have been put forward regarding the interaction mechanisms of IDPs and the potential advantages of disorder for protein-protein interactions. Whether disorder may increase, as proposed, e.g., in the “fly-casting” hypothesis, or decrease binding rates, increase or decrease binding specificity, or what role pre-formed structure might play in interactions involving IDPs (conformational selection vs. induced fit), are subjects of intense debate. Experimentally, these questions remain difficult to address. Here, we review experimental studies of binding mechanisms of IDPs using NMR spectroscopy and transient kinetic techniques, as well as the underlying theoretical concepts and numerical methods that can be applied to describe these interactions at the atomic level. The available literature suggests that the kinetic and thermodynamic parameters characterizing interactions involving IDPs can vary widely and that there may be no single common mechanism that can explain the different binding modes observed experimentally. Rather, disordered proteins appear to make combined use of features such as pre-formed structure and flexibility, depending on the individual system and the functional context.