The Alzheimer's peptides Abeta40 and 42 adopt distinct conformations in water: a combined MD / NMR study

Structures of beta-amyloid peptide 1-40, 1-42, and 1-55-the 672-726 fragment of APP-in a membrane environment with implications for interactions with gamma-secretase

Aggregation of Amyloid beta (Abeta) peptide has been linked to the neurodegenerative Alzheimer's Disease and implicated in other amyloid diseases including cerebral amyloid angiopathy. Abeta peptide is generated by cleavage of the amyloid precursor protein (APP) by transmembrane proteases. It is crucial to determine the structures of beta-amyloid peptides in a membrane to provide a molecular basis for the cleavage mechanism. We report the structures of amyloid beta peptide (Abeta(1-40) and Abeta(1-42)) as well as the 672-726 fragment of APP (referred to as Abeta(1-55)) in a membrane environment determined by replica-exchange molecular dynamics simulation. Abeta(1-40) is found to have two helical domains A (13-22) and B(30-35) and a type I beta-turn at 23-27. The peptide is localized at the interface between membrane and solvent. Substantial fluctuations in domain A are observed. The dominant simulated tertiary structure of Abeta(1-40) is observed to be similar to the simulated Abeta(1-42) structure. However, there are differences observed in the overall conformational ensemble, as characterized by the two-dimensional free energy surfaces. The fragment of APP (Abeta(1-55)) is observed to have a long transmembrane helix. The position of the transmembrane region and ensemble of membrane structures are elucidated. The conformational transition between the transmembrane Abeta(1-55) structure, prior to cleavage, and the Abeta(1-40) structure, following cleavage, is proposed.

In silico study of full-length amyloid beta 1-42 tri- and penta-oligomers in solution

Amyloid oligomers are considered to play causal roles in the pathogenesis of amyloid-related degenerative diseases including Alzheimer's disease. Using MD simulation techniques, we explored the contributions of the different structural elements of trimeric and pentameric full-length Abeta1-42 aggregates in solution to their stability and conformational dynamics. We found that our models are stable at a temperature of 310 K, and converge toward an interdigitated side-chain packing for intermolecular contacts within the two beta-sheet regions of the aggregates: beta1 (residues 18-26) and beta2 (residues 31-42). MD simulations reveal that the beta-strand twist is a characteristic element of Abeta-aggregates, permitting a compact, interdigitated packing of side chains from neighboring beta-sheets. The beta2 portion formed a tightly organized beta-helix, whereas the beta1 portion did not show such a firm structural organization, although it maintained its beta-sheet conformation. Our simulations indicate that the hydrophobic core comprising the beta2 portion of the aggregate is a crucial stabilizing element in the Abeta aggregation process. On the basis of these structure-stability findings, the beta2 portion emerges as an optimal target for further antiamyloid drug design.

Structure of the amyloid-beta (1-42) monomer absorbed to model phospholipid bilayers: a molecular dynamics study

The amyloid-beta (Abeta) peptide, the 39 to 43 amino acid peptide that plays a substantial role in Alzheimer's disease, has been shown to interact strongly with lipids both in vitro and in vivo. Abeta-lipid interactions have been proposed as a considerable factor in accelerating Abeta aggregation through the templating role of membranes in aggregation disorders. Previous work has shown that anionic lipids are able to significantly increase Abeta aggregation rate and induce a structural conversion in Abeta from a random coil to a beta-structure that is similar to the monomer structure observed in mature fibrils. However, it is unclear if this structural change occurs with the Abeta monomer because of direct interactions with the lipids or if the structural change results from protein-protein interactions during oligomerization. We use extensive replica exchange molecular dynamics simulations of an Abeta monomer bound to a homogeneous model zwitterionic or anionic lipid bilayer. From these simulations, we do not observe any significant beta-structure formation except for a small, unstable beta-hairpin formed on the anionic dioleylphosphatidylserine bilayer. Further, we see that the Asp23-Lys28 salt bridge that plays a role in beta-hairpin formation is not substantially formed on the bilayer surface and that Lys28 preferentially interacts with lipids when bound to the bilayer. These results suggest that the structural conversion seen in experiments are not due to the ordering of monomeric Abeta on the bilayer surface but are a result of protein-protein interactions enhanced by Abeta binding to the cell membrane.

Principles governing oligomer formation in amyloidogenic peptides

Identifying the principles that describe the formation of protein oligomers and fibrils with distinct morphologies is a daunting problem. Here we summarize general principles of oligomer formation gleaned from molecular dynamics simulations of Abeta-peptides. The spectra of high free energy structures sampled by the monomer provide insights into the plausible fibril structures, providing a rationale for the 'strain phenomenon.' Heterogeneous growth dynamics of small oligomers of Abeta(16-22), whose lowest free energy structures are like nematic droplets, can be broadly described using a two-stage dock-lock mechanism. In the growth process, water is found to play various roles depending on the oligomer size, and peptide length, and sequence. Water may be an explicit element of fibril structure linked to various fibril morphologies.

Poly-N-methylated amyloid beta-peptide (Abeta) C-terminal fragments reduce Abeta toxicity in vitro and in Drosophila melanogaster

Alzheimer's disease (AD), an age related neurodegenerative disorder, threatens to become a major health-economic problem. Assembly of 40- or 42-residue amyloid beta-peptides (Abeta) into neurotoxic oligo-/polymeric beta-sheet structures is an important pathogenic feature in AD, thus, inhibition of this process has been explored to prevent or treat AD. The C-terminal part plays an important role in Abeta aggregation, but most Abeta aggregation inhibitors have targeted the central region around residues 16-23. Herein, we synthesized hexapeptides with varying extents of N-methylation based on residues 32-37 of Abeta, to target its C-terminal region. We measured the peptides' abilities to retard beta-sheet and fibril formation of Abeta and to reduce Abeta neurotoxicity. A penta-N-methylated peptide was more efficient than peptides with 0, 2, or 3 N-methyl groups. This penta-N-methylated peptide moreover increased life span and locomotor activity in Drosophila melanogaster flies overexpressing human Abeta(1-42).

Molecular dynamics simulations of Ibuprofen binding to Abeta peptides

Using replica exchange molecular dynamics simulations and the implicit solvent model we probed binding of ibuprofen to Abeta(10-40) monomers and amyloid fibrils. We found that the concave (CV) fibril edge has significantly higher binding affinity for ibuprofen than the convex edge. Furthermore, binding of ibuprofen to Abeta monomers, as compared to fibrils, results in a smaller free energy gain. The difference in binding free energies is likely to be related to the presence of the groove on the CV fibril edge, in which ibuprofen tends to accumulate. The confinement effect of the groove promotes the formation of large low-energy ibuprofen clusters, which rarely occur on the surface of Abeta monomers. These observations led us to suggest that the ibuprofen binding mechanism for Abeta fibrils is different from that for monomers. In general, ibuprofen shows a preference to bind to those regions of Abeta monomers (amino terminal) and fibrils (the CV edge) that are also the primary aggregation interfaces. Based on our findings and on available experimental data, we propose a rationale for the ibuprofen antiaggregation effect.

Interpeptide interactions induce helix to strand structural transition in Abeta peptides

Replica exchange molecular dynamics and all-atom implicit solvent model are used to compute the structural propensities in Abeta monomers, dimers, and Abeta peptides bound to the edge of amyloid fibril. These systems represent, on an approximate level, different stages in Abeta aggregation. Abeta monomers are shown to form helical structure in the N-terminal (residues 13 to 21). Interpeptide interactions in Abeta dimers and, especially, in the peptides bound to the fibril induce a dramatic shift in the secondary structure, from helical states toward beta-strand conformations. The sequence region 10-23 in Abeta peptide is found to form most of interpeptide interactions upon aggregation. Simulation results are tested by comparing the chemical shifts in Abeta monomers computed from simulations and obtained experimentally. Possible implications of our simulations for designing aggregation-resistant variants of Abeta are discussed.

Structures and thermodynamics of Alzheimer's amyloid-beta Abeta(16-35) monomer and dimer by replica exchange molecular dynamics simulations: implication for full-length Abeta fibrillation

Many proteins display a strand-loop-strand motif in their amyloid fibrillar states. For instance, the amyloid beta-protein, Abeta1-40, associated with Alzheimer's disease, displays a loop at positions 22-28 in its amyloid fibril state. It has been suggested that this loop could appear early in the aggregation process, but quantitative information regarding its presence in small oligomers remains scant. Because residues 1-15 are disordered in Abeta1-42 fibrils and Abeta10-35 forms fibrils in vitro, we select the peptide Abeta16-35, centered on residues 22-28 and determine the structures and thermodynamics of the monomer and dimer using coarse-grained implicit solvent replica exchange molecular dynamics simulations. Our simulations totalling 5 mus for the monomer and 12 micros for the dimer show no sign of strong secondary structure signals in both instances and the significant impact of dimerization on the global structure of Abeta16-35. They reveal however that the loop 22-28 acts as a quasi-independent unit in both species. The loop structure ensemble we report in Abeta16-35 monomer and dimer has high similarity to the loop formed by the Abeta21-30 peptide in solution and, to a lesser extent, to the loop found in Abeta1-40 fibrils. We discuss the implications of our findings on the assembly of full-length Abeta.

Probing energetics of Abeta fibril elongation by molecular dynamics simulations

Using replica exchange molecular dynamics simulations and an all-atom implicit solvent model, we probed the energetics of Abeta(10-40) fibril growth. The analysis of the interactions between incoming Abeta peptides and the fibril led us to two conclusions. First, considerable variations in fibril binding propensities are observed along the Abeta sequence. The peptides in the fibril and those binding to its edge interact primarily through their N-termini. Therefore, the mutations affecting the Abeta positions 10-23 are expected to have the largest impact on fibril elongation compared with those occurring in the C-terminus and turn. Second, we performed weak perturbations of the binding free energy landscape by scanning partial deletions of side-chain interactions at various Abeta sequence positions. The results imply that strong side-chain interactions--in particular, hydrophobic contacts--impede fibril growth by favoring disordered docking of incoming peptides. Therefore, fibril elongation may be promoted by moderate reduction of Abeta hydrophobicity. The comparison with available experimental data is presented.

Modeling of peptides connecting the ligand-binding and transmembrane domains in the GluR2 glutamate receptor

Ligand-gated Glutamate receptors (GluR) mediate synaptic signals in the nervous system. Ionotropic GluRs of AMPA type, the subject of this study, are tetrameric assemblies of monomer subunits, each of which is constructed in a modular fashion from functional subdomains. The extracellular ligand-binding domain (LBD) changes its conformation upon binding of an agonist ligand followed by opening of a transmembrane (TM) ion channel. Peptides connecting the LBD and TM domains facilitate gating of the channel, and their structure and composition are important for the receptor functioning. In this study, we used replica exchange molecular dynamics (REMD) simulations to model S1M1 and S2M3 connecting peptides of the GluR2 receptor in two implicit solvents, water and interfacial water/lipid medium characterized by lower polarity. Propensity of these peptides to form helical structures was analyzed using helicity measure derived from the free energy of the simulated ensembles of structures. The S1M1 and S2M3 connecting peptides were not helical in our simulations in both dielectric environments in the absence of the rest of the protein. The structures of the LBD fragment with known high-resolution alpha-helical structure and of the TM3 helix were successfully predicted in the simulations, which in part validate our results. The S2M3 peptide, which is important in gating, formed a well-defined coil structure and salt-bridges with the S2 domain. The S1M1 peptide formed a loop structure via formation of internal salt-bridges. Potential implications of these structures on function of the receptor are discussed.

The interplay of catalysis and toxicity by amyloid intermediates on lipid bilayers: insights from type II diabetes

The dynamics, energies, and structures governing protein folding are critical to biological function. Amyloidoses are a class of disease defined, in part, by the misfolding and aggregation of functional protein precursors into fibrillar states. Amyloid fibers contribute to the pathology of many diseases, including type II diabetes, Alzheimer's, and Parkinson's. In these disorders, amyloid fibers are present in affected tissues. However, it has become clear that intermediate states, rather than mature fibers, represent the cytotoxic species. In this review, we focus particularly on lipid bilayer-bound intermediates. Remarkably, the precursors of these fibers are intrinsically disordered, and yet catalysis of beta-sheet formation appears to be mediated by the stabilization of alpha-helical states. On the lipid bilayer, these intermediate species have been implicated as cytotoxic through elimination of ionic homeostasis. Recent advances are enabling insights at a molecular level that promise to provide meaningful targets for the development of therapeutics.

Ca(2+), within the physiological concentrations, selectively accelerates Abeta42 fibril formation and not Abeta40 in vitro

Alzheimer's disease (AD) in humans is a common progressive neurodegenerative disease, associated with cognitive dysfunction, memory loss and neuronal loss. Alzheimer peptides Abeta40 and Abeta42 are precursors of the amyloid fibers that accumulate in the brain of patients. These peptides misfold and the monomers aggregate to neurotoxic oligomers and fibrils. Thus, the aggregation kinetics of these peptides is central to understanding the etiology of AD. Using size exclusion chromatography as well as filtration methods, we report here that Ca(2+) ions at physiological concentrations greatly accelerate the rate of aggregation of Abeta42 to form intermediate soluble associated species and fibrils. In the presence of 1 or 2 mM Ca(2+), CD spectra indicated that the secondary structure of Abeta42 changed from an unfolded to a predominantly beta-sheet conformation. These concentrations of Ca(2+) greatly decreased the lag time for Abeta42 fibril formation, measured with thioflavin T. However, the elongation rate was apparently unaffected. Ca(2+) appears to predominantly accelerate the nucleation stage of Abeta42 on pathway to the Alzheimer's fibril formation. Unlike Abeta42, Ca(2+) was not observed to trigger similar effect at any stage during the study of fibrillation kinetics of Abeta40 by any techniques. Abeta40 and Abeta42 seem to have distinct aggregation pathways.

Assembly dynamics of two-beta sheets revealed by molecular dynamics simulations

The assembly dynamics of two beta sheets with different initial separation distances are explored by multiple all-atom molecular dynamics simulations with the presence of explicit water solvent. The beta sheet is composed of seven identical peptides in an antiparallel fashion. The peptide sequence is the 20-29 segment of human Islet amyloid polypeptide. Our simulations show that the assembly occurs not only in the lateral direction but also along the longitudinal direction, which provides a new insight into the assembly pathway at the early stage of fibril elongation. Based on Poisson-Boltzmann free energy analysis and quasiharmonic configuration entropy estimation, the entropic contribution is found to play an important role in the longitudinal assembly. Moreover, a possible oligomeric state with cyclic form is suggested based on one assembly model found in the simulations, illustrating the polymorphic nature of aggregation of the amyloidogenic peptide.

Probing the effect of amino-terminal truncation for Abeta1-40 peptides

We examine the effect of deletion of the amino-terminal (residues 1-9) on the structure and energetics of Abeta1-40 peptides. To this end, we use replica exchange molecular dynamics to compare the conformational ensembles of Abeta1-40 and amino-truncated Abeta10-40 monomers and dimers. Overall, the deletion of the amino-terminal appears to cause minor structural and energetic changes in Abeta monomers and dimers. More specifically, our findings are as follows: (1) there is a small but discernible conversion of beta-strand structure into helix upon amino-terminal deletion, (2) secondary structure changes due to truncation are caused by missing side chain interactions formed by the amino-terminal, and (3) the amino-terminal together with the central sequence region (residues 10-23) represents the primary aggregation interface in Abeta1-40 dimers. The amino-truncated Abeta10-40 retains this aggregation interface, which is reduced to the central sequence region. We argue that the analysis of available experimental data supports our conclusions. Our findings also suggest that amino-truncated Abeta10-40 peptide is an adequate model for studying Abeta1-40 aggregation.

The structure of Abeta42 C-terminal fragments probed by a combined experimental and theoretical study

The C-terminus of amyloid beta-protein (Abeta) 42 plays an important role in this protein's oligomerization and may therefore be a good therapeutic target for the treatment of Alzheimer's disease. Certain C-terminal fragments (CTFs) of Abeta42 have been shown to disrupt oligomerization and to strongly inhibit Abeta42-induced neurotoxicity. Here we study the structures of selected CTFs [Abeta(x-42); x=29-31, 39] using replica exchange molecular dynamics simulations and ion mobility mass spectrometry. Our simulations in explicit solvent reveal that the CTFs adopt a metastable beta-structure: beta-hairpin for Abeta(x-42) (x=29-31) and extended beta-strand for Abeta(39-42). The beta-hairpin of Abeta(30-42) is converted into a turn-coil conformation when the last two hydrophobic residues are removed, suggesting that I41 and A42 are critical in stabilizing the beta-hairpin in Abeta42-derived CTFs. The importance of solvent in determining the structure of the CTFs is further highlighted in ion mobility mass spectrometry experiments and solvent-free replica exchange molecular dynamics simulations. A comparison between structures with solvent and structures without solvent reveals that hydrophobic interactions are critical for the formation of beta-hairpin. The possible role played by the CTFs in disrupting oligomerization is discussed.

Interaction between amyloid-beta (1-42) peptide and phospholipid bilayers: a molecular dynamics study

The amyloid-beta (Abeta) peptide is a key aggregate species in Alzheimer's disease. Although important aspects of Abeta peptide aggregation are understood, the initial stage of aggregation from monomer to oligomer is still not clear. One potential mediator of this early aggregation process is interactions of Abeta with anionic cell membranes. We used unconstrained and umbrella sampling molecular dynamics simulations to investigate interactions between the 42-amino acid Abeta peptide and model bilayers of zwitterionic dipalmitoylphosphatidylcholine (DPPC) lipids and anionic dioleoylphosphatidylserine (DOPS) lipids. Using these methods, we determined that Abeta is attracted to the surface of DPPC and DOPS bilayers over the small length scales used in these simulations. We also found supporting evidence that the charge on both the bilayer surface and the peptide affects the free energy of binding of the peptide to the bilayer surface and the distribution of the peptide on the bilayer surface. Our work demonstrates that interactions between the Abeta peptide and lipid bilayer promotes a peptide distribution on the bilayer surface that is prone to peptide-peptide interactions, which can influence the propensity of Abeta to aggregate into higher-order structures.

Influence of preformed Asp23-Lys28 salt bridge on the conformational fluctuations of monomers and dimers of Abeta peptides with implications for rates of fibril formation

Recent experiments have shown that the congener Abeta(1-40)[D23-K28], in which the side chains of charged residues Asp23 and Lys28 are linked by a lactam bridge, forms amyloid fibrils that are structurally similar to the wild type (WT) Abeta peptide, but at a rate that is nearly 1000 times faster. We used all atom molecular dynamics simulations in explicit water, and two force fields, of the WT dimer, a monomer with the lactam bridge (Abeta(10-35)-lactam[D23-K28]), and the monomer and dimers with harmonically constrained D23-K28 salt bridge (Abeta(10-35)[D23-K28]) to understand the origin of the enhanced fibril rate formation. The simulations show that the assembly competent fibril-like monomer (N*) structure, which is present among the conformations sampled by the isolated monomer, with strand conformations in the residues spanning the N and C termini and a bend involving residues D(23) VGSNKG(29), are populated to a much greater extent in Abeta(10-35)[D23-K28] and Abeta(10-35)-lactam[D23-K28] than in the WT, which has negligible probability of forming N*. The salt bridge in N* of Abeta(10-35)[D23-K28], whose topology is similar to that found in the fibril, is hydrated. The reduction in the free energy barrier to fibril formation in Abeta(10-35)[D23-K28] and in Abeta(10-35)-lactam[D23-K28], compared to the WT, arises largely due to entropic restriction which enables the bend formation. A decrease in the entropy of the unfolded state and the lesser penalty for conformational rearrangement including the formation of the salt bridge in Abeta peptides with D23-K28 constraint results in a reduction in the kinetic barrier in the Abeta(1-40)-lactam[D23-K28] congener compared to the WT. The decrease in the barrier, which is related to the free energy cost of forming a bend, is estimated to be in the range (4-7)k(B)T. Although a number of factors determine the growth of fibrils, the decrease in the free energy barrier, relative to the WT, to N* formation is a major factor in the rate enhancement in the fibril formation of Abeta(1-40)[D23-K28] congener. Qualitatively similar results were obtained using simulations of Abeta(9-40) peptides and various constructs related to the Abeta(10-35) systems that were probed using OPLS and CHARMM force fields. We hypothesize that mutations or other constraints that preferentially enhance the population of the N* species would speed up aggregation rates. Conversely, ligands that lock it in the fibril-like N* structure would prevent amyloid formation.

Conformation and dynamics of a rhodamine probe attached at two sites on a protein: implications for molecular structure determination in situ

Replica exchange molecular dynamics (REMD) calculations were used to determine the conformation and dynamics of bifunctional rhodamine probes attached to pairs of cysteines in three model systems: (a) a polyalanine helix, (b) the isolated C helix (residues 53-66) of troponin C, and (c) the C helix of the N-terminal region (residues 1-90) of troponin C (sNTnC). In each case, and for both diastereoisomers of each probe-protein complex, the hydrophobic face of the probe is close to the protein surface, and its carboxylate group is highly solvated. The visible-range fluorescence dipole of the probe is approximately parallel to the line joining the two cysteine residues, as assumed in previous in situ fluorescence polarization studies. The independent rotational motion of the probe with respect to the protein on the nanosecond time scale is highly restricted, in agreement with data from fluorescence polarization and NMR relaxation studies. The detailed interaction of the probe with the protein surface depends on steric factors, electrostatic and hydrophobic interactions, hydrogen bonds, and hydration effects. The interaction is markedly different between diastereoisomers, and multiple preferred conformations exist for a single diasteroisomer. These results show that the combination of the hydrophobic xanthylium moiety of bifunctional rhodamine with the carboxylate substitution in its pendant phenyl ring causes the probe to be immobilized on the protein surface, while the two-site cysteine attachment defines the orientation of its fluorescence dipole. These features allow the orientation of protein components to be accurately determined in situ by polarized fluorescence measurements from bifunctional rhodamine probes.

Amyloid beta-protein monomer folding: free-energy surfaces reveal alloform-specific differences

Alloform-specific differences in structural dynamics between amyloid beta-protein (Abeta) 40 and Abeta42 appear to underlie the pathogenesis of Alzheimer's disease. To elucidate these differences, we performed microsecond timescale replica-exchange molecular dynamics simulations to sample the conformational space of the Abeta monomer and constructed its free-energy surface. We find that neither peptide monomer is unstructured, but rather that each may be described as a unique statistical coil in which five relatively independent folding units exist, comprising residues 1-5, 10-13, 17-22, 28-37, and 39-42, which are connected by four turn structures. The free-energy surfaces of both peptides are characterized by two large basins, comprising conformers with either substantial alpha-helix or beta-sheet content. Conformational transitions within and between these basins are rapid. The two additional hydrophobic residues at the Abeta42 C-terminus, Ile41 and Ala42, significantly increase contacts within the C-terminus, and between the C-terminus and the central hydrophobic cluster (Leu17-Ala21). As a result, the beta-structure of Abeta42 is more stable than that of Abeta40, and the conformational equilibrium in Abeta42 shifts towards beta-structure. These results suggest that drugs stabilizing alpha-helical Abeta conformers (or destabilizing the beta-sheet state) would block formation of neurotoxic oligomers. The atomic-resolution conformer structures determined in our simulations may serve as useful targets for this purpose. The conformers also provide starting points for simulations of Abeta oligomerization-a process postulated to be the key pathogenetic event in Alzheimer's disease.

Simulating oligomerization at experimental concentrations and long timescales: A Markov state model approach

Here, we present a novel computational approach for describing the formation of oligomeric assemblies at experimental concentrations and timescales. We propose an extension to the Markovian state model approach, where one includes low concentration oligomeric states analytically. This allows simulation on long timescales (seconds timescale) and at arbitrarily low concentrations (e.g., the micromolar concentrations found in experiments), while still using an all-atom model for protein and solvent. As a proof of concept, we apply this methodology to the oligomerization of an Abeta peptide fragment (Abeta(21-43)). Abeta oligomers are now widely recognized as the primary neurotoxic structures leading to Alzheimer's disease. Our computational methods predict that Abeta trimers form at micromolar concentrations in 10 ms, while tetramers form 1000 times more slowly. Moreover, the simulation results predict specific intermonomer contacts present in the oligomer ensemble as well as putative structures for small molecular weight oligomers. Based on our simulations and statistical models, we propose a novel mutation to stabilize the trimeric form of Abeta in an experimentally verifiable manner.

C-terminal peptides coassemble into Abeta42 oligomers and protect neurons against Abeta42-induced neurotoxicity

Alzheimer's disease (AD) is an age-related disorder that threatens to become an epidemic as the world population ages. Neurotoxic oligomers of Abeta42 are believed to be the main cause of AD; therefore, disruption of Abeta oligomerization is a promising approach for developing therapeutics for AD. Formation of Abeta42 oligomers is mediated by intermolecular interactions in which the C terminus plays a central role. We hypothesized that peptides derived from the C terminus of Abeta42 may get incorporated into oligomers of Abeta42, disrupt their structure, and thereby inhibit their toxicity. We tested this hypothesis using Abeta fragments with the general formula Abeta(x-42) (x = 28-39). A cell viability screen identified Abeta(31-42) as the most potent inhibitor. In addition, the shortest peptide, Abeta(39-42), also had high activity. Both Abeta(31-42) and Abeta(39-42) inhibited Abeta-induced cell death and rescued disruption of synaptic activity by Abeta42 oligomers at micromolar concentrations. Biophysical characterization indicated that the action of these peptides likely involved stabilization of Abeta42 in nontoxic oligomers. Computer simulations suggested a mechanism by which the fragments coassembled with Abeta42 to form heterooligomers. Thus, Abeta(31-42) and Abeta(39-42) are leads for obtaining mechanism-based drugs for treatment of AD using a systematic structure-activity approach.

Effects of familial Alzheimer's disease mutations on the folding nucleation of the amyloid beta-protein

The effect of single amino acid substitutions associated with the Italian (E22K), Arctic (E22G), Dutch (E22Q) and Iowa (D23N) familial forms of Alzheimer's disease and cerebral amyloid angiopathy on the structure of the 21-30 fragment of the Alzheimer amyloid beta-protein (Abeta) is investigated by replica-exchange molecular dynamics simulations. The 21-30 segment has been shown in our earlier work to adopt a bend structure in solution that may serve as the folding nucleation site for Abeta. Our simulations reveal that the 24-28 bend motif is retained in all E22 mutants, suggesting that mutations involving residue E22 may not affect the structure of the folding nucleation site of Abeta. Enhanced aggregation in Abeta with familial Alzheimer's disease substitutions may result from the depletion of the E22-K28 salt bridge, which destabilizes the bend structure. Alternately, the E22 mutations may affect longer-range interactions outside the 21-30 segment that can impact the aggregation of Abeta. Substituting at residue D23, on the other hand, leads to the formation of a turn rather than a bend motif, implying that in contrast to E22 mutants, the D23N mutant may affect monomer Abeta folding and subsequent aggregation. Our simulations suggest that the mechanisms by which E22 and D23 mutations affect the folding and aggregation of Abeta are fundamentally different.

Contrasting disease and nondisease protein aggregation by molecular simulation

[Figurre: see text]. Protein aggregation can be defined as the sacrifice of stabilizing intrachain contacts of the functional state that are replaced with interchain contacts to form non-functional states. The resulting aggregate morphologies range from amorphous structures without long-range order typical of nondisease proteins involved in inclusion bodies to highly structured fibril assemblies typical of amyloid disease proteins. In this Account, we describe the development and application of computational models for the investigation of nondisease and disease protein aggregation as illustrated for the proteins L and G and the Alzheimer's Abeta systems. In each case, we validate the models against relevant experimental observables and then expand on the experimental window to better elucidate the link between molecular properties and aggregation outcomes. Our studies show that each class of protein exhibits distinct aggregation mechanisms that are dependent on protein sequence, protein concentration, and solution conditions. Nondisease proteins can have native structural elements in the denatured state ensemble or rapidly form early folding intermediates, which offers avenues of protection against aggregation even at relatively high concentrations. The possibility that early folding intermediates may be evolutionarily selected for their protective role against unwanted aggregation could be a useful strategy for reengineering sequences to slow aggregation and increase folding yield in industrial protein production. The observed oligomeric aggregates that we see for nondisease proteins L and G may represent the nuclei for larger aggregates, not just for large amorphous inclusion bodies, but potentially as the seeds of ordered fibrillar aggregates, since most nondisease proteins can form amyloid fibrils under conditions that destabilize the native state. By contrast, amyloidogenic protein sequences such as Abeta 1-40,42 and the familial Alzheimer's disease (FAD) mutants favor aggregation into ordered fibrils once the free-energy barrier for forming a critical nucleus is crossed. However, the structural characteristics and oligomer size of the soluble nucleation species have yet to be determined experimentally for any disease peptide sequence, and the molecular mechanism of polymerization that eventually delineates a mature fibril is unknown. This is in part due to the limited experimental access to very low peptide concentrations that are required to characterize these early aggregation events, providing an opportunity for theoretical studies to bridge the gap between the monomer and fibril end points and to develop testable hypotheses. Our model shows that Abeta 1-40 requires as few as 6-10 monomer chains (depending on sequence) to begin manifesting the cross-beta order that is a signature of formation of amyloid filaments or fibrils assessed in dye-binding kinetic assays. The richness of the oligomeric structures and viable filament and fibril polymorphs that we observe may offer structural clues to disease virulence variations that are seen for the WT and hereditary mutants.

Lipid membrane templates the ordering and induces the fibrillogenesis of Alzheimer's disease amyloid-beta peptide

The lipid membrane has been shown to mediate the fibrillogenesis and toxicity of Alzheimer's disease (AD) amyloid-beta (Abeta) peptide. Electrostatic interactions between Abeta40 and the phospholipid headgroup have been found to control the association and insertion of monomeric Abeta into lipid monolayers, where Abeta exhibited enhanced interactions with charged lipids compared with zwitterionic lipids. To elucidate the molecular-scale structural details of Abeta-membrane association, we have used complementary X-ray and neutron scattering techniques (grazing-incidence X-ray diffraction, X-ray reflectivity, and neutron reflectivity) in this study to investigate in situ the association of Abeta with lipid monolayers composed of either the anionic lipid 1,2-dipalmitoyl-sn-glycero-3-[phospho-rac-(1-glycerol)] (DPPG), the zwitterionic lipid 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), or the cationic lipid 1,2-dipalmitoyl 3-trimethylammonium propane (DPTAP) at the air-buffer interface. We found that the anionic lipid DPPG uniquely induced crystalline ordering of Abeta at the membrane surface that closely mimicked the beta-sheet structure in fibrils, revealing an intriguing templated ordering effect of DPPG on Abeta. Furthermore, incubating Abeta with lipid vesicles containing the anionic lipid 1-palmitoyl-2-oleoyl-sn-glycero-3-[phospho-rac-(1-glycerol)] (POPG) induced the formation of amyloid fibrils, confirming that the templated ordering of Abeta at the membrane surface seeded fibril formation. This study provides a detailed molecular-scale characterization of the early structural fluctuation and assembly events that may trigger the misfolding and aggregation of Abeta in vivo. Our results implicate that the adsorption of Abeta to anionic lipids, which could become exposed to the outer membrane leaflet by cell injury, may serve as an in vivo mechanism of templated-aggregation and drive the pathogenesis of AD.

Energy landscapes of the monomer and dimer of the Alzheimer's peptide Abeta(1-28)

The cytotoxicity of Alzheimer's disease has been linked to the self-assembly of the 4042 amino acid of the amyloid-beta (Abeta) peptide into oligomers. To understand the assembly process, it is important to characterize the very first steps of aggregation at an atomic level of detail. Here, we focus on the N-terminal fragment 1-28, known to form fibrils in vitro. Circular dichroism and NMR experiments indicate that the monomer of Abeta(1-28) is alpha-helical in a membranelike environment and random coil in aqueous solution. Using the activation-relaxation technique coupled with the OPEP coarse grained force field, we determine the structures of the monomer and of the dimer of Abeta(1-28). In agreement with experiments, we find that the monomer is predominantly random coil in character, but displays a non-negligible beta-strand probability in the N-terminal region. Dimerization impacts the structure of each chain and leads to an ensemble of intertwined conformations with little beta-strand content in the region Leu17-Ala21. All these structural characteristics are inconsistent with the amyloid fibril structure and indicate that the dimer has to undergo significant rearrangement en route to fibril formation.

Structure-function relationships of pre-fibrillar protein assemblies in Alzheimer's disease and related disorders

Several neurodegenerative diseases, including Alzheimer's, Parkinson's, Huntington's and prion diseases, are characterized pathognomonically by the presence of intra- and/or extracellular lesions containing proteinaceous aggregates, and by extensive neuronal loss in selective brain regions. Related non-neuropathic systemic diseases, e.g., light-chain and senile systemic amyloidoses, and other organ-specific diseases, such as dialysis-related amyloidosis and type-2 diabetes mellitus, also are characterized by deposition of aberrantly folded, insoluble proteins. It is debated whether the hallmark pathologic lesions are causative. Substantial evidence suggests that these aggregates are the end state of aberrant protein folding whereas the actual culprits likely are transient, pre-fibrillar assemblies preceding the aggregates. In the context of neurodegenerative amyloidoses, the proteinaceous aggregates may eventuate as potentially neuroprotective sinks for the neurotoxic, oligomeric protein assemblies. The pre-fibrillar, oligomeric assemblies are believed to initiate the pathogenic mechanisms that lead to synaptic dysfunction, neuronal loss, and disease-specific regional brain atrophy. The amyloid beta-protein (Abeta), which is believed to cause Alzheimer's disease (AD), is considered an archetypal amyloidogenic protein. Intense studies have led to nominal, functional, and structural descriptions of oligomeric Abeta assemblies. However, the dynamic and metastable nature of Abeta oligomers renders their study difficult. Different results generated using different methodologies under different experimental settings further complicate this complex area of research and identification of the exact pathogenic assemblies in vivo seems daunting. Here we review structural, functional, and biological experiments used to produce and study pre-fibrillar Abeta assemblies, and highlight similar studies of proteins involved in related diseases. We discuss challenges that contemporary researchers are facing and future research prospects in this demanding yet highly important field.

Structure and dynamics of the Abeta(21-30) peptide from the interplay of NMR experiments and molecular simulations

We combine molecular dynamics simulations and new high-field NMR experiments to describe the solution structure of the Abeta(21-30) peptide fragment that may be relevant for understanding structural mechanisms related to Alzheimer's disease. By using two different empirical force-field combinations, we provide predictions of the three-bond scalar coupling constants ((3)J(H(N)H(alpha))), chemical-shift values, (13)C relaxation parameters, and rotating-frame nuclear Overhauser effect spectroscopy (ROESY) crosspeaks that can then be compared directly to the same observables measured in the corresponding NMR experiment of Abeta(21-30). We find robust prediction of the (13)C relaxation parameters and medium-range ROESY crosspeaks by using new generation TIP4P-Ew water and Amber ff99SB protein force fields, in which the NMR validates that the simulation yields both a structurally and dynamically correct ensemble over the entire Abeta(21-30) peptide. Analysis of the simulated ensemble shows that all medium-range ROE restraints are not satisfied simultaneously and demonstrates the structural diversity of the Abeta(21-30) conformations more completely than when determined from the experimental medium-range ROE restraints alone. We find that the structural ensemble of the Abeta(21-30) peptide involves a majority population (approximately 60%) of unstructured conformers, lacking any secondary structure or persistent hydrogen-bonding networks. However, the remaining minority population contains a substantial percentage of conformers with a beta-turn centered at Val24 and Gly25, as well as evidence of the Asp23 to Lys28 salt bridge important to the fibril structure. This study sets the stage for robust theoretical work on Abeta(1-40) and Abeta(1-42), for which collection of detailed NMR data on the monomer will be more challenging because of aggregation and fibril formation on experimental timescales at physiological conditions. In addition, we believe that the interplay of modern molecular simulation and high-quality NMR experiments has reached a fruitful stage for characterizing structural ensembles of disordered peptides and proteins in general.

The beta-strand-loop-beta-strand conformation is marginally populated in beta2-microglobulin (20-41) peptide in solution as revealed by replica exchange molecular dynamics simulations

Solid-state NMR study shows that the 22-residue K3 peptide (Ser(20)-Lys(41)) from beta(2)-microglobulin (beta(2)m) adopts a beta-strand-loop-beta-strand conformation in its fibril state. Residue Pro(32) has a trans conformation in the fibril state of the peptide, while it adopts a cis conformation in the native state of full-length beta(2)m. To get insights into the structural properties of the K3 peptide, and determine whether the strand-loop-strand conformation is encoded at the monomeric level, we run all-atom explicit solvent replica exchange molecular dynamics on both the cis and trans variants. Our simulations show that the conformational space of the trans- and cis-K3 peptides is very different, with 1% of the sampled conformations in common at room temperature. In addition, both variants display only 0.3-0.5% of the conformations with beta-strand-loop-beta-strand character. This finding, compared to results on the Alzheimer's Abeta peptide, suggests that the biases toward aggregation leading to the beta-strand-loop-beta-strand conformation in fibrils are peptide-dependent.

Role of the familial Dutch mutation E22Q in the folding and aggregation of the 15-28 fragment of the Alzheimer amyloid-beta protein

Amyloid fibrils, large ordered aggregates of amyloid beta proteins (Abeta), are clinical hallmarks of Alzheimer's disease (AD). The aggregation properties of amyloid beta proteins can be strongly affected by single-point mutations at positions 22 and 23. The Dutch mutation involves a substitution at position 22 (E22Q) and leads to increased deposition rates of the AbetaE22Q peptide onto preseeded fibrils. We investigate the effect of the E22Q mutation on two key regions involved in the folding and aggregation of the Abeta peptide through replica exchange molecular dynamics simulations of the 15-28 fragment of the Abeta peptide. The Abeta15-28 peptide encompasses the 22-28 region that constitutes the most structured part of the Abeta peptide (the E22-K28 bend), as well as the central hydrophobic cluster (CHC) (segment 17-21), the primary docking site for Abeta monomers depositing onto fibrils. Our simulations show that the 22-28 bend is preserved in the Abeta(15-28) peptide and that the CHC, which is mostly unstructured, interacts with this bend region. The E22Q mutation does not affect the structure of the bend but weakens the interactions between the CHC and the bend. This leads to an increased population of beta-structure in the CHC. Our analysis of the fibril elongation reaction reveals that the CHC adopts a beta-strand conformation in the transition state ensemble, and that the E22Q mutation increases aggregation rates by lowering the barrier for Abeta monomer deposition onto a fibril. Thermodynamic signatures of this enhanced fibrillization process from our simulations are in good agreement with experimental observations.

Computational study on the structural diversity of amyloid Beta Peptide (abeta(10-35)) oligomers

We studied the oligomerization of Alzheimer amyloid beta peptide (Abeta) using a replica exchange molecular dynamics (REMD) simulation. The simulation was performed with Abeta(10-35) dimers, trimers, and tetramers. Extensive REMD simulations illustrated several possible oligomer conformations. As the size of the oligomer increased from a dimer to a tetramer, the number of possible configurations was reduced. We identified all the possible conformations for each oligomer and characterized their temperature dependence. It was found that the detailed structures of the oligomers, which may act as folding intermediates, are highly sensitive to the parameters of the simulation environment such as temperature and concentration. Structural diversities of Abeta oligomers suggest multiple pathways of the aggregation process.

Methyl dynamics of the amyloid-β peptides Aβ40 and Aβ42

To probe the role of side chain dynamics in Abeta aggregation, we studied the methyl dynamics of native Abeta40 and Abeta42 by measuring cross relaxation rates with interleaved data collection. The methyl groups in the C-terminus are in general more rigid in Abeta42 than in Abeta40, consistent with previous results from backbone (15)N dynamics. This lends support to the hypothesis that a rigid C-terminus in Abeta42 may serve as an internal aggregation seed. Interestingly, two methyl groups of V18 located in the central hydrophobic cluster are more mobile in Abeta42 than in Abeta40, most likely due to the paucity of V18 intra-molecular interactions in Abeta42. V18 may then be more available for inter-molecular interactions to form Abeta42 aggregates. Thus, the side chain mobility of the central hydrophobic cluster may play an important role in Abeta aggregation and may contribute to the difference in aggregation propensity between Abeta40 and Abeta42.