Mapping the Topography of a Protein Energy Landscape

Snapshots of urea-induced early structural changes and unfolding of an ankyrin repeat protein at atomic resolution

Protein folding and unfolding is a complex process, underscored by the many proteotoxic diseases associated with misfolded proteins. Mapping pathways from a native structure to an unfolded protein or vice versa, identifying the intermediates, and defining the role of sequence and structure en route remain outstanding problems in the field. It is even more challenging to capture the events at atomistic resolution. X-ray diffraction has so far been used to understand how urea interacts with and unfolds two stable globular proteins. Here, we present the case study on PSMD10Gankyrin , a prototype for a moderately stable, non-globular repeat protein, long and rigid, with its termini located at either end.   We define structural changes in the time dimension using low urea concentrations to arrive at the following conclusions. (a) Unfolding is rapidly initiated at the C-terminus, slowly at the N-terminus, and proceeds inwards from both ends. (b) C-terminus undergoes an α to 310 helix transition, representing the structure of a potential early unfolding intermediate before disorder sets in. (c) Distinct and progressive changes in the electrostatic landscape of PSMD10Gankyrin were observed, indicative of conformational changes in the seemingly inflexible motif involved in protein-protein interaction. We believe this unique study will open up the field for better and bolder queries and increase the choice of model proteins for a better understanding of the challenging problems of protein folding, protein interactions, protein degradation, and diseases associated with misfolding.

The Wako-Saitô-Muñoz-Eaton Model for Predicting Protein Folding and Dynamics

Despite the recent advances in the prediction of protein structures by deep neutral networks, the elucidation of protein-folding mechanisms remains challenging. A promising theory for describing protein folding is a coarse-grained statistical mechanical model called the Wako-Saitô-Muñoz-Eaton (WSME) model. The model can calculate the free-energy landscapes of proteins based on a three-dimensional structure with low computational complexity, thereby providing a comprehensive understanding of the folding pathways and the structure and stability of the intermediates and transition states involved in the folding reaction. In this review, we summarize previous and recent studies on protein folding and dynamics performed using the WSME model and discuss future challenges and prospects. The WSME model successfully predicted the folding mechanisms of small single-domain proteins and the effects of amino-acid substitutions on protein stability and folding in a manner that was consistent with experimental results. Furthermore, extended versions of the WSME model were applied to predict the folding mechanisms of multi-domain proteins and the conformational changes associated with protein function. Thus, the WSME model may contribute significantly to solving the protein-folding problem and is expected to be useful for predicting protein folding, stability, and dynamics in basic research and in industrial and medical applications.

Multivalent Interaction of Beta-Catenin With its Intrinsically Disordered Binding Partner Adenomatous Polyposis Coli

The Wnt signalling pathway plays key roles in cell proliferation, differentiation and fate decisions in embryonic development and maintenance of adult tissues, and the twelve Armadillo (ARM) repeat-containing protein β-catenin acts as the signal transducer in this pathway. Here we investigate the interaction between β-catenin's ARM repeat domain and the intrinsically disordered protein adenomatous polyposis coli (APC). APC is a giant multivalent scaffold that brings together the different components of the so-called "β-catenin destruction complex", which drives β-catenin degradation via the ubiquitin-proteasome pathway. Mutations and truncations in APC, resulting in loss of APC function and hence elevated β-catenin levels and upregulation of Wnt signalling, are associated with numerous cancers including colorectal carcinomas. APC has a long intrinsically disordered region (IDR) that contains a series of 15-residue and 20-residue binding regions for β-catenin. Here we explore the multivalent nature of the interaction of β-catenin with the highest affinity APC repeat, both at equilibrium and under kinetic conditions. We use a combination of single-site substitutions, deletions and insertions to dissect the mechanism of molecular recognition and the roles of the three β-catenin-binding subdomains of APC.

Heterogeneity in Protein Folding and Unfolding Reactions

Proteins have dynamic structures that undergo chain motions on time scales spanning from picoseconds to seconds. Resolving the resultant conformational heterogeneity is essential for gaining accurate insight into fundamental mechanistic aspects of the protein folding reaction. The use of high-resolution structural probes, sensitive to population distributions, has begun to enable the resolution of site-specific conformational heterogeneity at different stages of the folding reaction. Different states populated during protein folding, including the unfolded state, collapsed intermediate states, and even the native state, are found to possess significant conformational heterogeneity. Heterogeneity in protein folding and unfolding reactions originates from the reduced cooperativity of various kinds of physicochemical interactions between various structural elements of a protein, and between a protein and solvent. Heterogeneity may arise because of functional or evolutionary constraints. Conformational substates within the unfolded state and the collapsed intermediates that exchange at rates slower than the subsequent folding steps give rise to heterogeneity on the protein folding pathways. Multiple folding pathways are likely to represent distinct sequences of structure formation. Insight into the nature of the energy barriers separating different conformational states populated during (un)folding can also be obtained by resolving heterogeneity.

A hierarchy of coupling free energies underlie the thermodynamic and functional architecture of protein structures

Protein sequences and structures evolve by satisfying varied physical and biochemical constraints. This multi-level selection is enabled not just by the patterning of amino acids on the sequence, but also via coupling between residues in the native structure. Here, we employ an energetically detailed statistical mechanical model with millions of microstates to extract such long-range structural correlations, i.e. thermodynamic coupling free energies, from a diverse family of protein structures. We find that despite the intricate and anisotropic distribution of coupling patterns, the majority of residues (>70%) are only marginally coupled contributing to functional motions and catalysis. Physical origins of ‘sectors’, determinants of native ensemble heterogeneity in extant, ancient and designed proteins, and the basis for allostery emerge naturally from coupling free energies. The statistical framework highlights how evolutionary selection and optimization occur at the level of global interaction network for a given protein fold impacting folding, function, and allosteric outputs.

•

Evolution of protein structures occurs at the level of global interaction network.

•

More than 70% of the protein residues are weakly or marginally coupled.

•

Functional ‘sector’ regions are a manifestation of marginal coupling.

•

Coupling indices vary across the entire proteins in extant-ancient and natural-designed pairs.

•

The proposed methodology can be used to understand allostery and epistasis.

Emergence of effective temperatures in an out-of-equilibrium model of biopolymer folding

We investigate the possibility of extending the notion of temperature in a stochastic model for the RNA or protein folding driven out of equilibrium. We simulate the dynamics of a small RNA hairpin subject to an external pulling force, which is time-dependent. First, we consider a fluctuation-dissipation relation (FDR) whereby we verify that various effective temperatures can be obtained for different observables, only when the slowest intrinsic relaxation timescale of the system regulates the dynamics of the system. Then, we introduce a different nonequilibrium temperature, which is defined from the rate of heat exchanged with a weakly interacting thermal bath. Notably, this "kinetic" temperature can be defined for any frequency of the external switching force. We also discuss and compare the behavior of these two emerging parameters, by discriminating the time-delayed nature of the FDR temperature from the instantaneous character of the kinetic temperature. The validity of our numerics are corroborated by a simple four-state Markov model which describes the long-time behavior of the RNA molecule.

The Effect of Mutations in the TPR and Ankyrin Families of Alpha Solenoid Repeat Proteins

Protein repeats are short, highly similar peptide motifs that occur several times within a single protein, for example the TPR and Ankyrin repeats. Understanding the role of mutation in these proteins is complicated by the competing facts that 1) the repeats are much more restricted to a set sequence than non-repeat proteins, so mutations should be harmful much more often because there are more residues that are heavily restricted due to the need of the sequence to repeat and 2) the symmetry of the repeats in allows the distribution of functional contributions over a number of residues so that sometimes no specific site is singularly responsible for function (unlike enzymatic active site catalytic residues). To address this issue, we review the effects of mutations in a number of natural repeat proteins from the tetratricopeptide and Ankyrin repeat families. We find that mutations are context dependent. Some mutations are indeed highly disruptive to the function of the protein repeats while mutations in identical positions in other repeats in the same protein have little to no effect on structure or function.

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.

Parallel and Sequential Pathways of Molecular Recognition of a Tandem-Repeat Protein and Its Intrinsically Disordered Binding Partner

The Wnt signalling pathway plays an important role in cell proliferation, differentiation, and fate decisions in embryonic development and the maintenance of adult tissues. The twelve armadillo (ARM) repeat-containing protein β-catenin acts as the signal transducer in this pathway. Here, we investigated the interaction between β-catenin and the intrinsically disordered transcription factor TCF7L2, comprising a very long nanomolar-affinity interface of approximately 4800 Å2 that spans ten of the twelve ARM repeats of β-catenin. First, a fluorescence reporter system for the interaction was engineered and used to determine the kinetic rate constants for the association and dissociation. The association kinetics of TCF7L2 and β-catenin were monophasic and rapid (7.3 ± 0.1 × 107 M-1·s-1), whereas dissociation was biphasic and slow (5.7 ± 0.4 × 10-4 s-1, 15.2 ± 2.8 × 10-4 s-1). This reporter system was then combined with site-directed mutagenesis to investigate the striking variability in the conformation adopted by TCF7L2 in the three different crystal structures of the TCF7L2-β-catenin complex. We found that the mutation had very little effect on the association kinetics, indicating that most interactions form after the rate-limiting barrier for association. Mutations of the N- and C-terminal subdomains of TCF7L2 that adopt relatively fixed conformations in the crystal structures had large effects on the dissociation kinetics, whereas the mutation of the labile sub-domain connecting them had negligible effect. These results point to a two-site avidity mechanism of binding with the linker region forming a "fuzzy" complex involving transient contacts that are not site-specific. Strikingly, the two mutations in the N-terminal subdomain that had the largest effects on the dissociation kinetics showed two additional phases, indicating partial flux through an alternative dissociation pathway that is inaccessible to the wild type. The results presented here provide insights into the kinetics of the molecular recognition of a long intrinsically disordered region with an elongated repeat-protein surface, a process found to involve parallel routes with sequential steps in each.

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.

Application of a transparent window vibrational probe (azido probe) to the structural dynamics of model dipeptides and amyloid β-peptide

The azido asymmetric stretching motion is widely used for the elucidation of the intrinsic conformational preference and folding mechanism of protein since it has strong vibrational absorbance in the spectral transparent windows. However, the possible secondary structural disturbance induced by the insertion of azido group in the side chain of polypeptides should be carefully evaluated. Here, DFT calculation and enhanced sampling method were employed for model dipeptides with or without azido substitution, and the outcome results show that the lower potential energy basins of isolated model dipeptides are consistent with the preferred structural distributions of model dipeptides in aqueous solution. The azido asymmetric stretching frequency shows its sensitivity to the backbone configurations just like amide-I vibration does, and the azido vibration exhibits great potential as a structural reporter in the transparent window. For the evaluation of the application of azido group in biologically related system, the structural dynamics of Aβ_37-42 and N₃-Aβ_37-42 fragments and the self-assemble process of their protofiliments in aqueous solution were demonstrated. The outcome results show that the structural fluctuations of Aβ_37-42 and its protofilament in aqueous solution are quite similar with or without azido substitution, and the dewetting transitions of Aβ_37-42 and N₃-Aβ_37-42 β-sheet layers are both complete within 30 ns and assemble into stable protofilaments. Therefore, the azido asymmetric vibrational motion is a minimally invasive structural probe and would not introduce much disturbance to the structural dynamics of polypeptides.

Decoupling a tandem-repeat protein: Impact of multiple loop insertions on a modular scaffold

The simple topology and modular architecture of tandem-repeat proteins such as tetratricopeptide repeats (TPRs) and ankyrin repeats makes them straightforward to dissect and redesign. Repeat-protein stability can be manipulated in a predictable way using site-specific mutations. Here we explore a different type of modification - loop insertion - that will enable a simple route to functionalisation of this versatile scaffold. We previously showed that a single loop insertion has a dramatically different effect on stability depending on its location in the repeat array. Here we dissect this effect by a combination of multiple and alternated loop insertions to understand the origins of the context-dependent loss in stability. We find that the scaffold is remarkably robust in that its overall structure is maintained. However, adjacent repeats are now only weakly coupled, and consequently the increase in solvent protection, and thus stability, with increasing repeat number that defines the tandem-repeat protein class is lost. Our results also provide us with a rulebook with which we can apply these principles to the design of artificial repeat proteins with precisely tuned folding landscapes and functional capabilities, thereby paving the way for their exploitation as a versatile and truly modular platform in synthetic biology.

Thermodynamics and folding landscapes of large proteins from a statistical mechanical model

Statistical mechanical models that afford an intermediate resolution between macroscopic chemical models and all-atom simulations have been successful in capturing folding behaviors of many small single-domain proteins. However, the applicability of one such successful approach, the Wako-Saitô-Muñoz-Eaton (WSME) model, is limited by the size of the protein as the number of conformations grows exponentially with protein length. In this work, we surmount this size limitation by introducing a novel approximation that treats stretches of 3 or 4 residues as blocks, thus reducing the phase space by nearly three orders of magnitude. The performance of the 'bWSME' model is validated by comparing the predictions for a globular enzyme (RNase H) and a repeat protein (IκBα), against experimental observables and the model without block approximation. Finally, as a proof of concept, we predict the free-energy surface of the 370-residue, multi-domain maltose binding protein and identify an intermediate in good agreement with single-molecule force-spectroscopy measurements. The bWSME model can thus be employed as a quantitative predictive tool to explore the conformational landscapes of large proteins, extract the structural features of putative intermediates, identify parallel folding paths, and thus aid in the interpretation of both ensemble and single-molecule experiments.

Context-Dependent Energetics of Loop Extensions in a Family of Tandem-Repeat Proteins

Consensus-designed tetratricopeptide repeat proteins are highly stable, modular proteins that are strikingly amenable to rational engineering. They therefore have tremendous potential as building blocks for biomaterials and biomedicine. Here, we explore the possibility of extending the loops between repeats to enable further diversification, and we investigate how this modification affects stability and folding cooperativity. We find that extending a single loop by up to 25 residues does not disrupt the overall protein structure, but, strikingly, the effect on stability is highly context-dependent: in a two-repeat array, destabilization is relatively small and can be accounted for purely in entropic terms, whereas extending a loop in the middle of a large array is much more costly because of weakening of the interaction between the repeats. Our findings provide important and, to our knowledge, new insights that increase our understanding of the structure, folding, and function of natural repeat proteins and the design of artificial repeat proteins in biotechnology.

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.

Extracting the Hidden Distributions Underlying the Mean Transition State Structures in Protein Folding

The inherent conflict between noncovalent interactions and the large conformational entropy of the polypeptide chain forces folding reactions and their mechanisms to deviate significantly from chemical reactions. Accordingly, measures of structure in the transition state ensemble (TSE) are strongly influenced by the underlying distributions of microscopic folding pathways that are challenging to discern experimentally. Here, we present a detailed analysis of 150,000 folding transition paths of five proteins at three different thermodynamic conditions from an experimentally consistent statistical mechanical model. We find that the underlying TSE structural distributions are rarely unimodal, and the average experimental measures arise from complex underlying distributions. Unfolding pathways also exhibit subtle differences from folding counterparts due to a combination of Hammond behavior and native-state movements. Local interactions and topological complexity, to a lesser extent, are found to determine pathway heterogeneity, underscoring the importance of the balance between local and nonlocal energetics in protein folding.

Toward a quantitative description of microscopic pathway heterogeneity in protein folding

Experimentally consistent statistical modeling of protein folding thermodynamics reveals unprecedented complexity with numerous parallel folding routes in five different proteins.

How cooperative are protein folding and unfolding transitions?

A thermodynamically and kinetically simple picture of protein folding envisages only two states, native (N) and unfolded (U), separated by a single activation free energy barrier, and interconverting by cooperative two-state transitions. The folding/unfolding transitions of many proteins occur, however, in multiple discrete steps associated with the formation of intermediates, which is indicative of reduced cooperativity. Furthermore, much advancement in experimental and computational approaches has demonstrated entirely non-cooperative (gradual) transitions via a continuum of states and a multitude of small energetic barriers between the N and U states of some proteins. These findings have been instrumental towards providing a structural rationale for cooperative versus noncooperative transitions, based on the coupling between interaction networks in proteins. The cooperativity inherent in a folding/unfolding reaction appears to be context dependent, and can be tuned via experimental conditions which change the stabilities of N and U. The evolution of cooperativity in protein folding transitions is linked closely to the evolution of function as well as the aggregation propensity of the protein. A large activation energy barrier in a fully cooperative transition can provide the kinetic control required to prevent the accumulation of partially unfolded forms, which may promote aggregation. Nevertheless, increasing evidence for barrier-less "downhill" folding, as well as for continuous "uphill" unfolding transitions, indicate that gradual non-cooperative processes may be ubiquitous features on the free energy landscape of protein folding.