Quantification of Conformational Entropy Unravels Effect of Disordered Flanking Region in Coupled Folding and Binding

Advances in uncovering the mechanisms of macromolecular conformational entropy

During protein folding, proteins transition from a disordered polymer into a globular structure, markedly decreasing their conformational degrees of freedom, leading to a substantial reduction in entropy. Nonetheless, folded proteins retain substantial entropy as they fluctuate between the conformations that make up their native state. This residual entropy contributes to crucial functions like binding and catalysis, supported by growing evidence primarily from NMR and simulation studies. Here, we propose three major ways that macromolecules use conformational entropy to perform their functions; first, prepaying entropic cost through ordering of the ground state; second, spatially redistributing entropy, in which a decrease in entropy in one area is reciprocated by an increase in entropy elsewhere; third, populating catalytically competent ensembles, in which conformational entropy within the enzymatic scaffold aids in lowering transition state barriers. We also provide our perspective on how solving the current challenge of structurally defining the ensembles encoding conformational entropy will lead to new possibilities for controlling binding, catalysis and allostery.

Protein binding and folding through an evolutionary lens

Protein-protein associations are often mediated by an intrinsically disordered protein region interacting with a folded domain in a coupled binding and folding reaction. Classic physical organic chemistry approaches together with structural biology have shed light on mechanistic aspects of such reactions. Further insight into general principles may be obtained by interpreting the results through an evolutionary lens. This review attempts to provide an overview on how the analysis of binding and folding reactions can benefit from an evolutionary approach, and is aimed at protein scientists without a background in evolution. Evolution constantly reshapes existing proteins by sampling more or less fit variants. Most new variants are weeded out as generations and new species come and go over hundreds to hundreds of millions of years. The huge ongoing genome sequencing efforts have provided us with a snapshot of existing adapted fit-for-purpose protein homologs in thousands of different organisms. Comparison of present-day orthologs and paralogs highlights general principles of the evolution of coupled binding and folding reactions and demonstrate a great potential for evolution to operate on disordered regions and modulate affinity and specificity of the interactions.

Stereochemistry in the disorder-order continuum of protein interactions

Intrinsically disordered proteins can bind via the formation of highly disordered protein complexes without the formation of three-dimensional structure1. Most naturally occurring proteins are levorotatory (L)-that is, made up only of L-amino acids-imprinting molecular structure and communication with stereochemistry2. By contrast, their mirror-image dextrorotatory (D)-amino acids are rare in nature. Whether disordered protein complexes are truly independent of chiral constraints is not clear. Here, to investigate the chiral constraints of disordered protein-protein interactions, we chose as representative examples a set of five interacting protein pairs covering the disorder-order continuum. By observing the natural ligands and their stereochemical mirror images in free and bound states, we found that chirality was inconsequential in a fully disordered complex. However, if the interaction relied on the ligand undergoing extensive coupled folding and binding, correct stereochemistry was essential. Between these extremes, binding could be observed for the D-ligand with a strength that correlated with disorder in the final complex. These findings have important implications for our understanding of the molecular processes that lead to complex formation, the use of D-peptides in drug discovery and the chemistry of protein evolution of the first living entities on Earth.

Positive activation entropy of Bacillus circulans xylanase catalyzed ONPX2 hydrolysis: A mechanistic and engineering study

Transition state (TS) stabilization by enzymes greatly accelerates catalytic reactions. For some enzymes, the TS complex has entropy higher than enzyme substrate (ES) complex. But the origin of favorable entropy remains unclear. In this work, we studied the mechanism of Bacillus Circulans xylanase (BCX) 11 catalyzed o-nitrophenyl β-xylobioside (ONPX2) glycoside hydrolysis. The catalytic reaction exhibits a positive activation entropy, and an increase in ionic strength leads to a decrease in entropy without affecting the activation free energy, indicating that the entropy is predominantly influenced by electrostatic forces. Moreover, NMR measurements of electrostatic attractions within the active site demonstrate a positive entropy, aligning with molecular dynamics (MD) simulations showing that electrostatic interactions contribute to the entropic stabilization of the TS complex. These findings suggest that the positive entropy primarily originates from alterations in electrostatic interactions due to the formation of the oxocarbenium ion at C1 in the TS. Differences of electrostatic interactions in ES and TS modify hydrogen bonding of surrounding residues in the active site which causes their side chain dynamics and thus conformational entropy changes. Residues critical for the positive activation entropy are identified. A new BCX mutant with an increased activation entropy and catalytic activity is found.

The ribosome lowers the entropic penalty of protein folding

Most proteins fold during biosynthesis on the ribosome1, and co-translational folding energetics, pathways and outcomes of many proteins have been found to differ considerably from those in refolding studies2-10. The origin of this folding modulation by the ribosome has remained unknown. Here we have determined atomistic structures of the unfolded state of a model protein on and off the ribosome, which reveal that the ribosome structurally expands the unfolded nascent chain and increases its solvation, resulting in its entropic destabilization relative to the peptide chain in isolation. Quantitative 19F NMR experiments confirm that this destabilization reduces the entropic penalty of folding by up to 30 kcal mol-1 and promotes formation of partially folded intermediates on the ribosome, an observation that extends to other protein domains and is obligate for some proteins to acquire their active conformation. The thermodynamic effects also contribute to the ribosome protecting the nascent chain from mutation-induced unfolding, which suggests a crucial role of the ribosome in supporting protein evolution. By correlating nascent chain structure and dynamics to their folding energetics and post-translational outcomes, our findings establish the physical basis of the distinct thermodynamics of co-translational protein folding.

Quantification of thermodynamic effects of carbohydrate multivalency on avidity using synthetic discrete glycooligomers

A quantitative understanding of thermodynamic effects of avidity in biomolecular interactions is important. Herein, we synthesized discrete glycooligomers and evaluated their interactions with a model protein using isothermal titration calorimetry. The dimeric glycooligomer exhibited higher binding constants compared to the glycomonomer, attributed to the reduced conformational entropy loss through local presentation of multiple carbohydrate units. Conversely, divalent glycoligands with polyethylene glycol linkers, aiming for multivalent binding, showed enhanced interactions through increased enthalpy. These findings emphasize the importance of distinguishing between the "local avidity" and the "multipoint avidity".

Structural Adaptation of Secondary p53 Binding Sites on MDM2 and MDMX

The thermodynamics of secondary p53 binding sites on MDM2 and MDMX were evaluated using p53 peptides containing residues 16-29, 17-35, and 1-73. All the peptides had large, negative heat capacity (ΔCp), consistent with the burial of p53 residues F19, W23, and L26 in the primary binding sites of MDM2 and MDMX. MDMX has a higher affinity and more negative ΔCp than MDM2 for p5317-35, which is due to MDMX stabilization and not additional interactions with the secondary binding site. ΔCp measurements show binding to the secondary site is inhibited by the disordered tails of MDM2 for WT p53 but not a more helical mutant where proline 27 is changed to alanine. This result is supported by all-atom molecular dynamics simulations showing that p53 residues 30-35 turn away from the disordered tails of MDM2 in P27A17-35 and make direct contact with this region in p5317-35. Molecular dynamics simulations also suggest that an intramolecular methionine-aromatic motif found in both MDM2 and MDMX structurally adapts to support multiple p53 binding modes with the secondary site. ΔCp measurements also show that tighter binding of the P27A mutant to MDM2 and MDMX is due to increased helicity, which reduces the energetic penalty associated with coupled folding and binding. Our results will facilitate the design of selective p53 inhibitors for MDM2 and MDMX.

Assessment of machine-learning predictions for the Mediator complex subunit MED25 ACID domain interactions with transactivation domains

The human Mediator complex subunit MED25 binds transactivation domains (TADs) present in various cellular and viral proteins using two binding interfaces, named H1 and H2, which are found on opposite sides of its ACID domain. Here, we use and compare deep learning methods to characterize human MED25-TAD interfaces and assess the predicted models to published experimental data. For the H1 interface, AlphaFold produces predictions with high-reliability scores that agree well with experimental data, while the H2 interface predictions appear inconsistent, preventing reliable binding modes. Despite these limitations, we experimentally assess the validity of MED25 interface predictions with the viral transcriptional activators Lana-1 and IE62. AlphaFold predictions also suggest the existence of a unique hydrophobic pocket for the Arabidopsis MED25 ACID domain.

Leucine Motifs Stabilize Residual Helical Structure in Disordered Proteins

Many examples are known of regions of intrinsically disordered proteins that fold into α-helices upon binding to their targets. These helical binding motifs (HBMs) can be partially helical also in the unbound state, and this so-called residual structure can affect binding affinity and kinetics. To investigate the underlying mechanisms governing the formation of residual helical structure, we assembled a dataset of experimental helix contents of 65 peptides containing HBM that fold-upon-binding. The average residual helicity is 17% and increases to 60% upon target binding. The helix contents of residual and target-bound structures do not correlate, however the relative location of helix elements in both states shows a strong overlap. Compared to the general disordered regions, HBMs are enriched in amino acids with high helix preference and these residues are typically involved in target binding, explaining the overlap in helix positions. In particular, we find that leucine residues and leucine motifs in HBMs are the major contributors to helix stabilization and target-binding. For the two model peptides, we show that substitution of leucine motifs to other hydrophobic residues (valine or isoleucine) leads to reduction of residual helicity, supporting the role of leucine as helix stabilizer. From the three hydrophobic residues only leucine can efficiently stabilize residual helical structure. We suggest that the high occurrence of leucine motifs and a general preference for leucine at binding interfaces in HBMs can be explained by its unique ability to stabilize helical elements.

Conservation of Affinity Rather Than Sequence Underlies a Dynamic Evolution of the Motif-Mediated p53/MDM2 Interaction in Ray-Finned Fishes

The transcription factor and cell cycle regulator p53 is marked for degradation by the ubiquitin ligase MDM2. The interaction between these 2 proteins is mediated by a conserved binding motif in the disordered p53 transactivation domain (p53TAD) and the folded SWIB domain in MDM2. The conserved motif in p53TAD from zebrafish displays a 20-fold weaker interaction with MDM2, compared to the interaction in human and chicken. To investigate this apparent difference, we tracked the molecular evolution of the p53TAD/MDM2 interaction among ray-finned fishes (Actinopterygii), the largest vertebrate clade. Intriguingly, phylogenetic analyses, ancestral sequence reconstructions, and binding experiments showed that different loss-of-affinity changes in the canonical binding motif within p53TAD have occurred repeatedly and convergently in different fish lineages, resulting in relatively low extant affinities (KD = 0.5 to 5 μM). However, for 11 different fish p53TAD/MDM2 interactions, nonconserved regions flanking the canonical motif increased the affinity 4- to 73-fold to be on par with the human interaction. Our findings suggest that compensating changes at conserved and nonconserved positions within the motif, as well as in flanking regions of low conservation, underlie a stabilizing selection of "functional affinity" in the p53TAD/MDM2 interaction. Such interplay complicates bioinformatic prediction of binding and calls for experimental validation. Motif-mediated protein-protein interactions involving short binding motifs and folded interaction domains are very common across multicellular life. It is likely that the evolution of affinity in motif-mediated interactions often involves an interplay between specific interactions made by conserved motif residues and nonspecific interactions by nonconserved disordered regions.

Multivalency emerges as a common feature of intrinsically disordered protein interactions

Intrinsically disordered proteins (IDPs) use their unique molecular properties and conformational plasticity to interact with cellular partners in a wide variety of biological contexts. Multivalency is an important feature of IDPs that allows for utilization of an expanded toolkit for interactions with other macromolecules and confers additional complexity to molecular recognition processes. Recent studies have offered insights into how multivalent interactions of IDPs enable responsive and sensitive regulation in the context of transcription and cellular signaling. Multivalency is also widely recognized as an important feature of IDP interactions that mediate formation of biomolecular condensates. We highlight recent examples of multivalent interactions of IDPs across diverse contexts to illustrate the breadth of biological processes that utilize multivalency in molecular interactions.

Molecular switching in transcription through splicing and proline-isomerization regulates stress responses in plants

The Arabidopsis thaliana DREB2A transcription factor interacts with the negative regulator RCD1 and the ACID domain of subunit 25 of the transcriptional co-regulator mediator (Med25) to integrate stress signals for gene expression, with elusive molecular interplay. Using biophysical and structural analyses together with high-throughput screening, we reveal a bivalent binding switch in DREB2A containing an ACID-binding motif (ABS) and the known RCD1-binding motif (RIM). The RIM is lacking in a stress-induced DREB2A splice variant with retained transcriptional activity. ABS and RIM bind to separate sites on Med25-ACID, and NMR analyses show a structurally heterogeneous complex deriving from a DREB2A-ABS proline residue populating cis- and trans-isomers with remote impact on the RIM. The cis-isomer stabilizes an α-helix, while the trans-isomer may introduce energetic frustration facilitating rapid exchange between activators and repressors. Thus, DREB2A uses a post-transcriptionally and post-translationally modulated switch for transcriptional regulation.

Conformational entropy in molecular recognition of intrinsically disordered proteins

Broad conformational ensembles make intrinsically disordered proteins or regions entropically intriguing. Although methodologically challenging and understudied, emerging studies into their changes in conformational entropy (ΔS°conf) upon complex formation have provided both quantitative and qualitative insight. Recent work based on thermodynamics from isothermal titration calorimetry and NMR spectroscopy uncovers an expanded repertoire of regulatory mechanisms, where ΔS°conf plays roles in partner selection, state behavior, functional buffering, allosteric regulation, and drug design. We highlight these mechanisms to display the large entropic reservoir of IDPs for the regulation of molecular communication. We call upon the field to make efforts to contribute to this insight as more studies are needed for forwarding mechanistic decoding of intrinsically disordered proteins and their complexes.

Driving forces of the complex formation between highly charged disordered proteins

Highly disordered complexes between oppositely charged intrinsically disordered proteins present a new paradigm of biomolecular interactions. Here, we investigate the driving forces of such interactions for the example of the highly positively charged linker histone H1 and its highly negatively charged chaperone, prothymosin α (ProTα). Temperature-dependent single-molecule Förster resonance energy transfer (FRET) experiments and isothermal titration calorimetry reveal ProTα-H1 binding to be enthalpically unfavorable, and salt-dependent affinity measurements suggest counterion release entropy to be an important thermodynamic driving force. Using single-molecule FRET, we also identify ternary complexes between ProTα and H1 in addition to the heterodimer at equilibrium and show how they contribute to the thermodynamics observed in ensemble experiments. Finally, we explain the observed thermodynamics quantitatively with a mean-field polyelectrolyte theory that treats counterion release explicitly. ProTα-H1 complex formation resembles the interactions between synthetic polyelectrolytes, and the underlying principles are likely to be of broad relevance for interactions between charged biomolecules in general.

Intrinsically Disordered Flanking Regions Increase the Affinity of a Transcriptional Coactivator Interaction across Vertebrates

Interactions between two proteins are often mediated by a disordered region in one protein binding to a groove in a folded interaction domain in the other one. While the main determinants of a certain interaction are typically found within a well-defined binding interface involving the groove, recent studies show that nonspecific contacts by flanking regions may increase the affinity. One example is the coupled binding and folding underlying the interaction between the two transcriptional coactivators NCOA3 (ACTR) and CBP, where the flanking regions of an intrinsically disordered region in human NCOA3 increases the affinity for CBP. However, it is not clear whether this flanking region-mediated effect is a peculiarity of this single protein interaction or if it is of functional relevance in a broader context. To further assess the role of flanking regions in the interaction between NCOA3 and CBP, we analyzed the interaction across orthologs and paralogs (NCOA1, 2, and 3) in human, zebra fish, and ghost shark. We found that flanking regions increased the affinity 2- to 9-fold in the six interactions tested. Conservation of the amino acid sequence is a strong indicator of function. Analogously, the observed conservation of increased affinity provided by flanking regions, accompanied by moderate sequence conservation, suggests that flanking regions may be under selection to promote the affinity between NCOA transcriptional coregulators and CBP.

Interpreting Transient Interactions of Intrinsically Disordered Proteins

The flexible nature of intrinsically disordered proteins (IDPs) gives rise to a conformational ensemble with a diverse set of conformations. The simplest way to describe this ensemble is through a homopolymer model without any specific interactions. However, there has been growing evidence that the conformational properties of IDPs and their relevant functions can be affected by transient interactions between specific and even nonlocal pairs of amino acids. Interpreting these interactions from experimental methods, each of which is most sensitive to a different distance regime referred to as probing length, remains a challenging and unsolved problem. Here, we first show that transient interactions can be realized between short fragments of charged amino acids by generating conformational ensembles using model disordered peptides and coarse-grained simulations. Using these ensembles, we investigate how sensitive different types of experimental measurements are to the presence of transient interactions. We find methods with shorter probing lengths to be more appropriate for detecting these transient interactions, but one experimental method is not sufficient due to the existence of other weak interactions typically seen in IDPs. Finally, we develop an adjusted polymer model with an additional short-distance peak which can robustly reproduce the distance distribution function from two experimental measurements with complementary short and long probing lengths. This new model can suggest whether a homopolymer model is insufficient for describing a specific IDP and meets the challenge of quantitatively identifying specific, transient interactions from a background of nonspecific, weak interactions.

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.

Computationally Designed Anti-LuxP DNA Aptamer Suppressed Flagellar Assembly- and Quorum Sensing-Related Gene Expression in Vibrio parahaemolyticus

(1) Background: Quorum sensing (QS) is the chemical communication between bacteria that sense chemical signals in the bacterial population to control phenotypic changes through the regulation of gene expression. The inhibition of QS has various potential applications, particularly in the prevention of bacterial infection. QS can be inhibited by targeting the LuxP, a periplasmic receptor protein that is involved in the sensing of the QS signaling molecule known as the autoinducer 2 (AI-2). The sensing of AI-2 by LuxP transduces the chemical information through the inner membrane sensor kinase LuxQ protein and activates the QS cascade. (2) Methods: An in silico approach was applied to design DNA aptamers against LuxP in this study. A method combining molecular docking and molecular dynamics simulations was used to select the oligonucleotides that bind to LuxP, which were then further characterized using isothermal titration calorimetry. Subsequently, the bioactivity of the selected aptamer was examined through comparative transcriptome analysis. (3) Results: Two aptamer candidates were identified from the ITC, which have the lowest dissociation constants (Kd) of 0.2 and 0.5 micromolar. The aptamer with the lowest Kd demonstrated QS suppression and down-regulated the flagellar-assembly-related gene expression. (4) Conclusions: This study developed an in silico approach to design an aptamer that possesses anti-QS properties.

Intrinsic and extrinsic regulators of Aux/IAA protein degradation dynamics

The plant hormone auxin acts through regulated degradation of Auxin/INDOLE-3-ACETIC ACID (Aux/IAA) proteins to regulate transcriptional events. In this review, we examine the composition and function of each Aux/IAA structural motif. We then focus on recent characterization of Aux/IAA N-terminal disordered regions, formation of secondary structure within these disordered regions, and post-translational modifications (PTMs) that affect Aux/IAA function and stability. We propose how structural variations between Aux/IAA family members may be tuned for differential transcriptional repression and degradation dynamics.

Microsecond MD Simulations of the Plexin-B1 RBD: N-H Probability Density as Descriptor of Structural Dynamics, Dimerization-Related Conformational Entropy, and Transient Dimer Asymmetry

Amide-bond equilibrium probability density, P_eq = exp(-u) (u, local potential), and associated conformational entropy, S_k = -∫P_eq (ln P_eq) dΩ ─ln ∫dΩ, are derived for the Rho GTPase binding domain of Plexin-B1 (RBD) as monomer and dimer from 1 μs MD simulations. The objective is to elucidate the effect of dimerization on the dynamic structure of the RBD. Dispersed (peaked) P_eq functions indicate "flexibility" ("rigidity"; the respective concepts are used below in this context). The L1 and L3 loops are throughout highly flexible, the L2 loop and the secondary structure elements are generally rigid, and the L4 loop is flexible in the monomer and rigid in the dimer. Overall, many residues are more flexible in the dimer. These features, and their implications, are discussed. Unexpectedly, we find that monomer unit 1 of the dimer (in short, d1) is unusually flexible, whereas monomer unit 2 (in short, d2) is as rigid as the RBD monomer. This is revealed due to their engagement in slow-to-intermediate conformational exchange detected previously by ¹⁵N relaxation experiments. Such motions occur with rates on the order of 10³-10⁴ s^-1; hence, they cannot be completely sampled over the course of 1 μs simulation. However, the extent to which rigid d2 is affected is small enough to enable physically relevant analysis. The entropy difference between d2 and the monomer yields an entropic contribution of -7 ± 0.7 kJ/mol to the free energy of RBD dimerization. In previous work aimed at similar objectives we used 50-100 ns MD simulations. Those results and the present result differ considerably. In summary, bond-vector P_eq functions derived directly from long MD simulations are useful descriptors of protein structural dynamics and provide accurate conformational entropy. Within the scope of slow conformational exchange, they can be useful, even in the presence of incomplete sampling.

Conformational buffering underlies functional selection in intrinsically disordered protein regions

Many disordered proteins conserve essential functions in the face of extensive sequence variation, making it challenging to identify the mechanisms responsible for functional selection. Here we identify the molecular mechanism of functional selection for the disordered adenovirus early gene 1A (E1A) protein. E1A competes with host factors to bind the retinoblastoma (Rb) protein, subverting cell cycle regulation. We show that two binding motifs tethered by a hypervariable disordered linker drive picomolar affinity Rb binding and host factor displacement. Compensatory changes in amino acid sequence composition and sequence length lead to conservation of optimal tethering across a large family of E1A linkers. We refer to this compensatory mechanism as conformational buffering. We also detect coevolution of the motifs and linker, which can preserve or eliminate the tethering mechanism. Conformational buffering and motif-linker coevolution explain robust functional encoding within hypervariable disordered linkers and could underlie functional selection of many disordered protein regions.

The free energy folding penalty accompanying binding of intrinsically disordered α-helical motifs

Intrinsically disordered proteins (IDPs) are abundant in eukaryotic proteomes and preform critical roles in many cellular processes, most often through the association with globular proteins. Despite lacking a stable three-dimensional structure by themselves, they may acquire a defined conformation upon binding globular targets. The most common type of secondary structure acquired by these binding motifs entails formation of an α-helix. It has been hypothesized that such disorder-to-order transitions are associated with a significant free energy penalty due to IDP folding, which reduces the overall IDP-target affinity. However, the exact magnitude of IDP folding penalty in α-helical binding motifs has not been systematically estimated. Here, we report the folding penalty contributions for 30 IDPs undergoing folding-upon-binding and find that the average IDP folding penalty is +2.0 kcal/mol and ranges from 0.7 to 3.5 kcal/mol. We observe that the folding penalty scales approximately linearly with the change in IDP helicity upon binding, which provides a simple empirical way to estimate folding penalty. We analyze to what extent do pre-structuring and target-bound IDP dynamics (fuzziness) reduce the folding penalty and find that these effects combined, on average, reduce the folding cost by around half. Taken together, the presented analysis provides a quantitative basis for understanding the role of folding penalty in IDP-target interactions and introduces a method estimate this quantity. Estimation and reduction of IDP folding penalty may prove useful in the rational design of helix-stabilized inhibitors of IDP-target interactions. STATEMENT: The α-helical binding motifs are ubiquitous among the intrinsically disordered proteins (IDPs). Upon binding their targets, they undergo a disorder-to-order transition, which is accompanied by a significant folding penalty whose magnitude is generally not known. Here, we use recently developed statistical-thermodynamic model to estimate the folding penalties for 30 IDPs and clarify the roles of IDP pre-folding and bound-state dynamics in reducing the folding penalty.

αα-hub coregulator structure and flexibility determine transcription factor binding and selection in regulatory interactomes

Formation of transcription factor (TF)-coregulator complexes is a key step in transcriptional regulation, with coregulators having essential functions as hub nodes in molecular networks. How specificity and selectivity are maintained in these nodes remain open questions. In this work, we addressed specificity in transcriptional networks using complexes formed between TFs and αα-hubs, which are defined by a common αα-hairpin secondary structure motif, as a model. Using NMR spectroscopy and binding thermodynamics, we analyzed the structure, dynamics, stability, and ligand-binding properties of the Arabidopsis thaliana RST domains from TAF4 and known binding partner RCD1, and the TAFH domain from human TAF4, allowing comparison across species, functions, and architectural contexts. While these αα-hubs shared the αα-hairpin motif, they differed in length and orientation of accessory helices as well as in their thermodynamic profiles of ligand binding. Whereas biologically relevant RCD1-ligand pairs displayed high affinity driven by enthalpy, TAF4-ligand interactions were entropy driven and exhibited less binding-induced structuring. We in addition identified a thermal unfolding state with a structured core for all three domains, although the temperature sensitivity differed. Thermal stability studies suggested that initial unfolding of the RCD1-RST domain localized around helix 1, lending this region structural malleability, while effects in TAF4-RST were more stochastic, suggesting variability in structural adaptability upon binding. Collectively, our results support a model in which hub structure, flexibility, and binding thermodynamics contribute to αα-hub-TF binding specificity, a finding of general relevance to the understanding of coregulator-ligand interactions and interactome sizes.

Drugging Fuzzy Complexes in Transcription

Transcription factors (TFs) are one of the most promising but underutilized classes of drug targets. The high degree of intrinsic disorder in both the structure and the interactions (i.e., “fuzziness”) of TFs is one of the most important challenges to be addressed in this context. Here, we discuss the impacts of fuzziness on transcription factor drug discovery, describing how disorder poses fundamental problems to the typical drug design, and screening approaches used for other classes of proteins such as receptors or enzymes. We then speculate on ways modern biophysical and chemical biology approaches could synergize to overcome many of these challenges by directly addressing the challenges imposed by TF disorder and fuzziness.