Separating instability from aggregation propensity in γS-crystallin variants

Molecular dynamics (MD) simulations, circular dichroism (CD), and dynamic light scattering (DLS) measurements were used to investigate the aggregation propensity of the eye-lens protein γS-crystallin. The wild-type protein was investigated along with the cataract-related G18V variant and the symmetry-related G106V variant. The MD simulations suggest that local sequence differences result in dramatic differences in dynamics and hydration between these two apparently similar point mutations. This finding is supported by the experimental measurements, which show that although both variants appear to be mostly folded at room temperature, both display increased aggregation propensity. Although the disease-related G18V variant is not the most strongly destabilized, it aggregates more readily than either the wild-type or the G106V variant. These results indicate that γS-crystallin provides an excellent model system for investigating the role of dynamics and hydration in aggregation by locally unfolded proteins.     

Role of histidine interruption in mitigating the pathological effects of long polyglutamine stretches in SCA1: A molecular approach

Polyglutamine expansions, leading to aggregation, have been implicated in various neurodegenerative disorders. The range of repeats observed in normal individuals in most of these diseases is 19-36, whereas mutant proteins carry 40-81 repeats. In one such disorder, spinocerebellar ataxia (SCA1), it has been reported that certain individuals with expanded polyglutamine repeats in the disease range (Q(12)HQHQ(12)HQHQ(14/15)) but with histidine interruptions were found to be phenotypically normal. To establish the role of histidine, a comparative study of conformational properties of model peptide sequences with (Q(12)HQHQ(12)HQHQ(12)) and without (Q(42)) interruptions is presented here. Q(12)HQHQ(12)HQHQ(12) displays greater solubility and lesser aggregation propensity compared to uninterrupted Q(42) as well as much shorter Q(22). The solvent and temperature-driven conformational transitions (beta structure <--> random coil --> alpha helix) displayed by these model polyQ stretches is also discussed in the present report. The study strengthens our earlier hypothesis of the importance of histidine interruptions in mitigating the pathogenicity of expanded polyglutamine tract at the SCA1 locus. The relatively lower propensity for aggregation observed in case of histidine interrupted stretches even in the disease range suggests that at a very low concentration, the protein aggregation in normal cells, is possibly not initiated at all or the disease onset is significantly delayed. Our present study also reveals that besides histidine interruption, proline interruption in polyglutamine stretches can lower their aggregation propensity.     

Separating Instability from Aggregation Propensity in γS-Crystallin Variants

Molecular dynamics (MD) simulations, circular dichroism (CD), and dynamic light scattering (DLS) measurements were used to investigate the aggregation propensity of the eye-lens protein γS-crystallin. The wild-type protein was investigated along with the cataract-related G18V variant and the symmetry-related G106V variant. The MD simulations suggest that local sequence differences result in dramatic differences in dynamics and hydration between these two apparently similar point mutations. This finding is supported by the experimental measurements, which show that although both variants appear to be mostly folded at room temperature, both display increased aggregation propensity. Although the disease-related G18V variant is not the most strongly destabilized, it aggregates more readily than either the wild-type or the G106V variant. These results indicate that γS-crystallin provides an excellent model system for investigating the role of dynamics and hydration in aggregation by locally unfolded proteins.     

Targeting the early steps of Abeta16-22 protofibril disassembly by N-methylated inhibitors: a numerical study

Aggregation of the Abeta1-40/Abeta1-42 peptides is a key factor in Alzheimer's disease. Though the inhibitory effect of N-methylated Abeta16-22 (mAbeta16-22) peptides is well characterized in vitro, there is little information on how they disassemble full-length Abeta fibrils or block fibril formation. Here, we report coarse-grained implicit solvent molecular dynamics (MD) and replica exchange molecular dynamics (REMD) simulations on Abeta16-22 and mAbeta16-22 monomers, and then a preformed six-chain Abeta16-22 bilayer with either four copies of Abeta16-22 or four copies of mAbeta16-22. Our simulations show that the effect of N-methylation on mAbeta16-22 monomer is to reduce the density of compact forms. While 100 ns MD trajectories do not reveal any significant differences between the two ten-chain systems, the REMD simulations totaling 1 micros help understand the first steps of Abeta16-22 protofibril disassembly by N-methylated inhibitors. Notably, we find that mAbeta16-22 preferentially interacts with Abeta16-22 by blocking both beta-sheet extension and lateral association of layers, but also by intercalation of the inhibitors allowing sequestration of Abeta16-22 peptides. This third binding mode is particularly appealing for blocking Abeta fibrillogenesis.     

In vitro analysis of SpUre2p, a prion-related protein, exemplifies the relationship between amyloid and prion

The yeast Saccharomyces cerevisiae contains in its proteome at least three prion proteins. These proteins (Ure2p, Sup35p, and Rnq1p) share a set of remarkable properties. In vivo, they form aggregates that self-perpetuate their aggregation. This aggregation is controlled by Hsp104, which plays a major role in the growth and severing of these prions. In vitro, these prion proteins form amyloid fibrils spontaneously. The introduction of such fibrils made from Ure2p or Sup35p into yeast cells leads to the prion phenotypes [URE3] and [PSI], respectively. Previous studies on evolutionary biology of yeast prions have clearly established that [URE3] is not well conserved in the hemiascomycetous yeasts and particularly in S. paradoxus. Here we demonstrated that the S. paradoxus Ure2p is able to form infectious amyloid. These fibrils are more resistant than S. cerevisiae Ure2p fibrils to shear force. The observation, in vivo, of a distinct aggregation pattern for GFP fusions confirms the higher propensity of SpUre2p to form fibrillar structures. Our in vitro and in vivo analysis of aggregation propensity of the S. paradoxus Ure2p provides an explanation for its loss of infective properties and suggests that this protein belongs to the non-prion amyloid world.     

Amyloid of the Candida albicans Ure2p prion domain is infectious and has an in-register parallel β-sheet structure

Ure2p of Candida albicans (Ure2(albicans) or CaUre2p) can be a prion in Saccharomyces cerevisiae, but Ure2p of Candida glabrata (Ure2(glabrata)) cannot, even though the Ure2(glabrata) N-terminal domain is more similar to that of the S. cerevisiae Ure2p (Ure2(cerevisiae)) than Ure2(albicans) is. We show that the N-terminal N/Q-rich prion domain of Ure2(albicans) forms amyloid that is infectious, transmitting [URE3alb] to S. cerevisiae cells expressing only C. albicans Ure2p. Using solid-state nuclear magnetic resonance of selectively labeled C. albicans Ure2p(1-90), we show that this infectious amyloid has an in-register parallel β-sheet structure, like that of the S. cerevisiae Ure2p prion domain and other S. cerevisiae prion amyloids. In contrast, the N/Q-rich N-terminal domain of Ure2(glabrata) does not readily form amyloid, and that formed upon prolonged incubation is not infectious.     

Circular dichroism and aggregation studies of amyloid beta (11-8) fragment and its variants

Aggregation of Abeta peptides is a seminal event in Alzheimer's disease. Detailed understanding of Abeta assembly would facilitate the targeting and design of fibrillogenesis inhibitors. Here comparative conformational and aggregation studies using CD spectroscopy and thioflavine T fluorescence assay are presented. As a model peptide, the 11-28 fragment of Abeta was used. This model peptide is known to contain the core region responsible for Abeta aggregation. The structural and aggregational behaviour of the peptide was compared with the properties of its variants corresponding to natural, clinically relevant mutants at positions 21-23 (A21G, E22K, E22G, E22Q and D23N). In HFIP (hexafluoro-2-propanol), a strong alpha-helix inducer, the CD spectra revealed an unexpectedly high amount of beta-sheet conformation. The aggregation process of Abeta(11-28) variants provoked by water addition to HFIP was found to be consistent with a model of an alpha-helix-containing intermediate. The aggregation propensity of all Abeta(11-28) variants was also compared and discussed.     

Long time dynamics of Met-enkephalin: comparison of explicit and implicit solvent models

Met-enkephalin is one of the smallest opiate peptides. Yet, its dynamical structure and receptor docking mechanism are still not well understood. The conformational dynamics of this neuron peptide in liquid water are studied here by using all-atom molecular dynamics (MD) and implicit water Langevin dynamics (LD) simulations with AMBER potential functions and the three-site transferable intermolecular potential (TIP3P) model for water. To achieve the same simulation length in physical time, the full MD simulations require 200 times as much CPU time as the implicit water LD simulations. The solvent hydrophobicity and dielectric behavior are treated in the implicit solvent LD simulations by using a macroscopic solvation potential, a single dielectric constant, and atomic friction coefficients computed using the accessible surface area method with the TIP3P model water viscosity as determined here from MD simulations for pure TIP3P water. Both the local and the global dynamics obtained from the implicit solvent LD simulations agree very well with those from the explicit solvent MD simulations. The simulations provide insights into the conformational restrictions that are associated with the bioactivity of the opiate peptide dermorphin for the delta-receptor.     

A comparative study of the relationship between protein structure and beta-aggregation in globular and intrinsically disordered proteins

A growing number of proteins are being identified that are biologically active though intrinsically disordered, in sharp contrast with the classic notion that proteins require a well-defined globular structure in order to be functional. At the same time recent work showed that aggregation and amyloidosis are initiated in amino acid sequences that have specific physico-chemical properties in terms of secondary structure propensities, hydrophobicity and charge. In intrinsically disordered proteins (IDPs) such sequences would be almost exclusively solvent-exposed and therefore cause serious solubility problems. Further, some IDPs such as the human prion protein, synuclein and Tau protein are related to major protein conformational diseases. However, this scenario contrasts with the large number of unstructured proteins identified, especially in higher eukaryotes, and the fact that the solubility of these proteins is often particularly good. We have used the algorithm TANGO to compare the beta aggregation tendency of a set of globular proteins derived from SCOP and a set of 296 experimentally verified, non-redundant IDPs but also with a set of IDPs predicted by the algorithms DisEMBL and GlobPlot. Our analysis shows that the beta-aggregation propensity of all-alpha, all-beta and mixed alpha/beta globular proteins as well as membrane-associated proteins is fairly similar. This illustrates firstly that globular structures possess an appreciable amount of structural frustration and secondly that beta-aggregation is not determined by hydrophobicity and beta-sheet propensity alone. We also show that globular proteins contain almost three times as much aggregation nucleating regions as IDPs and that the formation of highly structured globular proteins comes at the cost of a higher beta-aggregation propensity because both structure and aggregation obey very similar physico-chemical constraints. Finally, we discuss the fact that although IDPs have a much lower aggregation propensity than globular proteins, this does not necessarily mean that they have a lower potential for amyloidosis.     

Assessing the effect of D59P mutation in the DE loop region in amyloid aggregation propensity of β2‐microglobulin: A molecular dynamics simulation study

Dialysis‐related amyloidosis (DRA) is a severe condition characterized by the accumulation of amyloidogenic β2‐microglobulin (β2m) protein around skeletal joints and bones. The recent studies highlighted a critical role of the DE loop region for β2m stability and amyloid aggregation propensity. Despite significant efforts, the molecular mechanism of enhanced aggregation due to D59P mutation in the DE loop region remain elusive. In the present study, explicit‐solvent molecular dynamics (MD) simulations were performed to examine the key changes in the structural and dynamic properties of wild type (wt) β2m upon D59P mutation. MD simulations reveal a decrease in the average number of hydrogen bonds in the loop regions on D59P mutation that enhances conformational flexibility, which lead to higher aggregation propensity of D59P as compare to wt β2m. The principal component analysis (PCA) highlight that D59P covers a larger region of phase space and display a higher trace value than wt β2m, which suggest an overall enhancement in the conformational flexibility. D59P display two minimum energy basins in the free energy landscape (FEL) that are associated with thermodynamically less stable conformational states as compare to single minimum energy basin in wt β2m. The present study provides theoretical insights into the molecular mechanism behind the higher aggregation propensity of D59P as compare to wt β2m.     

Molecular dynamics simulations of peptides and proteins with a continuum electrostatic model based on screened Coulomb potentials

A continuum electrostatics approach for molecular dynamics (MD) simulations of macromolecules is presented and analyzed for its performance on a peptide and a globular protein. The approach incorporates the screened Coulomb potential (SCP) continuum model of electrostatics, which was reported earlier. The model was validated in a broad set of tests some of which were based on Monte Carlo simulations that included single amino acids, peptides, and proteins. The implementation for large-scale MD simulations presented in this article is based on a pairwise potential that makes the electrostatic model suitable for fast analytical calculation of forces. To assess the suitability of the approach, a preliminary validation is conducted, which consists of (i) a 3-ns MD simulation of the immunoglobulin-binding domain of streptococcal protein G, a 56-residue globular protein and (ii) a 3-ns simulation of Dynorphin, a biological peptide of 17 amino acids. In both cases, the results are compared with those obtained from MD simulations using explicit water (EW) molecules in an all-atom representation. The initial structure of Dynorphin was assumed to be an alpha-helix between residues 1 and 9 as suggested from NMR measurements in micelles. The results obtained in the MD simulations show that the helical structure collapses early in the simulation, a behavior observed in the EW simulation and consistent with spectroscopic data that suggest that the peptide may adopt mainly an extended conformation in water. The dynamics of protein G calculated with the SCP implicit solvent model (SCP-ISM) reveals a stable structure that conserves all the elements of secondary structure throughout the entire simulation time. The average structures calculated from the trajectories with the implicit and explicit solvent models had a cRMSD of 1.1 A, whereas each average structure had a cRMSD of about 0.8A with respect to the X-ray structure. The main conformational differences of the average structures with respect to the crystal structure occur in the loop involving residues 8-14. Despite the overall similarity of the simulated dynamics with EW and SCP models, fluctuations of side-chains are larger when the implicit solvent is used, especially in solvent exposed side-chains. The MD simulation of Dynorphin was extended to 40 ns to study its behavior in an aqueous environment. This long simulation showed that the peptide has a tendency to form an alpha-helical structure in water, but the stabilization free energy is too weak, resulting in frequent interconversions between random and helical conformations during the simulation time. The results reported here suggest that the SCP implicit solvent model is adequate to describe electrostatic effects in MD simulation of both peptides and proteins using the same set of parameters. It is suggested that the present approach could form the basis for the development of a reliable and general continuum approach for use in molecular biology, and directions are outlined for attaining this long-term goal.     

Free energy landscape of protein folding in water: explicit vs. implicit solvent

The Generalized Born (GB) continuum solvent model is arguably the most widely used implicit solvent model in protein folding and protein structure prediction simulations; however, it still remains an open question on how well the model behaves in these large-scale simulations. The current study uses the beta-hairpin from C-terminus of protein G as an example to explore the folding free energy landscape with various GB models, and the results are compared to the explicit solvent simulations and experiments. All free energy landscapes are obtained from extensive conformation space sampling with a highly parallel replica exchange method. Because solvation model parameters are strongly coupled with force fields, five different force field/solvation model combinations are examined and compared in this study, namely the explicit solvent model: OPLSAA/SPC model, and the implicit solvent models: OPLSAA/SGB (Surface GB), AMBER94/GBSA (GB with Solvent Accessible Surface Area), AMBER96/GBSA, and AMBER99/GBSA. Surprisingly, we find that the free energy landscapes from implicit solvent models are quite different from that of the explicit solvent model. Except for AMBER96/GBSA, all other implicit solvent models find the lowest free energy state not the native state. All implicit solvent models show erroneous salt-bridge effects between charged residues, particularly in OPLSAA/SGB model, where the overly strong salt-bridge effect results in an overweighting of a non-native structure with one hydrophobic residue F52 expelled from the hydrophobic core in order to make better salt bridges. On the other hand, both AMBER94/GBSA and AMBER99/GBSA models turn the beta-hairpin in to an alpha-helix, and the alpha-helical content is much higher than the previously reported alpha-helices in an explicit solvent simulation with AMBER94 (AMBER94/TIP3P). Only AMBER96/GBSA shows a reasonable free energy landscape with the lowest free energy structure the native one despite an erroneous salt-bridge between D47 and K50. Detailed results on free energy contour maps, lowest free energy structures, distribution of native contacts, alpha-helical content during the folding process, NOE comparison with NMR, and temperature dependences are reported and discussed for all five models.     

Early events in protein aggregation: molecular flexibility and hydrophobicity/charge interaction in amyloid peptides as studied by molecular dynamics simulations

In a previous article (Zbilut et al., Biophys J 2003;85:3544-3557), we demonstrated how an aggregation versus folding choice could be approached considering hydrophobicity distribution and charge. In this work, our aim is highlighting the mutual interaction of charge and hydrophobicity distribution in the aggregation process. Use was made of two different peptides, both derived from a transmembrane protein (amyloid precursor protein; APP), namely, Abeta(1-28) and Abeta(1-40). Abeta(1-28) has a much lower aggregation propensity than Abeta(1-40). The results obtained by means of molecular dynamics simulations show that, when submitted to the most "aggregation-prone" environment, corresponding to the isoelectric point and consequently to zero net charge, both peptides acquire their maximum flexibility, but Abeta(1-40) has a definitely higher conformational mobility than Abeta(1-28). The absence of a hydrophobic "tail," which is the most mobile part of the molecule in Abeta(1-40), is the element lacking in Abeta(1-28) for obtaining a "fully aggregating" phenotype. Our results suggest that conformational flexibility, determined by both hydrophobicity and charge effect, is the main mechanistic determinant of aggregation propensity.     

Impact of the mutation A21G (Flemish variant) on Alzheimer's beta-amyloid dimers by molecular dynamics simulations

Soluble oligomers of the amyloid beta-protein (Abeta) are linked to Alzheimer's disease. Irrespective of the nature of the nucleus before fibril growth, dimers are essential species in Abeta assembly, but their transient character has precluded, thus far, high-resolution structure determination. We have investigated the effects of the point mutation A21G on Abeta dimers by performing high temperature all-atom molecular dynamics simulations of Abeta(40), Abeta(42), and their Flemish variants (A21G) starting from their fibrillar conformations. Abeta dimers are found in equilibrium between various topologies, and the absence of common structural features shared by the four species makes problematic the design of a unique inhibitor for blocking dimers. We also show that the impact of the point mutation A21G on Abeta structure and dynamics varies from Abeta(40) to Abeta(42). Finally, we provide a possible structural explanation for the reduced aggregation rate of Abeta fibrils containing the Flemish disease-causing mutation.     

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

The role of peptides Abeta40 and Abeta42 in the early pathogenesis of Alzheimer's disease (AD) is frequently emphasized in the literature. It is known that Abeta42 is more prone to aggregation than Abeta40, even though they differ in only two (IA) amino acid residues at the C-terminal end. A direct comparison of the ensembles of conformations adopted by the monomers in solution has been limited by the inherent flexibility of the unfolded peptides. Here, we characterize the conformations of Abeta40 and Abeta42 in water by using a combination of molecular dynamics (MD) and measured scalar (3)J(HNHalpha) data from NMR experiments. We perform replica exchange MD (REMD) simulations and find that classical forcefields reproduce the NMR data quantitatively when the sampling is extended to the microseconds time-scale. Using the quantitative agreement of the NMR data as a validation of the model, we proceed to compare the conformational ensembles of the Abeta40 and Abeta42 peptide monomers. Our analysis confirms the existence of structured regions within the otherwise flexible Abeta peptides. We find that the C terminus of Abeta42 is more structured than that of Abeta40. The formation of a beta-hairpin in the sequence (31)IIGLMVGGVVIA involving short strands at residues 31-34 and 38-41 (in bold) reduces the C-terminal flexibility of the Abeta42 peptide and may be responsible for the higher propensity of this peptide to form amyloids.     

Modeling loop reorganization free energies of acetylcholinesterase: a comparison of explicit and implicit solvent models

The treatment of hydration effects in protein dynamics simulations varies in model complexity and spans the range from the computationally intensive microscopic evaluation to simple dielectric screening of charge-charge interactions. This paper compares different solvent models applied to the problem of estimating the free-energy difference between two loop conformations in acetylcholinesterase. Molecular dynamics (MD) simulations were used to sample potential energy surfaces of the two basins with solvent treated by means of explicit and implicit methods. Implicit solvent methods studied include the generalized Born (GB) model, atomic solvation potential (ASP), and the distance-dependent dieletric constant. By using the linear response approximation (LRA), the explicit solvent calculations determined a free-energy difference that is in excellent agreement with the experimental estimate, while rescoring the protein conformations with GB or the Poisson equation showed inconsistent and inferior results. While the approach of rescoring conformations from explicit water simulations with implicit solvent models is popular among many applications, it perturbs the energy landscape by changing the solvent contribution to microstates without conformational relaxation, thus leading to non-optimal solvation free energies. Calculations applying MD with a GB solvent model produced results of comparable accuracy as observed with LRA, yet the electrostatic free-energy terms were significantly different due to optimization on a potential energy surface favored by an implicit solvent reaction field. The simpler methods of ASP and the distance-dependent scaling of the dielectric constant both produced considerable distortions in the protein internal free-energy terms and are consequently unreliable.     

Computational investigation on the effects of H50Q and G51D mutations on the α-Synuclein aggregation propensity

The aggregation of α-synuclein is linked directly to the histopathology of Parkinson’s disease (PD). However, several missense mutations present in the α-synuclein gene (SNCA) have been known to be associated with PD. Several studies have highlighted the effect of SNCA mutations on the α-synuclein aggregation, but their pathological roles are not completely established. In this study, we have focused on the effects of the recently discovered α-synuclein missense mutants (H50Q and G51D) on the aggregation using computational approaches. We performed all atom molecular dynamics (MD) simulation on these mutants and compared their conformational dynamics with Wild-Type (WT) α-synuclein. We noticed the solvent accessible surface area (SASA), radius of gyration, atomic fluctuations, and beta strand content to be higher in H50Q than G51D and WT. Using PDBSum online server; we analyzed the inter-molecular interactions that drive the association of monomeric units of H50Q, WT, and G51D in forming the respective homo-dimer. We noticed the interface area, number of interacting residues and binding free energy to be higher for H50Q homo-dimer than the WT and G51D homo-dimers. Our findings in this study suggest that in comparison to WT and G51D, H50Q mutation to have a positive effect on increasing the α-synuclein aggregation propensity. Hence, we see that H50Q and G51D mutation show conflicting effect on the aggregation propensity of α-synuclein.     

Prion Domain of Yeast Ure2 Protein Adopts a Completely Disordered Structure: A Solid-Support EPR Study

Amyloid fibril formation is associated with a range of neurodegenerative diseases in humans, including Alzheimer’s, Parkinson’s, and prion diseases. In yeast, amyloid underlies several non-Mendelian phenotypes referred to as yeast prions. Mechanism of amyloid formation is critical for a complete understanding of the yeast prion phenomenon and human amyloid-related diseases. Ure2 protein is the basis of yeast prion [URE3]. The Ure2p prion domain is largely disordered. Residual structures, if any, in the disordered region may play an important role in the aggregation process. Studies of Ure2p prion domain are complicated by its high aggregation propensity, which results in a mixture of monomer and aggregates in solution. Previously we have developed a solid-support electron paramagnetic resonance (EPR) approach to address this problem and have identified a structured state for the Alzheimer’s amyloid-β monomer. Here we use solid-support EPR to study the structure of Ure2p prion domain. EPR spectra of Ure2p prion domain with spin labels at every fifth residue from position 10 to position 75 show similar residue mobility profile for denaturing and native buffers after accounting for the effect of solution viscosity. These results suggest that Ure2p prion domain adopts a completely disordered structure in the native buffer. A completely disordered Ure2p prion domain implies that the amyloid formation of Ure2p, and likely other Q/N-rich yeast prion proteins, is primarily driven by inter-molecular interactions.     

Spontaneous beta-barrel formation: an all-atom Monte Carlo study of Abeta16-22 oligomerization

Using all-atom Monte Carlo simulations with implicit water, combined with a cluster size analysis, we study the aggregation of Abeta(16) (-22), a peptide capable of forming amyloid fibrils. We consider a system of six initially randomly oriented Abeta(16) (-22) peptides, and investigate the thermodynamics and structural properties of aggregates formed by this system. The system is unaggregated without ordered secondary structure at high temperature, and forms beta-sheet rich aggregates at low temperature. At the crossover between these two regimes, we find that clusters of all sizes occur, whereas the beta-strand content is low. In one of several runs, we observe the spontaneous formation of a beta-barrel with six antiparallel strands. The beta-barrel stands out as the by far most long-lived aggregate seen in our simulations.     

The structure of the Alzheimer amyloid beta 10-35 peptide probed through replica-exchange molecular dynamics simulations in explicit solvent

The conformational states sampled by the Alzheimer amyloid beta (10-35) (Abeta 10-35) peptide were probed using replica-exchange molecular dynamics (REMD) simulations in explicit solvent. The Abeta 10-35 peptide is a fragment of the full-length Abeta 40/42 peptide that possesses many of the amyloidogenic properties of its full-length counterpart. Under physiological temperature and pressure, our simulations reveal that the Abeta 10-35 peptide does not possess a single unique folded state. Rather, this peptide exists as a mixture of collapsed globular states that remain in rapid dynamic equilibrium with each other. This conformational ensemble is dominated by random coil and bend structures with insignificant presence of an alpha-helical or beta-sheet structure. The 3D structure of Abeta 10-35 is seen to be defined by a salt bridge formed between the side-chains of K28 and D23. This salt bridge is also observed in Abeta fibrils and our simulations suggest that monomeric conformations of Abeta 10-35 contain pre-folded structural motifs that promote rapid aggregation of this peptide.