Folding kinetics of the cooperatively folded subdomain of the IκBα ankyrin repeat domain

Site-specific incorporation of biophysical probes into NF-ĸB with non-canonical amino acids

The transcription factor NF-ĸB is a central mediator of immune and inflammatory responses. To understand the regulation of NF-ĸB, it is important to probe the underlying thermodynamics, kinetics, and conformational dynamics of the NF-ĸB/IĸBα/DNA interaction network. The development of genetic incorporation of non-canonical amino acids (ncAA) has enabled the installation of biophysical probes into proteins with site specificity. Recent single-molecule FRET (smFRET) studies of NF-ĸB with site-specific labeling via ncAA incorporation revealed the conformational dynamics for kinetic control of DNA-binding mediated by IĸBα. Here we report the design and protocols for incorporating the ncAA p-azidophenylalanine (pAzF) into NF-ĸB and site-specific fluorophore labeling with copper-free click chemistry for smFRET. We also expanded the ncAA toolbox of NF-ĸB to include p-benzoylphenylalanine (pBpa) for UV crosslinking mass spectrometry (XL-MS) and incorporated both pAzF and pBpa into the full-length NF-ĸB RelA subunit which includes the intrinsically disordered transactivation domain.

Evolution and folding of repeat proteins

Repeat proteins are made with tandem copies of similar amino acid stretches that fold into elongated architectures. These proteins constitute excellent model systems to investigate how evolution relates to structure, folding, and function. Here, we propose a scheme to map evolutionary information at the sequence level to a coarse-grained model for repeat-protein folding and use it to investigate the folding of thousands of repeat proteins. We model the energetics by a combination of an inverse Potts-model scheme with an explicit mechanistic model of duplications and deletions of repeats to calculate the evolutionary parameters of the system at the single-residue level. These parameters are used to inform an Ising-like model that allows for the generation of folding curves, apparent domain emergence, and occupation of intermediate states that are highly compatible with experimental data in specific case studies. We analyzed the folding of thousands of natural Ankyrin repeat proteins and found that a multiplicity of folding mechanisms are possible. Fully cooperative all-or-none transitions are obtained for arrays with enough sequence-similar elements and strong interactions between them, while noncooperative element-by-element intermittent folding arose if the elements are dissimilar and the interactions between them are energetically weak. Additionally, we characterized nucleation-propagation and multidomain folding mechanisms. We show that the global stability and cooperativity of the repeating arrays can be predicted from simple sequence scores.

Folding and Stability of Ankyrin Repeats Control Biological Protein Function

Ankyrin repeat proteins are found in all three kingdoms of life. Fundamentally, these proteins are involved in protein-protein interaction in order to activate or suppress biological processes. The basic architecture of these proteins comprises repeating modules forming elongated structures. Due to the lack of long-range interactions, a graded stability among the repeats is the generic properties of this protein family determining both protein folding and biological function. Protein folding intermediates were frequently found to be key for the biological functions of repeat proteins. In this review, we discuss most recent findings addressing this close relation for ankyrin repeat proteins including DARPins, Notch receptor ankyrin repeat domain, IκBα inhibitor of NFκB, and CDK inhibitor p19^INK4d. The role of local folding and unfolding and gradual stability of individual repeats will be discussed during protein folding, protein-protein interactions, and post-translational modifications. The conformational changes of these repeats function as molecular switches for biological regulation, a versatile property for modern drug discovery.

Large Ankyrin repeat proteins are formed with similar and energetically favorable units

Ankyrin containing proteins are one of the most abundant repeat protein families present in all extant organisms. They are made with tandem copies of similar amino acid stretches that fold into elongated architectures. Here, we built and curated a dataset of 200 thousand proteins that contain 1.2 million Ankyrin regions and characterize the abundance, structure and energetics of the repetitive regions in natural proteins. We found that there is a continuous roughly exponential variety of array lengths with an exceptional frequency at 24 repeats. We described that individual repeats are seldom interrupted with long insertions and accept few deletions, in line with the known tertiary structures. We found that longer arrays are made up of repeats that are more similar to each other than shorter arrays, and display more favourable folding energy, hinting at their evolutionary origin. The array distributions show that there is a physical upper limit to the size of an array of repeats of about 120 copies, consistent with the limit found in nature. The identity patterns within the arrays suggest that they may have originated by sequential copies of more than one Ankyrin unit.

Folding cooperativity and allosteric function in the tandem-repeat protein class

The term allostery was originally developed to describe structural changes in one binding site induced by the interaction of a partner molecule with a distant binding site, and it has been studied in depth in the field of enzymology. Here, we discuss the concept of action at a distance in relation to the folding and function of the solenoid class of tandem-repeat proteins such as tetratricopeptide repeats (TPRs) and ankyrin repeats. Distantly located repeats fold cooperatively, even though only nearest-neighbour interactions exist in these proteins. A number of repeat-protein scaffolds have been reported to display allosteric effects, transferred through the repeat array, that enable them to direct the activity of the multi-subunit enzymes within which they reside. We also highlight a recently identified group of tandem-repeat proteins, the RRPNN subclass of TPRs, recent crystal structures of which indicate that they function as allosteric switches to modulate multiple bacterial quorum-sensing mechanisms. We believe that the folding cooperativity of tandem-repeat proteins and the biophysical mechanisms that transform them into allosteric switches are intimately intertwined. This opinion piece aims to combine our understanding of the two areas and develop ideas on their common underlying principles.This article is part of a discussion meeting issue 'Allostery and molecular machines'.

PyFolding: Open-Source Graphing, Simulation, and Analysis of the Biophysical Properties of Proteins

For many years, curve-fitting software has been heavily utilized to fit simple models to various types of biophysical data. Although such software packages are easy to use for simple functions, they are often expensive and present substantial impediments to applying more complex models or for the analysis of large data sets. One field that is reliant on such data analysis is the thermodynamics and kinetics of protein folding. Over the past decade, increasingly sophisticated analytical models have been generated, but without simple tools to enable routine analysis. Consequently, users have needed to generate their own tools or otherwise find willing collaborators. Here we present PyFolding, a free, open-source, and extensible Python framework for graphing, analysis, and simulation of the biophysical properties of proteins. To demonstrate the utility of PyFolding, we have used it to analyze and model experimental protein folding and thermodynamic data. Examples include: 1) multiphase kinetic folding fitted to linked equations, 2) global fitting of multiple data sets, and 3) analysis of repeat protein thermodynamics with Ising model variants. Moreover, we demonstrate how PyFolding is easily extensible to novel functionality beyond applications in protein folding via the addition of new models. Example scripts to perform these and other operations are supplied with the software, and we encourage users to contribute notebooks and models to create a community resource. Finally, we show that PyFolding can be used in conjunction with Jupyter notebooks as an easy way to share methods and analysis for publication and among research teams.

PFDB: A standardized protein folding database with temperature correction

We constructed a standardized protein folding kinetics database (PFDB) in which the logarithmic rate constants of all listed proteins are calculated at the standard temperature (25 °C). A temperature correction based on the Eyring–Kramers equation was introduced for proteins whose folding kinetics were originally measured at temperatures other than 25 °C. We verified the temperature correction by comparing the logarithmic rate constants predicted and experimentally observed at 25 °C for 14 different proteins, and the results demonstrated improvement of the quality of the database. PFDB consists of 141 (89 two-state and 52 non-two-state) single-domain globular proteins, which has the largest number among the currently available databases of protein folding kinetics. PFDB is thus intended to be used as a standard for developing and testing future predictive and theoretical studies of protein folding. PFDB can be accessed from the following link: http://lee.kias.re.kr/~bala/PFDB.

Inferring repeat-protein energetics from evolutionary information

Natural protein sequences contain a record of their history. A common constraint in a given protein family is the ability to fold to specific structures, and it has been shown possible to infer the main native ensemble by analyzing covariations in extant sequences. Still, many natural proteins that fold into the same structural topology show different stabilization energies, and these are often related to their physiological behavior. We propose a description for the energetic variation given by sequence modifications in repeat proteins, systems for which the overall problem is simplified by their inherent symmetry. We explicitly account for single amino acid and pair-wise interactions and treat higher order correlations with a single term. We show that the resulting evolutionary field can be interpreted with structural detail. We trace the variations in the energetic scores of natural proteins and relate them to their experimental characterization. The resulting energetic evolutionary field allows the prediction of the folding free energy change for several mutants, and can be used to generate synthetic sequences that are statistically indistinguishable from the natural counterparts.

Functional importance of stripping in NFκB signaling revealed by a stripping-impaired IκBα mutant

Stress-response transcription factors such as NFκB turn on hundreds of genes and must have a mechanism for rapid cessation of transcriptional activation. We recently showed that the inhibitor of NFκB signaling, IκBα, dramatically accelerates the dissociation of NFκB from transcription sites, a process we have called "stripping." To test the role of the IκBα C-terminal PEST (rich in proline, glutamic acid, serine, and threonine residues) sequence in NFκB stripping, a mutant IκBα was generated in which five acidic PEST residues were mutated to their neutral analogs. This IκBα(5xPEST) mutant was impaired in stripping NFκB from DNA and formed a more stable intermediate ternary complex than that formed from IκBα(WT) because DNA dissociated more slowly. NMR and amide hydrogen-deuterium exchange mass spectrometry showed that the IκBα(5xPEST) appears to be "caught in the act of stripping" because it is not yet completely in the folded and NFκB-bound state. When the mutant was introduced into cells, the rate of postinduction IκBα-mediated export of NFκB from the nucleus decreased markedly.

Dissecting and reprogramming the folding and assembly of tandem-repeat proteins

Studying protein folding and protein design in globular proteins presents significant challenges because of the two related features, topological complexity and co-operativity. In contrast, tandem-repeat proteins have regular and modular structures composed of linearly arrayed motifs. This means that the biophysics of even giant repeat proteins is highly amenable to dissection and to rational design. Here we discuss what has been learnt about the folding mechanisms of tandem-repeat proteins. The defining features that have emerged are: (i) accessibility of multiple distinct routes between denatured and native states, both at equilibrium and under kinetic conditions; (ii) different routes are favoured for folding compared with unfolding; (iii) unfolding energy barriers are broad, reflecting stepwise unravelling of an array repeat by repeat; (iv) highly co-operative unfolding at equilibrium and the potential for exceptionally high thermodynamic stabilities by introducing consensus residues; (v) under force, helical-repeat structures are very weak with non-co-operative unfolding leading to elasticity and buffering effects. This level of understanding should enable us to create repeat proteins with made-to-measure folding mechanisms, in which one can dial into the sequence the order of repeat folding, number of pathways taken, step size (co-operativity) and fine-structure of the kinetic energy barriers.

Protein Folding Mechanism of the Dimeric AmphiphysinII/Bin1 N-BAR Domain

The human AmphyphisinII/Bin1 N-BAR domain belongs to the BAR domain superfamily, whose members sense and generate membrane curvatures. The N-BAR domain is a 57 kDa homodimeric protein comprising a six helix bundle. Here we report the protein folding mechanism of this protein as a representative of this protein superfamily. The concentration dependent thermodynamic stability was studied by urea equilibrium transition curves followed by fluorescence and far-UV CD spectroscopy. Kinetic unfolding and refolding experiments, including rapid double and triple mixing techniques, allowed to unravel the complex folding behavior of N-BAR. The equilibrium unfolding transition curve can be described by a two-state process, while the folding kinetics show four refolding phases, an additional burst reaction and two unfolding phases. All fast refolding phases show a rollover in the chevron plot but only one of these phases depends on the protein concentration reporting the dimerization step. Secondary structure formation occurs during the three fast refolding phases. The slowest phase can be assigned to a proline isomerization. All kinetic experiments were also followed by fluorescence anisotropy detection to verify the assignment of the dimerization step to the respective folding phase. Based on these experiments we propose for N-BAR two parallel folding pathways towards the homodimeric native state depending on the proline conformation in the unfolded state.

Direct observation of parallel folding pathways revealed using a symmetric repeat protein system

Although progress has been made to determine the native fold of a polypeptide from its primary structure, the diversity of pathways that connect the unfolded and folded states has not been adequately explored. Theoretical and computational studies predict that proteins fold through parallel pathways on funneled energy landscapes, although experimental detection of pathway diversity has been challenging. Here, we exploit the high translational symmetry and the direct length variation afforded by linear repeat proteins to directly detect folding through parallel pathways. By comparing folding rates of consensus ankyrin repeat proteins (CARPs), we find a clear increase in folding rates with increasing size and repeat number, although the size of the transition states (estimated from denaturant sensitivity) remains unchanged. The increase in folding rate with chain length, as opposed to a decrease expected from typical models for globular proteins, is a clear demonstration of parallel pathways. This conclusion is not dependent on extensive curve-fitting or structural perturbation of protein structure. By globally fitting a simple parallel-Ising pathway model, we have directly measured nucleation and propagation rates in protein folding, and have quantified the fluxes along each path, providing a detailed energy landscape for folding. This finding of parallel pathways differs from results from kinetic studies of repeat-proteins composed of sequence-variable repeats, where modest repeat-to-repeat energy variation coalesces folding into a single, dominant channel. Thus, for globular proteins, which have much higher variation in local structure and topology, parallel pathways are expected to be the exception rather than the rule.

Highly polarized C-terminal transition state of the leucine-rich repeat domain of PP32 is governed by local stability

The leucine-rich repeat domain of PP32 is composed of five β-strand-containing repeats anchored by terminal caps. These repeats differ in sequence but are similar in structure, providing a means to connect topology, sequence, and folding pathway selection. Through kinetic studies of PP32, we find folding to be rate-limited by the formation of an on-pathway intermediate. Destabilizing core substitutions reveal a transition state ensemble that is highly polarized toward the C-terminal repeat and cap. To determine if this nucleus for folding corresponds to the most stable region of PP32, we monitored amide hydrogen exchange by NMR spectroscopy. Indeed, we find the highest protection to be biased toward the C terminus. Sequence manipulations that destabilize the C terminus spread out the transition state toward the middle of the protein. Consistent with results for helical ankyrin repeat proteins, these results suggest that local stabilities determine folding pathways.

Falling down: landscape and kinetics of one-dimensional protein folding

In this issue of Structure, Tsytlonok and colleagues describe the folding landscape of the giant HEAT-repeat protein PR65/A (a molecular adaptor of protein phosphatase 2A) by using experimental and theoretical methods. Both approaches agree in suggesting the presence of parallel folding pathways with several intermediates.

Diffuse transition state structure for the unfolding of a leucine-rich repeat protein

Tandem-repeat proteins, such as leucine-rich repeats, comprise arrays of small structural motifs that pack in a linear fashion to produce elongated architectures. They lack contacts between residues that are distant in primary sequence, a feature that distinguishes them from the complex topologies of globular proteins. Here we have investigated the unfolding pathway of the leucine-rich repeat domain of the mRNA export protein TAP (TAPLRR) using Φ-value analysis. Whereas most of the tandem-repeat proteins studied to date have been found to unfold via a polarised mechanism in which only a small, localised number of repeats are structured in the transition state, the unfolding mechanism of TAPLRR is more diffuse in nature. In the transition state for unfolding of TAPLRR, three of the four LRRs are highly structured and non-native interactions are formed within the N-terminal α-helical cap and the first LRR. Thus, the α-helical cap plays an important role in which non-native interactions are required to provide a scaffold for the LRRs to pack against in the folding reaction.

A disorder-induced domino-like destabilization mechanism governs the folding and functional dynamics of the repeat protein IκBα

The stability of the repeat protein IκBα, a transcriptional inhibitor in mammalian cells, is critical in the functioning of the NF-κB signaling module implicated in an array of cellular processes, including cell growth, disease, immunity and apoptosis. Structurally, IκBα is complex, with both ordered and disordered regions, thus posing a challenge to the available computational protocols to model its conformational behavior. Here, we introduce a simple procedure to model disorder in systems that undergo binding-induced folding that involves modulation of the contact map guided by equilibrium experimental observables in combination with an Ising-like Wako-Saitô-Muñoz-Eaton model. This one-step procedure alone is able to reproduce a variety of experimental observables, including ensemble thermodynamics (scanning calorimetry, pre-transitions, m-values) and kinetics (roll-over in chevron plot, intermediates and their identity), and is consistent with hydrogen-deuterium exchange measurements. We further capture the intricate distance-dynamics between the domains as measured by single-molecule FRET by combining the model predictions with simple polymer physics arguments. Our results reveal a unique mechanism at work in IκBα folding, wherein disorder in one domain initiates a domino-like effect partially destabilizing neighboring domains, thus highlighting the effect of symmetry-breaking at the level of primary sequences. The offshoot is a multi-state and a dynamic conformational landscape that is populated by increasingly partially folded ensembles upon destabilization. Our results provide, in a straightforward fashion, a rationale to the promiscuous binding and short intracellular half-life of IκBα evolutionarily engineered into it through repeats with variable stabilities and expand the functional repertoire of disordered regions in proteins.

Detecting repetitions and periodicities in proteins by tiling the structural space

The notion of energy landscapes provides conceptual tools for understanding the complexities of protein folding and function. Energy landscape theory indicates that it is much easier to find sequences that satisfy the “Principle of Minimal Frustration” when the folded structure is symmetric (Wolynes, P. G. Symmetry and the Energy Landscapes of Biomolecules. Proc. Natl. Acad. Sci. U.S.A.1996, 93, 14249–14255). Similarly, repeats and structural mosaics may be fundamentally related to landscapes with multiple embedded funnels. Here we present analytical tools to detect and compare structural repetitions in protein molecules. By an exhaustive analysis of the distribution of structural repeats using a robust metric, we define those portions of a protein molecule that best describe the overall structure as a tessellation of basic units. The patterns produced by such tessellations provide intuitive representations of the repeating regions and their association toward higher order arrangements. We find that some protein architectures can be described as nearly periodic, while in others clear separations between repetitions exist. Since the method is independent of amino acid sequence information, we can identify structural units that can be encoded by a variety of distinct amino acid sequences.

Tandem-repeat proteins: regularity plus modularity equals design-ability

Researchers in the field of rational protein design face a significant challenge, which arises from the two defining and inter-related features of typical globular protein structures, namely topological complexity and cooperativity. In striking contrast to globular proteins, tandem repeat proteins, such as ankyrin, tetratricopeptide and leucine-rich repeats, have regular, modular, linearly arrayed structures which makes it especially straightforward to dissect and redesign their properties. Here we review what we have learnt about the biophysics of natural repeat proteins and recent progress in applying that knowledge to engineer the thermodynamics, folding pathways and molecular recognition properties of tandem repeat proteins, and we discuss the wealth of possibilities presented for the extension of this modular construction process to build new molecules for use in medicine and biotechnology.

Single-molecule FRET reveals the native-state dynamics of the IκBα ankyrin repeat domain

Previous single-molecule fluorescence resonance energy transfer (smFRET) studies in which the second and sixth ankyrin repeats (ARs) of IκBα were labeled with FRET pairs showed slow fluctuations as if the IκBα AR domain was unfolding in its native state. To systematically probe where these slow dynamic fluctuations occur, we now present data from smFRET studies wherein FRET labels were placed at ARs 1 and 4 (mutant named AR 1-4), at ARs 2 and 5 (AR 2-5), and at ARs 3 and 6 (AR 3-6). The results presented here reveal that AR 6 most readily detaches/unfolds from the AR domain, undergoing substantial fluctuations at room temperature. AR 6 has fewer stabilizing consensus residues than the other IκBα ARs, probably contributing to the ease with which AR 6 "loses grip". AR 5 shows almost no fluctuations at room temperature, but a significant fraction of molecules shows fluctuations at 37 °C. Introduction of stabilizing mutations that are known to fold AR 6 dampen the fluctuations of AR 5, indicating that the AR 5 fluctuations are likely due to weakened inter-repeat stabilization from AR 6. AR 1 also fluctuates somewhat at room temperature, suggesting that fluctuations are a general behavior of ARs at ends of AR domains. Remarkably, AR 1 still fluctuates in the bound state, but mainly between 0.6 and 0.9 FRET efficiency, whereas in the free IκBα, the fluctuations extend to <0.5 FRET efficiency. Overall, our results provide a more complete picture of the energy landscape of the native state dynamics of an AR domain.

Long-range effects and functional consequences of stabilizing mutations in the ankyrin repeat domain of IκBα

Protein domains containing three or more ankyrin repeats (ARs) are ubiquitous in all phyla. Sequence alignments previously identified certain conserved positions, which have been shown to stabilize AR domains and promote their folding. Consensus mutations [Y254L/T257A (YLTA) and C186P/A220P (CPAP)] stabilize the naturally occuring AR domain of human IκBα to denaturation; however, only the YLTA mutations stabilize the protein to proteasomal degradation. We present results from NMR experiments designed to probe the roles of these consensus mutations in IκBα. According to residual dipolar coupling analysis, the gross structures of the AR domains of both mutants appear to be similar to the wild type (WT). Comparison of chemical shifts of mutant and WT proteins reveals that the YLTA and CPAP consensus mutations cause unexpected long-range effects throughout the AR domains. Backbone dynamics experiments reveal that the YLTA mutations in the sixth AR order the C-terminal PEST sequence on the picosecond-to-nanosecond timescale, compared to either the WT or the CPAP mutant IκBαs. This property is likely the mechanism by which the half-life of YLTA IκBα is extended in vivo.

Discrete kinetic models from funneled energy landscape simulations

A general method for facilitating the interpretation of computer simulations of protein folding with minimally frustrated energy landscapes is detailed and applied to a designed ankyrin repeat protein (4ANK). In the method, groups of residues are assigned to foldons and these foldons are used to map the conformational space of the protein onto a set of discrete macrobasins. The free energies of the individual macrobasins are then calculated, informing practical kinetic analysis. Two simple assumptions about the universality of the rate for downhill transitions between macrobasins and the natural local connectivity between macrobasins lead to a scheme for predicting overall folding and unfolding rates, generating chevron plots under varying thermodynamic conditions, and inferring dominant kinetic folding pathways. To illustrate the approach, free energies of macrobasins were calculated from biased simulations of a non-additive structure-based model using two structurally motivated foldon definitions at the full and half ankyrin repeat resolutions. The calculated chevrons have features consistent with those measured in stopped flow chemical denaturation experiments. The dominant inferred folding pathway has an “inside-out”, nucleation-propagation like character.

Equilibrium and kinetic studies of protein cooperativity using urea-induced folding/unfolding of a Ubq-UIM fusion protein

Understanding the origins of cooperativity in proteins remains an important topic in protein folding. This study describes experimental folding/unfolding equilibrium and kinetic studies of the engineered protein Ubq-UIM, consisting of ubiquitin (Ubq) fused to the sequence of the ubiquitin interacting motif (UIM) via a short linker. Urea-induced folding/unfolding profiles of Ubq-UIM were monitored by far-UV circular dichroism and fluorescence spectroscopies and compared to those of the isolated Ubq domain. It was found that the equilibrium data for Ubq-UIM is inconsistent with a two-state model. Analysis of the kinetics of folding shows similarity in the folding transition state ensemble between Ubq and Ubq-UIM, suggesting that formation of Ubq domain is independent of UIM. The major contribution to the stabilization of Ubq-UIM, relative to Ubq, was found to be in the rates of unfolding. Moreover, it was found that the kinetic m-values for Ubq-UIM unfolding, monitored by different probes (far-UV circular dichroism and fluorescence spectroscopies), are different; thereby, further supporting deviations from a two-state behavior. A thermodynamic linkage model that involves four states was found to be applicable to the urea-induced unfolding of Ubq-UIM, which is in agreement with the previous temperature-induced unfolding study. The applicability of the model was further supported by site-directed variants of Ubq-UIM that have altered stabilities of Ubq/UIM interface and/or stabilities of individual Ubq- and UIM-domains. All variants show increased cooperativity and one variant, E43N_Ubq-UIM, appears to behave very close to an equilibrium two-state.