Probing the origins of increased activity of the E22Q "Dutch" mutant Alzheimer's beta-amyloid peptide

High-resolution probing of early events in amyloid-β aggregation related to Alzheimer's disease

In Alzheimer's disease (AD), soluble oligomers of amyloid-β (Aβ) are emerging as a crucial entity in driving disease progression as compared to insoluble amyloid deposits. The lacuna in establishing the structure to function relationship for Aβ oligomers prevents the development of an effective treatment for AD. While the transient and heterogeneous properties of Aβ oligomers impose many challenges for structural investigation, an effective use of a combination of NMR techniques has successfully identified and characterized them at atomic-resolution. Here, we review the successful utilization of solution and solid-state NMR techniques to probe the aggregation and structures of small and large oligomers of Aβ. Biophysical studies utilizing the commonly used solution and 19F based NMR experiments to identify the formation of small size early intermediates and to obtain their structures, and dock-lock mechanism of fiber growth at atomic-resolution are discussed. In addition, the use of proton-detected magic angle spinning (MAS) solid-state NMR experiments to obtain high-resolution insights into the aggregation pathways and structures of large oligomers and other aggregates is also presented. We expect these NMR based studies to be valuable for real-time monitoring of the depletion of monomers and the formation of toxic oligomers and high-order aggregates under a variety of conditions, and to solve the high-resolution structures of small and large size oligomers for most amyloid proteins, and therefore to develop inhibitors and drugs.

Flexibility of short DNA helices with finite-length effect: From base pairs to tens of base pairs

Flexibility of short DNA helices is important for the biological functions such as nucleosome formation and DNA-protein recognition. Recent experiments suggest that short DNAs of tens of base pairs (bps) may have apparently higher flexibility than those of kilo bps, while there is still the debate on such high flexibility. In the present work, we have studied the flexibility of short DNAs with finite-length of 5-50 bps by the all-atomistic molecular dynamics simulations and Monte Carlo simulations with the worm-like chain model. Our microscopic analyses reveal that short DNAs have apparently high flexibility which is attributed to the significantly strong bending and stretching flexibilities of ∼6 bps at each helix end. Correspondingly, the apparent persistence length lp of short DNAs increases gradually from ∼29 nm to ∼45 nm as DNA length increases from 10 to 50 bps, in accordance with the available experimental data. Our further analyses show that the short DNAs with excluding ∼6 bps at each helix end have the similar flexibility with those of kilo bps and can be described by the worm-like chain model with lp ∼ 50 nm.

Effect of the English familial disease mutation (H6R) on the monomers and dimers of Aβ40 and Aβ42

The self-assembly of the amyloid beta (Aβ) peptides into senile plaques is the hallmark of Alzheimer's disease. Recent experiments have shown that the English familial disease mutation (H6R) speeds up the fibril formation process of alloforms Aβ40 and Aβ42 peptides altering their toxicity to cells. We used all-atom molecular dynamics simulations at microsecond time scales with the OPLS-AA force field and TIP4P explicit water model to study the structural dynamics of the monomer and dimer of H6R sequences of both peptides. The reason behind the self-assembly acceleration is common that upon mutation the net charge is reduced leading to the weaker repulsive interaction between chains that facilitates the peptide association. In addition, our estimation of the solvation free energy shows that the mutation enhances the hydrophobicity of both peptides speeding up their aggregation. However, we can show that the acceleration mechanisms are different for different peptides: the rate of fibril formation of Aβ42 increases due to increased β-structure at the C-terminal in both monomer and dimer and enhanced stability of salt bridge Asp23-Lys28 in monomer, while the enhancement of turn at residues 25-29 and reduction of coil in regions 10-13, 26-19, and 30-34 would play the key role for Aβ40. Overall, our study provides a detailed atomistic picture of the H6R-mediated conformational changes that are consistent with the experimental findings and highlights the important role of the N-terminal in Aβ peptide aggregation.

Effect of the Tottori familial disease mutation (D7N) on the monomers and dimers of Aβ40 and Aβ42

Recent experiments have shown that the mutation Tottori (D7N) alters the toxicity, assembly and rate of fibril formation of the wild type (WT) amyloid beta (Aβ) Aβ40 and Aβ42 peptides. We used all-atom molecular dynamics simulations in explicit solvent of the monomer and dimer of both alloforms with their WT and D7N sequences. The monomer simulations starting from a random coil and totaling 3 μs show that the D7N mutation changes the fold and the network of salt bridges in both alloforms. The dimer simulations starting from the amyloid fibrillar states and totaling 4.4 μs also reveal noticeable changes in terms of secondary structure, salt bridge, and topology. Overall, this study provides physical insights into the enhanced rate of fibril formation upon D7N mutation and an atomic picture of the D7N-mediated conformational change on Aβ40 and Aβ42 peptides.

Effects of familial mutations on the monomer structure of Aβ₄₂

Amyloid beta (Aβ) peptide plays an important role in Alzheimer's disease. A number of mutations in the Aβ sequence lead to familial Alzheimer's disease, congophilic amyloid angiopathy, or hereditary cerebral hemorrhage with amyloid. Using molecular dynamics simulations of ∼200 μs for each system, we characterize and contrast the consequences of four pathogenic mutations (Italian, Dutch, Arctic, and Iowa) for the structural ensemble of the Aβ monomer. The four familial mutations are found to have distinct consequences for the monomer structure.

Role of cholesterol in APP metabolism and its significance in Alzheimer's disease pathogenesis

Alzheimer's disease (AD) is a complex multifactorial neurodegenerative disorder believed to be initiated by accumulation of amyloid β (Aβ)-related peptides derived from proteolytic processing of amyloid precursor protein (APP). Research over the past two decades provided a mechanistic link between cholesterol and AD pathogenesis. Genetic polymorphisms in genes regulating the pivotal points in cholesterol metabolism have been suggested to enhance the risk of developing AD. Altered neuronal membrane cholesterol level and/or subcellular distribution have been implicated in aberrant formation, aggregation, toxicity, and degradation of Aβ-related peptides. However, the results are somewhat contradictory and we still do not have a complete understanding on how cholesterol can influence AD pathogenesis. In this review, we summarize our current understanding on the role of cholesterol in regulating the production/function of Aβ-related peptides and also examine the therapeutic potential of regulating cholesterol homeostasis in the treatment of AD pathology.

Investigating how peptide length and a pathogenic mutation modify the structural ensemble of amyloid beta monomer

The aggregation of amyloid beta (Aβ) peptides plays an important role in the development of Alzheimer's disease. Despite extensive effort, it has been difficult to characterize the secondary and tertiary structure of the Aβ monomer, the starting point for aggregation, due to its hydrophobicity and high aggregation propensity. Here, we employ extensive molecular dynamics simulations with atomistic protein and water models to determine structural ensembles for Aβ(42), Aβ(40), and Aβ(42)-E22K (the Italian mutant) monomers in solution. Sampling of a total of >700 microseconds in all-atom detail with explicit solvent enables us to observe the effects of peptide length and a pathogenic mutation on the disordered Aβ monomer structural ensemble. Aβ(42) and Aβ(40) have crudely similar characteristics but reducing the peptide length from 42 to 40 residues reduces β-hairpin formation near the C-terminus. The pathogenic Italian E22K mutation induces helix formation in the region of residues 20-24. This structural alteration may increase helix-helix interactions between monomers, resulting in altered mechanism and kinetics of Aβ oligomerization.

Comparative fibril formation of analogs corresponding to the (12-24) segment of the β-amyloid peptide

The (1-42) β-amyloid peptide is a main component of the plaques found in the brain of patients suffering from the Alzheimer's disease. As the single substitution of Glu for Gln at position 22 of this peptide seems to be responsible for the manifestation of the more severe amyloidosis (Dutch-type), we decided to evaluate the aggregation characteristics of peptide analogs interchanging Glu and Gln residues at positions 22 and also 15 in the minor (12-24) (VHHQ(15)KLVFFAE(22)DV) fragment. The Q15Q22, E15E22, E15Q22 and the native Q15E22 were compared to the (1-42) β-amyloid peptide in terms of fibril or structured aggregates formation propensity. In contrast to a rather similar solubility data measured of all analogs, fluorescence and light scattering methods indicated that only Q15E22 and Q15Q22 displayed relevant fibril formation capacity. Conversely, E15E22 and E15Q22 were not capable of the formation of this type of structure thus suggesting a key role for the Q(15) residue in the unique aggregation characteristic of the β-amyloid peptide.

The role of β-amyloid peptide in neurodegenerative diseases

Studies of neurodegenerative disorders (NDDs) are drawing more attention of researchers worldwide due to the high incidence of Alzheimer's disease (AD). The pathophysiology of such disorders is, in part, characterized by the transition of a wild-type peptide from its native conformation into a very stable pathological isoform. Subsequently, these abnormal proteins form aggregates of amyloid fibrils that continuously increase in size. Changes in the metabolic processes of neurons (e.g. oxidative stress, hyperphosphorylation of the tau protein, and resulting secondary changes in the cell metabolism) ultimately lead to cell death. We hypothesize that extracellular deposition of β-amyloid peptide fibrils and neurofibrillary tangles represents the body's adaptation mechanism, aimed at preservation of autonomic functioning; while the cognitive decline is severe, the rest of the organ systems remain unaffected and continue to function. This hypothesis is supported by the fact that destruction of pathological plaques, fibrils, and tangles and the use of vaccines targeting β-amyloid result in undesirable side effects. To gain a better understanding of the pathophysiology of Alzheimer's disease and to develop novel therapies, continued studies of the sporadic form of disease and the mechanisms triggering conformational changes in β-amyloid peptide fragments are essential. This review is focused on studies investigating the formation of amyloid fibrils and their role in the pathogenesis of neurodegenerative diseases. In addition, we discuss a related disorder--amyloidosis--where formation of fibrils, tangles, and plaques leads to neuronal death which may occur as a result of a failed adaptation process. Further in-depth investigation and comprehensive analysis of alterations in the metabolism of APP, β-amyloid, and tau protein, which have a pathological effect on cell membrane, alter phosphate exchange, and impair other key metabolic functions of the cell long before the characteristic amyloid deposition takes place, is warranted. A better understanding of intraneuronal processes is crucial in identifying specific inhibitors of pathologic neuronal processes and, consequently, will allow for targeted therapy, thus maximizing efficacy of selected therapeutic regimens.

Characterizing the structural behavior of selected Aβ-42 monomers with different solubilities

The conformational behavior of the wild-type amyloid β-42 (Aβ-42) monomer and two of its mutants was explored via all-atom replica exchange molecular dynamics simulations in explicit solvent, to identify structural features that may promote or deter early-stage oligomerization. The markers used for this purpose indicate that while the three peptides are relatively flexible they have distinct preferential structures and degree of rigidity. In particular, we found that one mutant that remains in the monomeric state in experiments displays a characteristic N-terminal structure that significantly enhances its rigidity. This finding is consistent with various studies that have detected a reduction in oligomerization frequency and Aβ-related toxicity upon sequence-specific antibody or ligand binding to the N-terminal tail of wild-type monomers, likely leading to the stabilization of this region. In general, our results highlight a potential role of the N-terminal segment on Aβ oligomerization and give insights into specific interactions that may be responsible for promoting the pronounced structural changes observed upon introducing point mutations on the wild-type Aβ-42 peptide.

Association thermodynamics and conformational stability of beta-sheet amyloid beta(17-42) oligomers: effects of E22Q (Dutch) mutation and charge neutralization

Amyloid fibrils are associated with many neurodegenerative diseases. It was found that amyloidogenic oligomers, not mature fibrils, are neurotoxic agents related to these diseases. Molecular mechanisms of infectivity, pathways of aggregation, and molecular structure of these oligomers remain elusive. Here, we use all-atom molecular dynamics, molecular mechanics combined with solvation analysis by statistical-mechanical, three-dimensional molecular theory of solvation (also known as 3D-RISM-KH) in a new MM-3D-RISM-KH method to study conformational stability, and association thermodynamics of small wild-type Abeta(17-42) oligomers with different protonation states of Glu(22), as well the E22Q (Dutch) mutants. The association free energy of small beta-sheet oligomers shows near-linear trend with the dimers being thermodynamically more stable relative to the larger constructs. The linear (within statistical uncertainty) dependence of the association free energy on complex size is a consequence of the unilateral stacking of monomers in the beta-sheet oligomers. The charge reduction of the wild-type Abeta(17-42) oligomers upon protonation of the solvent-exposed Glu(22) at acidic conditions results in lowering the association free energy compared to the wild-type oligomers at neutral pH and the E22Q mutants. The neutralization of the peptides because of the E22Q mutation only marginally affects the association free energy, with the reduction of the direct electrostatic interactions mostly compensated by the unfavorable electrostatic solvation effects. For the wild-type oligomers at acidic conditions such compensation is not complete, and the electrostatic interactions, along with the gas-phase nonpolar energetic and the overall entropic effects, contribute to the lowering of the association free energy. The differences in the association thermodynamics between the wild-type Abeta(17-42) oligomers at neutral pH and the Dutch mutants, on the one hand, and the Abeta(17-42) oligomers with protonated Glu(22), on the other, may be explained by destabilization of the inter- and intrapeptide salt bridges between Asp(23) and Lys(28). Peculiarities in the conformational stability and the association thermodynamics for the different models of the Abeta(17-42) oligomers are rationalized based on the analysis of the local physical interactions and the microscopic solvation structure.

Elucidation of amyloid beta-protein oligomerization mechanisms: discrete molecular dynamics study

Oligomers of amyloid beta-protein (Abeta) play a central role in the pathology of Alzheimer's disease. Of the two predominant Abeta alloforms, Abeta(1-40) and Abeta(1-42), Abeta(1-42) is more strongly implicated in the disease. We elucidated the structural characteristics of oligomers of Abeta(1-40) and Abeta(1-42) and their Arctic mutants, [E22G]Abeta(1-40) and [E22G]Abeta(1-42). We simulated oligomer formation using discrete molecular dynamics (DMD) with a four-bead protein model, backbone hydrogen bonding, and residue-specific interactions due to effective hydropathy and charge. For all four peptides under study, we derived the characteristic oligomer size distributions that were in agreement with prior experimental findings. Unlike Abeta(1-40), Abeta(1-42) had a high propensity to form paranuclei (pentameric or hexameric) structures that could self-associate into higher-order oligomers. Neither of the Arctic mutants formed higher-order oligomers, but [E22G]Abeta(1-40) formed paranuclei with a similar propensity to that of Abeta(1-42). Whereas the best agreement with the experimental data was obtained when the charged residues were modeled as solely hydrophilic, further assembly from spherical oligomers into elongated protofibrils was induced by nonzero electrostatic interactions among the charged residues. Structural analysis revealed that the C-terminal region played a dominant role in Abeta(1-42) oligomer formation whereas Abeta(1-40) oligomerization was primarily driven by intermolecular interactions among the central hydrophobic regions. The N-terminal region A2-F4 played a prominent role in Abeta(1-40) oligomerization but did not contribute to the oligomerization of Abeta(1-42) or the Arctic mutants. The oligomer structure of both Arctic peptides resembled Abeta(1-42) more than Abeta(1-40), consistent with their potentially more toxic nature.

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.

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.

Benchmarking implicit solvent folding simulations of the amyloid beta(10-35) fragment

A pathogenetic feature of Alzhemier disease is the aggregation of monomeric beta-amyloid proteins (Abeta) to form oligomers. Usually these oligomers of long peptides aggregate on time scales of microseconds or longer, making computational studies using atomistic molecular dynamics models prohibitively expensive and making it essential to develop computational models that are cheaper and at the same time faithful to physical features of the process. We benchmark the ability of our implicit solvent model to describe equilibrium and dynamic properties of monomeric Abeta(10-35) using all-atom Langevin dynamics (LD) simulations, since Alphabeta(10-35) is the only fragment whose monomeric properties have been measured. The accuracy of the implicit solvent model is tested by comparing its predictions with experiment and with those from a new explicit water MD simulation, (performed using CHARMM and the TIP3P water model) which is approximately 200 times slower than the implicit water simulations. The dependence on force field is investigated by running multiple trajectories for Alphabeta(10-35) using the CHARMM, OPLS-aal, and GS-AMBER94 force fields, whereas the convergence to equilibrium is tested for each force field by beginning separate trajectories from the native NMR structure, a completely stretched structure, and from unfolded initial structures. The NMR order parameter, S2, is computed for each trajectory and is compared with experimental data to assess the best choice for treating aggregates of Alphabeta. The computed order parameters vary significantly with force field. Explicit and implicit solvent simulations using the CHARMM force fields display excellent agreement with each other and once again support the accuracy of the implicit solvent model. Alphabeta(10-35) exhibits great flexibility, consistent with experiment data for the monomer in solution, while maintaining a general strand-loop-strand motif with a solvent-exposed hydrophobic patch that is believed to be important for aggregation. Finally, equilibration of the peptide structure requires an implicit solvent LD simulation as long as 30 ns.

Structures and free-energy landscapes of the wild type and mutants of the Abeta(21-30) peptide are determined by an interplay between intrapeptide electrostatic and hydrophobic interactions

The initial events in protein aggregation involve fluctuations that populate monomer conformations, which lead to oligomerization and fibril assembly. The highly populated structures, driven by a balance between hydrophobic and electrostatic interactions in the protease-resistant wild-type Abeta(21-30) peptide and mutants E22Q (Dutch), D23N (Iowa), and K28N, are analyzed using molecular dynamics simulations. Intrapeptide electrostatic interactions were connected to calculated pK(a) values that compare well with the experimental estimates. The pK(a) values of the titratable residues show that E22 and D23 side chains form salt bridges only infrequently with the K28 side chain. Contacts between E22-K28 are more probable in "dried" salt bridges, whereas D23-K28 contacts are more probable in solvated salt bridges. The strength of the intrapeptide hydrophobic interactions increases as D23N<WT<E22Q<K28A. Free-energy profiles and disconnectivity representation of the energy landscapes show that the monomer structures partition into four distinct basins. The hydrophobic interactions cluster the Abeta(21-30) peptide into two basins, differentiated by the relative position of the DVG(23-25) and GSN(25-27) fragments about the G25 residue. The E22Q mutation increases the population with intact VGSN turn compared to the wild-type (WT) peptide. The increase in the population of the structures in the aggregation-prone Basin I in E22Q, which occurs solely due to the difference in charge states between the Dutch mutant and the WT, gives a structural explanation of the somewhat larger aggregation rate in the mutant. The D23N mutation dramatically reduces the intrapeptide interactions. The K28A mutation increases the intrapeptide hydrophobic interactions that promote population of structures in Basin I and Basin II whose structures are characterized by hydrophobic interaction between V24 and K28 side chains but with well-separated ends of the backbone atoms in the VGSN turn. The intrapeptide electrostatic interactions in the WT and E22Q peptides roughen the free-energy surface compared to the K28A peptide. The D23N mutation has a flat free-energy surface, corresponding to an increased population of random coil-like structures with weak hydrophobic and electrostatic interactions. We propose that mutations or sequences that enhance the probability of occupying Basin I would promote aggregation of Abeta peptides.

Familial Alzheimer's disease mutations alter the stability of the amyloid beta-protein monomer folding nucleus

Amyloid beta-protein (Abeta) oligomers may be the proximate neurotoxins in Alzheimer's disease (AD). Recently, to elucidate the oligomerization pathway, we studied Abeta monomer folding and identified a decapeptide segment of Abeta, (21)Ala-(22)Glu-(23)Asp-(24)Val-(25)Gly-(26)Ser-(27)Asn-(28)Lys-(29)Gly-(30)Ala, within which turn formation appears to nucleate monomer folding. The turn is stabilized by hydrophobic interactions between Val-24 and Lys-28 and by long-range electrostatic interactions between Lys-28 and either Glu-22 or Asp-23. We hypothesized that turn destabilization might explain the effects of amino acid substitutions at Glu-22 and Asp-23 that cause familial forms of AD and cerebral amyloid angiopathy. To test this hypothesis, limited proteolysis, mass spectrometry, and solution-state NMR spectroscopy were used here to determine and compare the structure and stability of the Abeta(21-30) turn within wild-type Abeta and seven clinically relevant homologues. In addition, we determined the relative differences in folding free energies (DeltaDeltaG(f)) among the mutant peptides. We observed that all of the disease-associated amino acid substitutions at Glu-22 or Asp-23 destabilized the turn and that the magnitude of the destabilization correlated with oligomerization propensity. The Ala21Gly (Flemish) substitution, outside the turn proper (Glu-22-Lys-28), displayed a stability similar to that of the wild-type peptide. The implications of these findings for understanding Abeta monomer folding and disease causation are discussed.

Computational simulations of the early steps of protein aggregation

There is strong evidence that the oligomers of key proteins, formed during the early steps of aggregation, could be the primary toxic species associated with human neuro-degenerative diseases, such as Alzheimer's and prion diseases. Here, we review recent progress in the development of computational approaches in order to understand the structures, dynamics and free energy surfaces of oligomers. We also discuss possible research directions for the coming years.

Effect of single-point sequence alterations on the aggregation propensity of a model protein

Sequences of contemporary proteins are believed to have evolved through a process that optimized their overall fitness, including their resistance to deleterious aggregation. Biotechnological processing may expose therapeutic proteins to conditions that are much more conducive to aggregation than those encountered in a cellular environment. An important task of protein engineering is to identify alternative sequences that would protect proteins when processed at high concentrations without altering their native structure associated with specific biological function. Our computational studies exploit parallel tempering simulations of coarse-grained model proteins to demonstrate that isolated amino acid residue substitutions can result in significant changes in the aggregation resistance of the protein in a crowded environment while retaining protein structure in isolation. A thermodynamic analysis of protein clusters subject to competing processes of folding and association shows that moderate mutations can produce effects similar to those caused by changes in system conditions, including temperature, concentration, and solvent composition, that affect the aggregation propensity. The range of conditions where a protein can resist aggregation can therefore be tuned by sequence alterations, although the protein generally may retain its generic ability for aggregation.

Spontaneous fibril formation by polyalanines; discontinuous molecular dynamics simulations

Fibrillary protein aggregates rich in beta-sheet structure have been implicated in the pathology of several neurodegenerative diseases. In this work, we investigate the formation of fibrils by performing discontinuous molecular dynamics simulations on systems containing 12 to 96 model Ac-KA(14)K-NH(2) peptides using our newly developed off-lattice, implicit-solvent, intermediate-resolution model, PRIME. We find that, at a low concentration, random-coil peptides assemble into alpha-helices at low temperatures. At intermediate concentrations, random-coil peptides assemble into alpha-helices at low temperatures and large beta-sheet structures at high temperatures. At high concentrations, the system forms beta-sheets over a wide range of temperatures. These assemble into fibrils above a critical temperature which decreases with concentration and exceeds the isolated peptide's folding temperature. At very high temperatures and all concentrations, the system is in a random-coil state. All of these results are in good qualitative agreement with those by Blondelle and co-workers on Ac-KA(14)K-NH(2) peptides. The fibrils observed in our simulations mimic the structural characteristics observed in experiments in terms of the number of sheets formed, the values of the intra- and intersheet separations, and the parallel peptide arrangement within each beta-sheet. Finally, we find that when the strength of the hydrophobic interaction between nonpolar side chains is high compared to the strength of hydrogen bonding, amorphous aggregates, rather than fibrillar aggregates, are formed.

Molecular simulation of protein aggregation

Computer simulation offers unique possibilities for investigating molecular-level phenomena difficult to probe experimentally. Drawing from a wealth of studies concerning protein folding, computational studies of protein aggregation are emerging. These studies have been successful in capturing aspects of aggregation known from experiment and are being used to refine experimental methods aimed at abating aggregation. Here we review molecular-simulation studies of protein aggregation conducted in our laboratory. Specific attention is devoted to issues with implications for biotechnology.

Solvent and mutation effects on the nucleation of amyloid beta-protein folding

Experimental evidence suggests that the folding and aggregation of the amyloid beta-protein (Abeta) into oligomers is a key pathogenetic event in Alzheimer's disease. Inhibiting the pathologic folding and oligomerization of Abeta could be effective in the prevention and treatment of Alzheimer's disease. Here, using all-atom molecular dynamics simulations in explicit solvent, we probe the initial stages of folding of a decapeptide segment of Abeta, Abeta(21-30), shown experimentally to nucleate the folding process. In addition, we examine the folding of a homologous decapeptide containing an amino acid substitution linked to hereditary cerebral hemorrhage with amyloidosis-Dutch type, [Gln-22]Abeta(21-30). We find that: (i) when the decapeptide is in water, hydrophobic interactions and transient salt bridges between Lys-28 and either Glu-22 or Asp-23 are important in the formation of a loop in the Val-24-Lys-28 region of the wild-type decapeptide; (ii) in the presence of salt ions, salt bridges play a more prominent role in the stabilization of the loop; (iii) in water with a reduced density, the decapeptide forms a helix, indicating the sensitivity of folding to different aqueous environments; and (iv) the "Dutch" peptide in water, in contrast to the wild-type peptide, fails to form a long-lived Val-24-Lys-28 loop, suggesting that loop stability is a critical factor in determining whether Abeta folds into pathologic structures.

Insulin and cholesterol pathways in neuronal function, memory and neurodegeneration

Insulin and cholesterol play important roles in basic metabolic processes in peripheral tissues. Both insulin and cholesterol can also act as signalling molecules in the central nervous system that participate in neuronal function, memory and neurodegenerative diseases. A high-cholesterol diet improves spatial memory in experimental animals. β-Amyloid, the toxic peptide in neurons of AD (Alzheimer's disease) patients, binds cholesterol and catalyses its oxidation to 7β-hydroxycholesterol, a highly toxic oxysterol that is a potent inhibitor of α-PKC (α-protein kinase C), an enzyme critical in memory consolidation and synaptic plasticity and implicated in AD. Oxidized cholesterol also can act as a second messenger for insulin. Oxidized low-density lipoprotein inhibits insulin-dependent phosphorylation of the signalling kinases ERK (extracellular-signal-regulated kinase) and PKB/Akt. In sporadic AD patients, insulin levels are decreased, suggesting links between AD and diabetes. Insulin signalling is also important in synaptic plasticity. Insulin receptors are up-regulated and undergo translocation after spatial learning. Insulin modulates the activity of excitatory and inhibitory receptors including the glutamate and γ-aminobutyric acid receptors and activates two biochemical pathways: the shc-ras-mitogen-activated protein kinase pathway and the PI3K (phosphoinositide 3-kinase)/PKC pathway, both of which are involved in memory processing. These findings point to a convergence at the biochemical level between pathways involved in AD and those important for normal memory.

A Strand-Loop-Strand Structure Is a Possible Intermediate in Fibril Elongation: Long Time Simulations of Amyloid-β Peptide (10−35)

A total of 6.2 micros molecular dynamics simulations of amyloid-beta (10-35) (Abeta) were performed in explicit water solvent. The results reveal that the collapsed-coil (cc) structure determined by experiments is stable at pH 5.6 for hundreds of nanoseconds, but it can exchange with a strand-loop-strand (SLS) structure on the microsecond time scale. The SLS structure has D23-K28 as a reverse loop and the central hydrophobic core and the C-terminal in hydrophobic contact. This SLS structure topologically resembles the proposed monomer conformation in fibrils. Since it has been suggested that a special conformation of Abeta is needed when the monomer binds to fibril ends to elongate fibrils, we propose that the SLS structure may be an important intermediate binding structure for Abeta fibril growth. Simulations at pH 2.0, which is used to mimic the mutation of E22Q and D23N, and at high temperature (400 K) indicate that the SLS structure is considerably populated under these conditions while the cc structure is disrupted. These results imply that the SLS structures may also be a binding intermediate in other conditions such as E22Q and/or D23N mutations and high temperature, which have been proved to promote fibril formation previously.

Conformational transition of amyloid beta-peptide

The amyloid beta-peptides (Abetas), containing 39-43 residues, are the key protein components of amyloid deposits in Alzheimer's disease. To structurally characterize the dynamic behavior of Abeta(40), 12 independent long-time molecular dynamics (MD) simulations for a total of 850 ns were performed on both the wide-type peptide and its mutant in both aqueous solution and a biomembrane environment. In aqueous solution, an alpha-helix to beta-sheet conformational transition for Abeta(40) was observed, and an entire unfolding process from helix to coil was traced by MD simulation. Structures with beta-sheet components were observed as intermediates in the unfolding pathway of Abeta(40). Four glycines (G(25), G(29), G(33), and G(37)) are important for Abeta(40) to form beta-sheet in aqueous solution; mutations of these glycines to alanines almost abolished the beta-sheet formation and increased the content of the helix component. In the dipalmitoyl phosphatidylcholine (DPPC) bilayer, the major secondary structure of Abeta(40) is a helix; however, the peptide tends to exit the membrane environment and lie down on the surface of the bilayer. The dynamic feature revealed by our MD simulations rationalized several experimental observations for Abeta(40) aggregation and amyloid fibril formation. The results of MD simulations are beneficial to understanding the mechanism of amyloid formation and designing the compounds for inhibiting the aggregation of Abeta and amyloid fibril formation.

Probing the initial stage of aggregation of the Abeta(10-35)-protein: assessing the propensity for peptide dimerization

Characterization of the early stages of peptide aggregation is of fundamental importance in elucidating the mechanism of the formation of deposits associated with amyloid disease. The initial step in the pathway of aggregation of the Abeta-protein, whose monomeric NMR structure is known, was studied through the simulation of the structure and stability of the peptide dimer in aqueous solution. A protocol based on shape complementarity was used to generate an assortment of possible dimer structures. The structures generated based on shape complementarity were evaluated using rapidly computed estimates of the desolvation and electrostatic interaction energies to identify a putative stable dimer structure. The potential of mean force associated with the dimerization of the peptides in aqueous solution was computed for both the hydrophobic and the electrostatic driven forces using umbrella sampling and classical molecular dynamics simulation at constant temperature and pressure with explicit solvent and periodic boundary conditions. The comparison of the two free energy profiles suggests that the structure of the peptide dimer is determined by the favorable desolvation of the hydrophobic residues at the interface. Molecular dynamics trajectories originating from two putative dimer structures indicate that the peptide dimer is stabilized primarily through hydrophobic interactions, while the conformations of the peptide monomers undergo substantial structural reorganization in the dimerization process. The finding that the phi-dimer may constitute the ensemble of stable Abeta(10-35) dimer has important implications for fibril formation. In particular, the expulsion of water molecules at the interface might be a key event, just as in the oligomerization of Abeta(16-22) fragments. We conjecture that events prior to the nucleation process themselves might involve crossing free energy barriers which depend on the peptide-peptide and peptide-water interactions. Consistent with existing experimental studies, the peptides within the ensemble of aggregated states show no signs of formation of secondary structure.

Kinetics of fibril formation by polyalanine peptides

Ordered beta-sheet complexes, termed amyloid fibrils, are the underlying structural components of the intra- and extracellular fibrillar protein deposits that are associated with a variety of human diseases, including Alzheimer's, Parkinson's, and the prion diseases. In this work, we investigated the kinetics of fibril formation using our newly developed off-lattice intermediate resolution model, PRIME. The model is simple enough to allow the treatment of large multichain systems while maintaining a fairly realistic description of protein dynamics without built-in bias toward any conformation when used in conjunction with constant temperature discontinuous molecular dynamics, a fast alternative to conventional molecular dynamics. Simulations were performed on systems containing 48-96 model Ac-KA14K-NH2 peptides. We found that fibril formation for polyalanines incorporate features that are characteristic of three models, the templated assembly, nucleated polymerization, and nucleated conformational conversion models, but that none of them gave a completely satisfactory description of the simulation kinetics. Fibril formation was nucleation-dependent, occurring after a lag time that decreased with increasing peptide concentration and increased with increasing temperature. Fibril formation appeared to be a conformational conversion process in which small amorphous aggregates --> beta-sheets --> ordered nucleus --> subsequent rapid growth of a small stable fibril or protofilament. Fibril growth in our simulations involved both beta-sheet elongation, in which the fibril grew by adding individual peptides to the end of each beta-sheet, and lateral addition, in which the fibril grew by adding already formed beta-sheets to its side. The initial rate of fibril formation increased with increasing concentration and decreased with increasing temperature.

Oxidation of cholesterol by amyloid precursor protein and beta-amyloid peptide

Alzheimer's disease (AD) is characterized by accumulation of the neurotoxic peptide beta-amyloid, which is produced by proteolysis of amyloid precursor protein (APP). APP is a large membrane-bound copper-binding protein that is essential in maintaining synaptic function and may play a role in synaptogenesis. beta-Amyloid has been shown to contribute to the oxidative stress that accompanies AD. Later stages of AD are characterized by neuronal apoptosis. However, the biochemical function of APP and the mechanism of the toxicity of beta-amyloid are still unclear. In this study, we show that both beta-amyloid and APP can oxidize cholesterol to form 7beta-hydroxycholesterol, a proapoptotic oxysterol that was neurotoxic at nanomolar concentrations. 7beta-Hydroxycholesterol inhibited secretion of soluble APP from cultured rat hippocampal H19-7/IGF-IR neuronal cells and inhibited tumor necrosis factor-alpha-converting enzyme alpha-secretase activity but had no effect on beta-site APP-cleaving enzyme 1 activity. 7beta-Hydroxycholesterol was also a potent inhibitor of alpha-protein kinase C, with a K(i) of approximately 0.2 nm. The rate of reaction between cholesterol and beta-amyloid was comparable to the rates of cholesterol-metabolizing enzymes (k(cat) = 0.211 min(-)1). The rate of production of 7beta-hydroxycholesterol by APP was approximately 200 times lower than by beta-amyloid. Oxidation of cholesterol was accompanied by stoichiometric production of hydrogen peroxide and required divalent copper. The results suggest that a function of APP may be to produce low levels of 7-hydroxycholesterol. Higher levels produced by beta-amyloid could contribute to the oxidative stress and cell loss observed in Alzheimer's disease.

Phase diagrams describing fibrillization by polyalanine peptides

Amyloid fibrils are the structural components underlying the intra- and extracellular protein deposits that are associated with a variety of human diseases, including Alzheimer's, Parkinson's, and the prion diseases. In this work, we examine the thermodynamics of fibril formation using our newly-developed off-lattice intermediate-resolution protein model, PRIME. The model is simple enough to allow the treatment of large multichain systems while maintaining a fairly realistic description of protein dynamics when used in conjunction with constant-temperature discontinuous molecular dynamics, a fast alternative to conventional molecular dynamics. We conduct equilibrium simulations on systems containing 96 Ac-KA14K-NH2 peptides over a wide range of temperatures and peptide concentrations using the replica-exchange method. Based on measured values of the heat capacity, radius of gyration, and percentage of peptides that form the various structures, a phase diagram in the temperature-concentration plane is constructed delineating the regions where each structure is stable. There are four distinct single-phase regions: alpha-helices, fibrils, nonfibrillar beta-sheets, and random coils; and four two-phase regions: random coils/nonfibrillar beta-sheets, random coils/fibrils, fibrils/nonfibrillar beta-sheets, and alpha-helices/nonfibrillar beta-sheets. The alpha-helical region is at low temperature and low concentration. The nonfibrillar beta-sheet region is at intermediate temperatures and low concentrations and expands to higher temperatures as concentration is increased. The fibril region occurs at intermediate temperatures and intermediate concentrations and expands to lower as the peptide concentration is increased. The random-coil region is at high temperatures and all concentrations; this region shifts to higher temperatures as the concentration is increased.

Differential degradation of amyloid beta genetic variants associated with hereditary dementia or stroke by insulin-degrading enzyme

Inherited amino acid substitutions at position 21, 22, or 23 of amyloid beta (Abeta) lead to presenile dementia or stroke. Insulin-degrading enzyme (IDE) can hydrolyze Abeta wild type, yet whether IDE is capable of degrading Abeta bearing pathogenic substitutions is not known. We studied the degradation of all of the published Abeta genetic variants by recombinant rat IDE (rIDE). Monomeric Abeta wild type, Flemish (A21G), Italian (E22K), and Iowa (D23N) variants were readily degraded by rIDE with a similar efficiency. However, proteolysis of Abeta Dutch (E22Q) and Arctic (E22G) was significantly lower as compared with Abeta wild type and the rest of the mutant peptides. In the case of Abeta Dutch, inefficient proteolysis was related to a high content of beta structure as assessed by circular dichroism. All of the Abeta variants were cleaved at Glu3-Phe4 and Phe4-Arg5 in addition to the previously described major sites within positions 13-15 and 18-21. SDS-stable Abeta dimers were highly resistant to proteolysis by rIDE regardless of the variant, suggesting that IDE recognizes a conformation that is available for interaction only in monomeric Abeta. These results raise the possibility that upregulation of IDE may promote the clearance of soluble Abeta in hereditary forms of Abeta diseases.

Structural and dynamical analysis of the hydration of the Alzheimer's beta-amyloid peptide

An analysis of the water molecules in the first solvation shell obtained from the molecular dynamics simulation of the amyloid beta(10-35)NH2 peptide and the amyloid beta(10-35)NH2E22Q "Dutch" mutant peptide is presented. The structure, energetics, and dynamics of water in the hydration shell have been investigated using a variety of measures, including the hydrogen bond network, the water residence times for all the peptide residues, the diffusion constant, experimentally determined HN amide proton exchange, and the transition probabilities for water to move from one residue to another or into the bulk. The results of the study indicate that: (1) the water molecules at the peptide-solvent interface are organized in an ordered structure similar for the two peptide systems but different from that of the bulk, (2) the peptide structure inhibits diffusion perpendicular to the peptide surface by a factor of 3 to 5 relative to diffusion parallel to the peptide surface, which is comparable to diffusion of bulk water, (3) water in the first solvation shell shows dynamical relaxation on fast (1-2 ps) and slow (10-40 ps) time scales, (4) a novel solvent relaxation master equation is shown to capture the details of the fast relaxation of water in the peptide's first solvation shell, (5) the interaction between the peptide and the solvent is stronger in the wild type than in the E22Q mutant peptide, in agreement with earlier results obtained from computer simulations [Massi, F.; Straub, J. E. Biophys J 2001, 81, 697] correlated with the observed enhanced activity of the E22Q mutant peptide.

Charge states rather than propensity for beta-structure determine enhanced fibrillogenesis in wild-type Alzheimer's beta-amyloid peptide compared to E22Q Dutch mutant

The activity of the Alzheimer's amyloid beta-peptide is a sensitive function of the peptide's sequence. Increased fibril elongation rate of the E22Q Dutch mutant of the Alzheimer's amyloid beta-peptide relative to that of the wild-type peptide has been observed. The increased activity has been attributed to a larger propensity for the formation of beta structure in the monomeric E22Q mutant peptide in solution relative to the WT peptide. That hypothesis is tested using four nanosecond timescale simulations of the WT and Dutch mutant forms of the Abeta(10-35)-peptide in aqueous solution. The simulation results indicate that the propensity for formation of beta-structure is no greater in the E22Q mutant peptide than in the WT peptide. A significant measure of "flickering" of helical structure in the central hydrophobic cluster region of both the WT and mutant peptides is observed. The simulation results argue against the hypothesis that the Dutch mutation leads to a higher probability of formation of beta-structure in the monomeric peptide in aqueous solution. We propose that the greater stability of the solvated WT peptide relative to the E22Q mutant peptide leads to decreased fibril elongation rate in the former. Stability difference is due to the differing charge state of the two peptides. The other proposal leads to the prediction that the fibril elongation rates for the WT and the mutant E22Q should be similar under acid conditions.