The mechanism of folding of Im7 reveals competition between functional and kinetic evolutionary constraints

Frustraevo: a web server to localize and quantify the conservation of local energetic frustration in protein families

According to the Principle of Minimal Frustration, folded proteins can only have a minimal number of strong energetic conflicts in their native states. However, not all interactions are energetically optimized for folding but some remain in energetic conflict, i.e. they are highly frustrated. This remaining local energetic frustration has been shown to be statistically correlated with distinct functional aspects such as protein-protein interaction sites, allosterism and catalysis. Fuelled by the recent breakthroughs in efficient protein structure prediction that have made available good quality models for most proteins, we have developed a strategy to calculate local energetic frustration within large protein families and quantify its conservation over evolutionary time. Based on this evolutionary information we can identify how stability and functional constraints have appeared at the common ancestor of the family and have been maintained over the course of evolution. Here, we present FrustraEvo, a web server tool to calculate and quantify the conservation of local energetic frustration in protein families.

Thermodynamic and kinetic analysis of the LAO binding protein and its isolated domains reveal non-additivity in stability, folding and function

Substrate-binding proteins (SBPs) are used by organisms from the three domains of life for transport and signalling. SBPs are composed of two domains that collectively trap ligands with high affinity and selectivity. To explore the role of the domains and the integrity of the hinge region between them in the function and conformation of SBPs, here, we describe the ligand binding, conformational stability and folding kinetics of the Lysine Arginine Ornithine (LAO) binding protein from Salmonella thiphimurium and constructs corresponding to its two independent domains. LAO is a class II SBP formed by a continuous and a discontinuous domain. Contrary to the expected behaviour based on their connectivity, the discontinuous domain shows a stable native-like structure that binds l-arginine with moderate affinity, whereas the continuous domain is barely stable and shows no detectable ligand binding. Regarding folding kinetics, studies of the entire protein revealed the presence of at least two intermediates. While the unfolding and refolding of the continuous domain exhibited only a single intermediate and simpler and faster kinetics than LAO, the folding mechanism of the discontinuous domain was complex and involved multiple intermediates. These findings suggest that in the complete protein the continuous domain nucleates folding and that its presence funnels the folding of the discontinuous domain avoiding nonproductive interactions. The strong dependence of the function, stability and folding pathway of the lobes on their covalent association is most likely the result of the coevolution of both domains as a single unit.

The Right-Handed Parallel β-Helix Topology of Erwinia chrysanthemi Pectin Methylesterase Is Intimately Associated with Both Sequential Folding and Resistance to High Pressure

The complex topologies of large multi-domain globular proteins make the study of their folding and assembly particularly demanding. It is often characterized by complex kinetics and undesired side reactions, such as aggregation. The structural simplicity of tandem-repeat proteins, which are characterized by the repetition of a basic structural motif and are stabilized exclusively by sequentially localized contacts, has provided opportunities for dissecting their folding landscapes. In this study, we focus on the Erwinia chrysanthemi pectin methylesterase (342 residues), an all-β pectinolytic enzyme with a right-handed parallel β-helix structure. Chemicals and pressure were chosen as denaturants and a variety of optical techniques were used in conjunction with stopped-flow equipment to investigate the folding mechanism of the enzyme at 25 °C. Under equilibrium conditions, both chemical- and pressure-induced unfolding show two-state transitions, with average conformational stability (ΔG° = 35 ± 5 kJ·mol⁻¹) but exceptionally high resistance to pressure (P_m = 800 ± 7 MPa). Stopped-flow kinetic experiments revealed a very rapid (τ < 1 ms) hydrophobic collapse accompanied by the formation of an extended secondary structure but did not reveal stable tertiary contacts. This is followed by three distinct cooperative phases and the significant population of two intermediate species. The kinetics followed by intrinsic fluorescence shows a lag phase, strongly indicating that these intermediates are productive species on a sequential folding pathway, for which we propose a plausible model. These combined data demonstrate that even a large repeat protein can fold in a highly cooperative manner.

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.

Effects of Topology and Sequence in Protein Folding Linked via Conformational Fluctuations

Experiments have compared the folding of proteins with different amino acid sequences but the same basic structure, or fold. Results indicate that folding is robust to sequence variations for proteins with some nonlocal folds, such as all-β, whereas the folding of more local, all-α proteins typically exhibits a stronger sequence dependence. Here, we use a coarse-grained model to systematically study how variations in sequence perturb the folding energy landscapes of three model sequences with 3α, 4β + α, and β-barrel folds, respectively. These three proteins exhibit folding features in line with experiments, including expected rank order in the cooperativity of the folding transition and stability-dependent shifts in the location of the free-energy barrier to folding. Using a generalized-ensemble simulation approach, we determine the thermodynamics of around 2000 sequence variants representing all possible hydrophobic or polar single- and double-point mutations. From an analysis of the subset of stability-neutral mutations, we find that folding is perturbed in a topology-dependent manner, with the β-barrel protein being the most robust. Our analysis shows, in particular, that the magnitude of mutational perturbations of the transition state is controlled in part by the size or "width" of the underlying conformational ensemble. This result suggests that the mutational robustness of the folding of the β-barrel protein is underpinned by its conformationally restricted transition state ensemble, revealing a link between sequence and topological effects in protein folding.

Conformational flexibility within the nascent polypeptide-associated complex enables its interactions with structurally diverse client proteins

As newly synthesized polypeptides emerge from the ribosome, it is crucial that they fold correctly. To prevent premature aggregation, nascent chains interact with chaperones that facilitate folding or prevent misfolding until protein synthesis is complete. Nascent polypeptide-associated complex (NAC) is a ribosome-associated chaperone that is important for protein homeostasis. However, how NAC binds its substrates remains unclear. Using native electrospray ionization MS (ESI-MS), limited proteolysis, NMR, and cross-linking, we analyzed the conformational properties of NAC from Caenorhabditis elegans and studied its ability to bind proteins in different conformational states. Our results revealed that NAC adopts an array of compact and expanded conformations and binds weakly to client proteins that are unfolded, folded, or intrinsically disordered, suggestive of broad substrate compatibility. Of note, we found that this weak binding retards aggregation of the intrinsically disordered protein α-synuclein both in vitro and in vivo These findings provide critical insights into the structure and function of NAC. Specifically, they reveal the ability of NAC to exploit its conformational plasticity to bind a repertoire of substrates with unrelated sequences and structures, independently of actively translating ribosomes.

Molecular recognition and packing frustration in a helical protein

Biomolecular recognition entails attractive forces for the functional native states and discrimination against potential nonnative interactions that favor alternate stable configurations. The challenge posed by the competition of nonnative stabilization against native-centric forces is conceptualized as frustration. Experiment indicates that frustration is often minimal in evolved biological systems although nonnative possibilities are intuitively abundant. Much of the physical basis of minimal frustration in protein folding thus remains to be elucidated. Here we make progress by studying the colicin immunity protein Im9. To assess the energetic favorability of nonnative versus native interactions, we compute free energies of association of various combinations of the four helices in Im9 (referred to as H1, H2, H3, and H4) by extensive explicit-water molecular dynamics simulations (total simulated time > 300 μs), focusing primarily on the pairs with the largest native contact surfaces, H1-H2 and H1-H4. Frustration is detected in H1-H2 packing in that a nonnative packing orientation is significantly stabilized relative to native, whereas such a prominent nonnative effect is not observed for H1-H4 packing. However, in contrast to the favored nonnative H1-H2 packing in isolation, the native H1-H2 packing orientation is stabilized by H3 and loop residues surrounding H4. Taken together, these results showcase the contextual nature of molecular recognition, and suggest further that nonnative effects in H1-H2 packing may be largely avoided by the experimentally inferred Im9 folding transition state with native packing most developed at the H1-H4 rather than the H1-H2 interface.

Biomolecules need to recognize one another with high specificity: promoting “native” functional intermolecular binding events while avoiding detrimental “nonnative” bound configurations; i.e., “frustration”—the tendency for nonnative interactions—has to be minimized. Folding of globular proteins entails a similar discrimination. To gain physical insight, we computed the binding affinities of helical structures of the protein Im9 in various native or nonnative configurations by atomic simulations, discovering that partial packing of the Im9 core is frustrated. This frustration is overcome when the entire core of the protein is assembled, consistent with experiment indicating no significant kinetic trapping in Im9 folding. Our systematic analysis thus reveals a subtle, contextual aspect of biomolecular recognition and provides a general approach to characterize folding frustration.

Monitoring protein folding through high pressure NMR spectroscopy

High-pressure is a well-known perturbation method used to destabilize globular proteins. It is perfectly reversible, which is essential for a proper thermodynamic characterization of a protein equilibrium. In contrast to other perturbation methods such as heat or chemical denaturant that destabilize protein structures uniformly, pressure exerts local effects on regions or domains of a protein containing internal cavities. When combined with NMR spectroscopy, hydrostatic pressure offers the possibility to monitor at a residue level the structural transitions occurring upon unfolding and to determine the kinetic properties of the process. High-pressure NMR experiments can now be routinely performed, owing to the recent development of commercially available high-pressure sample cells. This review summarizes recent advances and some future directions of high-pressure NMR techniques for the characterization at atomic resolution of the energy landscape of protein folding.

Folding while bound to chaperones

Chaperones are important in preventing protein aggregation and aiding protein folding. How chaperones aid protein folding remains a key question in understanding their mechanism. The possibility of proteins folding while bound to chaperones was reintroduced recently with the chaperone Spy, many years after the phenomenon was first reported with the chaperones GroEL and SecB. In this review, we discuss the salient features of folding while bound in the cases for which it has been observed and speculate about its biological importance and possible occurrence in other chaperones.

A molecular mechanism of chaperone-client recognition

Molecular chaperones are essential in aiding client proteins to fold into their native structure and in maintaining cellular protein homeostasis. However, mechanistic aspects of chaperone function are still not well understood at the atomic level. We use nuclear magnetic resonance spectroscopy to elucidate the mechanism underlying client recognition by the adenosine triphosphate-independent chaperone Spy at the atomic level and derive a structural model for the chaperone-client complex. Spy interacts with its partially folded client Im7 by selective recognition of flexible, locally frustrated regions in a dynamic fashion. The interaction with Spy destabilizes a partially folded client but spatially compacts an unfolded client conformational ensemble. By increasing client backbone dynamics, the chaperone facilitates the search for the native structure. A comparison of the interaction of Im7 with two other chaperones suggests that the underlying principle of recognizing frustrated segments is of a fundamental nature.

Folding Mechanism of Proteins Im7 and Im9: Insight from All-Atom Simulations in Implicit and Explicit Solvent

Im7 and Im9 are evolutionary related proteins with almost identical native structures. In spite of their structural similarity, experiments show that Im7 folds through a long-lived on-pathway intermediate, while Im9 folds according to two-state kinetics. In this work, we use a recently developed enhanced path sampling method to generate many folding trajectories for these proteins, using realistic atomistic force fields, in both implicit and explicit solvent. Overall, our results are in good agreement with the experimental ϕ values and with the result of ϕ-value-restrained molecular dynamics (MD) simulations. However, our implicit solvent simulations fail to predict a qualitative difference in the folding pathways of Im7 and Im9. In contrast, our simulations in explicit solvent correctly reproduce the fact that only protein Im7 folds through a on-pathway intermediate. By analyzing our atomistic trajectories, we provide a physical picture which explains the observed difference in the folding kinetics of these chains.

Visualizing chaperone-assisted protein folding

Challenges in determining the structures of heterogeneous and dynamic protein complexes have greatly hampered past efforts to obtain a mechanistic understanding of many important biological processes. One such process is chaperone-assisted protein folding. Obtaining structural ensembles of chaperone-substrate complexes would ultimately reveal how chaperones help proteins fold into their native state. To address this problem, we devised a new structural biology approach based on X-ray crystallography, termed residual electron and anomalous density (READ). READ enabled us to visualize even sparsely populated conformations of the substrate protein immunity protein 7 (Im7) in complex with the Escherichia coli chaperone Spy, and to capture a series of snapshots depicting the various folding states of Im7 bound to Spy. The ensemble shows that Spy-associated Im7 samples conformations ranging from unfolded to partially folded to native-like states and reveals how a substrate can explore its folding landscape while being bound to a chaperone.

Forces Driving Chaperone Action

It is still unclear what molecular forces drive chaperone-mediated protein folding. Here, we obtain a detailed mechanistic understanding of the forces that dictate the four key steps of chaperone-client interaction: initial binding, complex stabilization, folding, and release. Contrary to the common belief that chaperones recognize unfolding intermediates by their hydrophobic nature, we discover that the model chaperone Spy uses long-range electrostatic interactions to rapidly bind to its unfolded client protein Im7. Short-range hydrophobic interactions follow, which serve to stabilize the complex. Hydrophobic collapse of the client protein then drives its folding. By burying hydrophobic residues in its core, the client's affinity to Spy decreases, which causes client release. By allowing the client to fold itself, Spy circumvents the need for client-specific folding instructions. This mechanism might help explain how chaperones can facilitate the folding of various unrelated proteins.

Using the folding landscapes of proteins to understand protein function

Proteins fold on a biologically-relevant timescale because of a funnel-shaped energy landscape. This landscape is sculpted through evolution by selecting amino-acid sequences that stabilize native interactions while suppressing stable non-native interactions that occur during folding. However, there is strong evolutionary selection for functional residues and these cannot be chosen to optimize folding. Their presence impacts the folding energy landscape in a variety of ways. Here, we survey the effects of functional residues on folding by providing several examples. We then review how such effects can be detected computationally and be used as assays for protein function. Overall, an understanding of how functional residues modulate folding should provide insights into the design of natural proteins and their homeostasis.

Substrate protein folds while it is bound to the ATP-independent chaperone Spy

Chaperones assist the folding of many proteins in the cell. While the most well studied chaperones use cycles of ATP binding and hydrolysis to assist protein folding, a number of chaperones have been identified that promote protein folding in the absence of high-energy cofactors. Precisely how ATP-independent chaperones accomplish this feat is unclear. Here we have characterized the kinetic mechanism of substrate folding by the small, ATP-independent chaperone, Spy. Spy rapidly associates with its substrate, Immunity protein 7 (Im7), eliminating its potential for aggregation. Remarkably, Spy then allows Im7 to fully fold into its native state while remaining bound to the surface of the chaperone. These results establish a potentially widespread mechanism whereby ATP-independent chaperones can assist in protein refolding. They also provide compelling evidence that substrate proteins can fold while continuously bound to a chaperone.

Using hydroxyl radical footprinting to explore the free energy landscape of protein folding

Characterisation of the conformational states adopted during protein folding, including globally unfolded/disordered structures and partially folded intermediate species, is vital to gain fundamental insights into how a protein folds. In this work we employ fast photochemical oxidation of proteins (FPOP) to map the structural changes that occur in the folding of the four-helical bacterial immunity protein, Im7. Oxidative footprinting coupled with mass spectrometry (MS) is used to probe changes in the solvent accessibility of amino acid side-chains concurrent with the folding process, by quantifying the degree of oxidation experienced by the wild-type protein relative to a kinetically trapped, three-helical folding intermediate and an unfolded variant that lacks secondary structure. Analysis of the unfolded variant by FPOP-MS shows oxidative modifications consistent with the species adopting a solution conformation with a high degree of solvent accessibility. The folding intermediate, by contrast, experiences increased levels of oxidation relative to the wild-type, native protein only in regions destabilised by the amino acid substitutions introduced. The results demonstrate the utility of FPOP-MS to characterise protein variants in different conformational states and to provide insights into protein folding mechanisms that are complementary to measurements such as hydrogen/deuterium exchange labelling and Φ-value analysis.

Native contact density and nonnative hydrophobic effects in the folding of bacterial immunity proteins

The bacterial colicin-immunity proteins Im7 and Im9 fold by different mechanisms. Experimentally, at pH 7.0 and 10°C, Im7 folds in a three-state manner via an intermediate but Im9 folding is two-state-like. Accordingly, Im7 exhibits a chevron rollover, whereas the chevron arm for Im9 folding is linear. Here we address the biophysical basis of their different behaviors by using native-centric models with and without additional transferrable, sequence-dependent energies. The Im7 chevron rollover is not captured by either a pure native-centric model or a model augmented by nonnative hydrophobic interactions with a uniform strength irrespective of residue type. By contrast, a more realistic nonnative interaction scheme that accounts for the difference in hydrophobicity among residues leads simultaneously to a chevron rollover for Im7 and an essentially linear folding chevron arm for Im9. Hydrophobic residues identified by published experiments to be involved in nonnative interactions during Im7 folding are found to participate in the strongest nonnative contacts in this model. Thus our observations support the experimental perspective that the Im7 folding intermediate is largely underpinned by nonnative interactions involving large hydrophobics. Our simulation suggests further that nonnative effects in Im7 are facilitated by a lower local native contact density relative to that of Im9. In a one-dimensional diffusion picture of Im7 folding with a coordinate- and stability-dependent diffusion coefficient, a significant chevron rollover is consistent with a diffusion coefficient that depends strongly on native stability at the conformational position of the folding intermediate.

In order to fold correctly, a globular protein must avoid being trapped in wrong, i.e., nonnative conformations. Thus a biophysical account of how attractive nonnative interactions are bypassed by some amino acid sequences but not others is key to deciphering protein structure and function. We examine two closely related bacterial immunity proteins, Im7 and Im9, that are experimentally known to fold very differently: Whereas Im9 folds directly, Im7 folds through a mispacked conformational intermediate. A simple model we developed accounts for their intriguingly different folding kinetics in terms of a balance between the density of native-promoting contacts and the hydrophobicity of local amino acid sequences. This emergent principle is extensible to other biomolecular recognition processes.

Theoretical perspectives on nonnative interactions and intrinsic disorder in protein folding and binding

The diverse biological functions of intrinsically disordered proteins (IDPs) have markedly raised our appreciation of protein conformational versatility, whereas the existence of energetically favorable yet functional detrimental nonnative interactions underscores the physical limitations of evolutionary optimization. Here we survey recent advances in using biophysical modeling to gain insight into experimentally observed nonnative behaviors and IDP properties. Simulations of IDP interactions to date focus mostly on coupled folding-binding, which follows essentially the same organizing principle as the local-nonlocal coupling mechanism in cooperative folding of monomeric globular proteins. By contrast, more innovative theories of electrostatic and aromatic interactions are needed for the conceptually novel but less-explored 'fuzzy' complexes in which the functionally bound IDPs remain largely disordered.

Mechanistic studies of the biogenesis and folding of outer membrane proteins in vitro and in vivo: what have we learned to date?

Research into the mechanisms by which proteins fold into their native structures has been on-going since the work of Anfinsen in the 1960s. Since that time, the folding mechanisms of small, water-soluble proteins have been well characterised. By contrast, progress in understanding the biogenesis and folding mechanisms of integral membrane proteins has lagged significantly because of the need to create a membrane mimetic environment for folding studies in vitro and the difficulties in finding suitable conditions in which reversible folding can be achieved. Improved knowledge of the factors that promote membrane protein folding and disfavour aggregation now allows studies of folding into lipid bilayers in vitro to be performed. Consequently, mechanistic details and structural information about membrane protein folding are now emerging at an ever increasing pace. Using the panoply of methods developed for studies of the folding of water-soluble proteins. This review summarises current knowledge of the mechanisms of outer membrane protein biogenesis and folding into lipid bilayers in vivo and in vitro and discusses the experimental techniques utilised to gain this information. The emerging knowledge is beginning to allow comparisons to be made between the folding of membrane proteins with current understanding of the mechanisms of folding of water-soluble proteins.

Understanding the frustration arising from the competition between function, misfolding, and aggregation in a globular protein

Folding and function may impose different requirements on the amino acid sequences of proteins, thus potentially giving rise to conflict. Such a conflict, or frustration, can result in the formation of partially misfolded intermediates that can compromise folding and promote aggregation. We investigate this phenomenon by studying frataxin, a protein whose normal function is to facilitate the formation of iron-sulfur clusters but whose mutations are associated with Friedreich's ataxia. To characterize the folding pathway of this protein we carry out a Φ-value analysis and use the resulting structural information to determine the structure of the folding transition state, which we then validate by a second round of rationally designed mutagenesis. The analysis of the transition-state structure reveals that the regions involved in the folding process are highly aggregation-prone. By contrast, the regions that are functionally important are partially misfolded in the transition state but highly resistant to aggregation. Taken together, these results indicate that in frataxin the competition between folding and function creates the possibility of misfolding, and that to prevent aggregation the amino acid sequence of this protein is optimized to be highly resistant to aggregation in the regions involved in misfolding.

Frustration in biomolecules

Biomolecules are the prime information processing elements of living matter. Most of these inanimate systems are polymers that compute their own structures and dynamics using as input seemingly random character strings of their sequence, following which they coalesce and perform integrated cellular functions. In large computational systems with finite interaction-codes, the appearance of conflicting goals is inevitable. Simple conflicting forces can lead to quite complex structures and behaviors, leading to the concept of frustration in condensed matter. We present here some basic ideas about frustration in biomolecules and how the frustration concept leads to a better appreciation of many aspects of the architecture of biomolecules, and especially how biomolecular structure connects to function by means of localized frustration. These ideas are simultaneously both seductively simple and perilously subtle to grasp completely. The energy landscape theory of protein folding provides a framework for quantifying frustration in large systems and has been implemented at many levels of description. We first review the notion of frustration from the areas of abstract logic and its uses in simple condensed matter systems. We discuss then how the frustration concept applies specifically to heteropolymers, testing folding landscape theory in computer simulations of protein models and in experimentally accessible systems. Studying the aspects of frustration averaged over many proteins provides ways to infer energy functions useful for reliable structure prediction. We discuss how frustration affects folding mechanisms. We review here how the biological functions of proteins are related to subtle local physical frustration effects and how frustration influences the appearance of metastable states, the nature of binding processes, catalysis and allosteric transitions. In this review, we also emphasize that frustration, far from being always a bad thing, is an essential feature of biomolecules that allows dynamics to be harnessed for function. In this way, we hope to illustrate how Frustration is a fundamental concept in molecular biology.

Misfolding of luciferase at the single-molecule level

The folding of complex proteins can be dramatically affected by misfolding transitions. Directly observing misfolding and distinguishing it from aggregation is challenging. Experiments with optical tweezers revealed transitions between the folded states of a single protein in the absence of mechanical tension. Nonfolded chains of the multidomain protein luciferase folded within seconds to different partially folded states, one of which was stable over several minutes and was more resistant to forced unfolding than other partially folded states. Luciferase monomers can thus adopt a stable misfolded state and can do so without interacting with aggregation partners. This result supports the notion that luciferase misfolding is the cause of the low refolding yields and aggregation observed with this protein. This approach could be used to study misfolding transitions in other large proteins, as well as the factors that affect misfolding.

Conformational dynamics is more important than helical propensity for the folding of the all α-helical protein Im7

Im7 folds via an on-pathway intermediate that contains three of the four native α-helices. The missing helix, helix III, is the shortest and its failure to be formed until late in the pathway is related to frustration in the structure. Im7H3M3, a 94-residue variant of the 87-residue Im7 in which helix III is the longest of the four native helices, also folds via an intermediate. To investigate the structural basis for this we calculated the frustration in the structure of Im7H3M3 and used NMR to investigate its dynamics. We found that the native state of Im7H3M3 is highly frustrated and in equilibrium with an intermediate state that lacks helix III, similar to Im7. Model-free analysis identified residues with chemical exchange contributions to their relaxation that aligned with the residues predicted to have highly frustrated interactions, also like Im7. Finally, we determined properties of urea-denatured Im7H3M3 and identified four clusters of interacting residues that corresponded to the α-helices of the native protein. In Im7 the cluster sizes were related to the lengths of the α-helices with cluster III being the smallest but in Im7H3M3 cluster III was also the smallest, despite this region forming the longest helix in the native state. These results suggest that the conformational properties of the urea-denatured states promote formation of a three-helix intermediate in which the residues that form helix III remain non-helical. Thus it appears that features of the native structure are formed early in folding linked to collapse of the unfolded state.

Cold denaturation of a protein dimer monitored at atomic resolution

Protein folding and unfolding are crucial for a range of biological phenomena and human diseases. Defining the structural properties of the involved transient species is therefore of prime interest. Using a combination of cold denaturation with NMR spectroscopy, we reveal detailed insight into the unfolding of the homodimeric repressor protein CylR2. Seven three-dimensional structures of CylR2 at temperatures from 25 °C to -16 °C reveal a progressive dissociation of the dimeric protein into a native-like monomeric intermediate followed by transition into a highly dynamic, partially folded state. The core of the partially folded state seems critical for biological function and misfolding.

Reassessing the folding of the KIX domain: evidence for a two-state mechanism

The debate about the presence and role of intermediates in the folding of proteins has been a critical issue, especially for fast folders. One of the classical methodologies to identify such metastable species is the "burst-phase analysis," whereby the observed signal amplitude from stopped-flow traces is determined as a function of denaturant concentration. However, a complication may arise when folding is sufficiently fast to jeopardize the reliability of the stopped-flow technique. In this study, we reassessed the folding of the KIX domain from cAMP Response Element-Binding (CREB)-binding protein, which has been proposed to involve the formation of an intermediate that accumulates in the dead time of the stopped flow. By using an in-house-built capillary continuous flow with a 50-μs dead time, we demonstrate that this intermediate is not present; the problem arose because of the instrumental limitation of the standard stopped flow to assess very fast refolding rate constants (e.g., ≥ 500 s⁻¹).

Dynamic regulation on energy landscape evolution of single-molecule protein by conformational fluctuation

We formalize a theory to help explore the effect of conformational fluctuation on the energy landscape evolution of single-molecule protein. Using this formalization, we investigate the photon emission from single photoactivated fluorescent protein. A bimodal regulation on the energy landscape evolution was discovered, and its origin was attributed to slow conformational fluctuations of the protein matrix.

Nonnative interactions regulate folding and switching of myristoylated protein

We present an integrated experimental and computational study of the molecular mechanisms by which myristoylation affects protein folding and function, which has been little characterized to date. Myristoylation, the covalent linkage of a hydrophobic C14 fatty acyl chain to the N-terminal glycine in a protein, is a common modification that plays a critical role in vital regulated cellular processes by undergoing reversible energetic and conformational switching. Coarse-grained folding simulations for the model pH-dependent actin- and membrane-binding protein hisactophilin reveal that nonnative hydrophobic interactions of the myristoyl with the protein as well as nonnative electrostatic interactions have a pronounced effect on folding rates and thermodynamic stability. Folding measurements for hydrophobic residue mutations of hisactophilin and atomistic simulations indicate that the nonnative interactions of the myristoyl group in the folding transition state are nonspecific and robust, and so smooth the energy landscape for folding. In contrast, myristoyl interactions in the native state are highly specific and tuned for sensitive control of switching functionality. Simulations and amide hydrogen exchange measurements provide evidence for increases as well as decreases in stability localized on one side of the myristoyl binding pocket in the protein, implicating strain and altered dynamics in switching. The effects of folding and function arising from myristoylation are profoundly different from the effects of other post-translational modifications.

Visualizing transient protein-folding intermediates by tryptophan-scanning mutagenesis

To understand how proteins fold, assemble and function, it is necessary to characterize the structure and dynamics of each state they adopt during their lifetime. Experimental characterization of the transient states of proteins remains a major challenge because high-resolution structural techniques, including NMR and X-ray crystallography, cannot be directly applied to study short-lived protein states. To circumvent this limitation, we show that transient states during protein folding can be characterized by measuring the fluorescence of tryptophan residues, introduced at many solvent-exposed positions to determine whether each position is native-like, denatured-like or non-native-like in the intermediate state. We use this approach to characterize a late-folding-intermediate state of the small globular mammalian protein ubiquitin, and we show the presence of productive non-native interactions that suggest a 'flycatcher' mechanism of concerted binding and folding.

Understanding how small helical proteins fold: conformational dynamics of Im proteins relevant to their folding landscapes

Understanding the mechanism of folding of small proteins requires characterization of their starting unfolded states and any partially unfolded states populated during folding. Here, we review what is known from NMR about these states of Im7, a 4-helix bundle protein that folds via an on-pathway intermediate, and show that there is an alignment of non-native structure in urea-unfolded Im7 with the helices of native Im7 that is a consequence of hydrophobic helix-promoting residues also promoting cluster-formation in the unfolded protein. We suggest that this kind of alignment is present in other proteins and is relevant to how native state topology determines folding rates.

Conformational properties of the unfolded state of Im7 in nondenaturing conditions

The unfolded ensemble in aqueous solution represents the starting point of protein folding. Characterisation of this species is often difficult since the native state is usually predominantly populated at equilibrium. Previous work has shown that the four-helix protein, Im7 (immunity protein 7), folds via an on-pathway intermediate. While the transition states and folding intermediate have been characterised in atomistic detail, knowledge of the unfolded ensemble under the same ambient conditions remained sparse. Here, we introduce destabilising amino acid substitutions into the sequence of Im7, such that the unfolded state becomes predominantly populated at equilibrium in the absence of denaturant. Using far- and near-UV CD, fluorescence, urea titration and heteronuclear NMR experiments, we show that three amino acid substitutions (L18A-L19A-L37A) are sufficient to prevent Im7 folding, such that the unfolded state is predominantly populated at equilibrium. Using measurement of chemical shifts, (15)N transverse relaxation rates and sedimentation coefficients, we show that the unfolded species of L18A-L19A-L37A deviates significantly from random-coil behaviour. Specifically, we demonstrate that this unfolded species is compact (R(h)=25 Å) relative to the urea-denatured state (R(h)≥30 Å) and contains local clusters of hydrophobic residues in regions that correspond to the four helices in the native state. Despite these interactions, there is no evidence for long-range stabilising tertiary interactions or persistent helical structure. The results reveal an unfolded ensemble that is conformationally restricted in regions of the polypeptide chain that ultimately form helices I, II and IV in the native state.

Strand swapping regulates the iron-sulfur cluster in the diabetes drug target mitoNEET

MitoNEET is a recently identified diabetes drug target that coordinates a transferable 2Fe-2S cluster, and additionally contains an unusual strand swap. In this manuscript, we use a dual basin structure-based model to predict and characterize the folding and functionality of strand swapping in mitoNEET. We demonstrate that a strand unswapped conformation is kinetically accessible and that multiple levels of control are employed to regulate the conformational dynamics of the system. Environmental factors such as temperature can shift route preference toward the unswapped pathway. Additionally we see that a region recently identified as contributing to frustration in folding acts as a regulatory hinge loop that modulates conformational balance. Interestingly, strand unswapping transfers strain specifically to cluster-coordinating residues, opening the cluster-coordinating pocket. Strengthening contacts within the cluster-coordinating pocket opens a new pathway between the swapped and unswapped conformation that utilizes cracking to bypass the unfolded basin. These results suggest that local control within distinct regions affect motions important in regulating mitoNEET's 2Fe-2S clusters.

Knot formation in newly translated proteins is spontaneous and accelerated by chaperonins

Topological knots are found in a considerable number of protein structures, but it is not clear how they knot and fold within the cellular environment. We investigated the behavior of knotted protein molecules as they are first synthesized by the ribosome using a cell-free translation system. We found that newly translated knotted proteins can spontaneously self-tie and do not require the assistance of molecular chaperones to fold correctly to their trefoil-knotted structures. This process is slow but efficient, and we found no evidence of misfolded species. A kinetic analysis indicates that the knotting process is rate limiting, occurs post-translationally, and is specifically and significantly (P < 0.001) accelerated by the GroEL-GroES chaperonin complex. This demonstrates a new active mechanism for this molecular chaperone and suggests that chaperonin-catalyzed knotting probably dominates in vivo. These results explain how knotted protein structures have withstood evolutionary pressures despite their topological complexity.

Nuclease colicins and their immunity proteins

It is more than 80 years since Gratia first described 'a remarkable antagonism between two strains of Escherichia coli'. Shown subsequently to be due to the action of proteins (or peptides) produced by one bacterium to kill closely related species with which it might be cohabiting, such bacteriocins have since been shown to be commonplace in the internecine warfare between bacteria. Bacteriocins have been studied primarily from the twin perspectives of how they shape microbial communities and how they penetrate bacteria to kill them. Here, we review the modes of action of a family of bacteriocins that cleave nucleic acid substrates in E. coli, known collectively as nuclease colicins, and the specific immunity (inhibitor) proteins that colicin-producing organisms make in order to avoid committing suicide. In a process akin to targeting in mitochondria, nuclease colicins engage in a variety of cellular associations in order to translocate their cytotoxic domains through the cell envelope to the cytoplasm. As well as informing on the process itself, the study of nuclease colicin import has also illuminated functional aspects of the host proteins they parasitize. We also review recent studies where nuclease colicins and their immunity proteins have been used as model systems for addressing fundamental problems in protein folding and protein-protein interactions, areas of biophysics that are intimately linked to the role of colicins in bacterial competition and to the import process itself.

NMR characterisation of the relationship between frustration and the excited state of Im7

Previous work shows that Im9 folds in a two-state transition while its homologue Im7 folds in a three-state transition via an on-pathway kinetic intermediate state (KIS), with this difference being related to frustration in the structure of Im7. We have used NMR spectroscopy to study conformational dynamics connected to the frustration. A combination of equilibrium peptide N(1)H/N(2)H exchange, model-free analyses of backbone NH relaxation data and relaxation dispersion (RD)-NMR shows that the native state of Im7 is in equilibrium with an intermediate state that is lowly populated [equilibrium intermediate state (EIS)]. Comparison of kinetic and thermodynamic parameters describing the EIS native-state equilibrium obtained by RD-NMR with previously reported parameters describing the KIS native-state equilibrium obtained from stopped-flow fluorescence studies of refolding His-tagged Im7 shows that the KIS and the EIS are the same species. (15)N chemical shifts of the EIS obtained from the RD-NMR analysis show that residues forming helix III in the native state are unstructured in the EIS while other residues experiencing frustration in the native state are in structured regions of the EIS. We show that binding of Im7 and its L53A/I54A variant (which resembles the EIS as shown in previous work) to the cognate partner for Im7, the DNase domain of colicin E7, causes the dynamic processes associated with the frustration to be dampened.

A diversity of assembly mechanisms of a generic amyloid fold

Protein misfolding and amyloid assembly have long been recognized as being responsible for many devastating human diseases. Recent findings indicate that amyloid assemblies may facilitate crucial biological processes from bacteria to mammals. This review focuses on the mechanistic understanding of amyloid formation, including the transformation of initially innocuous proteins into oligomers and fibrils. The result is a competing folding and assembly energy landscape, which contains a number of routes by which the polypeptide chain can convert its primary sequence into functional structures, dysfunctional assemblies, or epigenetic entities that provide both threats and opportunities in the evolution of life.

The effect of glycosylation on the folding kinetics of erythropoietin

Glycosylation is a common posttranslational modification that generally increases protein solubility and thermodynamic stability. Less is known about how this modification influences protein folding, particularly folding processes involving intermediate species. In the present report, folding comparisons of a nonglycosylated erythropoietin (EPO) mutant are made with the fully glycosylated EPO, which was recently shown to fold by a three-state on-pathway mechanism. The absence of glycosylation did not alter the folding mechanism of EPO but did greatly decrease the stability of the intermediate species, change the rate-limiting step of the folding reaction, and accelerate the folding kinetics to both the intermediate state and the native state. Surprisingly, glycosylation stabilized the intermediate species to a greater extent than it increased the EPO equilibrium stability. These results suggest that glycosylation impedes the latter EPO folding steps rather than accelerating them by biasing particular folding pathways, as previously proposed for folding reactions initiated from unfolded ensembles with minimal residual structure. Due to the specific biological processes modulated by EPO glycosylation, however, there may be little evolutionary pressure to fold on a faster, more direct pathway at the expense of biological function, particularly given the protective role glycosylation has at preventing EPO aggregation. Lastly, evidence that is consistent with glycosylation destabilizing the unfolded state to some degree and contributing to the greater equilibrium stability of the glycosylated EPO is presented.

Cutting off functional loops from homodimeric enzyme superoxide dismutase 1 (SOD1) leaves monomeric β-barrels

Demetallation of the homodimeric enzyme Cu/Zn-superoxide dismutase (SOD1) is known to unleash pronounced dynamic motions in the long active-site loops that comprise almost a third of the folded structure. The resulting apo species, which shows increased propensity to aggregate, stands out as the prime disease precursor in amyotrophic lateral sclerosis (ALS). Even so, the detailed structural properties of the apoSOD1 framework have remained elusive and controversial. In this study, we examine the structural interplay between the central apoSOD1 barrel and the active-site loops by simply cutting them off; loops IV and VII were substituted with short Gly-Ala-Gly linkers. The results show that loop removal breaks the dimer interface and leads to soluble, monomeric β-barrels with high structural integrity. NMR-detected nuclear Overhauser effects are found between all of the constituent β-strands, confirming ordered interactions across the whole barrel. Moreover, the breathing motions of the SOD1 barrel are overall insensitive to loop removal and yield hydrogen/deuterium protection factors typical for cooperatively folded proteins (i.e. the active-site loops act as a "bolt-on" domain with little dynamic influence on its structural foundation). The sole exceptions are the relatively low protection factors in β-strand 5 and the turn around Gly-93, a hot spot for ALS-provoking mutations, which decrease even further upon loop removal. Taken together, these data suggest that the cytotoxic function of apoSOD1 does not emerge from its folded ground state but from a high energy intermediate or even from the denatured ensemble.

Malleability of folding intermediates in the homeodomain superfamily

Members of the homeodomain superfamily are three-helix bundle proteins whose second and third helices form a helix-turn-helix motif (HTH). Their folding mechanism slides from the ultrafast, three-state framework mechanism for the engrailed homeodomain (EnHD), in which the HTH motif is independently stable, to an apparent two-state nucleation-condensation model for family members with an unstable HTH motif. The folding intermediate of EnHD has nearly native HTH structure, but it is not docked with helix1. The determinant of whether two- or three-state folding was hypothesized to be the stability of the HTH substructure. Here, we describe a detailed Φ-value analysis of the folding of the Pit1 homeodomain, which has similar ultrafast kinetics to that of EnHD. Formation of helix1 was strongly coupled with formation of HTH, which was initially surprising because they are uncoupled in the EnHD folding intermediate. However, we found a key difference between Pit1 and EnHD: The isolated peptide corresponding to the HTH motif in Pit1 was not folded in the absence of H1. Independent molecular dynamics simulations of Pit1 unfolding found an intermediate with H1 misfolded onto the HTH motif. The Pit1 folding pathway is the connection between that of EnHD and the slower folding homeodomains and provides a link in the transition of mechanisms from two- to three-state folding in this superfamily. The malleability of folding intermediates can lead to unstable substructures being stabilized by a variety of nonnative interactions, adding to the continuum of folding mechanisms.

Dissecting key residues in folding and stability of the bacterial immunity protein 7

The small four-helix immunity protein, Im7, has previously been shown to fold via a compact intermediate containing three of the four native helices. The short, six-residue helix III only docks onto the developing Im7 structure after the rate-limiting second transition state has been traversed. Previous work demonstrated that mutation of the helix III sequence can be used to trap the protein in the on-pathway intermediate ensemble at equilibrium. Here the role played by individual residues in the native helix III sequence in locking Im7 into a stable native structure is further examined. This work commenced with an Im7 sequence trapped in the partially folded state by substitution of the six residues in helix III with a polyglycine sequence. Biophysical analysis of variants in which individual residues from the native helix III sequence, and combinations of these residues, were introduced into this background demonstrated a critical requirement for three residues, Leu 53, Ile 54 and Tyr 55, to lock Im7 into its unique native structure. The results demonstrate a stringent constraint on the evolution of the Im7 helix III sequence rationalizing its high-sequence identity in the fold family. Thus, Leu 53 and Ile 54 provide crucial stabilizing interactions in the hydrophobic core of native Im7, while Tyr 55 is required for both stability and function. In contrast, Tyr 56 is critical for colicin binding and has no role in maintaining a stable native fold.

Perturbing the folding energy landscape of the bacterial immunity protein Im7 by site-specific N-linked glycosylation

N-linked glycosylation modulates protein folding and stability through a variety of mechanisms. As such there is considerable interest in the development of general rules to predict the structural consequences of site-specific glycosylation and to understand how these effects can be exploited in the design and development of modified proteins with advantageous properties. In this study, expressed protein ligation is used to create site-specifically glycosylated variants of the bacterial immunity protein Im7 modified with the chitobiose disaccharide (GlcNAc-GlcNAc). Glycans were introduced at seven solvent exposed sites within the Im7 sequence and the kinetic and thermodynamic consequences of N-linked glycosylation analyzed. The ΔΔG° values for glycan incorporation were found to range from +5.2 to -3.8 kJ·mol(-1). In several cases, glycosylation influences folding by modulating the local conformational preferences of the glycosylated sequence. These locally mediated effects are most prominent in the center of α-helices where glycosylation negatively effects folding and in compact turn motifs between segments of ordered secondary structure where glycosylation promotes folding and enhances the overall stability of the native protein. The studies also provide insight into why glycosylation is commonly identified at the transition between different types of secondary structure and when glycosylation may be used to elaborate protein structure to protect disordered sequences from proteolysis or immune system recognition.

Dynameomics: a comprehensive database of protein dynamics

The dynamic behavior of proteins is important for an understanding of their function and folding. We have performed molecular dynamics simulations of the native state and unfolding pathways of over 2000 protein/peptide systems (approximately 11,000 independent simulations) representing the majority of folds in globular proteins. These data are stored and organized using an innovative database approach, which can be mined to obtain both general and specific information about the dynamics and folding/unfolding of proteins, relevant subsets thereof, and individual proteins. Here we describe the project in general terms and the type of information contained in the database. Then we provide examples of mining the database for information relevant to protein folding, structure building, the effect of single-nucleotide polymorphisms, and drug design. The native state simulation data and corresponding analyses for the 100 most populated metafolds, together with related resources, are publicly accessible through http://www.dynameomics.org.

Slow formation of aggregation-resistant beta-sheet folding intermediates

Protein folding has been studied extensively for decades, yet our ability to predict how proteins reach their native state from a mechanistic perspective is still rudimentary at best, limiting our understanding of folding-related processes in vivo and our ability to manipulate proteins in vitro. Here, we investigate the in vitro refolding mechanism of a large beta-helix protein, pertactin, which has an extended, elongated shape. At 55 kDa, this single domain, all-beta-sheet protein allows detailed analysis of the formation of beta-sheet structure in larger proteins. Using a combination of fluorescence and far-UV circular dichroism spectroscopy, we show that the pertactin beta-helix refolds remarkably slowly, with multiexponential kinetics. Surprisingly, despite the slow refolding rates, large size, and beta-sheet-rich topology, pertactin refolding is reversible and not complicated by off-pathway aggregation. The slow pertactin refolding rate is not limited by proline isomerization, and 30% of secondary structure formation occurs within the rate-limiting step. Furthermore, site-specific labeling experiments indicate that the beta-helix refolds in a multistep but concerted process involving the entire protein, rather than via initial formation of the stable core substructure observed in equilibrium titrations. Hence pertactin provides a valuable system for studying the refolding properties of larger, beta-sheet-rich proteins, and raises intriguing questions regarding the prevention of aggregation during the prolonged population of partially folded, beta-sheet-rich refolding intermediates. Proteins 2010. (c) 2009 Wiley-Liss, Inc.

Single-molecule studies of the Im7 folding landscape

Under appropriate conditions, the four-helical Im7 (immunity protein 7) folds from an ensemble of unfolded conformers to a highly compact native state via an on-pathway intermediate. Here, we investigate the unfolded, intermediate, and native states populated during folding using diffusion single-pair fluorescence resonance energy transfer by measuring the efficiency of energy transfer (or proximity or P ratio) between pairs of fluorophores introduced into the side chains of cysteine residues placed in the center of helices 1 and 4, 1 and 3, or 2 and 4. We show that while the native states of each variant give rise to a single narrow distribution with high P values, the distributions of the intermediates trapped at equilibrium (denoted I^eqm) are fitted by two Gaussian distributions. Modulation of the folding conditions from those that stabilize the intermediate to those that destabilize the intermediate enabled the distribution of lower P value to be assigned to the population of the unfolded ensemble in equilibrium with the intermediate state. The reduced stability of the I^eqm variants allowed analysis of the effect of denaturant concentration on the compaction and breadth of the unfolded state ensemble to be quantified from 0 to 6 M urea. Significant compaction is observed as the concentration of urea is decreased in both the presence and absence of sodium sulfate, as previously reported for a variety of proteins. In the presence of Na₂SO₄ in 0 M urea, the P value of the unfolded state ensemble approaches that of the native state. Concurrent with compaction, the ensemble displays increased peak width of P values, possibly reflecting a reduction in the rate of conformational exchange among iso-energetic unfolded, but compact conformations. The results provide new insights into the initial stages of folding of Im7 and suggest that the unfolded state is highly conformationally constrained at the outset of folding.

Fundamental link between folding states and functional states of proteins

Folding and function of proteins are two aspects of proteins which are usually considered as basically unrelated phenomena that are optimized by evolution independently. From the funnel model of folding/unfolding and the associated energy landscape, we infer the paradigm that the minimum number of folding intermediates is determined by the number of all functional states of a protein ("essential" folding intermediates). Here, we demonstrate the supposed fundamental link using the Ras protein complexed with the GTP analogue GppNHp that occurs in two structural states coexisting in solution. State 2 was shown earlier to represent the effector interacting state, and the function of state 1 was hitherto unknown. By (31)P NMR spectroscopy, we demonstrate that state 1 represents the conformation interacting with guanine nucleotide exchange factors (GEFs). Denaturation experiments of the protein with a chaotropic reagent show that both functional states coexist during folding and unfolding. Application of high pressure represents another perturbation of the energy landscape, leading to an increased population of the state 1 as observed by NMR spectroscopy. The specific volume difference between the two states DeltaV(12) is 17.2 +/- 0.5 mL mol(-1), indicating that state 1 represents a more open conformation of the protein. The free energies of stabilization for state 1 and state 2 at 278 K can be determined as 8.3 and 9.8 kJ mol(-1), respectively.