Discrete molecular dynamics study of oligomer formation by N-terminally truncated amyloid β-protein

Insights into the Mixture of Aβ24 and Aβ42 Peptides from Atomistic Simulations

Amyloid-β (Aβ) oligomers play a central role in Alzheimer's disease (AD). Plaques of AD patients consist of Aβ40 and Aβ42 peptides and truncated Aβ peptides. The Aβ24 peptide, identified in human AD brains, was found to impair Aβ42 clearance through the brain-blood barrier. The Aβ24 peptide was also shown to reduce Aβ42 aggregation kinetics in pure buffer, but the underlying mechanism is unknown at atomistic level. In this study, we explored the conformational ensemble of the equimolar mixture of Aβ24 and Aβ42 by replica exchange molecular dynamics simulations and compared it to our previous results on the pure Aβ42 dimer. Our simulations demonstrate that the truncation at residue 24 changes the secondary, tertiary, and quaternary structures of the dimer, offering an explanation of the slower aggregation kinetics of the mixture.

Binary Matrix Method to Enumerate, Hierarchically Order, and Structurally Classify Peptide Aggregation

Protein aggregation is a common and complex phenomenon in biological processes, yet a robust analysis of this aggregation process remains elusive. The commonly used methods such as center-of-mass to center-of-mass (COM-COM) distance, the radius of gyration (R_g), hydrogen bonding (HB), and solvent accessible surface area do not quantify the aggregation accurately. Herein, a new and robust method that uses an aggregation matrix (AM) approach to investigate peptide aggregation in a MD simulation trajectory is presented. An nxn two-dimensional AM is created by using the interpeptide C_α-C_α cutoff distances, which are binarily encoded (0 or 1). These aggregation matrices are analyzed to enumerate, hierarchically order, and structurally classify the aggregates. Comparison of the present AM method suggests that it is superior to the HB method since it can incorporate nonspecific interactions and the R_g and COM-COM methods since the cutoff distance is independent of the length of the peptide. More importantly, the present method can structurally classify the peptide aggregates, which the conventional R_g, COM-COM, and HB methods fail to do. The unique selling point of this method is its ability to structurally classify peptide aggregates using two-dimensional matrices.

Biological and methodological complexities of beta-amyloid peptide: Implications for Alzheimer's disease research

Although controversial, the amyloid cascade hypothesis remains central to the Alzheimer's disease (AD) field and posits amyloid-beta (Aβ) as the central factor initiating disease onset. In recent years, there has been an increase in emphasis on studying the role of low molecular weight aggregates, such as oligomers, which are suggested to be more neurotoxic than fibrillary Aβ. Other Aβ isoforms, such as truncated Aβ, have also been implicated in disease. However, developing a clear understanding of AD pathogenesis has been hampered by the complexity of Aβ biochemistry in vitro and in vivo. This review explores factors contributing to the lack of consistency in experimental approaches taken to model Aβ aggregation and toxicity and provides an overview of the different techniques available to analyse Aβ, such as electron and atomic force microscopy, nuclear magnetic resonance spectroscopy, dye-based assays, size exclusion chromatography, mass spectrometry and SDS-PAGE. The review also explores how different types of Aβ can influence Aβ aggregation and toxicity, leading to variation in experimental outcomes, further highlighting the need for standardisation in Aβ preparations and methods used in current research.

Soluble State of Villin Headpiece Protein as a Tool in the Assessment of MD Force Fields

Protein self-assembly plays an important role in cellular processes. Whereas molecular dynamics (MD) represents a powerful tool in studying assembly mechanisms, its predictions depend on the accuracy of underlying force fields, which are known to overly promote protein assembly. We here examine villin headpiece domain, HP36, which remains soluble at concentrations amenable to MD studies. The experimental characterization of soluble HP36 at concentrations of 0.05 to 1 mM reveals concentration-independent 90% monomeric and 10% dimeric populations. Extensive all-atom MD simulations at two protein concentrations, 0.9 and 8.5 mM, probe the HP36 dimer population, stability, and kinetics of dimer formation within two MD force fields, Amber ff14SB and CHARMM36m. MD results demonstrate that whereas CHARMM36m captures experimental HP36 monomer populations at the lower concentration, both force fields overly promote HP36 association at the higher concentration. Moreover, contacts stabilizing HP36 dimers are force-field-dependent. CHARMM36m produces consistently higher HP36 monomer populations, lower association rates, and weaker dependence of these quantities on the protein concentration than Amber ff14SB. Nonetheless, the highest monomer populations and dissociation constants are observed when the TIP3P water model in Amber ff14SB is replaced by TIP4P/2005, showcasing the critical role of the water model in addressing the protein solubility problem in MD.

C-Terminal Plays as the Possible Nucleation of the Self-Aggregation of the S-Shape Aβ11-42 Tetramer in Solution: Intensive MD Study

Amyloid beta (Aβ) peptides are characterized as the major factors associated with neuron death in Alzheimer's disease, which is listed as the most common form of neurodegeneration. Disordered Aβ peptides are released from proteolysis of the amyloid precursor protein. The Aβ self-assembly process roughly takes place via five steps: disordered forms → oligomers → photofibrils → mature fibrils → plaques. Although Aβ fibrils are often observed in patient brains, oligomers were recently indicated to be major neurotoxic elements. In this work, the neurotoxic compound S-shape Aβ11-42 tetramer (S4Aβ11-42) was investigated over 10 μs of unbiased MD simulations. In particular, the S4Aβ11-42 oligomer adopted a high dynamics structure, resulting in unsuccessful determination of their structures in experiments. The C-terminal was suggested as the possible nucleation of the Aβ42 aggregation. The sequences 27-35 and 39-40 formed rich β-content, whereas other residues mostly adopted coil structures. The mean value of the β-content over the equilibrium interval is ∼42 ± 3%. Furthermore, the dissociation free energy of the S4Aβ11-42 peptide was predicted using a biased sampling method. The obtained free energy is ΔG US = -58.44 kcal/mol which is roughly the same level as the corresponding value of the U-shape Aβ17-42 peptide. We anticipate that the obtained S4Aβ11-42 structures could be used as targets for AD inhibitor screening over the in silico study.

Seeding and Cross-Seeding Aggregations of Aβ40 and Its N-Terminal-Truncated Peptide Aβ11-40

In the amyloid plaques of Alzheimer's disease (AD) patients, a large number of N-terminal-truncated amyloid β (Aβ) peptides such as Aβ_11-40 have been identified in addition to the full-length Aβ peptides. However, little is known about the roles of the N-terminal-truncated peptides in AD pathological process. Herein, seeding and cross-seeding aggregations of Aβ₄₀ and its N-terminal-truncated Aβ_11-40 were investigated in the solution and on the surfaces of chips with immobilized seeds by extensive biophysical and biological analyses. The results showed that Aβ₄₀ and Aβ_11-40 aggregates could seed both homologous and heterologous aggregations of the two monomers. However, the capability and characteristics of the seeding (homologous aggregation) and cross-seeding (heterologous aggregation) were significantly different. Aβ₄₀ seeds showed stronger acceleration effects on the aggregations than Aβ_11-40 seeds and induced β-sheet-rich fibrous aggregates of similar cytotoxicities for the two monomers. This indicates that Aβ₄₀ and Aβ_11-40 had similar aggregation pathways in the seeding and cross-seeding on Aβ₄₀ seeds. By contrast, Aβ_11-40 seeds led to different aggregation pathways of Aβ₄₀ and Aβ_11-40. Pure Aβ_11-40 aggregates had higher toxicity than Aβ₄₀ aggregates, and as seeds, Aβ_11-40 seeds induced Aβ₄₀ to form aggregates of higher cytotoxicity. However, homologous Aβ_11-40 aggregates induced by Aβ_11-40 seeds showed lower cytotoxicity than pure Aβ_11-40 aggregates. The results suggest that Aβ_11-40 plays an important role in the pathological process of AD.

Antemortem-Postmortem Correlation of Florbetapir (18F) PET Amyloid Imaging with Quantitative Biochemical Measures of Aβ42 but not Aβ40

Amyloid imaging demonstrates the in vivo presence of amyloid-β (Aβ) deposits in the aging human brain but it is still unknown which structural forms and modifications of Aβ are detected. In Alzheimer's disease, most amyloid deposits are predominantly composed of Aβ ending at amino acid residues Val40 or Ala42. It has been reported that Aβ40 is largely restricted to neuritic plaques while Aβ42 may be deposited in amyloid plaques of all types, and is often the sole component of diffuse plaques. The distinction is important as it is mainly the neuritic plaques that correlate with cognitive impairment while diffuse plaques may be the initial type of Aβ deposited. Whether PET amyloid ligands such as florbetapir-18F (Amyvid) are partially or wholly selective for brain deposits of Aβ40 or Aβ42 is currently unknown. We compared antemortem florbetapir PET cortical/cerebellar signal intensity (SUVr) of 55 subjects with postmortem biochemical (ELISA) measurements employing specific antibodies against Aβ40 and Aβ42. Spearman's univariable correlations were significant for both Aβ40 and Aβ42, but were much stronger for Aβ42. Multiple linear regression showed significance only for Aβ42. These results suggest that florbetapir binds only weakly, if at all, to Aβ40. This may be in part due to the higher likelihood for Aβ42 to be present in a β-pleated sheet tertiary structure, or to differences between Aβ40 and Aβ42 in β-pleated sheet tertiary or quaternary structure.

Islet Amyloid Polypeptide Promotes Amyloid-Beta Aggregation by Binding-Induced Helix-Unfolding of the Amyloidogenic Core

Amyloid aggregation of amyloid-beta (Aβ) and islet amyloid polypeptide (IAPP) is associated with Alzheimer's disease (AD) and type-2 diabetes (T2D), respectively. With T2D being the risk factor for AD and the ability of IAPP to cross the blood-brain barrier, the coaggregation of Aβ and IAPP has been explored to understand the cross-talk between the two diseases. Recent studies demonstrated that soluble IAPP could significantly accelerate the aggregation of Aβ while preformed amyloids of IAPP were poor "seeds" for Aβ aggregation. Here, we apply all-atom discrete molecular dynamics simulations to investigate possible molecular mechanisms for the accelerated coaggregation of IAPP and Aβ42 comparing to Aβ42 aggregation alone, which was confirmed by the complementary thioflavin-T fluorescence assay. Our simulation results suggest that peptides in the mixture tend to form heterodimers as the first step toward their coaggregation. Strong interpeptide interactions with IAPP in the heterodimer shift the helical conformation of Aβ42 in its amyloidogenic central hydrophobic core, residues 16-22 (Aβ16-22), to the extended conformation ready to form β-sheets. Our study suggests that the unfolding of Aβ16-22 helix contributes to the aggregation free-energy barrier and corresponds to the rate-limiting conformational change for Aβ42 aggregation. Therefore, we propose that soluble IAPP promotes the aggregation of Aβ42 by binding-induced conformational change of Aβ42 in its amyloidogenic core and thus reduced aggregation free-energy barrier.

Initial Structural Models of the Aβ42 Dimer from Replica Exchange Molecular Dynamics Simulations

Experimental characterization of the molecular structure of small amyloid (A)β oligomers that are currently considered as toxic agents in Alzheimer's disease is a formidably difficult task due to their transient nature and tendency to aggregate. Such structural information is of importance because it can help in developing diagnostics and an effective therapy for the disease. In this study, molecular simulations and protein-protein docking are employed to explore a possible connection between the structure of Aβ monomers and the properties of the intermonomer interface in the Aβ42 dimer. A structurally diverse ensemble of conformations of the monomer was sampled in microsecond timescale implicit solvent replica exchange molecular dynamics simulations. Representative structures with different solvent exposure of hydrophobic residues and secondary structure content were selected to build structural models of the dimer. Analysis of these models reveals that formation of an intramonomer salt bridge (SB) between Asp23 and Lys28 residues can prevent the building of a hydrophobic interface between the central hydrophobic clusters (CHCs) of monomers upon dimerization. This structural feature of the Aβ42 dimer is related to the difference in packing of hydrophobic residues in monomers with the Asp23-Lys28 SB in on and off states, in particular, to a lower propensity to form hydrophobic contacts between the CHC domain and C-terminal residues in monomers with a formed SB. These findings could have important implications for understanding the difference between aggregation pathways of Aβ monomers leading to neurotoxic oligomers or inert fibrillar structures.

Microchip Circulation Drastically Accelerates Amyloid Aggregation of 1-42 β-amyloid Peptide from Felis catus

The amyloid aggregation process of amyloid β_1-42 peptide is responsible for Alzheimer's disease, affecting millions of elderly people worldwide. Although there has been a great deal of attention directed toward tackling this disease, still no medicine has been found for this fatal disorder. To address this challenge, it is vital to thoroughly understand the molecular mechanism underlying the amyloid peptide aggregation process, as well as seek substances that could hamper this aggregation. In order to shed light on mechanisms leading to amyloidogenesis, we employed a microfluidic system to determine the possible influence of in vivo-like flow in the microchip channel itself on feline Aβ_1-42 peptide amyloidogenesis. We have shown that shear forces occurring during such flow immensely accelerated peptide aggregation. We also tested the inhibitory influence of 3,3'-[1,6-(2,5-dioxahexane)]bis(1-dodecylimidazolium) dichloride gemini surfactant on peptide amyloidogenesis. Our results suggest that this surfactant may inhibit amyloid β_1-42 fibril formation.

Molecular Dynamics Simulations of Amyloid β-Peptide (1-42): Tetramer Formation and Membrane Interactions

The aggregation cascade and peptide-membrane interactions of the amyloid β-peptide (Aβ) have been implicated as toxic events in the development and progression of Alzheimer's disease. Aβ42 forms oligomers and ultimately plaques, and it has been hypothesized that these oligomeric species are the main toxic species contributing to neuronal cell death. To better understand oligomerization events and subsequent oligomer-membrane interactions of Aβ42, we performed atomistic molecular-dynamics (MD) simulations to characterize both interpeptide interactions and perturbation of model membranes by the peptides. MD simulations were utilized to first show the formation of a tetramer unit by four separate Aβ42 peptides. Aβ42 tetramers adopted an oblate ellipsoid shape and showed a significant increase in β-strand formation in the final tetramer unit relative to the monomers, indicative of on-pathway events for fibril formation. The Aβ42 tetramer unit that formed in the initial simulations was used in subsequent MD simulations in the presence of a pure POPC or cholesterol-rich raft model membrane. Tetramer-membrane simulations resulted in elongation of the tetramer in the presence of both model membranes, with tetramer-raft interactions giving rise to the rearrangement of key hydrophobic regions in the tetramer and the formation of a more rod-like structure indicative of a fibril-seeding aggregate. Membrane perturbation by the tetramer was manifested in the form of more ordered, rigid membranes, with the pure POPC being affected to a greater extent than the raft membrane. These results provide critical atomistic insight into the aggregation pathway of Aβ42 and a putative toxic mechanism in the pathogenesis of Alzheimer's disease.

Unmodified and pyroglutamylated amyloid β peptides form hypertoxic hetero-oligomers of unique secondary structure

Amyloid β (Aβ) peptide plays a major role in Alzheimer's disease (AD) and occurs in multiple forms, including pyroglutamylated Aβ (AβpE). Identification and characterization of the most cytotoxic Aβ species is necessary for advancement in AD diagnostics and therapeutics. While in brain tissue multiple Aβ species act in combination, structure/toxicity studies and immunotherapy trials have been focused on individual forms of Aβ. As a result, the molecular composition and the structural features of "toxic Aβ oligomers" have remained unresolved. Here, we have used a novel approach, hydration from gas phase coupled with isotope-edited Fourier transform infrared (FTIR) spectroscopy, to identify the prefibrillar assemblies formed by Aβ and AβpE and to resolve the structures of both peptides in combination. The peptides form unusual β-sheet oligomers stabilized by intramolecular H-bonding as opposed to intermolecular H-bonding in the fibrils. Time-dependent morphological changes in peptide assemblies have been visualized by atomic force microscopy. Aβ/AβpE hetero-oligomers exert unsurpassed cytotoxic effect on PC12 cells as compared to oligomers of individual peptides or fibrils. These findings lead to a novel concept that Aβ/AβpE hetero-oligomers, not just Aβ or AβpE oligomers, constitute the main neurotoxic conformation. The hetero-oligomers thus present a new biomarker that may be targeted for development of more efficient diagnostic and immunotherapeutic strategies to combat AD.

Alternative Selection of β-Site APP-Cleaving Enzyme 1 (BACE1) Cleavage Sites in Amyloid β-Protein Precursor (APP) Harboring Protective and Pathogenic Mutations within the Aβ Sequence

β-Site APP-cleaving enzyme 1 (BACE1) cleaves amyloid β-protein precursor (APP) at the bond between Met⁶⁷¹ and Asp⁶⁷² (β-site) to generate the carboxyl-terminal fragment (CTFβ/C99). BACE1 also cleaves APP at another bond between Thr⁶⁸¹ and Gln⁶⁸² (β'-site), yielding CTFβ'/C89. Cleavage of CTFβ/C99 by γ-secretase generates Aβ(1-XX), whereas cleavage of CTFβ'/C89 generates Aβ(11-XX). Thus, β'-site cleavage by BACE1 is amyloidolytic rather than amyloidogenic. β' cleavage of mouse APP is more common than the corresponding cleavage of human APP. We found that the H684R substitution within human Aβ, which replaces the histidine in the human protein with the arginine found at the corresponding position in mouse, facilitated β' cleavage irrespective of the species origin of BACE1, thereby significantly increasing the level of Aβ(11-XX) and decreasing the level of Aβ(1-XX). Thus, amino acid substitutions within the Aβ sequence influenced the selectivity of alternative β- or β'-site cleavage of APP by BACE1. In familial Alzheimer's disease (FAD), the APP gene harbors pathogenic variations such as the Swedish (K670N/M671L), Leuven (E682K), and A673V mutations, all of which decrease Aβ(11-40) generation, whereas the protective Icelandic mutation (A673T) increases generation of Aβ(11-40). Thus, A673T promotes β' cleavage of APP and protects subjects against AD. In addition, CTFβ/C99 was cleaved by excess BACE1 activity to generate CTFβ'/C89, followed by Aβ(11-40), even if APP harbored pathogenic mutations. The resultant Aβ(11-40) was more metabolically labile in vivo than Aβ(1-40). Our analysis suggests that some FAD mutations in APP are amyloidogenic and/or amyloidolytic via selection of alternative BACE1 cleavage sites.

Amyloidogenic amyloid-β-peptide variants induce microbial agglutination and exert antimicrobial activity

Amyloid-β (Aβ) peptides are the main components of the plaques found in the brains of patients with Alzheimer's disease. However, Aβ peptides are also detectable in secretory compartments and peripheral blood contains a complex mixture of more than 40 different modified and/or N- and C-terminally truncated Aβ peptides. Recently, anti-infective properties of Aβ peptides have been reported. Here, we investigated the interaction of Aβ peptides of different lengths with various bacterial strains and the yeast Candida albicans. The amyloidogenic peptides Aβ1-42, Aβ2-42, and Aβ3p-42 but not the non-amyloidogenic peptides Aβ1-40 and Aβ2-40 bound to microbial surfaces. As observed by immunocytochemistry, scanning electron microscopy and Gram staining, treatment of several bacterial strains and Candida albicans with Aβ peptide variants ending at position 42 (Aβx-42) caused the formation of large agglutinates. These aggregates were not detected after incubation with Aβx-40. Furthermore, Aβx-42 exerted an antimicrobial activity on all tested pathogens, killing up to 80% of microorganisms within 6 h. Aβ1-40 only had a moderate antimicrobial activity against C. albicans. Agglutination of Aβ1-42 was accelerated in the presence of microorganisms. These data demonstrate that the amyloidogenic Aβx-42 variants have antimicrobial activity and may therefore act as antimicrobial peptides in the immune system.

Assembly of Stefin B into Polymorphic Oligomers Probed by Discrete Molecular Dynamics

Assembly of an amyloidogenic protein stefin B into molten globule oligomers is studied by efficient discrete molecular dynamics. Consistent with in vitro findings, tetramers form primarily through dimer association, resulting in a decreased trimer abundance. Oligomers up to heptamers display elongated rod-like morphologies akin to protofibrils, whereas larger oligomers, decamers through dodecamers, form elongated, branched, as well as annular structures, providing structural insights into pore forming ability and toxicity of amyloidogenic proteins.

Amino acid substitutions [K16A] and [K28A] distinctly affect amyloid β-protein oligomerization

Amyloid β-protein (A β) assembles into oligomers that play a seminal role in Alzheimer's disease (AD), a leading cause of dementia among the elderly. Despite undisputed importance of A β oligomers, their structure and the basis of their toxicity remain elusive. Previous experimental studies revealed that the [K16A] substitution strongly inhibits toxicity of the two predominant A β alloforms in the brain, A β 40 and A β 42, whereas the [K28A] substitution exerts only a moderate effect. Here, folding and oligomerization of [A16]A β 40, [A28]A β 40, [A16]A β 42, and [A28]A β 42 are examined by discrete molecular dynamics (DMD) combined with a four-bead implicit solvent force field, DMD4B-HYDRA, and compared to A β 40 and A β 42 oligomer formation. Our results show that both substitutions promote A β 40 and A β 42 oligomerization and that structural changes to oligomers are substitution- and alloform-specific. The [K28A] substitution increases solvent-accessible surface area of hydrophobic residues and the intrapeptide N-to-C terminal distance within oligomers more than the [K16A] substitution. The [K16A] substitution decreases the overall β-strand content, whereas the [K28A] substitution exerts only a modest change. Substitution-specific tertiary and quaternary structure changes indicate that the [K16A] substitution induces formation of more compact oligomers than the [K28A] substitution. If the structure-function paradigm applies to A β oligomers, then the observed substitution-specific structural changes in A β 40 and A β 42 oligomers are critical for understanding the structural basis of A β oligomer toxicity and correct identification of therapeutic targets against AD.

An Account of Amyloid Oligomers: Facts and Figures Obtained from Experiments and Simulations

The deposition of amyloid in brain tissue in the context of neurodegenerative diseases involves the formation of intermediate species-termed oligomers-of lower molecular mass and with structures that deviate from those of mature amyloid fibrils. Because these oligomers are thought to be primarily responsible for the subsequent disease pathogenesis, the elucidation of their structure is of enormous interest. Nevertheless, because of the high aggregation propensity and the polydispersity of oligomeric species formed by the proteins or peptides in question, the preparation of appropriate samples for high-resolution structural methods has proven to be rather difficult. This is why theoretical approaches have been of particular importance in gaining insights into possible oligomeric structures for some time. Only recently has it been possible to achieve some progress with regard to the experimentally based structural characterization of defined oligomeric species. Here we discuss how theory and experiment are used to determine oligomer structures and what can be done to improve the integration of the two disciplines.

Soluble Prion Protein Binds Isolated Low Molecular Weight Amyloid-β Oligomers Causing Cytotoxicity Inhibition

A growing number of observations indicate that soluble amyloid-β (Aβ) oligomers play a major role in Alzheimer's disease. Recent studies strongly suggest that at least some of the neurotoxic effects of these oligomers are mediated by cellular, membrane-anchored prion protein and that Aβ neurotoxicity can be inhibited by soluble recombinant prion protein (rPrP) and its fragments. However, the mechanism by which rPrP interacts with Aβ oligomers and prevents their toxicity is largely unknown, and studies in this regard are hindered by the large structural heterogeneity of Aβ oligomers. To overcome this difficulty, here we used photoinduced cross-linking of unmodified proteins (PICUP) to isolate well-defined oligomers of Aβ42 and characterize these species with regard to their cytotoxicity and interaction with rPrP, as well the mechanism by which rPrP inhibits Aβ42 cytotoxicity. Our data shows that the addition of rPrP to the assembling Aβ42 results in a shift in oligomer size distribution, decreasing the population of toxic tetramers and higher order oligomers and increasing the population of nontoxic (and possibly neuroprotective) monomers. Isolated oligomeric species of Aβ42 are cytotoxic to primary neurons and cause permeation of model lipid bilayers. These toxic effects, which are oligomer size-dependent, can be inhibited by the addition of rPrP, and our data suggest potential mechanisms of this inhibitory action. This insight should help in current efforts to develop PrP-based therapeutics for Alzheimer's disease.

Applications of Discrete Molecular Dynamics in biology and medicine

Discrete Molecular Dynamics (DMD) is a physics-based simulation method using discrete energetic potentials rather than traditional continuous potentials, allowing microsecond time scale simulations of biomolecular systems to be performed on personal computers rather than supercomputers or specialized hardware. With the ongoing explosion in processing power even in personal computers, applications of DMD have similarly multiplied. In the past two years, researchers have used DMD to model structures of disease-implicated protein folding intermediates, study assembly of protein complexes, predict protein-protein binding conformations, engineer rescue mutations in disease-causative protein mutants, design a protein conformational switch to control cell signaling, and describe the behavior of polymeric dispersants for environmental cleanup of oil spills, among other innovative applications.

Baicalin attenuates alzheimer-like pathological changes and memory deficits induced by amyloid β1-42 protein

Baicalin is one bioactive flavone with anti-inflammatory and neuroprotective activities. The neuroprotective effects of baicalin on pathological changes and behavioral deficits were explored in a mouse model of amyloid β (Aβ)(1-42) protein-induced Alzheimer's disease (AD). Mice received a bilateral injection of Aβ(1-42) protein into the hippocampus, then they were treated with baicalin (30, 50 and 100 mg/kg body weight, orally) or Tween 80. The therapeutic effects of baicalin were monitored by Morris water maze trial and probe test. Then mice were sacrificed for immunohistochemistry and western blot analysis. After a relatively short-term treatment of 14 days, 100 mg/kg of baicalin significantly ameliorated memory impairment in the Morris water maze test and probe test, and also attenuated glial cell activations and increase of TNF-α and IL-6 expressions induced by Aβ(1-42) protein. These results suggest that baicalin ameliorated Aβ(1-42) protein-related pathology and cognitive dysfunction via its anti-neuroinflammatory activity, and may be a potential candidate for the treatment of AD.

Energetic contributions of residues to the formation of early amyloid-β oligomers

Low-weight amyloid-β (Aβ) oligomers formed at early stages of oligomerization rather than fibril assemblies seem to be the toxic components that drive neurodegeneration in Alzheimer's disease. Unfortunately, detailed knowledge of the structure of these early oligomers at the residue level is not yet available. In this study, we performed all-atom explicit solvent molecular dynamics simulations to examine the oligomerization process of Aβ10-35 monomers when forming dimers, trimers, tetramers and octamers, with four independent simulations of a total simulated time of 3 μs for each oligomer system. The decomposition of the stability free energy by MM-GBSA methodology allowed us to unravel the network of energetic interactions that stabilize such oligomers. The contribution of the intermonomeric van der Waals term is the most significant energy feature of the oligomerization process, consistent with the so-called hydrophobic effect. Furthermore, the decomposition of the stability free energy into residues and residue-pairwise terms revealed that it is mainly apolar interactions between the three specific hydrophobic fragments 31-35 (C-terminal region), 17-20 (central hydrophobic core) and 12-14 (N-terminal region) that are responsible for such a favourable effect. The conformation in which the hydrophobic cthr-chc interaction is oriented perpendicularly is particularly important. We propose three other model substructures that favour the oligomerization process and can thus be considered as molecular targets for future inhibitors. Understanding Aβ oligomerization at the residue level could lead to more efficient design of inhibitors of this process.

Disordered amyloidogenic peptides may insert into the membrane and assemble into common cyclic structural motifs

Aggregation of disordered amyloidogenic peptides into oligomers is the causative agent of amyloid-related diseases. In solution, disordered protein states are characterized by heterogeneous ensembles. Among these, β-rich conformers self-assemble via a conformational selection mechanism to form energetically-favored cross-β structures, regardless of their precise sequences. These disordered peptides can also penetrate the membrane, and electrophysiological data indicate that they form ion-conducting channels. Based on these and additional data, including imaging and molecular dynamics simulations of a range of amyloid peptides, Alzheimer's amyloid-β (Aβ) peptide, its disease-related variants with point mutations and N-terminal truncated species, other amyloidogenic peptides, as well as a cytolytic peptide and a synthetic gel-forming peptide, we suggest that disordered amyloidogenic peptides can also present a common motif in the membrane. The motif consists of curved, moon-like β-rich oligomers associated into annular organizations. The motif is favored in the lipid bilayer since it permits hydrophobic side chains to face and interact with the membrane and the charged/polar residues to face the solvated channel pores. Such channels are toxic since their pores allow uncontrolled leakage of ions into/out of the cell, destabilizing cellular ionic homeostasis. Here we detail Aβ, whose aggregation is associated with Alzheimer's disease (AD) and for which there are the most abundant data. AD is a protein misfolding disease characterized by a build-up of Aβ peptide as senile plaques, neurodegeneration, and memory loss. Excessively produced Aβ peptides may directly induce cellular toxicity, even without the involvement of membrane receptors through Aβ peptide-plasma membrane interactions.

N-truncation and pyroglutaminylation enhances the opsonizing capacity of Aβ-peptides and facilitates phagocytosis by macrophages and microglia

Abnormal accumulations of amyloid-β (Aβ)-peptides are one of the pathological hallmarks of Alzheimer's disease (AD). The precursor of the Aβ-peptides, the amyloid precursor protein (APP), is also found in peripheral blood cells, but its function in these cells remains elusive. We previously observed that mononuclear phagocytes release Aβ-peptides during activation and phagocytosis, suggesting a physiologic role in inflammatory processes. Here, we show that supplementing the media with soluble N-terminally truncated Aβ(2-40) and Aβ(2-42) as well as Aβ(1-42) induced the phagocytosis of polystyrene particles (PSPs) by primary human monocytes. If the PSPs were pre-incubated with Aβ-peptides, phagocytosis was induced by all tested Aβ-peptide species. N-terminally truncated Aβ(x-42) induced the phagocytosis of PSPs significantly more effectively than did Aβ(x-40). Similarly, the phagocytosis of Escherichia coli by GM-CSF- and M-CSF-elicited macrophages as well as microglia was particularly facilitated by pre-incubation with N-terminally truncated Aβ(x-42). The proinflammatory polarization of monocytes was indicated by the reduced MSRI expression and IL-10 secretion after phagocytosis of PSPs coated with Aβ(1-42), Aβ(2-42) and Aβ(3p-42). Polarization of the macrophages by GM-CSF reduced the phagocytic activity, but it did not affect the capabilities of Aβ-peptides to opsonize prey. Taken together, Aβ-peptides support phagocytosis as soluble factors and act as opsonins. Differential effects among the Aβ-peptide variants point to distinct mechanisms of interaction among monocytes/macrophages, prey and Aβ-peptides. A proinflammatory polarization induced by the phagocytosis of Aβ-peptide coated particles may provide a model for the chronic inflammatory reaction and sustained plaque deposition in AD.

Stability of amyloid oligomers

Molecular simulations are now commonly used to complement experimental techniques in investigating amyloids and their role in human diseases. In this chapter, we will summarize techniques and approaches often used in amyloid simulations and will present recent success stories. Our examples will be focused on lessons learned from molecular dynamics simulations in aqueous environments that start from preformed aggregates. These studies explore the limitations that arise from the choice of force field, the role of mutations in the growth of amyloid aggregates, segmental polymorphism, and the importance of cross-seeding. Furthermore, they give evidence for potential toxicity mechanisms. We finally discuss the role of molecular simulations in the search for aggregation inhibitors.

Minimal model of self-assembly: emergence of diversity and complexity

Molecular self-assembly is ubiquitous in nature, yet prediction of assembly pathways from fundamental interparticle interactions has yet to be achieved. Here, we introduce a minimal self-assembly model with two attractive and two repulsive beads bound into a tetrahedron. The model is associated with a single parameter η defined as the repulsive to attractive interaction ratio. We explore self-assembly pathways and resulting assembly morphologies for different η values by discrete molecular dynamics. Our results demonstrate that η governs the assembly dynamics and resulting assembly morphologies, revealing an unexpected diversity and complexity for 0.5 ≤ η < 1. One of the key processes that governs the assembly dynamics is assembly breakage, which emerges spontaneously at η > 0 with the breakage rate increasing with η. The observed assembly pathways display a broad variety of assembly structures characteristic of aggregation of amyloidogenic proteins, including quasi-spherical oligomers that coassemble into elongated protofibrils, followed by a conversion into ordered polymorphic fibril-like aggregates. We further demonstrate that η can be meaningfully mapped onto amyloidogenic protein sequences, with the majority of amyloidogenic proteins characterized by 0.5 ≤ η < 1. Prion proteins, which are known to form highly breakage-prone fibrils, are characterized by η > 1, consistent with the model predictions. Our model thus provides a theoretical basis for understanding the universal aspects of aggregation pathways of amyloidogenic proteins relevant to human disease. As the model is not specific to proteins, these findings represent an important step toward understanding and predicting assembly dynamics of not only proteins but also viruses, colloids, and nanoparticles.

Computational approaches to understanding protein aggregation in neurodegeneration

The generation of toxic non-native protein conformers has emerged as a unifying thread among disorders such as Alzheimer's disease, Parkinson's disease, and amyotrophic lateral sclerosis. Atomic-level detail regarding dynamical changes that facilitate protein aggregation, as well as the structural features of large-scale ordered aggregates and soluble non-native oligomers, would contribute significantly to current understanding of these complex phenomena and offer potential strategies for inhibiting formation of cytotoxic species. However, experimental limitations often preclude the acquisition of high-resolution structural and mechanistic information for aggregating systems. Computational methods, particularly those combine both all-atom and coarse-grained simulations to cover a wide range of time and length scales, have thus emerged as crucial tools for investigating protein aggregation. Here we review the current state of computational methodology for the study of protein self-assembly, with a focus on the application of these methods toward understanding of protein aggregates in human neurodegenerative disorders.

Atomic and dynamic insights into the beneficial effect of the 1,4-naphthoquinon-2-yl-L-tryptophan inhibitor on Alzheimer's Aβ1-42 dimer in terms of aggregation and toxicity

Aggregation of the amyloid β protein (Aβ) peptide with 40 or 42 residues is one key feature in Alzheimer's disease (AD). The 1,4-naphthoquinon-2-yl-L-tryptophan (NQTrp) molecule was reported to alter Aβ self-assembly and reduce toxicity. Though nuclear magnetic resonance experiments and various simulations provided atomic information about the interaction of NQTrp with Aβ peptides spanning the regions of residues 12-28 and 17-42, none of these studies were conducted on the full-length Aβ1-42 peptide. To this end, we performed extensive atomistic replica exchange molecular dynamics simulations of Aβ1-42 dimer with two NQTrp molecules in explicit solvent, by using a force field known to fold diverse proteins correctly. The interactions between NQTrp and Aβ1-42, which change the Aβ interface by reducing most of the intermolecular contacts, are found to be very dynamic and multiple, leading to many transient binding sites. The most favorable binding residues are Arg5, Asp7, Tyr10, His13, Lys16, Lys18, Phe19/Phe20, and Leu34/Met35, providing therefore a completely different picture from in vitro and in silico experiments with NQTrp with shorter Aβ fragments. Importantly, the 10 hot residues that we identified explain the beneficial effect of NQTrp in reducing both the level of Aβ1-42 aggregation and toxicity. Our results also indicate that there is room to design more efficient drugs targeting Aβ1-42 dimer against AD.

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.