Protein engineering as a strategy to avoid formation of amyloid fibrils

The activation domain of human procarboxypeptidase A2 (ADA2h) aggregates following thermal or chemical denaturation at acidic pH. The aggregated material contains well-defined ordered structures with all the characteristics of the fibrils associated with amyloidotic diseases. Variants of ADA2h containing a series of mutations designed to increase the local stability of each of the two helical regions of the protein have been found to have a substantially reduced propensity to form fibrils. This arises from a reduced tendency of the denatured species to aggregate rather than from a change in the overall stability of the native state. The reduction in aggregation propensity may result from an increase in the stability of local relative to longer range interactions within the polypeptide chain. These findings show that the intrinsic ability of a protein to form amyloid can be altered substantially by protein engineering methods without perturbing significantly its overall stability or activity. This suggests new strategies for combating diseases associated with the formation of aggregated proteins and for the design of novel protein or peptide therapeutics.     

A single mutation induces amyloid aggregation in the alpha-spectrin SH3 domain: analysis of the early stages of fibril formation

The Src-homology region 3 domain of chicken alpha-spectrin (Spc-SH3) is a small two-state folding protein, which has never been described to form amyloid fibrils under any condition investigated so far. We show here that the mutation of asparagine 47 to alanine at the distal loop, which destabilises similarly the native and folding transition states of the domain, induces the formation of amyloid fibrils under mild acid conditions. Amyloid aggregation of the mutant is enhanced by the increase in temperature, protein concentration and NaCl concentration. The early stages of amyloid formation have been monitored as a function of time and temperature using a variety of biophysical methods. Differential scanning calorimetry experiments under conditions of amyloid formation have allowed the identification of different thermal transitions corresponding to conformational and aggregation processes as well as to the high-temperature disaggregation and unfolding of the amyloid fibrils. Aggregation is preceded by a rapid conformational change in the monomeric domain involving about 40% of the global unfolding enthalpy, considerable change in secondary structure, large loss of tertiary structure and exposure of hydrophobic patches to the solvent. The conformational change is followed by formation of a majority of oligomeric species with apparent hydrodynamic radius between 2.5 nm and 10nm, depending on temperature, together with the appearance and progressive growth of protofibrillar aggregates. After these early aggregation stages, long and curved fibrils of up to several micrometers start to develop by elongation of the protofibrils. The calorimetric data indicate that the specific enthalpy of fibril disaggregation and unfolding is relatively low, suggesting a low density of interactions within the fibril structure as compared to the native protein and a main entropy contribution to the stability of the amyloid fibrils.     

The design of a hyperstable mutant of the Abp1p SH3 domain by sequence alignment analysis

We have characterized the thermodynamic stability of the SH3 domain from the Saccharomyces cerevisiae Abp1p protein and found it to be relatively low compared to most other SH3 domains, with a Tm of 60 degrees C and a deltaGu of 3.08 kcal/mol. Analysis of a large alignment of SH3 domains led to the identification of atypical residues at eight positions in the wild-type Abp1p SH3 domain sequence that were subsequently replaced by the residue seen most frequently at that position in the alignment. Three of the eight mutants constructed in this way displayed increases in Tm ranging from 8 to 15 degrees C with concomitant increases in deltaGu of up to 1.4 kcal/mol. The effects of these substitutions on folding thermodynamics and kinetics were entirely additive, and a mutant containing all three was dramatically stabilized with a Tm greater than 90 degrees C and a deltaGu more than double that of the wild-type domain. The folding rate of this hyperstable mutant was 10-fold faster than wild-type, while its unfolding rate was fivefold slower. All of the stabilized mutants were still able to bind a target peptide with wild-type affinity. We have analyzed the stabilizing amino acid substitutions isolated in this study and several other similar sequence alignment based studies. In approximately 25% of cases, increased stability can be explained by enhanced propensity of the substituted residue for the local backbone conformation at the mutagenized site.     

Protein aggregation and amyloid fibril formation by an SH3 domain probed by limited proteolysis

The SH3 domains are small protein modules of 60-85 amino acid residues that are found in many proteins involved in intracellular signal transduction. The SH3 domain of the p85alpha subunit of bovine phosphatidylinositol 3'-kinase (PI3-SH3) under acidic solution adopts a compact denatured state from which amyloid fibrils are readily formed. This aggregation process has been found to be modulated substantially by solution conditions. Here, we have analyzed the conformational features of the native and acid denatured states of PI3-SH3 by limited proteolysis experiments using proteinase K and pepsin, respectively. Moreover, we have analyzed the propensity of PI3-SH3 to be hydrolyzed by pepsin at different stages in the process of aggregation and amyloid formation at pH 1.2 and 2.0 and compared the sites of proteolysis under these conditions with the conformational features of both native and aggregated PI3-SH3. The results demonstrate that the denatured state of PI3-SH3 formed at low pH is relatively resistant to proteolysis, indicating that it is partially folded. The long loop connecting beta-strands b and c in the native protein is the region in this structure most susceptible to proteolysis. Remarkably, aggregates of PI3-SH3 that are formed initially from this denatured state in acid solution display enhanced susceptibility to proteolysis of the long loop, suggesting that the protein becomes more unfolded in the early stages of aggregation. By contrast, the more defined amyloid fibrils that are formed over longer periods of time are completely resistant to proteolysis. We suggest that the protein aggregates formed initially are relatively dynamic species that are able readily to reorganize their interactions to enable formation of very well ordered fibrillar structures. In addition, the disordered and non-native character of the polypeptide chains in the early aggregates could be important in determining the high cytotoxicity that has been revealed in previous studies of these species.     

Temperature dependence of the free energy landscape of the src-SH3 protein domain

The temperature dependence of the free energy landscape of the src-SH3 protein domain is investigated through fully atomic simulations in explicit solvent. Simulations are performed above and below the folding transition temperature, enabling an analysis of both protein folding and unfolding. The transition state for folding and unfolding, identified from the free energy surfaces, is found to be very similar, with structure in the central hydrophobic sheet and little structure throughout the rest of the protein. This is a result of a polarized folding (unfolding) mechanism involving early formation (late loss) of the central hydrophobic sheet at the transition state. Unfolding simulations map qualitatively well onto low-temperature free energy surfaces but appear, however, to miss important features observed in folding simulations. In particular, details of the folding mechanism involving the opening and closing of the hydrophobic core are not captured by unfolding simulations performed under strongly denaturing conditions. In addition, free energy surfaces at high temperatures do not display a desolvation barrier found at lower temperatures, involving the expulsion of water molecules from the hydrophobic core.     

Nucleation and the transition state of the SH3 domain

We present a verified computational model of the SH3 domain transition state (TS) ensemble. This model was built for three separate SH3 domains using experimental phi-values as structural constraints in all-atom protein folding simulations. While averaging over all conformations incorrectly considers non-TS conformations as transition states, quantifying structures as pre-TS, TS, and post-TS by measurement of their transmission coefficient ("probability to fold", or p(fold)) allows for rigorous conclusions regarding the structure of the folding nucleus and a full mechanistic analysis of the folding process. Through analysis of the TS, we observe a highly polarized nucleus in which many residues are solvent-exposed. Mechanistic analysis suggests the hydrophobic core forms largely after an early nucleation step. SH3 presents an ideal system for studying the nucleation-condensation mechanism and highlights the synergistic relationship between experiment and simulation in the study of protein folding.     

Thermodynamics and folding kinetics analysis of the SH3 domain form discrete molecular dynamics

We perform a detailed analysis of the thermodynamics and folding kinetics of the SH3 domain fold with discrete molecular dynamic simulations. We propose a protein model that reproduces some of the experimentally observed thermodynamic and folding kinetic properties of proteins. Specifically, we use our model to study the transition state ensemble of the SH3 fold family of proteins, a set of unstable conformations that fold to the protein native state with probability 1/2. We analyze the participation of each secondary structure element formed at the transition state ensemble. We also identify the folding nucleus of the SH3 fold and test extensively its importance for folding kinetics. We predict that a set of amino acid contacts between the RT-loop and the distal hairpin are the critical folding nucleus of the SH3 fold and propose a hypothesis that explains this result.     

FTIR reveals structural differences between native beta-sheet proteins and amyloid fibrils

The presence of beta-sheets in the core of amyloid fibrils raised questions as to whether or not beta-sheet-containing proteins, such as transthyretin, are predisposed to form such fibrils. However, we show here that the molecular structure of amyloid fibrils differs more generally from the beta-sheets in native proteins. This difference is evident from the amide I region of the infrared spectrum and relates to the distribution of the phi/psi dihedral angles within the Ramachandran plot, the average number of strands per sheet, and possibly, the beta-sheet twist. These data imply that amyloid fibril formation from native beta-sheet proteins can involve a substantial structural reorganization.     

Elongated oligomers assemble into mammalian PrP amyloid fibrils

In prion diseases, the mammalian prion protein PrP is converted from a monomeric, mainly alpha-helical state into beta-rich amyloid fibrils. To examine the structure of the misfolded state, amyloid fibrils were grown from a beta form of recombinant mouse PrP (residues 91-231). The beta-PrP precursors assembled slowly into amyloid fibrils with an overall helical twist. The fibrils exhibit immunological reactivity similar to that of ex vivo PrP Sc. Using electron microscopy and image processing, we obtained three-dimensional density maps of two forms of PrP fibrils with slightly different twists. They reveal two intertwined protofilaments with a subunit repeat of approximately 60 A. The repeating unit along each protofilament can be accounted for by elongated oligomers of PrP, suggesting a hierarchical assembly mechanism for the fibrils. The structure reveals flexible crossbridges between the two protofilaments, and subunit contacts along the protofilaments that are likely to reflect specific features of the PrP sequence, in addition to the generic, cross-beta amyloid fold.     

Src homology domains of v-Src stabilize an active conformation of the tyrosine kinase catalytic domain

To examine the interactions between Src homology,domains and the tyrosine kinase catalytic domain of v-Src, various combinations of domains have been expressed in bacteria as fusion proteins. Constructs containing the isolated catalytic domain, SH2 + catalytic domain, and SH3 + SH2 + catalytic domains were active in autophosphorylation assays. For the catalytic domain of v-Src, but not for v-Abl, addition of exogenous Src SH3-SH2 domains stimulated the autophosphorylation activity. In contrast to results for autophosphorylation, constructs containing Src homology domains were more active towards a synthetic peptide substrate than the isolated catalytic domain. The ability of the SH2 and SH3 domains of v-Src to stabilize an active enzyme conformation was also confirmed by refolding after denaturation in guanidinium hydrochloride. Collectively the data suggest that, in addition to their roles in intermolecular protein-protein interactions, the Src homology regions of v-Src exert a positive influence on tyrosine kinase function, potentially by maintaining an active conformation of the catalytic domain.     

Probing the mechanism of amyloidogenesis through a tandem repeat of the PI3-SH3 domain suggests a generic model for protein aggregation and fibril formation

Aggregation of the SH3 domain of the PI3 kinase, both as a single domain and as a tandem repeat in which the C terminus of one domain is linked to the N terminus of another by a flexible linker of ten glycine/serine residues, has been studied under a range of conditions in order to investigate the mechanism of protein aggregation and amyloid formation. The tandem repeat was found to form amyloid fibrils much more readily than the single domain under the acidic conditions used here, and the fibrils themselves have higher morphological homogeneity. The folding-unfolding transition of the PI3-SH3 domain shows two-state behaviour and is pH dependent; at pH 3.6, which is near the pH mid-point for folding and only slightly below the isoelectric point of the protein, both the single domain and the tandem repeat spontaneously form broad distributions of soluble oligomers without requirement for nucleation. Under prolonged incubation under these conditions, the oligomers convert into thin, curly fibrils that interact with thioflavin-T, suggesting that they contain an organised beta-sheet structure. Under more acidic conditions (pH 2.0) where the proteins are fully denatured and carry a positive net charge, long, straight fibrils are formed in a process having a pronounced lag phase. The latter was found to be reduced dramatically by the addition of oligomers exceeding a critical size of approximately 20 molecules. The results suggest that the process of aggregation of these SH3 domains can take place by a variety of mechanisms, ranging from downhill formation of relatively amorphous species to nucleated formation of highly organised structures, the relative importance of which varies greatly with solution conditions. Comparison with the behaviour of other amyloidogenic systems suggests that the general mechanistic features outlined here are likely to be common to at least a wide variety of peptides and proteins.     

Different propensity to form amyloid fibrils by two homologous proteins-Human stefins A and B: searching for an explanation

By using ThT fluorescence, X-ray diffraction, and atomic force microscopy (AFM), it has been shown that human stefins A and B (subfamily A of cystatins) form amyloid fibrils. Both protein fibrils show the 4.7 A and 10 A reflections characteristic for cross beta-structure. Similar height of approximately 3 nm and longitudinal repeat of 25-27 nm were observed by AFM for both protein fibrils. Fibrils with a double height of 5.6 nm were only observed with stefin A. The fibril's width for stefin A fibrils, as observed by transmission electron microscopy (TEM), was in the same range as previously reported for stefin B (Zerovnik et al., Biochem Biophys Acta 2002;1594:1-5). The conditions needed to undergo fibrillation differ, though. The amyloid fibrils start to form at pH 5 for stefin B, whereas in stefin A, preheated sample has to be acidified to pH < 2.5. In both cases, adding TFE, seeding, and alignment in a strong magnetic field accelerate the fibril growth. Visual analysis of the three-dimensional structures of monomers and domain-swapped dimers suggests that major differences in stability of both homologues stem from arrangement of specific salt bridges, which fix alpha-helix (and the alpha-loop) to beta-sheet in stefin A monomeric and dimeric forms.     

Protein folding kinetics beyond the phi value: using multiple amino acid substitutions to investigate the structure of the SH3 domain folding transition state

The SH3 domain folding transition state structure contains two well-ordered turn regions, known as the diverging turn and the distal loop. In the Src SH3 domain transition state, these regions are stabilized by a hydrogen bond between Glu30 in the diverging turn and Ser47 in the distal loop. We have examined the effects on folding kinetics of amino acid substitutions at the homologous positions (Glu24 and Ser41) in the Fyn SH3 domain. In contrast to most other folding kinetics studies which have focused primarily on non-disruptive substitutions with Ala or Gly, here we have examined the effects of substitutions with diverse amino acid residues. Using this approach, we demonstrate that the transition state structure is generally tolerant to amino acid substitutions. We also uncover a unique role for Ser at position 41 in facilitating folding of the distal loop, which can only be replicated by Asp at the same position. Both these residues appear to accelerate folding through the formation of short-range side-chain to backbone hydrogen bonds. The folding of the diverging turn region is shown to be driven primarily by local interactions. The diverging turn and distal loop regions are found to interact in the transition state structure, but only in the context of particular mutant backgrounds. This work demonstrates that studying the effects of a variety of amino acid substitutions on protein folding kinetics can provide unique insights into folding mechanisms which cannot be obtained by standard Phi value analysis.     

The role of protein sequence and amino acid composition in amyloid formation: scrambling and backward reading of IAPP amyloid fibrils

The specific functional structure of natural proteins is determined by the way in which amino acids are sequentially connected in the polypeptide. The tight sequence/structure relationship governing protein folding does not seem to apply to amyloid fibril formation because many proteins without any sequence relationship have been shown to assemble into very similar β-sheet-enriched structures. Here, we have characterized the aggregation kinetics, seeding ability, morphology, conformation, stability, and toxicity of amyloid fibrils formed by a 20-residue domain of the islet amyloid polypeptide (IAPP), as well as of a backward and scrambled version of this peptide. The three IAPP peptides readily aggregate into ordered, β-sheet-enriched, amyloid-like fibrils. However, the mechanism of formation and the structural and functional properties of aggregates formed from these three peptides are different in such a way that they do not cross-seed each other despite sharing a common amino acid composition. The results confirm that, as for globular proteins, highly specific polypeptide sequential traits govern the assembly pathway, final fine structure, and cytotoxic properties of amyloid conformations.     

A thermodynamic and kinetic analysis of the folding pathway of an SH3 domain entropically stabilised by a redesigned hydrophobic core

The folding thermodynamics and kinetics of the alpha-spectrin SH3 domain with a redesigned hydrophobic core have been studied. The introduction of five replacements, A11V, V23L, M25V, V44I and V58L, resulted in an increase of 16% in the overall volume of the side-chains forming the hydrophobic core but caused no remarkable changes to the positions of the backbone atoms. Judging by the scanning calorimetry data, the increased stability of the folded structure of the new SH3-variant is caused by entropic factors, since the changes in heat capacity and enthalpy upon the unfolding of the wild-type and mutant proteins were identical at 298 K. It appears that the design process resulted in an increase in burying both the hydrophobic and hydrophilic surfaces, which resulted in a compensatory effect upon the changes in heat capacity and enthalpy. Kinetic analysis shows that both the folding and unfolding rate constants are higher for the new variant, suggesting that its transition state becomes more stable compared to the folded and unfolded states. The phi(double dagger-U) values found for a number of side-chains are slightly lower than those of the wild-type protein, indicating that although the transition state ensemble (TSE) did not change overall, it has moved towards a more denatured conformation, in accordance with Hammond's postulate. Thus, the acceleration of the folding-unfolding reactions is caused mainly by an improvement in the specific and/or non-specific hydrophobic interactions within the TSE rather than by changes in the contact order. Experimental evidence showing that the TSE changes globally according to its hydrophobic content suggests that hydrophobicity may modulate the kinetic behaviour and also the folding pathway of a protein.     

Crystallographic structure of the SH3 domain of the human c-Yes tyrosine kinase: loop flexibility and amyloid aggregation

SH3 domains from the Src family of tyrosine kinases represent an interesting example of the delicate balance between promiscuity and specificity characteristic of proline-rich ligand recognition by SH3 domains. The development of inhibitors of therapeutic potential requires a good understanding of the molecular determinants of binding affinity and specificity and relies on the availability of high quality structural information. Here, we present the first high-resolution crystal structure of the SH3 domain of the c-Yes oncogen. Comparison with other SH3 domains from the Src family revealed significant deviations in the loop regions. In particular, the n-Src loop, highly flexible and partially disordered, is stabilized in an unusual conformation by the establishment of several intramolecular hydrogen bonds. Additionally, we present here the first report of amyloid aggregation by an SH3 domain from the Src family.     

The relationship between conservation, thermodynamic stability, and function in the SH3 domain hydrophobic core

To investigate the relationships between sequence conservation, protein stability, and protein function, we have measured the thermodynamic stability, folding kinetics, and in vitro peptide-binding activity of a large number of single-site substitutions in the hydrophobic core of the Fyn SH3 domain. Comparison of these data to that derived from an analysis of a large alignment of SH3 domain sequences revealed a very good correlation between the distinct pattern of conservation observed at each core position and the thermodynamic stability of mutants. Conservation was also found to correlate well with the unfolding rates of mutants, but not to the folding rates, suggesting that evolution selects more strongly for optimal native state packing interactions than for maximal folding rates. Structural analysis suggests that residue-residue core packing interactions are very similar in all SH3 domains, which provides an explanation for the correlation between conservation and mutant stability effects studied in a single SH3 domain. We also demonstrate a correlation between stability and the in vivo activity of mutants, and between conservation and activity. However, the relationship between conservation and activity was very strong only for the three most conserved hydrophobic core positions. The weaker correlation between activity and conservation seen at the other seven core positions indicates that maintenance of protein stability is the dominant selective pressure at these positions. In general, the pattern of conservation at hydrophobic core positions appears to arise from conserved packing constraints, and can be effectively utilized to predict the destabilizing effects of amino acid substitutions.     

Native state energetics of the Src SH2 domain: evidence for a partially structured state in the denatured ensemble

We have defined the free-energy profile of the Src SH2 domain using a variety of biophysical techniques. Equilibrium and kinetic experiments monitored by tryptophan fluorescence show that Src SH2 is quite stable and folds rapidly by a two-state mechanism, without populating any intermediates. Native state hydrogen-deuterium exchange confirms this two-state behavior; we detect no cooperative partially unfolded forms in equilibrium with the native conformation under any conditions. Interestingly, the apparent stability of the protein from hydrogen exchange is 2 kcal/mol greater than the stability determined by both equilibrium and kinetic studies followed by fluorescence. Native-state proteolysis demonstrates that this "super protection" does not result from a deviation from the linear extrapolation model used to fit the fluorescence data. Instead, it likely arises from a notable compaction in the unfolded state under native conditions, resulting in an ensemble of conformations with substantial solvent exposure of side chains and flexible regions sensitive to proteolysis, but backbone amides that exchange with solvent approximately 30-fold slower than would be expected for a random coil. The apparently simple behavior of Src SH2 in traditional unfolding studies masks the significant complexity present in the denatured-state ensemble.