Folding without charges

Beyond Misfolding: A New Paradigm for the Relationship Between Protein Folding and Aggregation

Aggregation is intricately linked to protein folding, necessitating a precise understanding of their relationship. Traditionally, aggregation has been viewed primarily as a sequential consequence of protein folding and misfolding. However, this conventional paradigm is inherently incomplete and can be deeply misleading. Remarkably, it fails to adequately explain how intrinsic and extrinsic factors, such as charges and cellular macromolecules, prevent intermolecular aggregation independently of intramolecular protein folding and structure. The pervasive inconsistencies between protein folding and aggregation call for a new framework. In all combined reactions of molecules, both intramolecular and intermolecular rate (or equilibrium) constants are mutually independent; accordingly, intrinsic and extrinsic factors independently affect both rate constants. This universal principle, when applied to protein folding and aggregation, indicates that they should be treated as two independent yet interconnected processes. Based on this principle, a new framework provides groundbreaking insights into misfolding, Anfinsen's thermodynamic hypothesis, molecular chaperones, intrinsic chaperone-like activities of cellular macromolecules, intermolecular repulsive force-driven aggregation inhibition, proteome solubility maintenance, and proteinopathies. Consequently, this paradigm shift not only refines our current understanding but also offers a more comprehensive view of how aggregation is coupled to protein folding in the complex cellular milieu.

Dissecting the stability determinants of a challenging de novo protein fold using massively parallel design and experimentation

Designing entirely new protein structures remains challenging because we do not fully understand the biophysical determinants of folding stability. Yet, some protein folds are easier to design than others. Previous work identified the 43-residue ɑββɑ fold as especially challenging: The best designs had only a 2% success rate, compared to 39 to 87% success for other simple folds [G. J. Rocklin et al., Science 357, 168-175 (2017)]. This suggested the ɑββɑ fold would be a useful model system for gaining a deeper understanding of folding stability determinants and for testing new protein design methods. Here, we designed over 10,000 new ɑββɑ proteins and found over 3,000 of them to fold into stable structures using a high-throughput protease-based assay. NMR, hydrogen-deuterium exchange, circular dichroism, deep mutational scanning, and scrambled sequence control experiments indicated that our stable designs fold into their designed ɑββɑ structures with exceptional stability for their small size. Our large dataset enabled us to quantify the influence of universal stability determinants including nonpolar burial, helix capping, and buried unsatisfied polar atoms, as well as stability determinants unique to the ɑββɑ topology. Our work demonstrates how large-scale design and test cycles can solve challenging design problems while illuminating the biophysical determinants of folding.

Physicochemical classification of organisms

The hypervariable residues that compose the major part of proteins’ surfaces are generally considered outside evolutionary control. Yet, these “nonconserved” residues determine the outcome of stochastic encounters in crowded cells. It has recently become apparent that these encounters are not as random as one might imagine, but carefully orchestrated by the intracellular electrostatics to optimize protein diffusion, interactivity, and partner search. The most influential factor here is the protein surface-charge density, which takes different optimal values across organisms with different intracellular conditions. In this study, we examine how far the net-charge density and other physicochemical properties of proteomes will take us in terms of distinguishing organisms in general. The results show that these global proteome properties not only follow the established taxonomical hierarchy, but also provide clues to functional adaptation. In many cases, the proteome–property divergence is even resolved at species level. Accordingly, the variable parts of the genes are not as free to drift as they seem in sequence alignment, but present a complementary tool for functional, taxonomic, and evolutionary assignment.

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.

Protein structure and aggregation: a marriage of necessity ruled by aggregation gatekeepers

Protein aggregation propensity is a pervasive and seemingly inescapable property of proteomes. Strikingly, a significant fraction of the proteome is supersaturated, meaning that, for these proteins, their native conformation is less stable than the aggregated state. Maintaining the integrity of a proteome under such conditions is precarious and requires energy-consuming proteostatic regulation. Why then is aggregation propensity maintained at such high levels over long evolutionary timescales? Here, we argue that the conformational stability of the native and aggregated states are correlated thermodynamically and that codon usage strengthens this correlation. As a result, the folding of stable proteins requires kinetic control to avoid aggregation, provided by aggregation gatekeepers. These unique residues are evolutionarily selected to kinetically favor native folding, either on their own or by coopting chaperones.

A Conceptual Framework for Integrating Cellular Protein Folding, Misfolding and Aggregation

How proteins properly fold and maintain solubility at the risk of misfolding and aggregation in the cellular environments still remains largely unknown. Aggregation has been traditionally treated as a consequence of protein folding (or misfolding). Notably, however, aggregation can be generally inhibited by affecting the intermolecular interactions leading to aggregation, independently of protein folding and conformation. We here point out that rigorous distinction between protein folding and aggregation as two independent processes is necessary to reconcile and underlie all observations regarding the combined cellular protein folding and aggregation. So far, the direct attractive interactions (e.g., hydrophobic interactions) between cellular macromolecules including chaperones and interacting polypeptides have been widely believed to mainly stabilize polypeptides against aggregation. However, the intermolecular repulsions by large excluded volume and surface charges of cellular macromolecules can play a key role in stabilizing their physically connected polypeptides against aggregation, irrespective of the connection types and induced conformational changes, underlying the generic intrinsic chaperone activity of cellular macromolecules. Such rigorous distinction and intermolecular repulsive force-driven aggregation inhibition by cellular macromolecules could give new insights into understanding the complex cellular protein landscapes that remain uncharted.

Surface charged amino acid-based strategy for rational engineering of kinetic stability and specific activity of enzymes: Linking experiments with computational modeling

A rational workflow for engineering kinetically stable enzymes with good specific activity by surface charged amino acids engineering was proposed based on systematically analyzing the results of mutating 44 negatively charged surface amino acids of a thermophilic β-mannanase (ManAK). Computational data, combined with experimental results indicated that percentage side-chain solvent accessibility (PSSA), changes in Gibbs free energy of unfolding (∆∆G_mut) and root-mean-square fluctuations (RMSF) could be suitable for screening kinetically stable mutants. A combinational standard (∆∆G_mut < -0.5 kJ/mol and RMSF >0.68 Å) resulted a decrease in the proportion of destabilizing mutants to 12.5%. The perturbations of substrate affinity and specific activity caused by mutation were weakened as the shortest distance from C_α of mutated site to C_α of catalytic sites (Ds_Cα-Cα) increased. Results indicated that hotspot zones contributing to the local stability and integrity of catalytic motif at elevated temperatures might be widely distributed across spatial structure of the protein, while the mutation perturbation on enzyme specific activity demonstrated a gradually weakening trend from the catalytic core to the protein surface. These findings further our understanding of the structural-functional relationships of protein and highlight a deduced workflow to engineering industrially useful enzymes.

Polyanions Cause Protein Destabilization Similar to That in Live Cells

The structural stability of proteins is found to markedly change upon their transfer to the crowded interior of live cells. For some proteins, the stability increases, while for others, it decreases, depending on both the sequence composition and the type of host cell. The mechanism seems to be linked to the strength and conformational bias of the diffusive in-cell interactions, where protein charge is found to play a decisive role. Because most proteins, nucleotides, and membranes carry a net-negative charge, the intracellular environment behaves like a polyanionic (Z:1) system with electrostatic interactions different from those of standard 1:1 ion solutes. To determine how such polyanion conditions influence protein stability, we use negatively charged polyacetate ions to mimic the net-negatively charged cellular environment. The results show that, per Na+ equivalent, polyacetate destabilizes the model protein SOD1barrel significantly more than monoacetate or NaCl. At an equivalent of 100 mM Na+, the polyacetate destabilization of SOD1barrel is similar to that observed in live cells. By the combined use of equilibrium thermal denaturation, folding kinetics, and high-resolution nuclear magnetic resonance, this destabilization is primarily assigned to preferential interaction between polyacetate and the globally unfolded protein. This interaction is relatively weak and involves mainly the outermost N-terminal region of unfolded SOD1barrel. Our findings point thus to a generic influence of polyanions on protein stability, which adds to the sequence-specific contributions and needs to be considered in the evaluation of in vivo data.

Engineering crowding sensitivity into protein linkers

The intracellular environment contains a high concentration of biomacromolecules that present steric barriers and ample surface area for weak chemical interactions. Consequently, these forces influence protein conformations and protein self-assembly, with an outcome that depends on the sum of the effects resulting from crowding. Linkers are disordered domains that lack tertiary structure, and this flexible nature would render them susceptible to compression or extension under crowded conditions, compared to the equilibrium conformation in a dilute buffer. The change in distance between the linked proteins can become essential where it attenuates protein activity. In this chapter, we first discuss the experimental findings in vitro and in the cell on how linkers and other relevant macromolecules are affected by crowding. We focus on the dependence on the linker's size, flexibility, and the intra- and intermolecular interactions. Although the experimental data on the systematic variation of proteins in a buffer and cells is limited, extrapolating the available insights allows us to propose a protocol on how to engineer the directionality of crowding effects in the linker. Finally, we describe a straightforward experimental protocol on the determination of crowding sensitivity in a buffer and cell.

Autonomous aggregation suppression by acidic residues explains why chaperones favour basic residues

Many chaperones favour binding to hydrophobic sequences that are flanked by basic residues while disfavouring acidic residues. However, the origin of this bias in protein quality control remains poorly understood. Here, we show that while acidic residues are the most efficient aggregation inhibitors, they are also less compatible with globular protein structure than basic amino acids. As a result, while acidic residues allow for chaperone-independent control of aggregation, their use is structurally limited. Conversely, we find that, while being more compatible with globular structure, basic residues are not sufficient to autonomously suppress protein aggregation. Using Hsp70, we show that chaperones with a bias towards basic residues are structurally adapted to prioritize aggregating sequences whose structural context forced the use of the less effective basic residues. The hypothesis that emerges from our analysis is that the bias of many chaperones for basic residues results from fundamental thermodynamic and kinetic constraints of globular structure. This also suggests the co-evolution of basic residues and chaperones allowed for an expansion of structural variety in the protein universe.

Behind closed gates - chaperones and charged residues determine protein fate

Charged residues flanking aggregation-prone regions play a role in protein folding and prevention of aggregation. In this issue of The EMBO Journal, Houben et al exploit the role of such charged gatekeepers in aggregation suppression and find that negative charges are more effective than positive ones. Strikingly, the prominent Hsp70 chaperone has a strong preference for the less effective, basic gate keepers. This implies co-adaptation of chaperone specificity and composition of protein sequences in evolution.

Proteins are Solitary! Pathways of Protein Folding and Aggregation in Protein Mixtures

We present a computational and experimental study on the folding and aggregation in solutions of multiple protein mixtures at different concentrations. We show how in protein mixtures each component is capable of maintaining its folded state at densities greater than the one at which they would precipitate in single-species solutions. We demonstrate the generality of our observation over many different proteins using computer simulations capable of fully characterizing the cross-aggregation phase diagram of all the mixtures. Dynamic light scattering experiments were performed to evaluate the aggregation of two proteins, bovine serum albumin (BSA) and consensus tetratricopeptide repeat (CTPR), in solutions of one or both proteins. The experiments confirm our hypothesis and the simulations. These findings elucidate critical aspects of the cross-regulation of expression and aggregation of proteins exerted by the cell and on the evolutionary selection of folding and non-aggregating protein sequences, paving the way for new experimental tests.

Why do eukaryotic proteins contain more intrinsically disordered regions?

Intrinsic disorder is more abundant in eukaryotic than prokaryotic proteins. Methods predicting intrinsic disorder are based on the amino acid sequence of a protein. Therefore, there must exist an underlying difference in the sequences between eukaryotic and prokaryotic proteins causing the (predicted) difference in intrinsic disorder. By comparing proteins, from complete eukaryotic and prokaryotic proteomes, we show that the difference in intrinsic disorder emerges from the linker regions connecting Pfam domains. Eukaryotic proteins have more extended linker regions, and in addition, the eukaryotic linkers are significantly more disordered, 38% vs. 12-16% disordered residues. Next, we examined the underlying reason for the increase in disorder in eukaryotic linkers, and we found that the changes in abundance of only three amino acids cause the increase. Eukaryotic proteins contain 8.6% serine; while prokaryotic proteins have 6.5%, eukaryotic proteins also contain 5.4% proline and 5.3% isoleucine compared with 4.0% proline and ≈ 7.5% isoleucine in the prokaryotes. All these three differences contribute to the increased disorder in eukaryotic proteins. It is tempting to speculate that the increase in serine frequencies in eukaryotes is related to regulation by kinases, but direct evidence for this is lacking. The differences are observed in all phyla, protein families, structural regions and type of protein but are most pronounced in disordered and linker regions. The observation that differences in the abundance of three amino acids cause the difference in disorder between eukaryotic and prokaryotic proteins raises the question: Are amino acid frequencies different in eukaryotic linkers because the linkers are more disordered or do the differences cause the increased disorder?

Intrinsic disorder is essential for various functions in eukaryotic cells and is a signature of eukaryotic proteins. Here, we try to understand the origin of the difference in disorder between eukaryotic and prokaryotic proteins. We show that eukaryotic proteins contain more extended linker regions and that these linker regions are significantly more disordered. Further, we show, for the first time, that the difference in disorder originates from a systematic difference in amino acid frequencies between eukaryotic and prokaryotic proteins. Three amino acids contribute to the difference in disorder; serine and proline are more abundant in eukaryotic linkers, while isoleucine is less frequent. These shifts in frequencies are observed in all phyla, protein families, structural regions and type of protein but are most pronounced in disordered and linker regions. It is tempting to speculate that the increase in serine frequencies in eukaryotes is related to regulation by kinases, but direct evidence for this is lacking. Anyhow the widespread of the shifts in abundance indicates that the differences are ancient and caused be some yet not fully understood selective difference acting on eukaryotic and prokaryotic proteins.

Conversion of a soluble protein into a potent chaperone in vivo

Molecular chaperones play an important role in cellular protein-folding assistance and aggregation inhibition. As a different but complementary model, we previously proposed that, in general, soluble cellular macromolecules with large excluded volume and surface charges exhibit intrinsic chaperone activity to prevent aggregation of their connected polypeptides irrespective of the connection type, thereby contributing to efficient protein folding. As a proof of concept, we here demonstrated that a model recombinant protein with a specific sequence-binding domain robustly exerted chaperone activity toward various proteins harbouring a short recognition tag of 7 residues in Escherichia coli. The chaperone activity of this protein was comparable to that of representative E. coli chaperones in vivo. Furthermore, in vitro refolding experiments confirmed the in vivo results. Our findings reveal that a soluble protein exhibits the intrinsic chaperone activity to prevent off-pathway aggregation of its interacting proteins, leading to more productive folding while allowing them to fold according to their intrinsic folding pathways. This study gives new insights into the plausible chaperoning role of soluble cellular macromolecules in terms of aggregation inhibition and indirect folding assistance.

Effects of protein size, thermodynamic stability, and net charge on cotranslational folding on the ribosome

During the last five decades, studies of protein folding in dilute buffer solutions have produced a rich picture of this complex process. In the cell, however, proteins can start to fold while still attached to the ribosome (cotranslational folding) and it is not yet clear how the ribosome affects the folding of protein domains of different sizes, thermodynamic stabilities, and net charges. Here, by using arrest peptides as force sensors and on-ribosome pulse proteolysis, we provide a comprehensive picture of how the distance from the peptidyl transferase center in the ribosome at which proteins fold correlates with protein size. Moreover, an analysis of a large collection of mutants of the Escherichia coli ribosomal protein S6 shows that the force exerted on the nascent chain by protein folding varies linearly with the thermodynamic stability of the folded state, and that the ribosome environment disfavors folding of domains of high net-negative charge.

Size and topology modulate the effects of frustration in protein folding

The presence of conflicting interactions, or frustration, determines how fast biomolecules can explore their configurational landscapes. Recent experiments have provided cases of systems with slow reconfiguration dynamics, perhaps arising from frustration. While it is well known that protein folding speed and mechanism are strongly affected by the protein native structure, it is still unknown how the response to frustration is modulated by the protein topology. We explore the effects of nonnative interactions in the reconfigurational and folding dynamics of proteins with different sizes and topologies. We find that structural correlations related to the folded state size and topology play an important role in determining the folding kinetics of proteins that otherwise have the same amount of nonnative interactions. In particular, we find that the reconfiguration dynamics of α-helical proteins are more susceptible to frustration than β-sheet proteins of the same size. Our results may explain recent experimental findings and suggest that attempts to measure the degree of frustration due to nonnative interactions might be more successful with α-helical proteins.

Electrostatic Interactions in Protein Structure, Folding, Binding, and Condensation

Charged and polar groups, through forming ion pairs, hydrogen bonds, and other less specific electrostatic interactions, impart important properties to proteins. Modulation of the charges on the amino acids, e.g., by pH and by phosphorylation and dephosphorylation, have significant effects such as protein denaturation and switch-like response of signal transduction networks. This review aims to present a unifying theme among the various effects of protein charges and polar groups. Simple models will be used to illustrate basic ideas about electrostatic interactions in proteins, and these ideas in turn will be used to elucidate the roles of electrostatic interactions in protein structure, folding, binding, condensation, and related biological functions. In particular, we will examine how charged side chains are spatially distributed in various types of proteins and how electrostatic interactions affect thermodynamic and kinetic properties of proteins. Our hope is to capture both important historical developments and recent experimental and theoretical advances in quantifying electrostatic contributions of proteins.

Unified understanding of folding and binding mechanisms of globular and intrinsically disordered proteins

Extensive experimental and theoretical studies have advanced our understanding of the mechanisms of folding and binding of globular proteins, and coupled folding and binding of intrinsically disordered proteins (IDPs). The forces responsible for conformational changes and binding are common in both proteins; however, these mechanisms have been separately discussed. Here, we attempt to integrate the mechanisms of coupled folding and binding of IDPs, folding of small and multi-subdomain proteins, folding of multimeric proteins, and ligand binding of globular proteins in terms of conformational selection and induced-fit mechanisms as well as the nucleation–condensation mechanism that is intermediate between them. Accumulating evidence has shown that both the rate of conformational change and apparent rate of binding between interacting elements can determine reaction mechanisms. Coupled folding and binding of IDPs occurs mainly by induced-fit because of the slow folding in the free form, while ligand binding of globular proteins occurs mainly by conformational selection because of rapid conformational change. Protein folding can be regarded as the binding of intramolecular segments accompanied by secondary structure formation. Multi-subdomain proteins fold mainly by the induced-fit (hydrophobic collapse) mechanism, as the connection of interacting segments enhances the binding (compaction) rate. Fewer hydrophobic residues in small proteins reduce the intramolecular binding rate, resulting in the nucleation–condensation mechanism. Thus, the folding and binding of globular proteins and IDPs obey the same general principle, suggesting that the coarse-grained, statistical mechanical model of protein folding is promising for a unified theoretical description of all mechanisms.

Modulation of electrostatic interactions to reveal a reaction network unifying the aggregation behaviour of the Aβ42 peptide and its variants

The aggregation of the amyloid β peptide (Aβ42), which is linked to Alzheimer's disease, can be altered significantly by modulations of the peptide's intermolecular electrostatic interactions. Variations in sequence and solution conditions have been found to lead to highly variable aggregation behaviour. Here we modulate systematically the electrostatic interactions governing the aggregation kinetics by varying the ionic strength of the solution. We find that changes in the solution ionic strength induce a switch in the reaction pathway, altering the dominant mechanisms of aggregate multiplication. This strategy thereby allows us to continuously sample a large space of different reaction mechanisms and develop a minimal reaction network that unifies the experimental kinetics under a wide range of different conditions. More generally, this universal reaction network connects previously separate systems, such as charge mutants of the Aβ42 peptide, on a continuous mechanistic landscape, providing a unified picture of the aggregation mechanism of Aβ42.

Physicochemical code for quinary protein interactions in Escherichia coli

How proteins sense and navigate the cellular interior to find their functional partners remains poorly understood. An intriguing aspect of this search is that it relies on diffusive encounters with the crowded cellular background, made up of protein surfaces that are largely nonconserved. The question is then if/how this protein search is amenable to selection and biological control. To shed light on this issue, we examined the motions of three evolutionary divergent proteins in the Escherichia coli cytoplasm by in-cell NMR. The results show that the diffusive in-cell motions, after all, follow simplistic physical-chemical rules: The proteins reveal a common dependence on (i) net charge density, (ii) surface hydrophobicity, and (iii) the electric dipole moment. The bacterial protein is here biased to move relatively freely in the bacterial interior, whereas the human counterparts more easily stick. Even so, the in-cell motions respond predictably to surface mutation, allowing us to tune and intermix the protein's behavior at will. The findings show how evolution can swiftly optimize the diffuse background of protein encounter complexes by just single-point mutations, and provide a rational framework for adjusting the cytoplasmic motions of individual proteins, e.g., for rescuing poor in-cell NMR signals and for optimizing protein therapeutics.

The Role of Evolutionary Selection in the Dynamics of Protein Structure Evolution

Homology modeling is a powerful tool for predicting a protein's structure. This approach is successful because proteins whose sequences are only 30% identical still adopt the same structure, while structure similarity rapidly deteriorates beyond the 30% threshold. By studying the divergence of protein structure as sequence evolves in real proteins and in evolutionary simulations, we show that this nonlinear sequence-structure relationship emerges as a result of selection for protein folding stability in divergent evolution. Fitness constraints prevent the emergence of unstable protein evolutionary intermediates, thereby enforcing evolutionary paths that preserve protein structure despite broad sequence divergence. However, on longer timescales, evolution is punctuated by rare events where the fitness barriers obstructing structure evolution are overcome and discovery of new structures occurs. We outline biophysical and evolutionary rationale for broad variation in protein family sizes, prevalence of compact structures among ancient proteins, and more rapid structure evolution of proteins with lower packing density.

Linking computation and experiments to study the role of charge-charge interactions in protein folding and stability

Over the past two decades there has been an increase in appreciation for the role of surface charge-charge interactions in protein folding and stability. The perception shifted from the belief that charge-charge interactions are not important for protein folding and stability to the near quantitative understanding of how these interactions shape the folding energy landscape. This led to the ability of computational approaches to rationally redesign surface charge-charge interactions to modulate thermodynamic properties of proteins. Here we summarize our progress in understanding the role of charge-charge interactions for protein stability using examples drawn from my own laboratory and touch upon unanswered questions.

Atomic Resolution Structure of Monomorphic Aβ42 Amyloid Fibrils

Amyloid-β (Aβ) is a 39-42 residue protein produced by the cleavage of the amyloid precursor protein (APP), which subsequently aggregates to form cross-β amyloid fibrils that are a hallmark of Alzheimer's disease (AD). The most prominent forms of Aβ are Aβ1-40 and Aβ1-42, which differ by two amino acids (I and A) at the C-terminus. However, Aβ42 is more neurotoxic and essential to the etiology of AD. Here, we present an atomic resolution structure of a monomorphic form of AβM01-42 amyloid fibrils derived from over 500 (13)C-(13)C, (13)C-(15)N distance and backbone angle structural constraints obtained from high field magic angle spinning NMR spectra. The structure (PDB ID: 5KK3 ) shows that the fibril core consists of a dimer of Aβ42 molecules, each containing four β-strands in a S-shaped amyloid fold, and arranged in a manner that generates two hydrophobic cores that are capped at the end of the chain by a salt bridge. The outer surface of the monomers presents hydrophilic side chains to the solvent. The interface between the monomers of the dimer shows clear contacts between M35 of one molecule and L17 and Q15 of the second. Intermolecular (13)C-(15)N constraints demonstrate that the amyloid fibrils are parallel in register. The RMSD of the backbone structure (Q15-A42) is 0.71 ± 0.12 Å and of all heavy atoms is 1.07 ± 0.08 Å. The structure provides a point of departure for the design of drugs that bind to the fibril surface and therefore interfere with secondary nucleation and for other therapeutic approaches to mitigate Aβ42 aggregation.

A Soluble, Folded Protein without Charged Amino Acid Residues

Charges are considered an integral part of protein structure and function, enhancing solubility and providing specificity in molecular interactions. We wished to investigate whether charged amino acids are indeed required for protein biogenesis and whether a protein completely free of titratable side chains can maintain solubility, stability, and function. As a model, we used a cellulose-binding domain from Cellulomonas fimi, which, among proteins of more than 100 amino acids, presently is the least charged in the Protein Data Bank, with a total of only four titratable residues. We find that the protein shows a surprising resilience toward extremes of pH, demonstrating stability and function (cellulose binding) in the pH range from 2 to 11. To ask whether the four charged residues present were required for these properties of this protein, we altered them to nontitratable ones. Remarkably, this chargeless protein is produced reasonably well in Escherichia coli, retains its stable three-dimensional structure, and is still capable of strong cellulose binding. To further deprive this protein of charges, we removed the N-terminal charge by acetylation and studied the protein at pH 2, where the C-terminus is effectively protonated. Under these conditions, the protein retains its function and proved to be both soluble and have a reversible folding-unfolding transition. To the best of our knowledge, this is the first time a soluble, functional protein with no titratable side chains has been produced.

Membrane-Induced Folding of the Plant Stress Dehydrin Lti30

Dehydrins are disordered proteins that are expressed in plants as a response to embryogenesis and water-related stress. The molecular function and structural action of the dehydrins are yet elusive, but increasing evidence points to a role in protecting the structure and functional dynamics of cell membranes. An intriguing example is the cold-induced dehydrin Lti30 that binds to membranes by its conserved K segments. Moreover, this binding can be regulated by pH and phosphorylation and shifts the membrane phase transition to lower temperatures, consistent with the protein's postulated function in cold stress. In this study, we reveal how the Lti30-membrane interplay works structurally at atomic level resolution in Arabidopsis (Arabidopsis thaliana). Nuclear magnetic resonance analysis suggests that negatively charged lipid head groups electrostatically capture the protein's disordered K segments, which locally fold up into α-helical segments on the membrane surface. Thus, Lti30 conforms to the general theme of structure-function relationships by folding upon binding, in spite of its disordered, atypically hydrophilic and repetitive sequence signatures. Moreover, the fixed and well-defined structure of the membrane-bound K segments suggests that dehydrins have the molecular prerequisites for higher level binding specificity and regulation, raising new questions about the complexity of their biological function.

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.

Thermodynamics of protein destabilization in live cells

Although protein folding and stability have been well explored under simplified conditions in vitro, it is yet unclear how these basic self-organization events are modulated by the crowded interior of live cells. To find out, we use here in-cell NMR to follow at atomic resolution the thermal unfolding of a β-barrel protein inside mammalian and bacterial cells. Challenging the view from in vitro crowding effects, we find that the cells destabilize the protein at 37 °C but with a conspicuous twist: While the melting temperature goes down the cold unfolding moves into the physiological regime, coupled to an augmented heat-capacity change. The effect seems induced by transient, sequence-specific, interactions with the cellular components, acting preferentially on the unfolded ensemble. This points to a model where the in vivo influence on protein behavior is case specific, determined by the individual protein's interplay with the functionally optimized "interaction landscape" of the cellular interior.

Increased Aggregation Is More Frequently Associated to Human Disease-Associated Mutations Than to Neutral Polymorphisms

Protein aggregation is a hallmark of over 30 human pathologies. In these diseases, the aggregation of one or a few specific proteins is often toxic, leading to cellular degeneration and/or organ disruption in addition to the loss-of-function resulting from protein misfolding. Although the pathophysiological consequences of these diseases are overt, the molecular dysregulations leading to aggregate toxicity are still unclear and appear to be diverse and multifactorial. The molecular mechanisms of protein aggregation and therefore the biophysical parameters favoring protein aggregation are better understood. Here we perform an in silico survey of the impact of human sequence variation on the aggregation propensity of human proteins. We find that disease-associated variations are statistically significantly enriched in mutations that increase the aggregation potential of human proteins when compared to neutral sequence variations. These findings suggest that protein aggregation might have a broader impact on human disease than generally assumed and that beyond loss-of-function, the aggregation of mutant proteins involved in cancer, immune disorders or inflammation could potentially further contribute to disease by additional burden on cellular protein homeostasis.

Protein aggregation has been recognized to contribute to the development of more than 30 human diseases such as Alzheimer and Parkinson disease. Here we have performed an in silico survey of human sequence variations to evaluate whether protein aggregation might impact human disease beyond the above-mentioned aggregation diseases. We find that human disease mutations are more likely to increase the aggregation potential of proteins than non-disease associated mutations. This survey therefore suggests the possibility that protein aggregation is a more widespread disease modifier than previously expected.

Electrostatics, structure prediction, and the energy landscapes for protein folding and binding

While being long in range and therefore weakly specific, electrostatic interactions are able to modulate the stability and folding landscapes of some proteins. The relevance of electrostatic forces for steering the docking of proteins to each other is widely acknowledged, however, the role of electrostatics in establishing specifically funneled landscapes and their relevance for protein structure prediction are still not clear. By introducing Debye-Hückel potentials that mimic long-range electrostatic forces into the Associative memory, Water mediated, Structure, and Energy Model (AWSEM), a transferable protein model capable of predicting tertiary structures, we assess the effects of electrostatics on the landscapes of thirteen monomeric proteins and four dimers. For the monomers, we find that adding electrostatic interactions does not improve structure prediction. Simulations of ribosomal protein S6 show, however, that folding stability depends monotonically on electrostatic strength. The trend in predicted melting temperatures of the S6 variants agrees with experimental observations. Electrostatic effects can play a range of roles in binding. The binding of the protein complex KIX-pKID is largely assisted by electrostatic interactions, which provide direct charge-charge stabilization of the native state and contribute to the funneling of the binding landscape. In contrast, for several other proteins, including the DNA-binding protein FIS, electrostatics causes frustration in the DNA-binding region, which favors its binding with DNA but not with its protein partner. This study highlights the importance of long-range electrostatics in functional responses to problems where proteins interact with their charged partners, such as DNA, RNA, as well as membranes.

Modulation of folding energy landscape by charge-charge interactions: linking experiments with computational modeling

The kinetics of folding-unfolding of a structurally diverse set of four proteins optimized for thermodynamic stability by rational redesign of surface charge-charge interactions is characterized experimentally. The folding rates are faster for designed variants compared with their wild-type proteins, whereas the unfolding rates are largely unaffected. A simple structure-based computational model, which incorporates the Debye-Hückel formalism for the electrostatics, was used and found to qualitatively recapitulate the experimental results. Analysis of the energy landscapes of the designed versus wild-type proteins indicates the differences in refolding rates may be correlated with the degree of frustration of their respective energy landscapes. Our simulations indicate that naturally occurring wild-type proteins have frustrated folding landscapes due to the surface electrostatics. Optimization of the surface electrostatics seems to remove some of that frustration, leading to enhanced formation of native-like contacts in the transition-state ensembles (TSE) and providing a less frustrated energy landscape between the unfolded and TS ensembles. Macroscopically, this results in faster folding rates. Furthermore, analyses of pairwise distances and radii of gyration suggest that the less frustrated energy landscapes for optimized variants are a result of more compact unfolded and TS ensembles. These findings from our modeling demonstrates that this simple model may be used to: (i) gain a detailed understanding of charge-charge interactions and their effects on modulating the energy landscape of protein folding and (ii) qualitatively predict the kinetic behavior of protein surface electrostatic interactions.

Calcium binding to gatekeeper residues flanking aggregation-prone segments underlies non-fibrillar amyloid traits in superoxide dismutase 1 (SOD1)

Calcium deregulation is a central feature among neurodegenerative diseases, including amyotrophic lateral sclerosis (ALS). Calcium accumulates in the spinal and brain stem motor neurons of ALS patients triggering multiple pathophysiological processes which have been recently shown to include direct effects on the aggregation cascade of superoxide dismutase 1 (SOD1). SOD1 is a Cu/Zn enzyme whose demetallated form is implicated in ALS protein deposits, contributing to toxic gain of function phenotypes. Here we undertake a combined experimental and computational study aimed at establishing the molecular details underlying the regulatory effects of Ca(2+) over SOD1 aggregation potential. Isothermal titration calorimetry indicates entropy driven low affinity association of Ca(2+) ions to apo SOD1, at pH7.5 and 37°C. Molecular dynamics simulations denote a noticeable loss of native structure upon Ca(2+) association that is especially prominent at the zinc-binding and electrostatic loops, whose decoupling is known to expose the central SOD1 β-barrel triggering aggregation. Structural mapping of the preferential apo SOD1 Ca(2+) binding locations reveals that among the most frequent ligands for Ca(2+) are negatively-charged gatekeeper residues located in boundary positions with respect to segments highly prone to edge-to-edge aggregation. Calcium interactions thus diminish gatekeeping roles of these residues, by shielding repulsive interactions via stacking between aggregating β-sheets, partly blocking fibril formation and promoting amyloidogenic oligomers such as those found in ALS inclusions. Interestingly, many fALS mutations occur at these positions, disclosing how Ca(2+) interactions recreate effects similar to those of genetic defects, a finding with relevance to understand sporadic ALS pathomechanisms.

Visualization of protein folding funnels in lattice models

Protein folding occurs in a very high dimensional phase space with an exponentially large number of states, and according to the energy landscape theory it exhibits a topology resembling a funnel. In this statistical approach, the folding mechanism is unveiled by describing the local minima in an effective one-dimensional representation. Other approaches based on potential energy landscapes address the hierarchical structure of local energy minima through disconnectivity graphs. In this paper, we introduce a metric to describe the distance between any two conformations, which also allows us to go beyond the one-dimensional representation and visualize the folding funnel in 2D and 3D. In this way it is possible to assess the folding process in detail, e.g., by identifying the connectivity between conformations and establishing the paths to reach the native state, in addition to regions where trapping may occur. Unlike the disconnectivity maps method, which is based on the kinetic connections between states, our methodology is based on structural similarities inferred from the new metric. The method was developed in a 27-mer protein lattice model, folded into a 3×3×3 cube. Five sequences were studied and distinct funnels were generated in an analysis restricted to conformations from the transition-state to the native configuration. Consistent with the expected results from the energy landscape theory, folding routes can be visualized to probe different regions of the phase space, as well as determine the difficulty in folding of the distinct sequences. Changes in the landscape due to mutations were visualized, with the comparison between wild and mutated local minima in a single map, which serves to identify different trapping regions. The extension of this approach to more realistic models and its use in combination with other approaches are discussed.

Water repellency in hydrophobic nanocapsules--molecular view on dewetting

The hydrophobic effect plays a major role in a variety of important phenomena in chemistry, materials science and biology, for instance in protein folding and protein-ligand interactions. Studies--performed within cavities of the unique metal oxide based porous capsules of the type {(pentagon)12(linker)30}≡{(W)W5}12{Mo2(ligand)}30 with different acetate/water ligand ratios--have provided unprecedented results revealing segregation/repellency of the encapsulated "water" from the internal hydrophobic ligand walls of the capsules, while the disordered water molecules, interacting strongly with each other via hydrogen bonding, form in all investigated cases the same type of spherical shell. The present results can be (formally) compared--but only regarding the repellency effect--with the amazing "action" of the (super)hydrophobic Lotus (Nelumbo) leaves, which are self-cleaning based on water repellency resulting in the formation of water droplets picking up dirt. The present results were obtained by constructing deliberately suitable hydrophobic interiors within the mentioned capsules.

A genome-wide sequence-structure analysis suggests aggregation gatekeepers constitute an evolutionary constrained functional class

Protein aggregation is geared by aggregation-prone regions that self-associate by β-strand interactions. Charged residues and prolines are enriched at the flanks of aggregation-prone regions resulting in decreased aggregation. It is still unclear what drives the overrepresentation of these "aggregation gatekeepers", that is, whether their presence results from structural constraints determining protein stability or whether they constitute a bona fide functional class selectively maintained to control protein aggregation. As functional residues are typically conserved regardless of their cost to protein stability, we compared sequence conservation and thermodynamic cost of these residues in 2659 protein families in Escherichia coli. Across protein families, we find gatekeepers to be under strong selective conservation while at the same time representing a significant thermodynamic cost to protein structure. This finding supports the notion that aggregation gatekeepers are not structurally determined but evolutionary selected to control protein aggregation.

Effect of salt on the formation of salt-bridges in β-hairpin peptides

The formation of salt-bridges in β-hairpin peptides is measured in increasing salt concentrations, indicating a decrease in the salt-bridged population due to charge–charge screening, as well as non-cooperative salt-bridge triads.

Large extent of disorder in Adenomatous Polyposis Coli offers a strategy to guard Wnt signalling against point mutations

Mutations in the central region of the signalling hub Adenomatous Polyposis Coli (APC) cause colorectal tumourigenesis. The structure of this region remained unknown. Here, we characterise the Mutation Cluster Region in APC (APC-MCR) as intrinsically disordered and propose a model how this structural feature may contribute to regulation of Wnt signalling by phosphorylation. APC-MCR was susceptible to proteolysis, lacked α-helical secondary structure and did not display thermal unfolding transition. It displayed an extended conformation in size exclusion chromatography and was accessible for phosphorylation by CK1ε in vitro. The length of disordered regions in APC increases with species complexity, from C. elegans to H. sapiens. We speculate that the large disordered region harbouring phosphorylation sites could be a successful strategy to stabilise tight regulation of Wnt signalling against single missense mutations.

Subdomain architecture and stability of a giant repeat protein

Tandem repeat proteins, which are widespread in the human genome, tend to exhibit high stability and favorable expression, and hence, they are emerging as promising protein scaffolds in alternative to antibodies in biotechnology. In order to investigate the origin of the stability of these proteins, we dissect the subdomain architecture of the giant repeat protein PR65/A, which comprises 15 α-helical HEAT repeats, using a series of truncations and deletions. We find that the N (HEAT 1-2) and the C (HEAT 14-15) subdomains are not capable of independent folding, but the addition of HEAT 13 to HEAT 14-15 results in an independently stable C-terminal subdomain (HEAT 13-15), which is in turn further stabilized by the inclusion of HEAT 12 (HEAT 12-15). We also further show that the stability of HEAT 13-15 is enhanced by its fusion to HEAT 1-2, and the artificial 5-HEAT-repeat protein thereby created (HEAT NC) behaves like a cooperative multidomain protein. We construct further variants, lacking one or both of the terminal subdomains, and find that such subdomains function as stabilizing caps within full-length PR65/A as in their absence, the central subdomain of the protein unfolds to form non-native β-sheet-like oligomers. Taken together, our results suggest that in full-length PR65/A, the more unstable regions within the central repeats are protected by the adjacent folded repeats, which thus act as gatekeepers by virtue of their greater stability.

Evolutionary selection for protein aggregation

Protein aggregation is being found to be associated with an increasing number of human diseases. Aggregation can lead to a loss of function (lack of active protein) or to a toxic gain of function (cytotoxicity associated with protein aggregates). Although potentially harmful, protein sequences predisposed to aggregation seem to be ubiquitous in all kingdoms of life, which suggests an evolutionary advantage to having such segments in polypeptide sequences. In fact, aggregation-prone segments are essential for protein folding and for mediating certain protein-protein interactions. Moreover, cells use protein aggregates for a wide range of functions. Against this background, life has adapted to tolerate the presence of potentially dangerous aggregation-prone sequences by constraining and counteracting the aggregation process. In the present review, we summarize the current knowledge of the advantages associated with aggregation-prone stretches in proteomes and the strategies that cellular systems have developed to control the aggregation process.

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

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