De Novo Evolutionary Emergence of a Symmetrical Protein Is Shaped by Folding Constraints.

Molecular evolution has focused on the divergence of molecular functions, yet we know little about how structurally distinct protein folds emerge de novo. We characterized the evolutionary trajectories and selection forces underlying emergence of ß-propeller proteins, a globular and symmetric fold group with diverse functions. The identification of short propeller-like motifs (<50 amino acids) in natural genomes indicated that they expanded via tandem duplications to form extant propellers. We phylogenetically reconstructed 47-residue ancestral motifs that form five-bladed lectin propellers via oligomeric assembly. We demonstrate a functional trajectory of tandem duplications of these motifs leading to monomeric lectins. Foldability, i.e., higher efficiency of folding, was the main parameter leading to improved functionality along the entire evolutionary trajectory. However, folding constraints changed along the trajectory: initially, conflicts between monomer folding and oligomer assembly dominated, whereas subsequently, upon tandem duplication, tradeoffs between monomer stability and foldability took precedence.

Plasticity of an ultrafast interaction between nucleoporins and nuclear transport receptors.

The mechanisms by which intrinsically disordered proteins engage in rapid and highly selective binding is a subject of considerable interest and represents a central paradigm to nuclear pore complex (NPC) function, where nuclear transport receptors (NTRs) move through the NPC by binding disordered phenylalanine-glycine-rich nucleoporins (FG-Nups). Combining single-molecule fluorescence, molecular simulations, and nuclear magnetic resonance, we show that a rapidly fluctuating FG-Nup populates an ensemble of conformations that are prone to bind NTRs with near diffusion-limited on rates, as shown by stopped-flow kinetic measurements. This is achieved using multiple, minimalistic, low-affinity binding motifs that are in rapid exchange when engaging with the NTR, allowing the FG-Nup to maintain an unexpectedly high plasticity in its bound state. We propose that these exceptional physical characteristics enable a rapid and specific transport mechanism in the physiological context, a notion supported by single molecule in-cell assays on intact NPCs.

Disorder in a two-domain neuronal Ca2+-binding protein regulates domain stability and dynamics using ligand mimicry.

Understanding the interplay between sequence, structure and function of proteins has been complicated in recent years by the discovery of intrinsically disordered proteins (IDPs), which perform biological functions in the absence of a well-defined three-dimensional fold. Disordered protein sequences account for roughly 30% of the human proteome and in many proteins, disordered and ordered domains coexist. However, few studies have assessed how either feature affects the properties of the other. In this study, we examine the role of a disordered tail in the overall properties of the two-domain, calcium-sensing protein neuronal calcium sensor 1 (NCS-1). We show that loss of just six of the 190 residues at the flexible C-terminus is sufficient to severely affect stability, dynamics, and folding behavior of both ordered domains. We identify specific hydrophobic contacts mediated by the disordered tail that may be responsible for stabilizing the distal N-terminal domain. Moreover, sequence analyses indicate the presence of an LSL-motif in the tail that acts as a mimic of native ligands critical to the observed order-disorder communication. Removing the disordered tail leads to a shorter life-time of the ligand-bound complex likely originating from the observed destabilization. This close relationship between order and disorder may have important implications for how investigations into mixed systems are designed and opens up a novel avenue of drug targeting exploiting this type of behavior.

Extrinsic conditions influence the self-association and structure of IF1, the regulatory protein of mitochondrial ATP synthase.

The endogenous inhibitor of ATP synthase in mitochondria, called IF1, conserves cellular energy when the proton-motive force collapses by inhibiting ATP hydrolysis. Around neutrality, the 84-amino-acid bovine IF1 is thought to self-assemble into active dimers and, under alkaline conditions, into inactive tetramers and higher oligomers. Dimerization is mediated by formation of an antiparallel α-helical coiled-coil involving residues 44-84. The inhibitory region of each monomer from residues 1-46 is largely α-helical in crystals, but disordered in solution. The formation of the inhibited enzyme complex requires the hydrolysis of two ATP molecules, and in the complex the disordered region from residues 8-13 is extended and is followed by an α-helix from residues 14-18 and a longer α-helix from residue 21, which continues unbroken into the coiled-coil region. From residues 21-46, the long α-helix binds to other α-helices in the C-terminal region of predominantly one of the ß-subunits in the most closed of the three catalytic interfaces. The definition of the factors that influence the self-association of IF1 is a key to understanding the regulation of its inhibitory properties. Therefore, we investigated the influence of pH and salt-types on the self-association of bovine IF1 and the folding of its unfolded region. We identified the equilibrium between dimers and tetramers as a potential central factor in the in vivo modulation of the inhibitory activity and suggest that the intrinsically disordered region makes its inhibitory potency exquisitely sensitive and responsive to physiological changes that influence the capability of mitochondria to make ATP.

The folding and unfolding behavior of ribonuclease H on the ribosome.

The health of a cell depends on accurate translation and proper protein folding, whereas misfolding can lead to aggregation and disease. The first opportunity for a protein to fold occurs during translation, when the ribosome and surrounding environment can affect the nascent chain energy landscape. However, quantifying these environmental effects is challenging because ribosomal proteins and rRNA preclude most spectroscopic measurements of protein energetics. Here, we have applied two gel-based approaches, pulse proteolysis and force-profile analysis, to probe the folding and unfolding pathways of RNase H (RNH) nascent chains stalled on the prokaryotic ribosome in vitro We found that ribosome-stalled RNH has an increased unfolding rate compared with free RNH. Because protein stability is related to the ratio of the unfolding and folding rates, this increase completely accounts for the observed change in protein stability and indicates that the folding rate is unchanged. Using arrest peptide-based force-profile analysis, we assayed the force generated during the folding of RNH on the ribosome. Surprisingly, we found that population of the RNH folding intermediate is required to generate sufficient force to release a stall induced by the SecM stalling sequence and that readthrough of SecM directly correlates with the stability of the RNH folding intermediate. Together, these results imply that the folding pathway of RNH is unchanged on the ribosome. Furthermore, our findings indicate that the ribosome promotes RNH unfolding while the nascent chain is proximal to the ribosome, which may limit the deleterious effects of RNH misfolding and assist in folding fidelity.

Selective Affimers Recognise the BCL-2 Family Proteins BCL-xL and MCL-1 through Noncanonical Structural Motifs*.

The BCL-2 family is a challenging group of proteins to target selectively due to sequence and structural homologies across the family. Selective ligands for the BCL-2 family regulators of apoptosis are useful as probes to understand cell biology and apoptotic signalling pathways, and as starting points for inhibitor design. We have used phage display to isolate Affimer reagents (non-antibody-binding proteins based on a conserved scaffold) to identify ligands for MCL-1, BCL-xL , BCL-2, BAK and BAX, then used multiple biophysical characterisation methods to probe the interactions. We established that purified Affimers elicit selective recognition of their target BCL-2 protein. For anti-apoptotic targets BCL-xL and MCL-1, competitive inhibition of their canonical protein-protein interactions is demonstrated. Co-crystal structures reveal an unprecedented mode of molecular recognition; where a BH3 helix is normally bound, flexible loops from the Affimer dock into the BH3 binding cleft. Moreover, the Affimers induce a change in the target proteins towards a desirable drug-bound-like conformation. These proof-of-concept studies indicate that Affimers could be used as alternative templates to inspire the design of selective BCL-2 family modulators and more generally other protein-protein interaction inhibitors.

Folding pathway of an Ig domain is conserved on and off the ribosome.

Proteins that fold cotranslationally may do so in a restricted configurational space, due to the volume occupied by the ribosome. How does this environment, coupled with the close proximity of the ribosome, affect the folding pathway of a protein? Previous studies have shown that the cotranslational folding process for many proteins, including small, single domains, is directly affected by the ribosome. Here, we investigate the cotranslational folding of an all-ß Ig domain, titin I27. Using an arrest peptide-based assay and structural studies by cryo-EM, we show that I27 folds in the mouth of the ribosome exit tunnel. Simulations that use a kinetic model for the force dependence of escape from arrest accurately predict the fraction of folded protein as a function of length. We used these simulations to probe the folding pathway on and off the ribosome. Our simulations-which also reproduce experiments on mutant forms of I27-show that I27 folds, while still sequestered in the mouth of the ribosome exit tunnel, by essentially the same pathway as free I27, with only subtle shifts of critical contacts from the C to the N terminus.

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.

Folding and binding pathways of BH3-only proteins are encoded within their intrinsically disordered sequence, not templated by partner proteins.

Intrinsically disordered regions are present in one-third of eukaryotic proteins and are overrepresented in cellular processes such as signaling, suggesting that intrinsically disordered proteins (IDPs) may have a functional advantage over folded proteins. Upon interacting with a partner macromolecule, a subset of IDPs can fold and bind to form a well-defined three-dimensional conformation. For example, disordered BH3-only proteins bind promiscuously to a large number of homologous BCL-2 family proteins, where they fold to a helical structure in a groove on the BCL-2-like protein surface. As two protein chains are involved in the folding reaction, and the structure is only formed in the presence of the partner macromolecule, this raises the question of where the folding information is encoded. Here, we examine these coupled folding and binding reactions to determine which component determines the folding and binding pathway. Using Φ value analysis to compare transition state interactions between the disordered BH3-only proteins PUMA and BID and the folded BCL-2-like proteins A1 and MCL-1, we found that, even though the BH3-only protein is disordered in isolation and requires a stabilizing partner to fold, its folding and binding pathway is encoded in the IDP itself; the reaction is not templated by the folded partner. We suggest that, by encoding both its transition state and level of residual structure, an IDP can evolve a specific kinetic profile, which could be a crucial functional advantage of disorder.

Disorder drives cooperative folding in a multidomain protein.

Many human proteins contain intrinsically disordered regions, and disorder in these proteins can be fundamental to their function-for example, facilitating transient but specific binding, promoting allostery, or allowing efficient posttranslational modification. SasG, a multidomain protein implicated in host colonization and biofilm formation in Staphylococcus aureus, provides another example of how disorder can play an important role. Approximately one-half of the domains in the extracellular repetitive region of SasG are intrinsically unfolded in isolation, but these E domains fold in the context of their neighboring folded G5 domains. We have previously shown that the intrinsic disorder of the E domains mediates long-range cooperativity between nonneighboring G5 domains, allowing SasG to form a long, rod-like, mechanically strong structure. Here, we show that the disorder of the E domains coupled with the remarkable stability of the interdomain interface result in cooperative folding kinetics across long distances. Formation of a small structural nucleus at one end of the molecule results in rapid structure formation over a distance of 10 nm, which is likely to be important for the maintenance of the structural integrity of SasG. Moreover, if this normal folding nucleus is disrupted by mutation, the interdomain interface is sufficiently stable to drive the folding of adjacent E and G5 domains along a parallel folding pathway, thus maintaining cooperative folding.

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.

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.

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.

Conserved Helix-Flanking Prolines Modulate Intrinsically Disordered Protein:Target Affinity by Altering the Lifetime of the Bound Complex.

Appropriate integration of cellular signals requires a delicate balance of ligand-target binding affinities. Increasing the level of residual structure in intrinsically disordered proteins (IDPs), which are overrepresented in these cellular processes, has been shown previously to enhance binding affinities and alter cellular function. Conserved proline residues are commonly found flanking regions of IDPs that become helical upon interacting with a partner protein. Here, we mutate these helix-flanking prolines in p53 and MLL and find opposite effects on binding affinity upon an increase in free IDP helicity. In both cases, changes in affinity were due to alterations in dissociation, not association, rate constants, which is inconsistent with conformational selection mechanisms. We conclude that, contrary to previous suggestions, helix-flanking prolines do not regulate affinity by modulating the rate of complex formation. Instead, they influence binding affinities by controlling the lifetime of the bound complex.

Insights into Coupled Folding and Binding Mechanisms from Kinetic Studies.

Intrinsically disordered proteins (IDPs) are characterized by a lack of persistent structure. Since their identification more than a decade ago, many questions regarding their functional relevance and interaction mechanisms remain unanswered. Although most experiments have taken equilibrium and structural perspectives, fewer studies have investigated the kinetics of their interactions. Here we review and highlight the type of information that can be gained from kinetic studies. In particular, we show how kinetic studies of coupled folding and binding reactions, an important class of signaling event, are needed to determine mechanisms.

Single-molecule fluorescence reveals sequence-specific misfolding in multidomain proteins.

A large range of debilitating medical conditions is linked to protein misfolding, which may compete with productive folding particularly in proteins containing multiple domains. Seventy-five per cent of the eukaryotic proteome consists of multidomain proteins, yet it is not understood how interdomain misfolding is avoided. It has been proposed that maintaining low sequence identity between covalently linked domains is a mechanism to avoid misfolding. Here we use single-molecule Förster resonance energy transfer to detect and quantify rare misfolding events in tandem immunoglobulin domains from the I band of titin under native conditions. About 5.5 per cent of molecules with identical domains misfold during refolding in vitro and form an unexpectedly stable state with an unfolding half-time of several days. Tandem arrays of immunoglobulin-like domains in humans show significantly lower sequence identity between neighbouring domains than between non-adjacent domains. In particular, the sequence identity of neighbouring domains has been found to be preferentially below 40 per cent. We observe no misfolding for a tandem of naturally neighbouring domains with low sequence identity (24 per cent), whereas misfolding occurs between domains that are 42 per cent identical. Coarse-grained molecular simulations predict the formation of domain-swapped structures that are in excellent agreement with the observed transfer efficiency of the misfolded species. We infer that the interactions underlying misfolding are very specific and result in a sequence-specific domain-swapping mechanism. Diversifying the sequence between neighbouring domains seems to be a successful evolutionary strategy to avoid misfolding in multidomain proteins.

Allostery within a transcription coactivator is predominantly mediated through dissociation rate constants.

The kinase-inducible domain interacting (KIX) domain of CREB binding protein binds to multiple intrinsically disordered transcription factors in vivo at two distinct sites on its surface. Several reports have been made of allosteric communication between these two sites in this well-characterized model system. In this work, we have performed fluorescence stopped-flow measurements to investigate the kinetics of binding of five KIX binding proteins. We find that they all have similar association and dissociation rate constants for complex formation, despite their wide range of intrinsic helical propensities. Furthermore, by careful arrangement of pseudofirst-order conditions, we have been able to show that both association and dissociation rate constants are decreased when a partner is bound at the alternative site. These decreases suggest that positive allosteric effects are not mediated by structural changes in binding sites but rather, through a more general mechanism, largely mediated through dissociation, which we propose is largely related to changes in the flexibility of the KIX domain itself.

Interplay between partner and ligand facilitates the folding and binding of an intrinsically disordered protein.

Protein-protein interactions are at the heart of regulatory and signaling processes in the cell. In many interactions, one or both proteins are disordered before association. However, this disorder in the unbound state does not prevent many of these proteins folding to a well-defined, ordered structure in the bound state. Here we examine a typical system, where a small disordered protein (PUMA, p53 upregulated modulator of apoptosis) folds to an α-helix when bound to a groove on the surface of a folded protein (MCL-1, induced myeloid leukemia cell differentiation protein). We follow the association of these proteins using rapid-mixing stopped flow, and examine how the kinetic behavior is perturbed by denaturant and carefully chosen mutations. We demonstrate the utility of methods developed for the study of monomeric protein folding, including ß-Tanford values, Leffler α, Φ-value analysis, and coarse-grained simulations, and propose a self-consistent mechanism for binding. Folding of the disordered protein before binding does not appear to be required and few, if any, specific interactions are required to commit to association. The majority of PUMA folding occurs after the transition state, in the presence of MCL-1. We also examine the role of the side chains of folded MCL-1 that make up the binding groove and find that many favor equilibrium binding but, surprisingly, inhibit the association process.

Experimental evidence for a frustrated energy landscape in a three-helix-bundle protein family.

Energy landscape theory is a powerful tool for understanding the structure and dynamics of complex molecular systems, in particular biological macromolecules. The primary sequence of a protein defines its free-energy landscape and thus determines the folding pathway and the rate constants of folding and unfolding, as well as the protein's native structure. Theory has shown that roughness in the energy landscape will lead to slower folding, but derivation of detailed experimental descriptions of this landscape is challenging. Simple folding models show that folding is significantly influenced by chain entropy; proteins in which the contacts are local fold quickly, owing to the low entropy cost of forming stabilizing, native contacts during folding. For some protein families, stability is also a determinant of folding rate constants. Where these simple metrics fail to predict folding behaviour, it is probable that there are features in the energy landscape that are unusual. Such general observations cannot explain the folding behaviour of the R15, R16 and R17 domains of alpha-spectrin. R15 folds approximately 3,000 times faster than its homologues, although they have similar structures, stabilities and, as far as can be determined, transition-state stabilities. Here we show that landscape roughness (internal friction) is responsible for the slower folding and unfolding of R16 and R17. We use chimaeric domains to demonstrate that this internal friction is a property of the core, and suggest that frustration in the landscape of the slow-folding spectrin domains may be due to misdocking of the long helices during folding. Theoretical studies have suggested that rugged landscapes will result in slower folding; here we show experimentally that such a phenomenon directly influences the folding kinetics of a 'normal' protein, that is, one with a significant energy barrier that folds on a relatively slow, millisecond-second, timescale.

Separating the effects of internal friction and transition state energy to explain the slow, frustrated folding of spectrin domains.

The elongated three-helix bundle domains spectrin R16 and R17 fold some two to three orders of magnitude more slowly than their homologue R15. We have shown that this slow folding is due, at least in part, to roughness in the free-energy landscape of R16 and R17. We have proposed that this roughness is due to a frustrated search for the correct docking of partly preformed helices. However, this accounts for only a small part of the slowing of folding and unfolding. Five residues on the A helix of R15, when inserted together into R16 or R17, increase the folding rate constants, reduce landscape roughness, and alter the folding mechanism to one resembling R15. The effect of each of these mutations individually is investigated here. No one mutation causes the behavior seen for the five in combination. However, two mutations, E18F and K25V, significantly increase the folding and unfolding rates of both R16 and R17 but without a concomitant loss in landscape roughness. E18F has the greatest effect on the kinetics, and a Φ-value analysis of the C helix reveals that the folding mechanism is unchanged. For both E18F and K25V the removal of the charge and resultant transition state stabilization is the main origin of the faster folding. Consequently, the major cause of the unusually slow folding of R16 and R17 is the non-native burial of the two charged residues in the transition state. The slowing due to landscape roughness is only about fivefold.