Native-like Interactions Research Articles

Mapping functional regions of essential bacterial proteins with dominant-negative protein fragments.

Fragments of native proteins can act as dominant negative inhibitors, representing a substantial stress for cells. Massively-parallel measurements of growth inhibition by protein fragments have been used to discover peptide inhibitors and map protein interaction sites. However, the underlying principles governing fragment-based inhibition have thus far remained unclear. Here, we adapt a high-throughput inhibitory fragment assay for use in Escherichia coli, applying it to a set of ten essential proteins. This approach yielded single amino acid resolution maps of inhibitory activity, with peaks localized to functionally important interaction sites, including oligomerization interfaces and folding contacts. Leveraging these data, we perform a systematic analysis to uncover principles of fragment-based inhibition. We determine a robust negative correlation between susceptibility to inhibition and cellular protein concentration, demonstrating that inhibitory fragments likely act primarily by titrating native protein interactions. We also characterize a series of trade-offs related to fragment length, showing that shorter peptides allow higher-resolution mapping but suffer from lower activity. We employ an unsupervised statistical analysis to show that the inhibitory activities of protein fragments are largely driven not by generic properties such as charge, hydrophobicity, and secondary structure, but by the more specific characteristics of their bespoke macromolecular interactions. AlphaFold computational modeling of peptide complexes with one protein shows that the inhibitory activity of peptides is associated with their predicted ability to form native-like interactions. Overall, this work demonstrates fundamental characteristics of inhibitory protein fragment function and provides a foundation for understanding and controlling protein interactions in vivo.

A Glimpse into the Structural Properties of the Intermediate and Transition State in the Folding of Bromodomain 2 Domain 2 by Φ Value Analysis.

Bromodomains (BRDs) are small protein interaction modules of about 110 amino acids that selectively recognize acetylated lysine in histones and other proteins. These domains have been identified in a variety of multi-domain proteins involved in transcriptional regulation or chromatin remodeling in eukaryotic cells. BRD inhibition is considered an attractive therapeutic approach in epigenetic disorders, particularly in oncology. Here, we present a Φ value analysis to investigate the folding pathway of the second domain of BRD2 (BRD2(2)). Using an extensive mutational analysis based on 25 site-directed mutants, we provide structural information on both the intermediate and late transition state of BRD2(2). The data reveal that the C-terminal region represents part of the initial folding nucleus, while the N-terminal region of the domain consolidates its structure only later in the folding process. Furthermore, only a small number of native-like interactions have been identified, suggesting the presence of a non-compact, partially folded state with scarce native-like characteristics. Taken together, these results indicate that, in BRD2(2), a hierarchical mechanism of protein folding can be described with non-native interactions that play a significant role in folding.

Structural characterization of an on-pathway intermediate and transition state in the folding of the N-terminal SH2 domain from SHP2.

Src Homology 2 (SH2) domains are a class of protein domains that present a conserved three-dimensional structure and possess a crucial role in mediating protein-protein interactions. Despite their importance and abundance in the proteome, knowledge about the folding properties of SH2 domain is limited. Here we present an extensive mutational analysis (Φ value analysis) of the folding pathway of the N-SH2 domain of the Src homology region 2 domain-containing phosphatase-2 (SHP2) protein, a 104 residues domain that presents the classical SH2 domain fold (two α-helices flanking a central β-sheet composed of 3-5 antiparallel β-strands), with a fundamental role in mediating the interaction of SHP2 with its substrates and triggering key metabolic pathways in the cell. By analysing folding kinetic data we demonstrated that the folding pathway of N-SH2 presents an obligatory on-pathway intermediate that accumulates during the folding reaction. The production of 24 conservative site-directed variants allowed us to perform a Φ value analysis, by which we could fully characterize the intermediate and the transition state native-like interactions, providing a detailed quantitative analysis of the folding pathway of N-SH2. Results highlight the presence of a hydrophobic nucleus that stabilizes the intermediate, leading to a higher degree of native-like interactions in the transition state. Data are discussed and compared with previous works on SH2 domains.

Local and non‐local topological information in the denatured state ensemble of a β‐barrel protein

The folding of predominantly β-sheet proteins is complicated by the presence of a large number of non-local interactions in their native states, which increase the ruggedness of their folding energy landscapes. However, forming non-local contacts early in folding or even in the unfolded state can smooth the energy landscape and facilitate productive folding. We report that several sequence regions of a β-barrel protein, cellular retinoic acid-binding protein 1 (CRABP1), populate native-like secondary structure to a significant extent in the denatured state in 8 M urea. In addition, we provide evidence for both local and non-local interactions in the denatured state of CRABP1. NMR chemical shift perturbations (CSPs) under denaturing conditions upon substitution of single residues by mutation support the presence of several non-local interactions in topologically key sites, arguing that the denatured state is conformationally restricted and contains topological information for the native fold. Among the most striking non-local interactions are those between the N- and C-terminal regions, which are involved in closure of the native β-barrel. In addition, CSPs support the presence of two features in the denatured state: a major hydrophobic cluster involving residues from various parts of the sequence and a native-like interaction similar to one identified in previous studies as forming early in folding (Budyak et al., Structure 21, 476 [2013]). Taken together, our data support a model in which transient structures involving nonlocal interactions prime early folding interactions in CRABP1, determine its barrel topology, and may protect this predominantly β-sheet protein against aggregation.

The Folding Pathway of the KIX Domain.

The KIX domain is an 89-residues globular domain with an important role in mediating protein-protein interactions. The presence of two distinct binding sites in such a small domain makes KIX a suitable candidate to investigate the effect of the potentially divergent demands between folding and function. Here, we report an extensive mutational analysis of the folding pathway of the KIX domain, based on 30 site-directed mutants, which allow us to assess the structures of both the transition and denatured states. Data reveal that, while the transition state presents mostly native-like interactions, the denatured state is somewhat misfolded. We mapped some of the non-native contacts in the denatured state using a second round of mutagenesis, based on double mutant cycles on 15 double mutants. Interestingly, such a misfolding arises from non-native interactions involving the residues critical for the function of the protein. The results described in this work appear to highlight the diverging demands between folding and function that may lead to misfolding, which may be observed in the early stages of folding.

Structural Consequences of Calmodulin EF Hand Mutations.

Calmodulin (CaM) is a cytosolic Ca2+-binding protein that serves as a control element for many enzymes. It consists of two globular domains, each containing two EF hand pairs capable of binding Ca2+, joined by a flexible central linker region. CaM is able to bind and activate its target proteins in the Ca2+-replete and Ca2+-deplete forms. To study the Ca2+-dependent/independent properties of binding and activation of target proteins by CaM, CaM constructs with Ca2+-binding disrupting mutations of Asp to Ala at position one of each EF hand have been used. These CaM mutant proteins are deficient in binding Ca2+ in either the N-lobe EF hands (CaM12), C-lobe EF hands (CaM34), or all four EF hands (CaM1234). To investigate potential structural changes these mutations may cause, we performed detailed NMR studies of CaM12, CaM34, and CaM1234 including determining the solution structure of CaM1234. We then investigated if these CaM mutants affected the interaction of CaM with a target protein known to interact with apoCaM by determining the solution structure of CaM34 bound to the iNOS CaM binding domain peptide. The structures provide direct structural evidence of changes that are present in these Ca2+-deficient CaM mutants and show these mutations increase the hydrophobic exposed surface and decrease the electronegative surface potential throughout each lobe of CaM. These Ca2+-deficient CaM mutants may not be a true representation of apoCaM and may not allow for native-like interactions of apoCaM with its target proteins.

Limited Proteolysis Reveals That Amyloids from the 3D Domain-Swapping Cystatin B Have a Non-Native β-Sheet Topology

3D domain-swapping proteins form multimers by unfolding and then sharing of secondary structure elements, often with native-like interactions. Runaway domain swapping is proposed as a mechanism for folded proteins to form amyloid fibres, with examples including serpins and cystatins. Cystatin C amyloids cause a hereditary form of cerebral amyloid angiopathy whilst cystatin B aggregates are found in cases of Unverricht-Lundborg Syndrome, a progressive form of myoclonic epilepsy. Under conditions that favour fibrillisation, cystatins populate stable 3D domain-swapped dimers both in vitro and in vivo that represent intermediates on route to the formation of fibrils. Previous work on cystatin B amyloid fibrils revealed that the α-helical region of the protein becomes disordered and identified the conservation of a continuous 20-residue elongated β-strand (residues 39–58), the latter being a salient feature of the dimeric 3D domain-swapped structure. Here we apply limited proteolysis to cystatin B amyloid fibrils and show that not only the α-helical N-terminal of the protein (residues 1–35) but also the C-terminal of the protein (residues 80–98) can be removed without disturbing the underlying fibril structure. This observation is incompatible with previous models of cystatin amyloid fibrils where the β-sheet is assumed to retain its native antiparallel arrangement. We conclude that our data favour a more generic, at least partially parallel, arrangement for cystatin β-sheet structure in mature amyloids and propose a model that remains consistent with available data for amyloids from either cystatin B or cystatin C.

Fab-based bispecific antibody formats with robust biophysical properties and biological activity

A myriad of innovative bispecific antibody (BsAb) platforms have been reported. Most require significant protein engineering to be viable from a development and manufacturing perspective. Single-chain variable fragments (scFvs) and diabodies that consist only of antibody variable domains have been used as building blocks for making BsAbs for decades. The drawback with Fv-only moieties is that they lack the native-like interactions with CH1/CL domains that make antibody Fab regions stable and soluble. Here, we utilize a redesigned Fab interface to explore 2 novel Fab-based BsAbs platforms. The redesigned Fab interface designs limit heavy and light chain mixing when 2 Fabs are co-expressed simultaneously, thus allowing the use of 2 different Fabs within a BsAb construct without the requirement of one or more scFvs. We describe the stability and activity of a HER2×HER2 IgG-Fab BsAb, and compare its biophysical and activity properties with those of an IgG-scFv that utilizes the variable domains of the same parental antibodies. We also generated an EGFR × CD3 tandem Fab protein with a similar format to a tandem scFv (otherwise known as a bispecific T cell engager or BiTE). We show that the Fab-based BsAbs have superior biophysical properties compared to the scFv-based BsAbs. Additionally, the Fab-based BsAbs do not simply recapitulate the activity of their scFv counterparts, but are shown to possess unique biological activity.

Surprising Abundance of Misfolding during Refolding of Multidomain Proteins

The covalently linked domains constituting a multidomain protein typically share a similar fold and sequence, are generally stable in isolation, and are largely capable of independent folding. However, the elegant interplay of forces leading an amino acid sequence to its native state is more complex in proteins with multiple domains because of inter-domain interactions, which complicate the folding energy landscape.In previous experiments, employing immunoglobulin-like (Ig-like) domains of the human multidomain protein titin, we showed that tandem repeats of domains with high sequence identity can form a stable misfolded state upon refolding in physiological solution conditions. Conversely, tandem constructs of natural neighbouring domains, which display a low sequence identity, did not misfold. This supported the hypothesis that sequence identity between neighbouring domains in multidomain proteins is reduced as a result of evolutionary pressure to avoid misfolding.With a combination of single-molecule Forster resonance energy transfer, microfluidic mixing, stopped-flow kinetics and molecular dynamics simulations, we now demonstrate that Ig-like domains can transiently populate a surprisingly broad range of misfolded conformations on the sub-second timescale. Using tandem repeats of Ig-like domains, we can resolve both strand-swapped misfolds dominated by native-like interactions and, remarkably, a non-native-like, largely disordered type of misfolded state which so far was never observed experimentally, characterized by promiscuous interactions. Even more surprisingly, both types of misfolding are detected also for the naturally occurring tandem repeat, showing how finely the propensities of folding and misfolding have been balanced by co-evolution of adjacent domains to avoid stable misfolded states formation. On longer timescales, however, all or most of the protein molecules are able to reach the native state, demonstrating that the overall free energy surface is still sufficiently optimized for the protein to efficiently reach its correctly folded state.

How Well Does a Funneled Energy Landscape Capture the Folding Mechanism of Spectrin Domains?

Three structurally similar domains from α-spectrin have been shown to fold very differently. First, there is a contrast in the folding mechanism, as probed by Φ-value analysis, between the R15 domain and the R16 and R17 domains. Second, there are very different contributions from internal friction to folding: the folding rate of the R15 domain was found to be inversely proportional to solvent viscosity, showing no apparent frictional contribution from the protein, but in the other two domains, a large internal friction component was evident. Non-native misdocking of helices has been suggested to be responsible for this phenomenon. Here, I study the folding of these three proteins with minimalist coarse-grained models based on a funneled energy landscape. Remarkably, I find that, despite the absence of non-native interactions, the differences in folding mechanism of the domains are well captured by the model, and the agreement of the Φ-values with experiment is fairly good. On the other hand, within the context of this model, there are no significant differences in diffusion coefficient along the chosen folding coordinate, and the model cannot explain the large differences in folding rates between the proteins found experimentally. These results are nonetheless consistent with the expectations from the energy landscape perspective of protein folding, namely, that the folding mechanism is primarily determined by the native-like interactions present in the Gō-like model, with missing non-native interactions being required to explain the differences in "internal friction" seen in experiment.

A residue in helical conformation in the native state adopts a β‐strand conformation in the folding transition state despite its high and canonical ϕ‐value

Delineating structures of the transition states in protein folding reactions has provided great insight into the mechanisms by which proteins fold. The most common method for obtaining this information is Φ-value analysis, which is carried out by measuring the changes in the folding and unfolding rates caused by single amino acid substitutions at various positions within a given protein. Canonical Φ-values range between 0 and 1, and residues displaying high values within this range are interpreted to be important in stabilizing the transition state structure, and to elicit this stabilization through native-like interactions. Although very successful in defining the general features of transition state structures, Φ-value analysis can be confounded when non-native interactions stabilize this state. In addition, direct information on backbone conformation within the transition state is not provided. In the work described here, we have investigated structure formation at a conserved β-bulge (with helical conformation) in the Fyn SH3 domain by characterizing the effects of substituting all natural amino acids at one position within this structural motif. By comparing the effects on folding rates of these substitutions with database-derived local structure propensity values, we have determined that this position adopts a non-native backbone conformation in the folding transition state. This result is surprising because this position displays a high and canonical Φ-value of 0.7. This work emphasizes the potential role of non-native conformations in folding pathways and demonstrates that even positions displaying high and canonical Φ-values may, nevertheless, adopt a non-native conformation in the transition state.

The H2A–H2B Dimeric Kinetic Intermediate Is Stabilized by Widespread Hydrophobic Burial with Few Fully Native Interactions

The H2A–H2B histone heterodimer folds via monomeric and dimeric kinetic intermediates. Within ∼5 ms, the H2A and H2B polypeptides associate in a nearly diffusion limited reaction to form a dimeric ensemble, denoted I2 and I2⁎, the latter being a subpopulation characterized by a higher content of nonnative structure (NNS). The I2 ensemble folds to the native heterodimer, N2, through an observable, first-order kinetic phase. To determine the regions of structure in the I2 ensemble, we characterized 26 Ala mutants of buried hydrophobic residues, spanning the three helices of the canonical histone folds of H2A and H2B and the H2B C-terminal helix. All but one targeted residue contributed significantly to the stability of I2, the transition state and N2; however, only residues in the hydrophobic core of the dimer interface perturbed the I2⁎ population. Destabilization of I2⁎ correlated with slower folding rates, implying that NNS is not a kinetic trap but rather accelerates folding. The pattern of Φ values indicated that residues forming intramolecular interactions in the peripheral helices contributed similar stability to I2 and N2, but residues involved in intermolecular interactions in the hydrophobic core are only partially folded in I2. These findings suggest a dimerize-then-rearrange model. Residues throughout the histone fold contribute to the stability of I2, but after the rapid dimerization reaction, the hydrophobic core of the dimer interface has few fully native interactions. In the transition state leading to N2, more native-like interactions are developed and nonnative interactions are rearranged.

Amyloid-like Fibrils from a Domain-swapping Protein Feature a Parallel, in-Register Conformation without Native-like Interactions

The formation of amyloid-like fibrils is characteristic of various diseases, but the underlying mechanism and the factors that determine whether, when, and how proteins form amyloid, remain uncertain. Certain mechanisms have been proposed based on the three-dimensional or runaway domain swapping, inspired by the fact that some proteins show an apparent correlation between the ability to form domain-swapped dimers and a tendency to form fibrillar aggregates. Intramolecular β-sheet contacts present in the monomeric state could constitute intermolecular β-sheets in the dimeric and fibrillar states. One example is an amyloid-forming mutant of the immunoglobulin binding domain B1 of streptococcal protein G, which in its native conformation consists of a four-stranded β-sheet and one α-helix. Under native conditions this mutant adopts a domain-swapped dimer, and it also forms amyloid-like fibrils, seemingly in correlation to its domain-swapping ability. We employ magic angle spinning solid-state NMR and other methods to examine key structural features of these fibrils. Our results reveal a highly rigid fibril structure that lacks mobile domains and indicate a parallel in-register β-sheet structure and a general loss of native conformation within the mature fibrils. This observation contrasts with predictions that native structure, and in particular intermolecular β-strand interactions seen in the dimeric state, may be preserved in "domain-swapping" fibrils. We discuss these observations in light of recent work on related amyloid-forming proteins that have been argued to follow similar mechanisms and how this may have implications for the role of domain-swapping propensities for amyloid formation.

Investigating the Refolding Pathway of Human Acidic Fibroblast Growth Factor (hFGF-1) from the Residual Structure(s) Obtained by Denatured-State Hydrogen/Deuterium Exchange

Human fibroblast growth factor 1 (hFGF-1) consists of 12 anti-parallel β-strands arranged into a β-trefoil architecture. We directly measured hydrogen/deuterium exchange rates on the urea-denatured hFGF-1 to obtain the information with regard to the persistent residual interaction(s) in the unfolded hFGF-1. Thirty-eight residues whose heteronuclear single quantum coherence cross-peaks can be observed after exchange show higher protections than those predicted for the same residues in a random coil conformation, suggesting the existence of residual structure(s). The urea-denaturation of hFGF-1 tested by both circular dichroism and fluorescence spectroscopy indicated that the unfolding process is a cooperative two-state process and that the residual structures observed did not originate from the existence of a partially structured intermediate. The coincident disappearance of the native heteronuclear single quantum coherence cross-peaks during the urea-denaturation process suggests that the residual structures observed contain no nativelike interactions. The protected residues (fold ons) in the urea-denatured state are mostly those that exchange slowly in the native state H/D exchange. The distribution of these fold ons in the native structure of hFGF-1 suggests that the refolding starts by collisions between the residual structures (microdomains) between the β-strands VI and VII, and between the β-strands II and III, which appear to be two independent refolding coordinates during the refolding process.

Experimental snapshots of a protein-DNA binding landscape

Protein recognition of DNA sites is a primary event for gene function. Its ultimate mechanistic understanding requires an integrated structural, dynamic, kinetic, and thermodynamic dissection that is currently limited considering the hundreds of structures of protein-DNA complexes available. We describe a protein-DNA-binding pathway in which an initial, diffuse, transition state ensemble with some nonnative contacts is followed by formation of extensive nonnative interactions that drive the system into a kinetic trap. Finally, nonnative contacts are slowly rearranged into native-like interactions with the DNA backbone. Dissimilar protein-DNA interfaces that populate along the DNA-binding route are explained by a temporary degeneracy of protein-DNA interactions, centered on "dual-role" residues. The nonnative species slow down the reaction allowing for extended functionality.

Intra- and Interdomain Effects Due to Mutation of Calcium-binding Sites in Calmodulin

The IQ-motif protein PEP-19, binds to the C-domain of calmodulin (CaM) with significantly different k(on) and k(off) rates in the presence and absence of Ca(2+), which could play a role in defining the levels of free CaM during Ca(2+) transients. The initial goal of the current study was to determine whether Ca(2+) binding to sites III or IV in the C-domain of CaM was responsible for affecting the kinetics of binding PEP-19. EF-hand Ca(2+)-binding sites were selectively inactivated by the common strategy of changing Asp to Ala at the X-coordination position. Although Ca(2+) binding to both sites III and IV appeared necessary for native-like interactions with PEP-19, the data also indicated that the mutations caused undesirable structural alterations as evidenced by significant changes in amide chemical shifts for apoCaM. Mutations in the C-domain also affected chemical shifts in the unmodified N-domain, and altered the Ca(2+) binding properties of the N-domain. Conversion of Asp(93) to Ala caused the greatest structural perturbations, possibly due to the loss of stabilizing hydrogen bonds between the side chain of Asp(93) and backbone amides in apo loop III. Thus, although these mutations inhibit binding of Ca(2+), the mutated CaM may not be able to support potentially important native-like activity of the apoprotein. This should be taken into account when designing CaM mutants for expression in cell culture.

Energetic Frustration of Apomyoglobin Folding: Role of the B Helix

Apomyoglobin folds by a sequential mechanism in which the A, G, and H helix regions undergo rapid collapse to form a compact intermediate onto which the central portion of the B helix subsequently docks. To investigate the factors that frustrate folding, we have made mutations in the N-terminus of the B helix to stabilize helical structure (in the mutant G23A/G25A) and to promote native-like hydrophobic packing interactions with helix G (in the mutant H24L/H119F). The kinetic and equilibrium intermediates of G23A/G25A and H24L/H119F were studied by hydrogen exchange pulse labeling and interrupted hydrogen/deuterium exchange combined with NMR. For both mutants, stabilization of helical structure in the N-terminal region of the B helix is confirmed by increased exchange protection in the equilibrium molten globule states near pH 4. Increased protection is also observed in the GH turn region in the G23A/G25A mutant, suggesting that stabilization of the B helix facilitates native-like interactions with the C-terminal region of helix G. These interactions are further enhanced in H24L/H119F. The kinetic burst phase intermediates of both mutants show increased protection, relative to wild-type protein, of amides in the N-terminus of the B helix and in part of the E helix. Stabilization of the E helix in the intermediate is attributed to direct interactions between E helix residues and the newly stabilized N-terminus of helix B. Stabilization of native packing between the B and G helices in H24L/H119F also favors formation of native-like interactions in the GH turn and between the G and H helices in the ensemble of burst phase intermediates. We conclude that instability at the N-terminus of the B helix of apomyoglobin contributes to the energetic frustration of folding by preventing docking and stabilization of the E helix.

Residue Specific Analysis of Frustration in Folding Landscape of Repeat Alpha/Beta Protein Apoflavodoxin

Topological frustrated proteins can give rise to complex folding pathways and transition state profiles. To better understand such systems, we combine experimental and computational methods to study Desulfovibrio desulfuricans apoflavodoxin by producing several point mutation variants. By equilibrium unfolding experiments, we first revealed how different secondary-structure elements contribute to overall protein resistance towards heat and urea. Next using stopped-flow mixing coupled to far-UV circular dichroism (CD), we probed how individual residues affect the amount of structure formed in the experimentally-detected burst-phase intermediate. Together with in silico folding route analysis of the same point-mutated variants and computation of the growth in nucleation size during early folding, computer simulations suggested the presence of two competing folding nuclei at opposite sides of the central β-strand 3 (i.e. at β-strands 1 and 4), which cause early topological frustration (i.e., misfolding) in the folding landscape. Particularly, the extent of heterogeneity in the folding nucleigrowth correlates with the in vitro burst phase CD amplitude. In addition, Φ-value analysis (in vitro and in silico) of the overall folding barrier to apoflavodoxin's native state revealed that native-like interactions in most of the β-strands must form in transition state. Our study reveals that an imbalanced competition between the two sides of apoflavodoxin's central β-sheet directs initial misfolding while proper alignment on both sides of β-strand 3 is necessary for productive folding.

Residue-Specific Analysis of Frustration in the Folding Landscape of Repeat β/α Protein Apoflavodoxin

Flavodoxin adopts the common repeat β/α topology and folds in a complex kinetic reaction with intermediates. To better understand this reaction, we analyzed a set of Desulfovibrio desulfuricans apoflavodoxin variants with point mutations in most secondary structure elements by in vitro and in silico methods. By equilibrium unfolding experiments, we first revealed how different secondary structure elements contribute to overall protein resistance to heat and urea. Next, using stopped-flow mixing coupled with far-UV circular dichroism, we probed how individual residues affect the amount of structure formed in the experimentally detected burst-phase intermediate. Together with in silico folding route analysis of the same point-mutated variants and computation of growth in nucleation size during early folding, computer simulations suggested the presence of two competing folding nuclei at opposite sides of the central β-strand 3 (i.e., at β-strands 1 and 4), which cause early topological frustration (i.e., misfolding) in the folding landscape. Particularly, the extent of heterogeneity in folding nuclei growth correlates with the in vitro burst-phase circular dichroism amplitude. In addition, ϕ-value analysis (in vitro and in silico) of the overall folding barrier to apoflavodoxin's native state revealed that native-like interactions in most of the β-strands must form in transition state. Our study reveals that an imbalanced competition between the two sides of apoflavodoxin's central β-sheet directs initial misfolding, while proper alignment on both sides of β-strand 3 is necessary for productive folding.

Thermodynamic Aspects of Coupled Binding and Folding of an Intrinsically Disordered Protein: A Computational Alanine Scanning Study

We have performed computational alanine scanning to gain insight into thermodynamic aspects of coupled folding and binding of the phosphorylated kinase-inducible activation domain (pKID) of CREB, and intrinsically disorder protein (IDP), to KIX. Residues of pKID directly involved in nativelike interactions at the interface were found to steeply lower their contribution to the free energy of binding. The results suggest that, with a steep lowering of the free energy of binding, it is ensured that pKID has full control of binding specificity, is not trapped in a non-native conformational state, has full plasticity for folding, and forms a thermodynamically stable complex with KIX without compromising binding kinetics.