Diffusive intracellular interactions: On the role of protein net charge and functional adaptation

A striking feature of nucleic acids and lipid membranes is that they all carry net negative charge and so is true for the majority of intracellular proteins. It is suggested that the role of this negative charge is to assure a basal intermolecular repulsion that keeps the cytosolic content suitably 'fluid' for function. We focus in this review on the experimental, theoretical and genetic findings which serve to underpin this idea and the new questions they raise. Unlike the situation in test tubes, any functional protein-protein interaction in the cytosol is subject to competition from the densely crowded background, i.e. surrounding stickiness. At the nonspecific limit of this stickiness is the 'random' protein-protein association, maintaining profuse populations of transient and constantly interconverting complexes at physiological protein concentrations. The phenomenon is readily quantified in studies of the protein rotational diffusion, showing that the more net negatively charged a protein is the less it is retarded by clustering. It is further evident that this dynamic protein-protein interplay is under evolutionary control and finely tuned across organisms to maintain optimal physicochemical conditions for the cellular processes. The emerging picture is then that specific cellular function relies on close competition between numerous weak and strong interactions, and where all parts of the protein surfaces are involved. The outstanding challenge is now to decipher the very basics of this many-body system: how the detailed patterns of charged, polar and hydrophobic side chains not only control protein-protein interactions at close- and long-range but also the collective properties of the cellular interior as a whole.

Mutant SOD1 aggregates formed in vitro and in cultured cells are polymorphic and differ from those arising in the CNS

Mutations in the human Superoxide dismutase 1 (hSOD1) gene are well-established cause of the motor neuron disease ALS. Patients and transgenic (Tg) ALS model mice carrying mutant variants develop hSOD1 aggregates in the CNS. We have identified two hSOD1 aggregate strains, which both transmit spreading template-directed aggregation and premature fatal paralysis when inoculated into adult transgenic mice. This prion-like spread of aggregation could be a primary disease mechanism in SOD1-induced ALS. Human SOD1 aggregation has been studied extensively both in cultured cells and under various conditions in vitro. To determine how the structure of aggregates formed in these model systems related to disease-associated aggregates in the CNS, we used a binary epitope-mapping assay to examine aggregates of hSOD1 variants G93A, G85R, A4V, D90A, and G127X formed in vitro, in four different cell lines and in the CNS of Tg mice. We found considerable variability between replicate sets of in vitro-generated aggregates. In contrast, there was a high similarity between replicates of a given hSOD1 mutant in a given cell line, but pronounced variations between different hSOD1 mutants and different cell lines in both structures and amounts of aggregates formed. The aggregates formed in vitro or in cultured cells did not replicate the aggregate strains that arise in the CNS. Our findings suggest that the distinct aggregate morphologies in the CNS could result from a micro-environment with stringent quality control combined with second-order selection by spreading ability. Explorations of pathogenesis and development of therapeutics should be conducted in models that replicate aggregate structures forming in the CNS.

Charge Engineering Reveals the Roles of Ionizable Side Chains in Electrospray Ionization Mass Spectrometry

In solution, the charge of a protein is intricately linked to its stability, but electrospray ionization distorts this connection, potentially limiting the ability of native mass spectrometry to inform about protein structure and dynamics. How the behavior of intact proteins in the gas phase depends on the presence and distribution of ionizable surface residues has been difficult to answer because multiple chargeable sites are present in virtually all proteins. Turning to protein engineering, we show that ionizable side chains are completely dispensable for charging under native conditions, but if present, they are preferential protonation sites. The absence of ionizable side chains results in identical charge state distributions under native-like and denaturing conditions, while coexisting conformers can be distinguished using ion mobility separation. An excess of ionizable side chains, on the other hand, effectively modulates protein ion stability. In fact, moving a single ionizable group can dramatically alter the gas-phase conformation of a protein ion. We conclude that although the sum of the charges is governed solely by Coulombic terms, their locations affect the stability of the protein in the gas phase.

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.

Connecting Longitudinal and Transverse Relaxation Rates in Live-Cell NMR

In the cytosolic environment, protein crowding and Brownian motions result in numerous transient encounters. Each such encounter event increases the apparent size of the interacting molecules, leading to slower rotational tumbling. The extent of transient protein complexes formed in live cells can conveniently be quantified by an apparent viscosity, based on NMR-detected spin-relaxation measurements, that is, the longitudinal (T1) and transverse (T2) relaxation. From combined analysis of three different proteins and surface mutations thereof, we find that T2 implies significantly higher apparent viscosity than T1. At first sight, the effect on T1 and T2 seems thus nonunifiable, consistent with previous reports on other proteins. We show here that the T1 and T2 deviation is actually not a inconsistency but an expected feature of a system with fast exchange between free monomers and transient complexes. In this case, the deviation is basically reconciled by a model with fast exchange between the free-tumbling reporter protein and a transient complex with a uniform 143 kDa partner. The analysis is then taken one step further by accounting for the fact that the cytosolic content is by no means uniform but comprises a wide range of molecular sizes. Integrating over the complete size distribution of the cytosolic interaction ensemble enables us to predict both T1 and T2 from a single binding model. The result yields a bound population for each protein variant and provides a quantification of the transient interactions. We finally extend the approach to obtain a correction term for the shape of a database-derived mass distribution of the interactome in the mammalian cytosol, in good accord with the existing data of the cellular composition.

Exposure of a cryptic Hsp70 binding site determines the cytotoxicity of the ALS-associated SOD1-mutant A4V

The accumulation of toxic protein aggregates is thought to play a key role in a range of degenerative pathologies, but it remains unclear why aggregation of polypeptides into non-native assemblies is toxic and why cellular clearance pathways offer ineffective protection. We here study the A4V mutant of SOD1, which forms toxic aggregates in motor neurons of patients with familial amyotrophic lateral sclerosis (ALS). A comparison of the location of aggregation prone regions (APRs) and Hsp70 binding sites in the denatured state of SOD1 reveals that ALS-associated mutations promote exposure of the APRs more than the strongest Hsc/Hsp70 binding site that we could detect. Mutations designed to increase the exposure of this Hsp70 interaction site in the denatured state promote aggregation but also display an increased interaction with Hsp70 chaperones. Depending on the cell type, in vitro this resulted in cellular inclusion body formation or increased clearance, accompanied with a suppression of cytotoxicity. The latter was also observed in a zebrafish model in vivo. Our results suggest that the uncontrolled accumulation of toxic SOD1A4V aggregates results from insufficient detection by the cellular surveillance network.

Diffusive protein interactions in human versus bacterial cells

Random encounters between proteins in crowded cells are by no means passive, but found to be under selective control. This control enables proteome solubility, helps to optimise the diffusive search for interaction partners, and allows for adaptation to environmental extremes. Interestingly, the residues that modulate the encounters act mesoscopically through protein surface hydrophobicity and net charge, meaning that their detailed signatures vary across organisms with different intracellular constraints. To examine such variations, we use in-cell NMR relaxation to compare the diffusive behaviour of bacterial and human proteins in both human and Escherichia coli cytosols. We find that proteins that ‘stick’ in E. coli are generally less restricted in mammalian cells. Furthermore, the rotational diffusion in the mammalian cytosol is less sensitive to surface-charge mutations. This implies that, in terms of protein motions, the mammalian cytosol is more forgiving to surface alterations than E. coli cells. The cellular differences seem not linked to the proteome properties per se, but rather to a 6-fold difference in protein concentrations. Our results outline a scenario in which the tolerant cytosol of mammalian cells, found in long-lived multicellular organisms, provides an enlarged evolutionary playground, where random protein-surface mutations are less deleterious than in short-generational bacteria.

•

Random protein encounters and diffusibility in cells are controlled by surface charge.

•

Protein rotational diffusion is less restricted in human cells than in E. coli.

•

Human cells are less sensitive to alterations of protein charge than Escherichia coli cells.

The Pierced Lasso Topology Leptin has a Bolt on Dynamic Domain Composed by the Disordered Loops I and III

Leptin is an important signaling hormone, mostly known for its role in energy expenditure and satiety. Furthermore, leptin plays a major role in other proteinopathies, such as cancer, marked hyperphagia, impaired immune function, and inflammation. In spite of its biological relevance in human health, there are no NMR resonance assignments of the human protein available, obscuring high-resolution characterization of the soluble protein and/or its conformational dynamics, suggested as being important for receptor interaction and biological activity. Here, we report the nearly complete backbone resonance assignments of human leptin. Chemical shift-based secondary structure prediction confirms that in solution leptin forms a four-helix bundle including a pierced lasso topology. The conformational dynamics, determined on several timescales, show that leptin is monomeric, has a rigid four-helix scaffold, and a dynamic domain, including a transiently formed helix. The dynamic domain is anchored to the helical scaffold by a secondary hydrophobic core, pinning down the long loops of leptin to the protein body, inducing motional restriction without a well-defined secondary or tertiary hydrogen bond stabilized structure. This dynamic region is well suited for and may be involved in functional allosteric dynamics upon receptor binding.

The Cost of Long Catalytic Loops in Folding and Stability of the ALS-Associated Protein SOD1

A conspicuous feature of the amyotrophic lateral sclerosis (ALS)-associated protein SOD1 is that its maturation into a functional enzyme relies on local folding of two disordered loops into a catalytic subdomain. To drive the disorder-to-order transition, the protein employs a single Zn²⁺ ion. The question is then if the entropic penalty of maintaining such disordered loops in the immature apoSOD1 monomer is large enough to explain its unusually low stability, slow folding, and pathological aggregation in ALS. To find out, we determined the effects of systematically altering the SOD1-loop lengths by protein redesign. The results show that the loops destabilize the apoSOD1 monomer by ∼3 kcal/mol, rendering the protein marginally stable and accounting for its aggregation behavior. Yet the effect on the global folding kinetics remains much smaller with a transition-state destabilization of <1 kcal/mol. Notably, this 1/3 transition-state to folded-state stability ratio provides a clear-cut example of the enigmatic disagreement between the Leffler α value from loop-length alterations (typically 1/3) and the "standard" reaction coordinates based on solvent perturbations (typically >2/3). Reconciling the issue, we demonstrate that the disagreement disappears when accounting for the progressive loop shortening that occurs along the folding pathway. The approach assumes a consistent Flory loop entropy scaling factor of c = 1.48 for both equilibrium and kinetic data and has the added benefit of verifying the tertiary interactions of the folding nucleus as determined by phi-value analysis. Thus, SOD1 not only represents a case where evolution of key catalytic function has come with the drawback of a destabilized apo state but also stands out as a well-suited model system for exploring the physicochemical details of protein self-organization.

Diffuse binding of Zn(2+) to the denatured ensemble of Cu/Zn superoxide dismutase 1

•

Nearly complete backbone assignments for monomeric SOD1 were obtained.

•

Monomeric SOD1 is unstructured in 9 M urea.

•

Zn²⁺-binding to monomeric SOD1 occurs through diffuse coordination to most His residues.

•

The binding geometry of Zn²⁺ is different than in the native, folded SOD1.

The stability and structural properties of the metalloprotein superoxide dismutase 1 (SOD1) are found to depend critically on metal ions. Native SOD1 monomers coordinate one structural Zn²⁺ and one redox-active Cu^2+/1+ to the active site. To do this, the Zn²⁺ ions need to interact with the SOD1 protein on the denatured side of the folding barrier, prior to the formation of the folding nucleus. In this study, we have examined at residue level the nature of this early Zn²⁺ binding by NMR studies on the urea denatured-state of SOD1. Nearly complete backbone chemical shift assignments were obtained in 9 M urea at physiological pH, conditions at which NMR studies are scarce. Our results demonstrate that SOD1 is predominantly unstructured under these conditions. Chemical-shift changes upon Zn²⁺ titration show that denatured SOD1 retains a significant affinity to Zn²⁺ ions, even in 9 M urea. However, the Zn²⁺ interactions are not limited to the native metal-binding ligands in the two binding sites, but are seen for all His residues. Moreover, the native Cu^2+/1+ ligand H46 seems not to bind as well as the other His residues, while the nearby non-native H43 does bind, indicating that the binding geometry is relaxed. The result suggests that the Zn²⁺-binding observed to catalyze folding of SOD1 in physiological buffer is initiated by diffuse, non-specific coordination to the coil, which subsequently funnels by ligand exchange into the native coordination geometry of the folded monomer. Altogether, this diffuse binding is a result with fundamental implications for folding of metalloproteins in general.

Direct detection of neuropeptide dynorphin A binding to the second extracellular loop of the κ opioid receptor using a soluble protein scaffold

The molecular determinants for selectivity of ligand binding to membrane receptors are of key importance for the understanding of cellular signalling, as well as for rational therapeutic intervention. In the present study, we target the interaction between the κ opioid receptor (KOR) and its native peptide ligand dynorphin A (DynA) using solution state NMR spectroscopy, which is generally made difficult by the sheer size of membrane bound receptors. Our method is based on 'transplantation' of an extracellular loop of KOR into a 'surrogate' scaffold; in this case, a soluble β-barrel. Our results corroborate the general feasibility of the method, showing that the inserted receptor segment has negligible effects on the properties of the scaffold protein, at the same time as maintaining an ability to bind its native DynA ligand. Upon DynA binding, only small induced chemical shift changes of the KOR loop were observed, whereas chemical shift changes of DynA and NMR paramagnetic relaxation data show conclusively that the peptide interacts with the inserted loop. The binding interface is composed of a disordered part of the KOR loop and involves both electrostatic and hydrophobic interactions. Even so, simultaneous effects along the DynA sequence upon binding show that control of the recognition is a concerted event.

Pruning the ALS-associated protein SOD1 for in-cell NMR

To efficiently deliver isotope-labeled proteins into mammalian cells poses a main challenge for structural and functional analysis by in-cell NMR. In this study we have employed cell-penetrating peptides (CPPs) to deliver the ALS-associated protein superoxide dismutase (SOD1) into HeLa cells. Our results show that, although full-length SOD1 cannot be efficiently internalized, a variant in which the active-site loops IV and VII have been truncated (SOD1(ΔIVΔVII)) yields high cytosolic delivery. The reason for the enhanced delivery of SOD1(ΔIVΔVII) seems to be the elimination of negatively charged side chains, which alters the net charge of the CPP-SOD1 complex from neutral to +4. The internalized SOD1(ΔIVΔVII) protein displays high-resolution in-cell NMR spectra similar to, but not identical to, those of the lysate of the cells. Spectral differences are found mainly in the dynamic β strands 4, 5, and 7, triggered by partial protonation of the His moieties of the Cu-binding site. Accordingly, SOD1(ΔIVΔVII) doubles here as an internal pH probe, revealing cytosolic acidification under the experimental treatment. Taken together, these observations show that CPP delivery, albeit inefficient at first trials, can be tuned by protein engineering to allow atomic-resolution NMR studies of specific protein structures that have evaded other in-cell NMR approaches: in this case, the structurally elusive apoSOD1 barrel implicated as precursor for misfolding in ALS.

Trimming down a protein structure to its bare foldons: spatial organization of the cooperative unit

Folding of the ribosomal protein S6 is a malleable process controlled by two competing, and partly overlapping, folding nuclei. Together, these nuclei extend over most of the S6 structure, except the edge strand β2, which is consistently missing in the folding transition states; despite being part of the S6 four-stranded sheet, β2 seems not to be part of the cooperative unit of the protein. The question is then whether β2 can be removed from the S6 structure without compromising folding cooperativity or native state integrity. To investigate this, we constructed a truncated variant of S6 lacking β2, reducing the size of the protein from 96 to 76 residues (S6(Δβ2)). The new S6 variant expresses well in Escherichia coli and has a well dispersed heteronuclear single quantum correlation spectrum and a perfectly wild-type-like crystal structure, but with a smaller three-stranded β-sheet. Moreover, S6(Δβ2) displays an archetypical v-shaped chevron plot with decreased slope of the unfolding limb, as expected from a protein with maintained folding cooperativity and reduced size. The results support the notion that foldons, as defined by the structural distribution of the folding nuclei, represent a property-based level of hierarchy in the build-up of larger protein structures and suggest that the role of β2 in S6 is mainly in intermolecular binding, consistent with the position of this strand in the ribosomal assembly.

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.

Folding catalysis by transient coordination of Zn2+ to the Cu ligands of the ALS-associated enzyme Cu/Zn superoxide dismutase 1

How coordination of metal ions modulates protein structures is not only important for elucidating biological function but has also emerged as a key determinant in protein turnover and protein-misfolding diseases. In this study, we show that the coordination of Zn(2+) to the ALS-associated enzyme Cu/Zn superoxide dismutase (SOD1) is directly controlled by the protein's folding pathway. Zn(2+) first catalyzes the folding reaction by coordinating transiently to the Cu ligands of SOD1, which are all contained within the folding nucleus. Then, after the global folding transition has commenced, the Zn(2+) ion transfers to the higher affinity Zn site, which structures only very late in the folding process. Here it remains dynamically coordinated with an off rate of ∼10(-5) s(-1). This relatively rapid equilibration of metals in and out of the SOD1 structure provides a simple explanation for how the exceptionally long lifetime, >100 years, of holoSOD1 is still compatible with cellular turnover: if a dissociated Zn(2+) ion is prevented from rebinding to the SOD1 structure then the lifetime of the protein is reduced to a just a few hours.

SOD1 mutations targeting surface hydrogen bonds promote amyotrophic lateral sclerosis without reducing apo-state stability

In good accord with the protein aggregation hypothesis for neurodegenerative diseases, ALS-associated SOD1 mutations are found to reduce structural stability or net repulsive charge. Moreover there are weak indications that the ALS disease progression rate is correlated with the degree of mutational impact on the apoSOD1 structure. A bottleneck for obtaining more conclusive information about these structure-disease relationships, however, is the large intrinsic variability in patient survival times and insufficient disease statistics for the majority of ALS-provoking mutations. As an alternative test of the structure-disease relationship we focus here on the SOD1 mutations that appear to be outliers in the data set. The results identify several ALS-provoking mutations whose only effect on apoSOD1 is the elimination or introduction of a single charge, i.e. D76V/Y, D101N, and N139D/K. The thermodynamic stability and folding behavior of these mutants are indistinguishable from the wild-type control. Moreover, D101N is an outlier in the plot of stability loss versus patient survival time by having rapid disease progression. Common to the identified mutations is that they truncate conserved salt-links and/or H-bond networks in the functional loops IV or VII. The results show that the local impact of ALS-associated mutations on the SOD1 molecule can sometimes overrun their global effects on apo-state stability and net repulsive charge, and point at the analysis of property outliers as an efficient strategy for mapping out new ALS-provoking features.

Solution structures and backbone dynamics of the ribosomal protein S6 and its permutant P(54-55)

The ribosomal protein S6 from Thermus thermophilus has served as a model system for the study of protein folding, especially for understanding the effects of circular permutations of secondary structure elements. This study presents the structure of a permutant protein, the 96-residue P(54-55), and the structure of its 101-residue parent protein S6(wt) in solution. The data also characterizes the effects of circular permutation on the backbone dynamics of S6. Consistent with crystallographic data on S6(wt), the overall solution structures of both P(54-55) and S6(wt) show a beta-sheet of four antiparallel beta-strands with two alpha-helices packed on one side of the sheet. In clear contrast to the crystal data, however, the solution structure of S6(wt) reveals a disordered loop in the region between beta-strands 2 and 3 (Leu43-Phe60) instead of a well-ordered stretch and associated hydrophobic mini-core observed in the crystal structure. Moreover, the data for P(54-55) show that the joined wild-type N- and C-terminals form a dynamically robust stretch with a hairpin structure that complies with the in silico design. Taken together, the results explain why the loop region of the S6(wt) structure is relatively insensitive to mutational perturbations, and why P(54-55) is more stable than S6(wt): the permutant incision at Lys54-Asp55 is energetically neutral by being located in an already disordered loop whereas the new hairpin between the wild-type N- and C-termini is stabilizing.

Distinguishing between smooth and rough free energy barriers in protein folding

Analysis of curved chevron plots is a powerful tool in investigating protein folding pathways, as the curvatures can be used to gain information about both early and late folding events. When and if accumulation of low-energy intermediates can be ruled out, two different models have classically been applied to describe curved chevron plots, namely , (i) Hammond effects along smooth barrier profiles and (ii) changes in the rate-limiting step between two discrete transition states. The two models lead to very similar numerical solutions, which are generally indistinguishable. This is not surprising, since the smooth barrier assumption approximates barrier profiles with a more complex topology involving multiple local maxima that are too close, or too broad, to yield clear-cut kinks in the chevron data. In this work, we have reconstructed the transition state shifts as a function of protein stability over a wide stability range for three small globular proteins, to screen for fingerprints more sensitive for different barrier profiles. We show that such an analysis represents a valuable test for the discrimination between the two different scenarios.

Transient structural distortion of metal-free Cu/Zn superoxide dismutase triggers aberrant oligomerization

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease linked to the misfolding of Cu/Zn superoxide dismutase (SOD1). ALS-related defects in SOD1 result in a gain of toxic function that coincides with aberrant oligomerization. The structural events triggering oligomerization have remained enigmatic, however, as is the case in other protein-misfolding diseases. Here, we target the critical conformational change that defines the earliest step toward aggregation. Using nuclear spin relaxation dispersion experiments, we identified a short-lived (0.4 ms) and weakly populated (0.7%) conformation of metal-depleted SOD1 that triggers aberrant oligomerization. This excited state emanates from the folded ground state and is suppressed by metal binding, but is present in both the disulfide-oxidized and disulfide-reduced forms of the protein. Our results pinpoint a perturbed region of the excited-state structure that forms intermolecular contacts in the earliest nonnative dimer/oligomer. The conformational transition that triggers oligomerization is a common feature of WT SOD1 and ALS-associated mutants that have widely different physicochemical properties. But compared with WT SOD1, the mutants have enhanced structural distortions in their excited states, and in some cases slightly higher excited-state populations and lower kinetic barriers, implying increased susceptibility to oligomerization. Our results provide a unified picture that highlights both (i) a common denominator among different SOD1 variants that may explain why diverse mutations cause the same disease, and (ii) a structural basis that may aid in understanding how different mutations affect disease propensity and progression.

SOD1-associated ALS: a promising system for elucidating the origin of protein-misfolding disease

Amyotropic lateral sclerosis (ALS) is a neurodegenerative disease linked to misfolding and aggregation of the homodimeric enzyme superoxide dismutase (SOD1). In contrast to the precursors of other neurodegenerative diseases, SOD1 is a soluble and simple-to-study protein with immunoglobulin-like structure. Also, there are more than 120 ALS-provoking SOD1 mutations at the disposal for detailed elucidation of the disease-triggering factors at molecular level. In this article, we review recent progress in the characterization of the folding and assembly pathway of the SOD1 dimer and how this is affected by ALS-provoking mutations. Despite the diverse nature of these mutations, the results offer so far a surprising simplicity. The ALS-provoking mutations decrease either protein stability or net repulsive charge: the classical hallmarks for a disease mechanism triggered by association of non-native protein. In addition, the mutant data identifies immature SOD1 monomers as the species from which the cytotoxic pathway emerges, and point at compromised folding cooperativity as a key disease determinant. The relative ease by which these data can be obtained makes SOD1 a promising model for elucidating also the origin of other neurodegenerative diseases where the precursor proteins are structurally more elusive.

Soluble misfolded subfractions of mutant superoxide dismutase-1s are enriched in spinal cords throughout life in murine ALS models

Mutants of superoxide dismutase-1 (SOD1) cause ALS by an unidentified cytotoxic mechanism. We have previously shown that the stable SOD1 mutants D90A and G93A are abundant and show the highest levels in liver and kidney in transgenic murine ALS models, whereas the unstable G85R and G127X mutants are scarce but enriched in the CNS. These data indicated that minute amounts of misfolded SOD1 enriched in the motor areas might exert the ALS-causing cytotoxicity. A hydrophobic interaction chromatography (HIC) protocol was developed with the aim to determine the abundance of soluble misfolded SOD1 in tissues in vivo. Most G85R and G127X mutant SOD1s bound in the assay, but only minute subfractions of the D90A and G93A mutants. The absolute levels of HIC-binding SOD1 were, however, similar and broadly inversely related to lifespans in the models. They were generally enriched in the susceptible spinal cord. The HIC-binding SOD1 was composed of disulfide-reduced subunits lacking metal ions and also subunits that apparently carried nonnative intrasubunit disulfide bonds. The levels were high from birth until death and were comparable to the amounts of SOD1 that become sequestered in aggregates in the terminal stage. The HIC-binding SOD1 species ranged from monomeric to trimeric in size. These species form a least common denominator amongst SOD1 mutants with widely different molecular characteristics and might be involved in the cytotoxicity that causes ALS.

Amyotrophic lateral sclerosis-associated copper/zinc superoxide dismutase mutations preferentially reduce the repulsive charge of the proteins

We provide bioinformatical evidence that protein charge plays a key role in the disease mechanism of amyotrophic lateral sclerosis (ALS). Analysis of 100 ALS-associated mutations in copper/zinc superoxide dismutase (SOD1) shows that these are site-selective with a preference to decrease the proteins' net repulsive charge. For each SOD1 monomer this charge is normally -6. Because biomolecules as a rule maintain net negative charge to assure solubility in the cellular interior, the result lends support to the hypothesis of protein aggregation as an initiating event in the ALS pathogenesis. The strength of the preferential reduction of repulsive charge is higher in SOD1-associated ALS than in other inherited protein disorders.

Malleability of protein folding pathways: a simple reason for complex behaviour

Although the structures of native proteins are generally unique, the pathways by which they form are often free to vary. Some proteins fold by a multitude of different pathways, whereas others seem restricted to only one choice. An explanation for this variation in folding behaviour has recently emerged from studies of transition state changes: the number of accessible pathways is linked to the number of nucleation motifs contained within the native topology. We refer to these nucleation motifs as 'foldons', as they approach the size of an independent cooperative unit. Thus, with respect to pathway malleability and the composition of the folding funnel, proteins can be seen as modular assemblies of competing foldons. For the split beta-alpha-beta fold, these foldons are two-strand-helix motifs coupled by spatial overlap.

The coupling between disulphide status, metallation and dimer interface strength in Cu/Zn superoxide dismutase

The gain of neurotoxic function in amyotrophic lateral sclerosis (ALS) has been linked to misfolding of the homodimeric enzyme Cu/Zn superoxide dismutase (SOD). Here, we present the crystal structure of fully cysteine-depleted human SOD (SOD(CallA)), representing a reduced, marginally stable intermediate on the folding pathway in vivo that has also been implicated as neurotoxic precursor state. A hallmark of this species is that it fails to dimerize and becomes trapped as a monomer in the absence of the active-site metals. The crystallographic data show that removal of the C57-C146 disulphide bond sets free the interface loop IV in the apo protein, whereas the same loop remains unaffected in the holo protein. Thus, the low dimerisation propensity of disulphide-reduced apoSOD seems to be of entropic origin due to increased loop flexibility in the monomeric state: in the disulphide-reduced holo protein this gain in configurational entropy upon splitting of the dimer interface is reduced by the metal coordination.

Folding of S6 structures with divergent amino acid composition: pathway flexibility within partly overlapping foldons

Studies of circular permutants have demonstrated that the folding reaction of S6 from Thermus thermophilus (S6(T)) is malleable and responds in an ordered manner to changes of the sequence separation between interacting residues: the S6(T) permutants retain a common nucleation pattern in the form of a two-strand-helix motif that can be recruited from different parts of the structure. To further test the robustness of the two-strand-helix nucleus we have here determined the crystal structure and folding reaction of an evolutionary divergent S6 protein from the hyperthermophilic bacterium Aquifex aeolicus (S6(A)). Although the overall topology of S6(A) is very similar to that of S6(T) the architecture of the hydrophobic core is radically different by containing a large proportion of stacked Phe side-chains. Despite this disparate core composition, the folding rate constant and the kinetic m values of S6(A) are identical to those of S6(T). The folding nucleus of S6(A) is also found to retain the characteristic two-strand-helix motif of the S6(T) permutants, but with a new structural emphasis. The results suggest that the protein folding reaction is linked to topology only in the sense that the native-state topology determines the repertoire of accessible nucleation motifs. If the native structure allows several equivalent ways of recruiting a productive nucleus the folding reaction is free to redistribute within these topological constraints.

Folding of Cu/Zn superoxide dismutase suggests structural hotspots for gain of neurotoxic function in ALS: parallels to precursors in amyloid disease

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease linked to misfolding of the ubiquitous enzyme Cu/Zn superoxide dismutase (SOD). In contrast to other protein-misfolding disorders with similar neuropathogenesis, ALS is not always associated with the in vivo deposition of protein aggregates. Thus, under the assumption that all protein-misfolding disorders share at primary level a similar disease mechanism, ALS constitutes an interesting disease model for identifying the yet-mysterious precursor states from which the cytotoxic pathway emerges. In this study, we have mapped out the conformational repertoire of the apoSOD monomer through analysis of its folding behavior. The results allow us to target the regions of the SOD structure that are most susceptible to unfolding locally under physiological conditions, leading to the exposure of structurally promiscuous interfaces that are normally hidden in the protein's interior. The structure of this putative ALS precursor is strikingly similar to those implicated in amyloid disease.

The experimental survey of protein-folding energy landscapes

We review what has been learned about the protein-folding problem from experimental kinetic studies. These studies reveal patterns of both great richness and surprising simplicity. The patterns can be interpreted in terms of proteins possessing an energy landscape which is largely, but not completely, funnel-like. Issues such as speed limitations of folding, the robustness of folding, the origin of barriers and cooperativity and the ensemble nature of transition states, intermediate and traps are assessed using the results from several experimental groups highlighting energy-landscape ideas as an interpretive framework.

Common Motifs and Topological Effects in the Protein Folding Transition State

Through extensive experiment, simulation, and analysis of protein S6 (1RIS), we find that variations in nucleation and folding pathway between circular permutations are determined principally by the restraints of topology and specific nucleation, and affected by changes in chain entropy. Simulations also relate topological features to experimentally measured stabilities. Despite many sizable changes in phi values and the structure of the transition state ensemble that result from permutation, we observe a common theme: the critical nucleus in each of the mutants share a subset of residues that can be mapped to the critical nucleus residues of the wild-type. Circular permutations create new N and C termini, which are the location of the largest disruption of the folding nucleus, leading to a decrease in both phi values and the role in nucleation. Mutant nuclei are built around the wild-type nucleus but are biased towards different parts of the S6 structure depending on the topological and entropic changes induced by the location of the new N and C termini.

Tryptophan-BODIPY: a versatile donor-acceptor pair for probing generic changes of intraprotein distances

We demonstrate that Tryptophan (Trp) and N-(4,4-difluoro-5,7-dimethyl-4-bora-3a,4a-diaza-s-indacene-3-yl)methyl iodoacetamide (BODIPY) is a suitable donor-acceptor (D-A) pair for intraprotein distance measurements, applicable to the study of protein folding. The suitability of the Trp-BODIPY electronic energy transfer is exemplified on the extensively-characterised two-state protein, S6, from Thermus thermophilus. This protein has proved to be useful for the elucidation of folding cooperativity and nucleation, as well as the changes upon induction of structural transitions. For a comprehensive structural coverage, BODIPY molecules were anchored by Cys insertions at four different positions on the S6 surface. Trp residues at position 33 or 62 acted as donors of electronic energy to the BODIPY groups. None of the D-A pairs show any detectable difference in the folding kinetics (or protein stability), which supports the notion that the two-state transition of S6 is a highly concerted process. Similar results are obtained for mutants affecting the N- and C-terminus. The kinetic analyses indicate that changes of the transition state occur through local unfolding of the native state, rather than by a decrease of the folding cooperativity. The distances obtained from the analysis of the time-resolved fluorescence experiments in the native state were compared to those calculated from X-ray structure. As an additional measure, molecular dynamics simulations of the different protein constructs were performed to account for variability in the BODIPY location on the protein surface. The agreement between fluorescence and X-ray data is quite convincing, and shows that energy transfer measurements between Trp and BODIPY can probe distances between ca. 17 to 34 A, with an error better than 10%.

Identification of the minimal protein-folding nucleus through loop-entropy perturbations

To explore the plasticity and structural constraints of the protein-folding nucleus we have constructed through circular permutation four topological variants of the ribosomal protein S6. In effect, these topological variants represent entropy mutants with maintained spatial contacts. The proteins were characterized at two complementary levels of detail: by phi-value analysis estimating the extent of contact formation in the transition-state ensemble and by Hammond analysis measuring the site-specific growth of the folding nucleus. The results show that, although the loop-entropy alterations markedly influence the appearance and structural location of the folding nucleus, it retains a common motif of one helix docking against two strands. This nucleation motif is built around a shared subset of side chains in the center of the hydrophobic core but extends in different directions of the S6 structure following the permutant-specific differences in local loop entropies. The adjustment of the critical folding nucleus to alterations in loop entropies is reflected by a direct correlation between the phi-value change and the accompanying change in local sequence separation.

Disulphide-reduced superoxide dismutase-1 in CNS of transgenic amyotrophic lateral sclerosis models

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease afflicting the voluntary motor system. More than 100 different mutations in the ubiquitously expressed enzyme superoxide dismutase-1 (SOD1) have been associated with the disease. To search for the nature of the cytotoxicity of mutant SOD1s, amounts, enzymic activities and structural properties of the protein as well as the CNS histopathology were examined in multiple transgenic murine models. In order to generate the ALS phenotype within the short lifespan of the mouse, more than 20-fold increased rates of synthesis of mutant SOD1s appear to be required. The organs of transgenic mice expressing human wild-type SOD1 or either of the G93A and D90A mutant proteins showed high steady-state protein levels. The major proportion of these SOD1s in the CNS were inactive due to insufficient Cu charging and all contained subfractions with a reduced C57-C146 intrasubunit disulphide bond. Both G85R and the truncated G127insTGGG mutant showed low steady-state protein levels, lacked enzyme activity and had no C57-C146 disulphide bond. These mutants were also enriched in the CNS relative to other organs, suggesting inefficient recognition and degradation of misfolded disulphide-reduced SOD1 in susceptible tissues. In end-stage disease, despite 35-fold differences in levels of mutant SOD1s, similar amounts of detergent-resistant aggregates accumulated in the spinal cord. Small granular as well as larger more diffuse human SOD1 (hSOD1)-inclusions developed in all strains, the latter more pronounced in those with high hSOD1 levels. Widespread vacuolizations were seen in the strains with high levels of hSOD1 but not those with low, suggesting these alterations to be artefacts related to high hSOD1 levels and not to the ALS-causing cytotoxicity. The findings suggest that the motoneuron degeneration could be due to long-term exposure to misfolded aggregation-prone disulphide-reduced SOD1, which constitutes minute subfractions of the stable mutants and larger proportions of the unstable mutants.

Systematically perturbed folding patterns of amyotrophic lateral sclerosis (ALS)-associated SOD1 mutants

Amyotrophic lateral sclerosis is a neurodegenerative syndrome associated with 114 mutations in the gene encoding the cytosolic homodimeric enzyme Cu/Zn superoxide dismutase (SOD). In this article, we report that amyotrophic lateral sclerosis-associated SOD mutations with distinctly different disease progression can be rationalized in terms of their folding patterns. The mutations are found to perturb the protein in multiple ways; they destabilize the precursor monomers (class 1), weaken the dimer interface (class 2), or both at the same time (class 1 + 2). A shared feature of the mutational perturbations is a shift of the folding equilibrium toward poorly structured SOD monomers. We observed a link, coupled to the altered folding patterns, between protein stability, net charge, and survival time for the patients carrying the mutations.

Scanning malleable transition state ensembles: comparing theory and experiment for folding protein U1A

Using a variational free energy functional, we calculate the characteristics of the transition state ensembles (TSE) for the folding of protein U1A and investigate how they respond to thermal and mutational changes. The functional directly yields predicted chevron plots both for the wild-type protein and for various mutants. The detailed variations of the TSE and changes in chevron plots predicted by the theory agree reasonably well with the results of the experiments. We also show how to visualize the folding nuclei using 3D isodensity plots.

Correspondence between anomalous m- and DeltaCp-values in protein folding

Proteins folding according to a classical two-state system characteristically show V-shaped chevron plots. We have previously interpreted the symmetrically curved chevron plot of the protein U1A as denaturant-dependent movements in the position of the transition state ensemble (TSE). S6, a structural analog of U1A, shows a classical V-shaped chevron plot indicative of straightforward two-state kinetics, but the mutant LA30 has a curved unfolding limb, which is most consistent with TSE mobility. The kinetic m-values (derivatives of the rate constants with respect to denaturant concentration) in themselves depend on denaturant concentration. To obtain complementary information about putative mobile TSEs, we have carried out a thermodynamic analysis of the three proteins, based on data for refolding and unfolding over the range 10 degrees C to 70 degrees C. The data at all temperatures can be fitted to two-state model systems. Importantly, for all three proteins the activation heat capacities are, within error, identical to the heat capacities measured in independent experiments under equilibrium conditions. Although the equilibrium heat capacities are essentially invariant with regard to denaturant concentration, the activation heat capacities, similar to the structurally equivalent kinetic m-values, show marked denaturant dependence. Furthermore, the values of beta++ at different denaturant concentrations measured by m-values and by heat capacity values are very similar. These observations are consistent with significant transition state movements within the framework of two-state folding. The basis for TSE movement appears to be enthalpic rather than entropic, suggesting that the binding energy of denaturant-protein interactions is a major determinant of the response of energy landscape contours to changing environments.

Protein folding: defining a "standard" set of experimental conditions and a preliminary kinetic data set of two-state proteins

Recent years have seen the publication of both empirical and theoretical relationships predicting the rates with which proteins fold. Our ability to test and refine these relationships has been limited, however, by a variety of difficulties associated with the comparison of folding and unfolding rates, thermodynamics, and structure across diverse sets of proteins. These difficulties include the wide, potentially confounding range of experimental conditions and methods employed to date and the difficulty of obtaining correct and complete sequence and structural details for the characterized constructs. The lack of a single approach to data analysis and error estimation, or even of a common set of units and reporting standards, further hinders comparative studies of folding. In an effort to overcome these problems, we define here a "consensus" set of experimental conditions (25 degrees C at pH 7.0, 50 mM buffer), data analysis methods, and data reporting standards that we hope will provide a benchmark for experimental studies. We take the first step in this initiative by describing the folding kinetics of 30 apparently two-state proteins or protein domains under the consensus conditions. The goal of our efforts is to set uniform standards for the experimental community and to initiate an accumulating, self-consistent data set that will aid ongoing efforts to understand the folding process.

Transient formation of nano-crystalline structures during fibrillation of an Abeta-like peptide

During the first few minutes of fibrillation of a 14-residue peptide homologous to the hydrophobic C-terminal part of the Abeta-peptide, EM micrographs reveal small crystalline areas (100 to 150 nm, repeating unit 47 A) scattered in more amorphous material. On a longer time scale, these crystalline areas disappear and are replaced by tangled clusters resembling protofilaments (hours), and eventually by more regular amyloid fibrils of 60 A to 120 A diameter (days). The transient population of the crystalline areas indicates the presence of ordered substructures in the early fibrillation process, the diameter of which matches the length of the 14-mer peptide in an extended beta-strand conformation.

Folding of human superoxide dismutase: disulfide reduction prevents dimerization and produces marginally stable monomers

The molecular mechanism by which the homodimeric enzyme Cu/Zn superoxide dismutase (SOD) causes neural damage in amytrophic lateral sclerosis is yet poorly understood. A striking, as well as an unusual, feature of SOD is that it maintains intrasubunit disulfide bonds in the reducing environment of the cytosol. Here, we investigate the role of these disulfide bonds in folding and assembly of the SOD apo protein (apoSOD) homodimer through extensive protein engineering. The results show that apoSOD folds in a simple three-state process by means of two kinetic barriers: 2D<==>2M<==>M(2). The early predominant barrier represents folding of the monomers (M), and the late barrier the assembly of the dimer (M(2)). Unique for this mechanism is a dependence of protein concentration on the unfolding rate constant under physiological conditions, which disappears above 6 M Urea where the transition state for unfolding shifts to first-order dissociation of the dimer in accordance with Hammond-postulate behavior. Although reduction of the intrasubunit disulfide bond C57-C146 is not critical for folding of the apoSOD monomer, it has a pronounced effect on its stability and abolishes subsequent dimerization. Thus, impaired ability to form, or retain, the C57-C146 bond in vivo is predicted to increase the cellular load of marginally stable apoSOD monomers, which may have implications for the amytrophic lateral sclerosis neuropathology.

Transient Aggregation and Stable Dimerization Induced by Introducing an Alzheimer Sequence into a Water-Soluble Protein

Transient contacts between denatured polypeptide chains are likely to play an important part in the initial stages of protein aggregation and fibrillation. To analyze the nature of such contacts, we have carried out a protein engineering study of the 102-residue protein U1A, which aggregates transiently in the wild-type form during refolding from the guanidinium chloride-denatured state. We have prepared a series of mutants with increased aggregation tendencies by increasing the homology between two beta-strands of U1A and the Alzheimer peptide (beta-AP). These mutants undergo transient aggregation during refolding, as measured by concentration dependence, double-jump experiments, and binding of ANS, a probe for exposed hydrophobic patches on protein surfaces. The propensity to aggregate increases with increasing homology to beta-AP. Further, the degree of transient ANS binding correlates reasonably well with the structural parameters recently shown to play a role in the fibrillation of natively unfolded proteins. Two mutants highly prone to transient aggregation, U1A-J and U1A-G, were also studied by NMR. Secondary structural elements of the U1A-J construct (with lower beta-AP homology) are very similar to those observed in U1A-wt. In contrast, the high-homology construct U1A-G exhibits local unfolding of the C-terminal helix, which packs against the beta-sheet in the wild-type protein. U1A-G is mainly dimeric according to (15)N spin relaxation data, and the dimer interface most likely involves the beta-sheet. Our data suggest that the transient aggregate relies on specific intermolecular interactions mediated by structurally flexible regions and that contacts may be formed in different beta-strand registers.

Scattered Hammond plots reveal second level of site-specific information in protein folding: phi' (beta++)

Site-specific information about structural heterogeneities of the protein-folding transition-state ensemble is commonly derived from the scatter of the Brønsted plot through the individual values of = Delta logk(f)/Delta logK(D-N). Here, we provide a second level of site-specific detail in the transition-state analysis by demonstrating that the scatter of the Hammond plot is related to heterogeneities in the -value growth. That is, the extent of transition-state movement (Delta beta(++)) is proportional to the free-energy gradient of the mutational perturbation across the top of the activation barrier, '(beta(++)) proportional, variant Delta logK(D-N). The analysis is applied to the two-state protein L23 where the site-specific free-energy gradients are used to identify the interactions that show the highest degree of consolidation after crossing the barrier top. These interactions are distributed as a shell around the high- initiation point and denote the side-chain contacts that add criticality to the folding nucleus.

Simulation, experiment, and evolution: understanding nucleation in protein S6 folding

In this study, we explore nucleation and the transition state ensemble of the ribosomal protein S6 using a Monte Carlo (MC) Go model in conjunction with restraints from experiment. The results are analyzed in the context of extensive experimental and evolutionary data. The roles of individual residues in the folding nucleus are identified, and the order of events in the S6 folding mechanism is explored in detail. Interpretation of our results agrees with, and extends the utility of, experiments that shift phi-values by modulating denaturant concentration and presents strong evidence for the realism of the mechanistic details in our MC Go model and the structural interpretation of experimental phi-values. We also observe plasticity in the contacts of the hydrophobic core that support the specific nucleus. For S6, which binds to RNA and protein after folding, this plasticity may result from the conformational flexibility required to achieve biological function. These results present a theoretical and conceptual picture that is relevant in understanding the mechanism of nucleation in protein folding.

High concentrations of viscogens decrease the protein folding rate constant by prematurely collapsing the coil

In several studies, viscogenic osmolytes have been suggested to decrease the folding rate constant of polypeptides by slowing their motion through the solvent. Here, we show that osmolytes may slow protein folding by prematurely collapsing the coil. At low or moderate concentrations of osmolytes (<30%), folding of the two-state protein CI2 becomes faster with increasing osmolyte concentrations, suggesting that the kinetics are governed by protein stability. However, at higher concentrations of osmolyte, the coil collapses in the dead-time of the refolding experiment, causing a dramatic drop in the folding rate. The collapsed state is non-native and appears to be different for different osmolytes.

Common denominator of Cu/Zn superoxide dismutase mutants associated with amyotrophic lateral sclerosis: decreased stability of the apo state

More than 100 point mutations of the superoxide scavenger Cu/Zn superoxide dismutase (SOD; EC ) have been associated with the neurodegenerative disease amyotrophic lateral sclerosis (ALS). However, these mutations are scattered throughout the protein and provide no clear functional or structural clues to the underlying disease mechanism. Therefore, we undertook to look for folding-related defects by comparing the unfolding behavior of five ALS-associated mutants with distinct structural characteristics: A4V at the interface between the N and C termini, C6F in the hydrophobic core, D90A at the protein surface, and G93A and G93C, which decrease backbone flexibility. With the exception of the disruptive replacements A4V and C6F, the mutations only marginally affect the stability of the native protein, yet all mutants share a pronounced destabilization of the metal-free apo state: the higher the stability loss, the lower the mean survival time for ALS patients carrying the mutation. Thus organism-level pathology may be directly related to the properties of the immature state of a protein rather than to those of the native species.