Amyloid Oligomers: A Joint Experimental/Computational Perspective on Alzheimer's Disease, Parkinson's Disease, Type II Diabetes, and Amyotrophic Lateral Sclerosis

Protein misfolding and aggregation is observed in many amyloidogenic diseases affecting either the central nervous system or a variety of peripheral tissues. Structural and dynamic characterization of all species along the pathways from monomers to fibrils is challenging by experimental and computational means because they involve intrinsically disordered proteins in most diseases. Yet understanding how amyloid species become toxic is the challenge in developing a treatment for these diseases. Here we review what computer, in vitro, in vivo, and pharmacological experiments tell us about the accumulation and deposition of the oligomers of the (Aβ, tau), α-synuclein, IAPP, and superoxide dismutase 1 proteins, which have been the mainstream concept underlying Alzheimer's disease (AD), Parkinson's disease (PD), type II diabetes (T2D), and amyotrophic lateral sclerosis (ALS) research, respectively, for many years.

N- and C-capping preferences for all 20 amino acids in alpha-helical peptides

We have determined the N- and C-capping preferences of all 20 amino acids by substituting residue X in the peptides NH2-XAKAAAAKAAAAKAAGY-CONH2 and in Ac-YGAAKAAAAKAAAAKAX-CO2H. Helix contents were measured by CD spectroscopy to obtain rank orders of capping preferences. The data were further analyzed by our modified Lifson-Roig helix-coil theory, which includes capping parameters (n and c), to find free energies of capping (-RT ln n and -RT ln c), relative to Ala. Results were obtained for charged and uncharged termini and for different charged states of titratable side chains. N-cap preferences varied from Asn (best) to Gln (worst). We find, as expected, that amino acids that can accept hydrogen bonds from otherwise free backbone NH groups, such as Asn, Asp, Ser, Thr, and Cys generally have the highest N-cap preference. Gly and acetyl group are favored, as are negative charges in side chains and at the N-terminus. Our N-cap preference scale agrees well with preferences in proteins. In contrast, we find little variation when changing the identity of the C-cap residue. We find no preference for Gly at the C-cap in contrast to the situation in proteins. Both N-cap and C-cap results for Tyr and Trp are inaccurate because their aromatic groups affect the CD spectrum. The data presented here are of value in rationalizing mutations at capping sites in proteins and in predicting the helix contents of peptides.

Side-chain conformational entropy in protein folding

An important, but often neglected, contribution to the thermodynamics of protein folding is the loss of entropy that results from restricting the number of accessible side-chain conformers in the native structure. Conformational entropy changes can be found by comparing the number of accessible rotamers in the unfolded and folded states, or by estimating fusion entropies. Comparison of several sets of results using different techniques shows that the mean conformational free energy change (T delta S) is 1 kcal.mol-1 per side chain or 0.5 kcal.mol-1 per bond. Changes in vibrational entropy appear to be negligible compared to the entropy change resulting from the loss of accessible rotamers. Side-chain entropies can help rationalize alpha-helix propensities, predict protein/inhibitor complex structures, and account for the distribution of side chains on the protein surface or interior.

Inhibition of toxicity in the beta-amyloid peptide fragment beta -(25-35) using N-methylated derivatives: a general strategy to prevent amyloid formation

beta-(25-35) is a synthetic derivative of beta-amyloid, the peptide that is believed to cause Alzheimer's disease. As it is highly toxic and forms fibrillar aggregates typical of beta-amyloid, it is suitable as a model for testing inhibitors of aggregation and toxicity. We demonstrate that N-methylated derivatives of beta-(25-35), which in isolation are soluble and non-toxic, can prevent the aggregation and inhibit the resulting toxicity of the wild type peptide. N-Methylation can block hydrogen bonding on the outer edge of the assembling amyloid. The peptides are assayed by Congo red and thioflavin T binding, electron microscopy, and a 3-(4, 5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) toxicity assay on PC12 cells. One peptide (Gly(25) N-methylated) has properties similar to the wild type, whereas five have varying effects on prefolded fibrils and fibril assembly. In particular, beta-(25-35) with Gly(33) N-methylated is able to completely prevent fibril assembly and to reduce the toxicity of prefolded amyloid. With Leu(34) N-methylated, the fibril morphology is altered and the toxicity reduced. We suggest that the use of N-methylated derivatives of amyloidogenic peptides and proteins could provide a general solution to the problem of amyloid deposition and toxicity.

Is polyproline II helix the killer conformation? A Raman optical activity study of the amyloidogenic prefibrillar intermediate of human lysozyme

The amyloidogenic prefibrillar partially denatured intermediate of human lysozyme, prepared by heating the native protein to 57 degrees C at pH 2.0, was studied using Raman optical activity (ROA). A positive band in the room temperature ROA spectrum of the native protein at approximately 1345 cm(-1), assigned to a hydrated form of alpha-helix, is not present in that of the prefibrillar intermediate, where a new strong positive band at approximately 1318 cm(-1) appears instead that is assigned to the poly(l-proline) II (PPII)-helical conformation. A sharp negative band at approximately 1241 cm(-1) in the native protein, assigned to beta-strand, shows little change in the ROA spectrum of the prefibrillar intermediate. The disappearance of a positive ROA band at approximately 1551 cm(-1) assigned to vibrations of tryptophan side-chains indicates that major conformational changes have occurred among the five tryptophan residues present in human lysozyme, four of which are located in the alpha-domain. The various ROA data suggest that a substantial loss of tertiary structure has occurred in the prefibrillar intermediate and that this is located more in the alpha-domain than in the beta-domain. There is no evidence for any increase in beta-structure. The ROA spectrum of hen lysozyme, which does not form amyloid fibrils so readily, remains much more native-like on heating to 57 degrees C at pH 2.0. The thermal behaviour of the alanine-rich alpha-helical peptide AK21 in aqueous solution was found to be similar to that of human lysozyme. Hydrated alpha-helix therefore appears to readily undergo a conformational change to PPII structure on heating, which may be a key step in the conversion of alpha-helix into beta-sheet in the formation of amyloid fibrils in human lysozyme. Since it is extended, flexible, lacks intrachain hydrogen bonds and is fully hydrated in aqueous solution, PPII helix has the appropriate characteristics to be implicated as a critical conformational element in many conformational diseases. Disorder of the PPII type may be a sine qua non for the formation of regular fibrils; whereas the more dynamic disorder of the random coil may lead only to amorphous aggregates.

Helix capping propensities in peptides parallel those in proteins

Helix content of peptides with various uncharged nonaromatic amino acids at either the N-terminal or C-terminal position has been determined. The choice of N-terminal amino acid has a major effect on helix stability: asparagine is the best, glycine is very good, and glutamine is the worst helix-stabilizing amino acid at this position. The rank order of helix stabilization parallels the frequencies of these amino acids at the N-terminal boundary (N-cap) position of helices in proteins found by Richardson and Richardson [Richardson, J. S. & Richardson, D. C. (1988) Science 240, 1648-1652], and the N-terminal amino acid in a peptide composed of helix-forming amino acids may be considered as the N-cap residue. The choice of C-terminal amino acid has only a minor effect on helix stability. N-capping interactions may be responsible for the asymmetric distribution of helix content within a given peptide found by various workers. An acetyl group on the N-terminal alpha-amino function cancels the N-cap effect and the acetyl group is equivalent to N-terminal asparagine in an unacetylated peptide. Our results demonstrate a close relationship between the mechanisms of alpha-helix formation in peptides and in proteins.

Why Is Research on Amyloid-β Failing to Give New Drugs for Alzheimer's Disease?

The two hallmarks of Alzheimer's disease (AD) are the presence of neurofibrillary tangles (NFT) made of aggregates of the hyperphosphorylated tau protein and of amyloid plaques composed of amyloid-β (Aβ) peptides, primarily Aβ1-40 and Aβ1-42. Targeting the production, aggregation, and toxicity of Aβ with small molecule drugs or antibodies is an active area of AD research due to the general acceptance of the amyloid cascade hypothesis, but thus far all drugs targeting Aβ have failed. From a review of the recent literature and our own experience based on in vitro, in silico, and in vivo studies, we present some reasons to explain this repetitive failure.

N-Methylated peptide inhibitors of beta-amyloid aggregation and toxicity. Optimization of the inhibitor structure

The key pathogenic event in the onset of Alzheimer's disease (AD) is believed to be the aggregation of the beta-amyloid (Abeta) peptide into toxic oligomers. Molecules that interfere with this process may therefore act as therapeutic agents for the treatment of AD. N-Methylated peptides (meptides) are a general class of peptide aggregation inhibitors that act by binding to one face of the aggregating peptide but are unable to hydrogen bond on the other face, because of the N-methyl group replacing a backbone NH group. Here, we optimize the structure of meptide inhibitors of Abeta aggregation, starting with the KLVFF sequence that is known to bind to Abeta. We varied the meptide length, N-methylation sites, acetylation, and amidation of the N and C termini, side-chain identity, and chirality, via five compound libraries. Inhibitor activity was tested by thioflavin T binding, affinity chromatography, electron microscopy, and an 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide toxicity assay. We found that inhibitors should have all d chirality, have a free N terminus but an amidated C terminus, and have large, branched hydrophobic side chains at positions 1-4, while the side chain at position 5 was less important. A single N-methyl group was necessary and sufficient. The most active compound, d-[(chGly)-(Tyr)-(chGly)-(chGly)-(mLeu)]-NH(2), was more active than all previously reported peptide inhibitors. Its related non-N-methylated analogues were insoluble and toxic.

Distinguishing enzyme structures from non-enzymes without alignments

The ability to predict protein function from structure is becoming increasingly important as the number of structures resolved is growing more rapidly than our capacity to study function. Current methods for predicting protein function are mostly reliant on identifying a similar protein of known function. For proteins that are highly dissimilar or are only similar to proteins also lacking functional annotations, these methods fail. Here, we show that protein function can be predicted as enzymatic or not without resorting to alignments. We describe 1178 high-resolution proteins in a structurally non-redundant subset of the Protein Data Bank using simple features such as secondary-structure content, amino acid propensities, surface properties and ligands. The subset is split into two functional groupings, enzymes and non-enzymes. We use the support vector machine-learning algorithm to develop models that are capable of assigning the protein class. Validation of the method shows that the function can be predicted to an accuracy of 77% using 52 features to describe each protein. An adaptive search of possible subsets of features produces a simplified model based on 36 features that predicts at an accuracy of 80%. We compare the method to sequence-based methods that also avoid calculating alignments and predict a recently released set of unrelated proteins. The most useful features for distinguishing enzymes from non-enzymes are secondary-structure content, amino acid frequencies, number of disulphide bonds and size of the largest cleft. This method is applicable to any structure as it does not require the identification of sequence or structural similarity to a protein of known function.

Determination of free energies of N-capping in alpha-helices by modification of the Lifson-Roig helix-coil therapy to include N- and C-capping

We have previously shown that varying the N-terminal amino acid in alpha-helical peptides can cause large variations in helix content (Chakrabartty et al., 1993a). The Lifson-Roig theory for the helix-coil transition predicts, however, that substitutions at the N-terminus in an unacetylated peptide should have no effect on alpha-helix stability. We have therefore modified the theory to include these N-capping effects by assigning a statistical weight (the "n-value") to the amino acid immediately preceding a stretch of helical residues. The n-value measures the N-capping propensity of an amino acid, and like the helix propensity (w-value), it is independent of neighboring residues or positions in sequence. The new theory was used, with the experimental data for these substitutions, to calculate n-values and, hence, free energies for N-capping for the amino acids Gln, Ala, Val, Met, Pro, Ile, Leu, Thr, Gly, Ser, and Asn as well as for the acetyl group, which is commonly used to cap peptides. The free energies vary by approximately 1 kcal mol-1 from Gln (worst) to Asn (best), and the acetyl group is nearly as effective as Asn. N-Capping free energies were also found for Leu, Thr, Gly, Ser, and Asn when the N-terminus is charged at pH 5. The unfavorable effect of protonation of the N-terminus in an alpha-helix was found to be approximately 0.5 kcal mol-1. Our results agree well with a survey of N-capping preferences from protein crystal structures and are compared to results from site-directed mutagenesis of N-caps in proteins.

Is the hydrophobic effect stabilizing or destabilizing in proteins? The contribution of disulphide bonds to protein stability

It has been recently concluded that the hydrophobic effect, hitherto regarded as a major driving force in the folding of proteins, destabilizes the folded state relative to the unfolded state. We summarize the properties of the hydrophobic effect obtained from solvent transfer experiments and show that the recent conclusion is an artifact of crosslinking in the unfolded state, caused by disulphide bonds, metals or cofactors. We show that, for the proteins in the data set, crosslinks surprisingly destabilize folded structures entropically, but stabilize them enthalpically to a greater extent. We also calculate non-polar surface areas of these unfolded proteins. These surface areas are decreased by crosslinks. The unfolded state of proteins lacking constraints, such as myoglobin, is well approximated by a mixture of residues containing alpha-helical and beta-sheet dihedral angles. Surface areas of unfolded proteins cannot be obtained by summing the surface areas of individual residues, since this ignores any unavoidable side-chain-side-chain interactions.

Effect of phosphorylation on alpha-helix stability as a function of position

We have investigated the effect of placing phosphoserine at the N-cap, N1, N2, N3, and interior position in alanine-based alpha-helical peptides. Helix contents of each peptide were measured by CD spectroscopy and titrations performed to determine pK(a) values. Data were analyzed with modified Lifson-Roig theory to determine helix-coil parameters (n, n(1), n(2), n(3), and w) and free energy changes for phosphoserine at each helical position. Results are given for a -1 and -2 phosphoserine charge state. Results show that phosphoserine stabilizes at the N-terminal positions by as much as 2.3 kcal.mol(-1), while destabilizes in the helix interior by 1.2 kcal.mol(-1), relative to serine. The rank order of free energies relative to serine at each position is N2 > N3 > N1 > N-cap > interior. Moreover, -2 phosphoserine is the most preferred residue known at each of these N-terminal positions. Experimental pK(a) values for the -1 to -2 phosphoserine transition are in the order N2 < N-cap < N1 < N3 < interior. This order agrees well with electrostatics calculations carried out with phosphoserine at the N-terminal positions and interior positions. Combining these with calculations at the C3, C2, C1, and C-cap positions gives results for phosphoserine along the length of the helix. We see a transition from phosphoserine stabilization at the N-terminus to destabilization at the C-terminus and can explain this in terms of the balance of protein solvation, favorable interactions, and dehydration. These results give insight into the phosphorylatable control of biological systems through positive or negative changes in stability.

Design strategies for anti-amyloid agents

Numerous diseases have been linked to a common pathogenic process called amyloidosis, whereby proteins or peptides clump together in the brain or body to form toxic soluble oligomers and/or insoluble fibres. An attractive strategy to develop therapies for these diseases is therefore to inhibit or reverse protein/peptide aggregation. A diverse range of small organic ligands have been found to act as aggregation inhibitors. Alternatively, the wild-type peptide can be derivatised so that it still binds to the amyloid target, but prevents further aggregation. This can be achieved by adding a bulky group or charged amino acid to either end of the peptide, or by incorporating proline residues or N-methylated amide groups.

The alpha-helix folds on the millisecond time scale

It has long been believed that nucleation of the alpha-helix is a very fast reaction, occurring in around 10(-7) s. We show here that helix nucleation, in fact, takes place on the millisecond time scale. The rate of alpha-helix nucleation in two polyalanine-based peptides and in lysine and glutamic acid homopolymers was measured directly by stopped-flow deep UV CD with synchrotron radiation as the light source. Synchrotron radiation CD gives far superior signal to noise than a conventional instrument. The 16-aa AK peptide folds with first-order kinetics and a rate constant of 15 s-1 at 0 degrees C. The rate-determining step is presumably the initiation of a new helix, which occurs at least 10(5) times slower than expected. Helix folding occurs in at least two steps on the millisecond time scale for the longer peptides, with a transient overshoot of helix content significantly greater than at equilibrium, similar to that seen in the folding of several proteins. We suggest that the overshoot is caused by the formation of a single long helix followed by its breakage into the two or more helices present at equilibrium.

Structures of N-termini of helices in proteins

We have surveyed 393 N-termini of alpha-helices and 156 N-termini of 3(10)-helices in 85 high resolution, non-homologous protein crystal structures for N-cap side-chain rotamer preferences, hydrogen bonding patterns, and solvent accessibilities. We find very strong rotamer preferences that are unique to N-cap sites. The following rules are generally observed for N-capping in alpha-helices: Thr and Ser N-cap side chains adopt the gauche - rotamer, hydrogen bond to the N3 NH and have psi restricted to 164 +/- 8 degrees. Asp and Asn N-cap side chains either adopt the gauche - rotamer and hydrogen bond to the N3 NH with psi = 172 +/- 10 degrees, or adopt the trans rotamer and hydrogen bond to both the N2 and N3 NH groups with psi = 1-7 +/- 19 degrees. With all other N-caps, the side chain is found in the gauche + rotamer so that the side chain does not interact unfavorably with the N-terminus by blocking solvation and psi is unrestricted. An i, i + 3 hydrogen bond from N3 NH to the N-cap backbone C = O in more likely to form at the N-terminus when an unfavorable N-cap is present. In the 3(10)-helix Asn and Asp remain favorable N-caps as they can hydrogen bond to the N2 NH while in the trans rotamer; in contrast, Ser and Thr are disfavored as their preferred hydrogen bonding partner (N3 NH) is inaccessible. This suggests that Ser is the optimum choice of N-cap when alpha-helix formation is to be encouraged while 3(10)-helix formation discouraged. The strong energetic and structural preferences found for N-caps, which differ greatly from positions within helix interiors, suggest that N-caps should be treated explicitly in any consideration of helical structure in peptides or proteins.

Recent advances in helix-coil theory

Peptide helices in solution form a complex mixture of all helix, all coil or, most frequently, central helices with frayed coil ends. In order to interpret experiments on helical peptides and make theoretical predictions on helices, it is therefore essential to use a helix-coil theory that takes account of this equilibrium. The original Zimm-Bragg and Lifson-Roig helix-coil theories have been greatly extended in the last 10 years to include additional interactions. These include preferences for the N-cap, N1, N2, N3 and C-cap positions, capping motifs, helix dipoles, side chain interactions and 3(10)-helix formation. These have been applied to determine energies for these preferences from experimental data and to predict the helix contents of peptides. This review discusses these newly recognised structural features of helices and how they have been included in helix-coil models.

Side-chain structures in the first turn of the alpha-helix

The first three residues at the N terminus of the alpha-helix are called N1, N2 and N3. We surveyed 2102 alpha-helix N termini in 298 high-resolution, non-homologous protein crystal structures for N1, N2 and N3 amino acid and side-chain rotamer propensities and hydrogen-bonding patterns. We find strong structural preferences that are unique to these sites. The rotamer distributions as a function of amino acid identity and position in the helix are often explained in terms of hydrogen-bonding interactions to the free N1, N2 and N3 backbone NH groups. Notably, the "good N2" amino acid residues Gln, Glu, Asp, Asn, Ser, Thr and His preferentially form i, i or i,i+1 hydrogen bonds to the backbone, though this is hindered by good N-caps (Asp, Asn, Ser, Thr and Cys) that compete for these hydrogen bond donors. We find a number of specific side-chain to side-chain interactions between N1 and N2 or between the N-cap and N2 or N3, such as Arg(N-cap) to Asp(N2). The strong energetic and structural preferences found for N1, N2 and N3, which differ greatly from positions within helix interiors, suggest that these sites should be treated explicitly in any consideration of helical structure in peptides or proteins.

Predicting Enzyme Class From Protein Structure Without Alignments

Methods for predicting protein function from structure are becoming more important as the rate at which structures are solved increases more rapidly than experimental knowledge. As a result, protein structures now frequently lack functional annotations. The majority of methods for predicting protein function are reliant upon identifying a similar protein and transferring its annotations to the query protein. This method fails when a similar protein cannot be identified, or when any similar proteins identified also lack reliable annotations. Here, we describe a method that can assign function from structure without the use of algorithms reliant upon alignments. Using simple attributes that can be calculated from any crystal structure, such as secondary structure content, amino acid propensities, surface properties and ligands, we describe each enzyme in a non-redundant set. The set is split according to Enzyme Classification (EC) number. We combine the predictions of one-class versus one-class support vector machine models to make overall assignments of EC number to an accuracy of 35% with the top-ranked prediction, rising to 60% accuracy with the top two ranks. In doing so we demonstrate the utility of simple structural attributes in protein function prediction and shed light on the link between structure and function. We apply our methods to predict the function of every currently unclassified protein in the Protein Data Bank.

Design of an N-methylated peptide inhibitor of alpha-synuclein aggregation guided by solid-state NMR

Many neurodegenerative diseases are associated with the aggregation of misfolded proteins into amyloid oligomers or fibrils that are deposited as pathological lesions within areas of the brain. An attractive therapeutic strategy for preventing or ameliorating amyloid formation is to identify agents that inhibit the onset or propagation of protein aggregation. Here we demonstrate how solid-state nuclear magnetic resonance (ssNMR) may be used to identify key residues within amyloidogenic protein sequences that may be targeted to inhibit the aggregation of the host protein. For alpha-synuclein, the major protein component of Lewy bodies associated with Parkinson's disease, we have used a combination of ssNMR and biochemical data to identify the key region for self-aggregation of the protein as residues 77-82 (VAQKTV). We used our new structural information to design a peptide derived from residues 77 to 82 of alpha-synuclein with an N-methyl group at the C-terminal residue, which was able to disrupt the aggregation of alpha-synuclein. Thus, we have shown how structural data obtained from ssNMR can guide the design of modified peptides for use as amyloid inhibitors, as a primary step toward developing therapeutic compounds for prevention and/or treatment of amyloid diseases.

Addition of side chain interactions to modified Lifson-Roig helix-coil theory: application to energetics of phenylalanine-methionine interactions

We introduce here i, i + 3 and i, i + 4 side chain interactions into the modified Lifson-Roig helix-coil theory of Doig et al. (1994, Biochemistry 33:3396-3403). The helix/coil equilibrium is a function of initiation, propagation, capping, and side chain interaction parameters. If each of these parameters is known, the helix content of any isolated peptide can be predicted. The model considers every possible conformation of a peptide, is not limited to peptides with only a single helical segment, and has physically meaningful parameters. We apply the theory to measure the i, i + 4 interaction energies between Phe and Met side chains. Peptides with these residues spaced i, i + 4 are significantly more helical than controls where they are spaced i, i + 5. Application of the model yields delta G for the Phe-Met orientation to be -0.75 kcal.mol-1, whereas that for the Met-Phe orientation is -0.54 kcal.mol-1. These orientational preferences can be explained, in part, by rotamer preferences for the interacting side chains. We place Phe-Met i, i + 4 at the N-terminus, the C-terminus, and in the center of the host peptide. The model quantitatively predicts the observed helix contents using a single parameter for the side chain-side chain interaction energy. This result indicates that the model works well even when the interaction is at different locations in the helix.

Stabilizing interactions between aromatic and basic side chains in alpha-helical peptides and proteins. Tyrosine effects on helix circular dichroism

Here we investigate the structures and energetics of interactions between aromatic (Phe or Tyr) and basic (Lys or Arg) amino acids in alpha-helices. Side chain interaction energies are measured using helical peptides, by quantifying their helicities with circular dichroism at 222 nm and interpreting the results with Lifson-Roig-based helix/coil theory. A difficulty in working with Tyr is that the aromatic ring perturbs the CD spectrum, giving an incorrect helicity. We calculated the effect of Tyr on the CD at 222 nm by deriving the intensities of the bands directly from the electronic and magnetic transition dipole moments through the rotational strengths corresponding to each excited state of the polypeptide. This gives an improved value of the helix preference of Tyr (from 0.48 to 0.35) and a correction to the helicity for the peptides containing Tyr. We find that Phe-Lys, Lys-Phe, Phe-Arg, Arg-Phe, and Tyr-Lys are all stabilizing by -0.10 to -0.18 kcal.mol-1 when placed i, i + 4 on the surface of a helix in aqueous solution, despite the great difference in polarity between these residues. Interactions between these side chains have previously been attributed to cation-pi bonds. A survey of protein structures shows that they are in fact predominantly hydrophobic interactions between the CH2 groups of Lys or Arg and the aromatic rings.