Controlling the aggregation and rate of release in order to improve insulin formulation: molecular dynamics study of full-length insulin amyloid oligomer models

Surface-Induced Protein Aggregation and Particle Formation in Biologics: Current Understanding of Mechanisms, Detection and Mitigation Strategies

Protein stability against aggregation is a major quality concern for the production of safe and effective biopharmaceuticals. Amongst the different drivers of protein aggregation, increasing evidence indicates that interactions between proteins and interfaces represent a major risk factor for the formation of protein aggregates in aqueous solutions. Potentially harmful surfaces relevant to biologics manufacturing and storage include air-water and silicone oil-water interfaces as well as materials from different processing units, storage containers, and delivery devices. The impact of some of these surfaces, for instance originating from impurities, can be difficult to predict and control. Moreover, aggregate formation may additionally be complicated by the simultaneous presence of interfacial, hydrodynamic and mechanical stresses, whose contributions may be difficult to deconvolute. As a consequence, it remains difficult to identify the key chemical and physical determinants and define appropriate analytical methods to monitor and predict protein instability at these interfaces. In this review, we first discuss the main mechanisms of surface-induced protein aggregation. We then review the types of contact materials identified as potentially harmful or detected as potential triggers of proteinaceous particle formation in formulations and discuss proposed mitigation strategies. Finally, we present current methods to probe surface-induced instabilities, which represent a starting point towards assays that can be implemented in early-stage screening and formulation development of biologics.

Progress in Simulation Studies of Insulin Structure and Function

Insulin is a peptide hormone known for chiefly regulating glucose level in blood among several other metabolic processes. Insulin remains the most effective drug for treating diabetes mellitus. Insulin is synthesized in the pancreatic β-cells where it exists in a compact hexameric architecture although its biologically active form is monomeric. Insulin exhibits a sequence of conformational variations during the transition from the hexamer state to its biologically-active monomer state. The structural transitions and the mechanism of action of insulin have been investigated using several experimental and computational methods. This review primarily highlights the contributions of molecular dynamics (MD) simulations in elucidating the atomic-level details of conformational dynamics in insulin, where the structure of the hormone has been probed as a monomer, dimer, and hexamer. The effect of solvent, pH, temperature, and pressure have been probed at the microscopic scale. Given the focus of this review on the structure of the hormone, simulation studies involving interactions between the hormone and its receptor are only briefly highlighted, and studies on other related peptides (e.g., insulin-like growth factors) are not discussed. However, the review highlights conformational dynamics underlying the activities of reported insulin analogs and mimetics. The future prospects for computational methods in developing promising synthetic insulin analogs are also briefly highlighted.

Molecular Docking, Dynamics, and WaterSwap Analysis to Identify Anti-aggregating Agents of Insulin and IFN-β

There are several challenges in the development, and formulation of biologics, particularly concerning their physical stabilities. The self-assembly of peptides like human insulin and interferon beta (IFN-β) has potential to form aggregates in pharmaceutical formulation. Therefore, it is a significant problem in the manufacturing, storage, and delivery of insulin and IFN-β formulations. Amino acids as aggregation suppressing additives have been used to stabilize proteins during manufacturing and storage. Several changes to the B chain's C-terminus have been proposed in an attempt to improve insulin formulation. The core segments of the A and B chains (SLYQLENY and LVEALYLV) have recently been identified as sheet-forming areas, and their microcrystalline structures have been exploited to construct a high-resolution insulin amyloid fibril model. Here, we have chosen twenty-one amino acids to develop as additives in rendering the insulin and IFN-β aggregations. Thereafter, integrated molecular docking studies of single layer monomers of full-length insulin and IFN-β have been performed to identify structural elements (amino acids) that can act as disaggregating agents. The stability of the best-docked amino acid complexes was judged using molecular dynamics studies. Finally, phenylalanine was identified as a disaggregation agent for insulin, and lysine, tyrosine, phenylalanine, and tryptophan were identified as disaggregation agents for IFN-β from the molecular dynamics study. These findings may open a novel proposal to explore further in vitro studies to increase the stability of the insulin and IFN-β formulation.

Study on the Kinetics and Mechanism of Ferrocene-Tripeptide Inhibiting Insulin Aggregation

Peptides are gaining popularity as neurodegenerative disease-targeted drugs due to their medicinal value and the simplicity in the biomedicine and pharmaceutical industry field. In this study, based on previously studied...

Determination of amyloid core regions of insulin analogues fibrils

A rapid-acting insulin lispro and long-acting insulin glargine are commonly used for the treatment of diabetes. Clinical cases have described the formation of injectable amyloidosis with these insulin analogues, but their amyloid core regions of fibrils were unknown. To reveal these regions, we have analysed the hydrolyzates of insulin fibrils and its analogues using high-performance liquid chromatography and mass spectrometry methods and found that insulin and its analogues have almost identical amyloid core regions that intersect with the predicted amyloidogenic regions. The obtained results can be used to create new insulin analogues with a low ability to form fibrils.

Abbreviations

a.a., amino acid residues; HPLC-MS, high-performance liquid chromatography/mass spectrometry; m/z, mass-to-charge ratio; TEM, transmission electron microscopy.

Computational Investigation of the Binding Dynamics of Oligo p-Phenylene Ethynylene Fluorescence Sensors and Aβ Oligomers

Amyloid protein aggregates are pathological hallmarks of neurodegenerative disorders such as Alzheimer's (AD) and Parkinson's (PD) diseases and are believed to be formed well before the onset of neurodegeneration and cognitive impairment. Monitoring the course of protein aggregation is thus vital to understanding and combating these diseases. We have recently demonstrated that a novel class of fluorescence sensors, oligomeric p-phenylene ethynylene (PE)-based electrolytes (OPEs) selectively bind to and detect prefibrillar and fibrillar aggregates of AD-related amyloid-β (Aβ) peptides over monomeric Aβ. In this study, we investigated the binding between two OPEs, anionic OPE12- and cationic OPE24+, and to two different β-sheet rich Aβ oligomers using classical all-atom molecular dynamics simulations. Our simulations have revealed a number of OPE binding sites on Aβ oligomer surfaces, and these sites feature hydrophobic amino acids as well as oppositely charged amino acids. Binding energy calculations show energetically favorable interactions between both anionic and cationic OPEs with Aβ oligomers. Moreover, OPEs bind as complexes as well as single molecules. Compared to free OPEs, Aβ protofibril bound OPEs show backbone planarization with restricted rotations and reduced hydration of the ethyl ester end groups. These characteristics, along with OPE complexation, align with known mechanisms of binding induced OPE fluorescence turn-on and spectral shifts from a quenched, unbound state in aqueous solutions. This study thus sheds light on the molecular-level details of OPE-Aβ protofibril interactions and provides a structural basis for fluorescence turn-on sensing modes of OPEs.

Human Serum Albumin Binds Native Insulin and Aggregable Insulin Fragments and Inhibits Their Aggregation

The purpose of this study was to investigate whether Human Serum Albumin (HSA) can bind native human insulin and its A13–A19 and B12–B17 fragments, which are responsible for the aggregation of the whole hormone. To label the hormone and both hot spots, so that their binding positions within the HSA could be identified, 4-(1-pyrenyl)butyric acid was used as a fluorophore. Triazine coupling reagent was used to attach the 4-(1-pyrenyl)butyric acid to the N-terminus of the peptides. When attached to the peptides, the fluorophore showed extended fluorescence lifetimes in the excited state in the presence of HSA, compared to the samples in buffer solution. We also analyzed the interactions of unlabeled native insulin and its hot spots with HSA, using circular dichroism (CD), the microscale thermophoresis technique (MST), and three independent methods recommended for aggregating peptides. The CD spectra indicated increased amounts of the α-helical secondary structure in all analyzed samples after incubation. Moreover, for each of the two unlabeled hot spots, it was possible to determine the dissociation constant in the presence of HSA, as 14.4 µM (A13–A19) and 246 nM (B12–B17). Congo Red, Thioflavin T, and microscopy assays revealed significant differences between typical amyloids formed by the native hormone or its hot-spots and the secondary structures formed by the complexes of HSA with insulin and A13–A19 and B12–B17 fragments. All results show that the tested peptide-probe conjugates and their unlabeled analogues interact with HSA, which inhibits their aggregation.

Search for New Aggregable Fragments of Human Insulin

In this study, three independent methods were used to identify short fragment of both chains of human insulin which are prone for aggregation. In addition, circular dichroism (CD) research was conducted to understand the progress of aggregation over time. The insulin fragments (deca- and pepta-peptides) were obtained by solid-phase synthesis using 4-(4,6-dimethoxy-1,3,5-triazin-2-yl)-4-methylmorpholinium toluene-4-sulfonate (DMT/NMM/TosO^-) as a coupling reagent. Systematic studies allowed identification of the new fragments, expected to be engaged in triggering aggregation of the entire structure of human insulin under physiological conditions. It was found that the aggregation process occurs through various structural conformers and may favor the formation of a fibrous structure of aggregate.

How accurate are your simulations? Effects of confined aqueous volume and AMBER FF99SB and CHARMM22/CMAP force field parameters on structural ensembles of intrinsically disordered proteins: Amyloid-β42 in water

Amyloid-β₄₂ (Aβ₄₂) is an intrinsically disordered peptide intimately related to the pathogenesis of several neurodegenerative diseases. Molecular dynamics (MD) simulations are extensively utilized in the characterization of the structures and conformational dynamics of intrinsically disordered proteins (IDPs) including Aβ₄₂, with AMBER and CHARMM parameters being commonly used in these studies. Recently, comparison of the effects of force field parameters on the Aβ₄₂ structures has started to gain significant attention. In this study, the structures of Aβ₄₂ are simulated using AMBER FF99SB and CHARMM22/CMAP parameters via replica exchange MD simulations utilizing a widely used clustering algorithm. These analyses show that the structural properties (extent and positioning of the elements of secondary and tertiary structure), radius of gyration values, number and position of salt bridges are extremely dependent on the chosen force field parameters notably with the usage of clustering algorithms. For example, predicted secondary structure elements, which are of the great importance for better understanding of the molecular mechanisms of neurodegenerative diseases, deviate enormously in models generated using currently available force field parameters for proteins. Based on the derived models, chemical shift values are calculated and compared to the experimentally determined data. This comparison revealed that although both force field parameters yield results in agreement with experiments, the obtained structural properties were rather different using a clustering algorithm. In other words, these results show that the predicted structures depend heavily on the force field parameters. Importantly, since none of the force field parameters currently utilized in MD studies were developed specifically taking into account the disordered nature of IDPs, these findings clearly indicate that new force field parameters have to be developed for IDPs considering their rapid flexibility and dynamics with high amplitude. Furthermore, molecular simulations of IDPs are typically conducted using one water volume. We show that the confined aqueous volume impacts the predicted structural properties of Aβ₄₂ in water. Although up to date, confined aqueous volume effects have been ignored in the MD simulations of IDPs in water, our data indicate that these effects have to be taken into account in predicting the structural and thermodynamic properties of disordered proteins in solution.

The quest for the shortest fragments of A (13-19) and B (12-17) responsible for the aggregation of human insulin

AIM

To identify the shortest components of A13-A19, B12-B17 fragments capable for fibrillation and to validate the dependability of aggregation on the presence of hydroxyl group engaged in the 'tyrosine kissing'.

MATERIALS & METHODS

Fragments A13-A19 and B12-B17 of insulin and all shortened analogues were obtained by using DMT/NMM/TosO(-) as a coupling reagent. The aggregation was studied by three independent tests.

RESULTS

Studies on the susceptibility to aggregation of truncated analogs of insulin amyloidogenic core show three groups of peptides.

CONCLUSION

Truncation of A13-A419 fragment shows that fibrous structures are formed by all peptides bearing (13)H-LeuTyr-OH(14). Propensity to aggregation was found for (16)H-TyrLeu-OH(17) B12-B17 fragment. Tyrosine residue modification by incorporation of tert-butyl group on hydroxyl function gave analogues still predisposed to aggregation.

Conformational changes of Aβ (1-42) monomers in different solvents

Amyloid proteins are known to be the main cause of numerous degenerative and neurodegenerative diseases. In general, amyloids are misfolded from monomers and they tend to have β-strand formations. These misfolded monomers are then transformed into oligomers, fibrils, and plaques. It is important to understand the forming mechanism of amyloids in order to prevent degenerative diseases to occur. Aβ protein is a highly noticeable protein which causes Alzheimer's disease. It is reported that solvents affect the forming mechanism of Aβ amyloids. In this research, Aβ1-42 was analyzed using an all-atom MD simulation with the consideration of effects induced by two disparate solvents: water and DMSO. As a result, two different conformation changes of Aβ1-42 were exhibited in each solvent. It was found that salt-bridge of Asp23 and Lys28 in Aβ1-42 was the key for amyloid folding based on the various analysis including hydrogen bond, electrostatic interaction energy and salt-bridge distance. Since this salt-bridge region plays a crucial role in initiating the misfolding of Aβ1-42, this research may shed a light for studies related in amyloid folding and misfolding.

The effect of structural heterogeneity on the conformation and stability of Aβ–tau mixtures

Oligomeric and fibrillar amyloids, which cause neurodegenerative diseases, are typically formed through repetitive fracture and elongation processes involving single homogeneous amyloid monomers.

Understanding structural characteristics of out-of-register hIAPP amyloid proteins via molecular dynamics

We investigated characteristics of out-of-register (OOR) hIAPP amyloids. By varying the length size of OOR hIAPP, we found 8 layers is most stable. In addition, OOR hIAPP has relative structural instability than in-register hAIPP.

Stability of Iowa mutant and wild type Aβ-peptide aggregates

Recent experiments indicate a connection between the structure of amyloid aggregates and their cytotoxicity as related to neurodegenerative diseases. Of particular interest is the Iowa Mutant, which causes early-onset of Alzheimer's disease. While wild-type Amyloid β-peptides form only parallel beta-sheet aggregates, the mutant also forms meta-stable antiparallel beta sheets. Since these structural variations may cause the difference in the pathological effects of the two Aβ-peptides, we have studied in silico the relative stability of the wild type and Iowa mutant in both parallel and antiparallel forms. We compare regular molecular dynamics simulations with such where the viscosity of the samples is reduced, which, we show, leads to higher sampling efficiency. By analyzing and comparing these four sets of all-atom molecular dynamics simulations, we probe the role of the various factors that could lead to the structural differences. Our analysis indicates that the parallel forms of both wild type and Iowa mutant aggregates are stable, while the antiparallel aggregates are meta-stable for the Iowa mutant and not stable for the wild type. The differences result from the direct alignment of hydrophobic interactions in the in-register parallel oligomers, making them more stable than the antiparallel aggregates. The slightly higher thermodynamic stability of the Iowa mutant fibril-like oligomers in its parallel organization over that in antiparallel form is supported by previous experimental measurements showing slow inter-conversion of antiparallel aggregates into parallel ones. Knowledge of the mechanism that selects between parallel and antiparallel conformations and determines their relative stability may open new avenues for the development of therapies targeting familial forms of early-onset Alzheimer's disease.

Atomistic mechanism of polyphenol amyloid aggregation inhibitors: molecular dynamics study of Curcumin, Exifone, and Myricetin interaction with the segment of tau peptide oligomer

Amyloid fibrils are highly ordered protein aggregates associated with many diseases affecting millions of people worldwide. Polyphenols such as Curcumin, Exifone, and Myricetin exhibit modest inhibition toward fibril formation of tau peptide which is associated with Alzheimer's disease. However, the molecular mechanisms of this inhibition remain elusive. We investigated the binding of three polyphenol molecules to the protofibrils of an amyloidogenic fragment VQIVYK of tau peptide by molecular dynamics simulations in explicit solvent. We find that polyphenols induce conformational changes in the oligomer aggregate. These changes disrupt the amyloid H bonding, perturbing the aggregate. While the structural evolution of the control oligomer with no ligand is limited to the twisting of the β-sheets without their disassembly, the presence of polyphenol molecule pushes the β-sheets apart, and leads to a loosely packed structure where two of four β-sheets dissociate in each of the three cases considered here. The H-bonding capacity of polyphenols is responsible for the observed behavior. The calculated binding free energies and its individual components enabled better understanding of the binding. Results indicated that the contribution from Van der Waals interactions is more significant than electrostatic contribution to the binding. The findings from this study are expected to assist in the development of aggregation inhibitors. Significant binding between polyphenols and aggregate oligomer identified in our simulations confirms the previous experimental observations in which polyphenols refold the tau peptide without forming covalent bonds.

Interplay of sequence, topology and termini charge in determining the stability of the aggregates of GNNQQNY mutants: a molecular dynamics study

This study explores the stabilities of single sheet parallel systems of three sequence variants of ¹GNNQQNY⁷, N2D, N2S and N6D, with variations in aggregate size (5–8) and termini charge (charged or neutral). The aggregates were simulated at 300 and 330 K. These mutations decrease amyloid formation in the yeast prion protein Sup35. The present study finds that these mutations cause instability even in the peptide context. The protonation status of termini is found to be a key determinant of stabilities; other determinants are sequence, position of mutation and aggregate size. All systems with charged termini are unstable, whereas both stable and unstable systems are found when the termini are neutral. When termini are charged, the largest stable aggregate for the N2S and N6D systems has 3 to 4 peptides whereas N2D mutation supports oligomers of larger size (5-and 6-mers) as well. Mutation at 2^nd position (N2S and N2D) results in fewer H-bonds at the mutated as well as neighboring (Gly1/Gln4) positions. However, no such effect is found if mutation is at 6^th position (N6D). The effect of Asn→Asp mutation depends on the position and termini charge: it is more destabilizing at the 2^nd position than at the 6^th in case of neutral termini, however, the opposite is true in case of charged termini. Appearance of twist in stable systems and in smaller aggregates formed in unstable systems suggests that twist is integral to amyloid arrangement. Disorder, dissociation or rearrangement of peptides, disintegration or collapse of aggregates and formation of amorphous aggregates observed in these simulations are likely to occur during the early stages of aggregation also. The smaller aggregates formed due to such events have a variety of arrangements of peptides. This suggests polymorphic nature of oligomers and presence of a heterogeneous mixture of oligomers during early stages of aggregation.

Inter-species cross-seeding: stability and assembly of rat-human amylin aggregates

Diseases such as type 2 diabetes, Alzheimer's and Parkinson's share as common feature the accumulation of mis-folded disease-specific protein aggregates into fibrillar structures, or plaques. These fibrils may either be toxic by themselves, or act as reservoirs for smaller cytotoxic oligomers. This suggests to investigate molecules as potential therapeutics that either reduce fibril formation or increase fibril stability. One example is rat amylin, which can inhibit aggregation of human amylin, a hallmark of type 2 diabetes. In the present paper, we use molecular dynamics to compare the stability of various preformed aggregates, built out of either human amylin, rat amylin, or mixtures of both. We considered two types of fibril-like oligomers: a single-layer in-register conformation, and a double-layer conformation in which the first U-shaped layer consists of rat amylin and the second layer of human amylin. Our results explain the weak amyloid-inhibiting properties of rat amylin and suggest that membrane leakage due to pore formation is responsible for the toxicity of rat amylin observed in a recent experiment. Together, our results put in question the use of rat amylin or the similar FDA approved drug pramlintide as an inhibitor of human amylin aggregation. They also point to mixed human-rat amylin fibril-like oligomers as possible model-systems for studies of amyloid formation that involve cross-species transmission.

In silico cross seeding of Aβ and amylin fibril-like oligomers

Recent epidemiological data have shown that patients suffering from Type 2 Diabetes Mellitus have an increased risk to develop Alzheimer's disease and vice versa. A possible explanation is the cross-sequence interaction between Aβ and amylin. Because the resulting amyloid oligomers are difficult to probe in experiments, we investigate stability and conformational changes of Aβ-amylin heteroassemblies through molecular dynamics simulations. We find that Aβ is a good template for the growth of amylin and vice versa. We see water molecules permeate the β-strand-turn-β-strand motif pore of the oligomers, supporting a commonly accepted mechanism for toxicity of β-rich amyloid oligomers. Aiming for a better understanding of the physical mechanisms of cross-seeding and cell toxicity of amylin and Aβ aggregates, our simulations also allow us to identify targets for the rational design of inhibitors against toxic fibril-like oligomers of Aβ and amylin oligomers.

Full length amylin oligomer aggregation: insights from molecular dynamics simulations and implications for design of aggregation inhibitors

Amyloid oligomers are considered to play essential roles in the pathogenesis of amyloid-related degenerative diseases including type 2 diabetes. Using an explicit solvent all atomic MD simulation, we explored the stability, conformational dynamics and association force of different single-layer models of the full-length wild-type and glycine mutants of amylin (pentamer) obtained from a recent high resolution fibril model. The RMSF profile shows enhanced flexibility in the disorder (Lys1-Cys7) and turn region (Ser19-Gly23), along with smallest fluctuation at the residues (Asn14-Phe15-Leu16-Val17-His18) of β1 region and (Ala25-Ile26-Leu27-Ser28-Ser29) of the β2 region. We obtained a significant difference in backbone RMSD between the wild-type and the mutants, indicating that mutations affected the stability of the peptide. The RMSD and RMSF profiles indicate the edge and loop residues are the primary contributors to the overall conformational changes. The degree of structural similarity between the oligomers in the simulation and the fibril conformation is proposed as the possible explanation for experimentally observed shortening of the nucleation lag phase of amylin with oligomer seeding. On the basis of structure-stability findings, the β1 and β2 portions are optimal target for further anti-amyloid drug design. The MM-PBSA binding energy calculation reveals the binding of amylin: amylin strands in single layer is dominated by contributions from van der Waals interactions. The non-polar solvation term is also found to be favorable. While the electrostatic interactions and polar solvation energy was found to be favorable for the interaction for the larger aggregate and unfavorable for the smaller aggregates. A per-residue decomposition of the binding free energy has been performed to identify the residues contributing most to the self-association free energy. Residues found in the β-sheet regions were found to be key residue making the largest favorable contributions to the single-layer association. The result from our simulation could be used in rational design of new amylinomimetic agent, amylin aggregation inhibitors and amylin-specific biomarkers.

The stability of cylindrin β-barrel amyloid oligomer models-a molecular dynamics study

Small-soluble amyloid oligomers are believed to play a significant role in the pathology of amyloid diseases. Recently, the atomic structure of a toxic oligomer formed by an 11 residue and its tandem repeat was found to have an out-off register antiparallel β-strands in the shape of a β-barrel. In the present article we investigate the effect of mutations in the hydrophobic cores on the structure and dynamic of the β-barrels using all atom multiple molecular dynamics simulations with an explicit solvent. Extending previous experiments with molecular dynamics simulations we systematically test how stability and formation of cylindrin depends on the interplay between hydrophobicity and steric effects of the core residues. We find that strong hydrophobic interactions between geometrically fitting residues keep the strands (both in register and out-off-register interface) in close proximity, which in turn stabilizes the side-chain and main-chain hydrogen bonds, and the salt bridges on the outer surface along the weak out-of-register interface. Our simulations also indicate presence of water molecules in the hydrophobic interior of the cylindrin β-barrel.Proteins 2013.

Side-chain hydrophobicity and the stability of Aβ₁₆₋₂₂ aggregates

Recent mutagenesis studies using the hydrophobic segment of Aβ suggest that aromatic π-stacking interactions may not be critical for fibril formation. We have tested this conjecture by probing the effect of Leu, Ile, and Ala mutation of the aromatic Phe residues at positions 19 and 20, on the double-layer hexametric chains of Aβ fragment Aβ₁₆₋₂₂ using explicit solvent all-atom molecular dynamics. As these simulations rely on the accuracy of the utilized force fields, we first evaluated the dynamic and stability dependence on various force fields of small amyloid aggregates. These initial investigations led us to choose AMBER99SB-ILDN as force field in multiple long molecular dynamics simulations of 100 ns that probe the stability of the wild-type and mutants oligomers. Single-point and double-point mutants confirm that size and hydrophobicity are key for the aggregation and stability of the hydrophobic core region (Aβ₁₆₋₂₂). This suggests as a venue for designing Aβ aggregation inhibitors the substitution of residues (especially, Phe 19 and 20) in the hydrophobic region (Aβ₁₆₋₂₂) with natural and non-natural amino acids of similar size and hydrophobicity.

Structure and dynamics of amyloid-β segmental polymorphisms

It is believed that amyloid-beta (Aβ) aggregates play a role in the pathogenesis of Alzheimer's disease. Aβ molecules form β-sheet structures with multiple interaction sites. This polymorphism gives rise to differences in morphology, physico-chemical property and level of cellular toxicity. We have investigated the conformational stability of various segmental polymorphisms using molecular dynamics simulations and find that the segmental polymorphic models of Aβ retain a U-shaped architecture. Our results demonstrate the importance of inter-sheet side chain-side chain contacts, hydrophobic contacts among the strands (β1 and β2) and of salt bridges in stabilizing the aggregates. Residues in β-sheet regions have smaller fluctuation while those at the edge and loop region are more mobile. The inter-peptide salt bridges between Asp23 and Lys28 are strong compared to intra-chain salt bridge and there is an exchange of the inter-chain salt-bridge with intra-chain salt bridge. As our results suggest that Aβ exists under physiological conditions as an ensemble of distinct segmental polymorphs, it may be necessary to account in the development of therapeutics for Alzheimer's disease the differences in structural stability and aggregation behavior of the various Aβ polymorphic forms.