The conformations of the amyloid-beta (21-30) fragment can be described by three families in solution

Secondary structure dependence of amyloid-β(1-40) on simulation techniques and force field parameters

Our recent studies revealed that none of the selected widely used force field parameters and molecular dynamics simulation techniques yield structural properties for the intrinsically disordered α-synuclein that are in agreement with various experiments via testing different force field parameters. Here, we extend our studies on the secondary structure properties of the disordered amyloid-β(1-40) peptide in aqueous solution. For these purposes, we conducted extensive replica exchange molecular dynamics simulations and obtained extensive molecular dynamics simulation trajectories from David E. Shaw group. Specifically, these molecular dynamics simulations were conducted using various force field parameters and obtained results are compared to our replica exchange molecular dynamics simulations and experiments. In this study, we calculated the secondary structure abundances and radius of gyration values for amyloid-β(1-40) that were simulated using varying force field parameter sets and different simulation techniques. In addition, the intrinsic disorder propensity, as well as sequence-based secondary structure predisposition of amyloid-β(1-40) and compared the findings with the results obtained from molecular simulations using various force field parameters and different simulation techniques. Our studies clearly show that the epitope region identification of amyloid-β(1-40) depends on the chosen simulation technique and chosen force field parameters. Based on comparison with experiments, we find that best computational results in agreement with experiments are obtained using the a99sb*-ildn, charmm36m, and a99sb-disp parameters for the amyloid-β(1-40) peptide in molecular dynamics simulations without parallel tempering or via replica exchange molecular dynamics simulations.

Experimentally Consistent Simulation of Aβ21-30 Peptides with a Minimal NMR Bias

Misfolded amyloid peptides are neurotoxic molecules associated with Alzheimer's disease. The Aβ_21-30 peptide fragment is a decapeptide fragment of the complete Aβ42 peptide which is a hypothesized cause of Alzheimer's disease via amyloid fibrillogenesis. Aβ_21-30 is investigated here with a combination of NMR (nuclear magnetic resonance) spectroscopy experiments and molecular dynamics simulations with experiment directed simulation (EDS). EDS is a maximum entropy biasing method that augments a molecular dynamics simulation with experimental data (NMR chemical shifts) to improve agreement with experiments and thus accuracy. EDS molecular dynamics shows that the Aβ_21-30 monomer has a β turn stabilized by the following interactions: S26-K28, D23-S26, and D23-K28. NMR, total correlation spectroscopy, and rotating frame Overhauser effect spectroscopy experiments provide independent agreement. Subsequent two- and four-monomer EDS simulations show aggregation. Diffusion coefficients calculated from molecular simulation also agreed with experimentally measured values only after using EDS, providing independent assessment of accuracy. This work demonstrates how accuracy can be improved by directly using experimental data in molecular dynamics of complex processes like self-assembly.

β-amyloid model core peptides: Effects of hydrophobes and disulfides

The mechanism by which a disordered peptide nucleates and forms amyloid is incompletely understood. A central domain of β-amyloid (Aβ21-30) has been proposed to have intrinsic structural propensities that guide the limited formation of structure in the process of fibrillization. In order to test this hypothesis, we examine several internal fragments of Aβ, and variants of these either cyclized or with an N-terminal Cys. While Aβ21-30 and variants were always monomeric and unstructured (circular dichroism (CD) and nuclear magnetic resonance spectroscopy (NMRS)), we found that the addition of flanking hydrophobic residues in Aβ16-34 led to formation of typical amyloid fibrils. NMR showed no long-range nuclear overhauser effect (nOes) in Aβ21-30, Aβ16-34, or their variants, however. Serial ¹ H-¹⁵ N-heteronuclear single quantum coherence spectroscopy, ¹ H-¹ H nuclear overhauser effect spectroscopy, and ¹ H-¹ H total correlational spectroscopy spectra were used to follow aggregation of Aβ16-34 and Cys-Aβ16-34 at a site-specific level. The addition of an N-terminal Cys residue (in Cys-Aβ16-34) increased the rate of fibrillization which was attributable to disulfide bond formation. We propose a scheme comparing the aggregation pathways for Aβ16-34 and Cys-Aβ16-34, according to which Cys-Aβ16-34 dimerizes, which accelerates fibril formation. In this context, cysteine residues form a focal point that guides fibrillization, a role which, in native peptides, can be assumed by heterogeneous nucleators of aggregation.

Insights into the Molecular Mechanisms of Alzheimer's and Parkinson's Diseases with Molecular Simulations: Understanding the Roles of Artificial and Pathological Missense Mutations in Intrinsically Disordered Proteins Related to Pathology

Amyloid-β and α-synuclein are intrinsically disordered proteins (IDPs), which are at the center of Alzheimer's and Parkinson's disease pathologies, respectively. These IDPs are extremely flexible and do not adopt stable structures. Furthermore, both amyloid-β and α-synuclein can form toxic oligomers, amyloid fibrils and other type of aggregates in Alzheimer's and Parkinson's diseases. Experimentalists face challenges in investigating the structures and thermodynamic properties of these IDPs in their monomeric and oligomeric forms due to the rapid conformational changes, fast aggregation processes and strong solvent effects. Classical molecular dynamics simulations complement experiments and provide structural information at the atomic level with dynamics without facing the same experimental limitations. Artificial missense mutations are employed experimentally and computationally for providing insights into the structure-function relationships of amyloid-β and α-synuclein in relation to the pathologies of Alzheimer's and Parkinson's diseases. Furthermore, there are several natural genetic variations that play a role in the pathogenesis of familial cases of Alzheimer's and Parkinson's diseases, which are related to specific genetic defects inherited in dominant or recessive patterns. The present review summarizes the current understanding of monomeric and oligomeric forms of amyloid-β and α-synuclein, as well as the impacts of artificial and pathological missense mutations on the structural ensembles of these IDPs using molecular dynamics simulations. We also emphasize the recent investigations on residual secondary structure formation in dynamic conformational ensembles of amyloid-β and α-synuclein, such as β-structure linked to the oligomerization and fibrillation mechanisms related to the pathologies of Alzheimer's and Parkinson's diseases. This information represents an important foundation for the successful and efficient drug design studies.

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.

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.

Force-Field Induced Bias in the Structure of Aβ21-30: A Comparison of OPLS, AMBER, CHARMM, and GROMOS Force Fields

In this work we examine the dynamics of an intrinsically disordered protein fragment of the amyloid β, the Aβ21-30, under seven commonly used molecular dynamics force fields (OPLS-AA, CHARMM27-CMAP, AMBER99, AMBER99SB, AMBER99SB-ILDN, AMBER03, and GROMOS53A6), and three water models (TIP3P, TIP4P, and SPC/E). We find that the tested force fields and water models have little effect on the measures of radii of gyration and solvent accessible surface area (SASA); however, secondary structure measures and intrapeptide hydrogen-bonding are significantly modified, with AMBER (99, 99SB, 99SB-ILDN, and 03) and CHARMM22/27 force-fields readily increasing helical content and the variety of intrapeptide hydrogen bonds. On the basis of a comparison between the population of helical and β structures found in experiments, our data suggest that force fields that suppress the formation of helical structure might be a better choice to model the Aβ21-30 peptide.

The peculiar role of the A2V mutation in amyloid-β (Aβ) 1-42 molecular assembly

We recently reported a novel Aβ precursor protein mutation (A673V), corresponding to position 2 of Aβ1-42 peptides (Aβ1-42A2V), that caused an early onset AD-type dementia in a homozygous individual. The heterozygous relatives were not affected as an indication of autosomal recessive inheritance of this mutation. We investigated the folding kinetics of native unfolded Aβ1-42A2V in comparison with the wild type sequence (Aβ1-42WT) and the equimolar solution of both peptides (Aβ1-42MIX) to characterize the oligomers that are produced in the early phases. We carried out the structural characterization of the three preparations using electron and atomic force microscopy, fluorescence emission, and x-ray diffraction and described the soluble oligomer formation kinetics by laser light scattering. The mutation promoted a peculiar pathway of oligomerization, forming a connected system similar to a polymer network with hydrophobic residues on the external surface. Aβ1-42MIX generated assemblies very similar to those produced by Aβ1-42WT, albeit with slower kinetics due to the difficulties of Aβ1-42WT and Aβ1-42A2V peptides in building up of stable intermolecular interaction.

The OPEP protein model: from single molecules, amyloid formation, crowding and hydrodynamics to DNA/RNA systems

The OPEP coarse-grained protein model has been applied to a wide range of applications since its first release 15 years ago. The model, which combines energetic and structural accuracy and chemical specificity, allows the study of single protein properties, DNA-RNA complexes, amyloid fibril formation and protein suspensions in a crowded environment. Here we first review the current state of the model and the most exciting applications using advanced conformational sampling methods. We then present the current limitations and a perspective on the ongoing developments.

Gly25-Ser26 amyloid β-protein structural isomorphs produce distinct Aβ42 conformational dynamics and assembly characteristics

One of the earliest events in amyloid β-protein (Aβ) self-association is nucleation of Aβ monomer folding through formation of a turn at Gly25-Lys28. We report here the effects of structural changes at the center of the turn, Gly25-Ser26, on Aβ42 conformational dynamics and assembly. We used "click peptide" chemistry to quasi-synchronously create Aβ42 from 26-O-acyliso-Aβ42 (iAβ42) through a pH jump from 3 to 7.4. We also synthesized Nα-acetyl-Ser26-iAβ42 (Ac-iAβ42), which cannot undergo O→N acyl chemistry, to study the behavior of this ester form of Aβ42 itself at neutral pH. Data from experiments monitoring increases in β-sheet formation (thioflavin T, CD), hydrodynamic radius (RH), scattering intensity (quasielastic light scattering spectroscopy), and extent of oligomerization (ion mobility spectroscopy-mass spectrometry) were quite consistent. A rank order of Ac-iAβ42>iAβ42>Aβ42 was observed. Photochemically cross-linked iAβ42 displayed an oligomer distribution with a prominent dimer band that was not present with Aβ42. These dimers also were observed selectively in iAβ42 in ion mobility spectrometry experiments. The distinct biophysical behaviors of iAβ42 and Aβ42 appear to be due to the conversion of iAβ42 into "pure" Aβ42 monomer, a nascent form of Aβ42 that does not comprise the variety of oligomeric and aggregated states present in pre-existent Aβ42. These results emphasize the importance of the Gly25-Ser26 dipeptide in organizing Aβ42 monomer structure and thus suggest that drugs altering the interactions of this dipeptide with neighboring side-chain atoms or with the peptide backbone could be useful in therapeutic strategies targeting formation of Aβ oligomers and higher-order assemblies.

Spontaneous dimer states of the Aβ21–30decapeptide

Computational examination of the spontaneous dimerization of Aβ_21–30and stability measures of the resulting parallel and anti-parallel aligned dimers.

Aβ monomers transiently sample oligomer and fibril-like configurations: ensemble characterization using a combined MD/NMR approach

Amyloid β (Aβ) peptides are a primary component of fibrils and oligomers implicated in the etiology of Alzheimer's disease (AD). However, the intrinsic flexibility of these peptides has frustrated efforts to investigate the secondary and tertiary structure of Aβ monomers, whose conformational landscapes directly contribute to the kinetics and thermodynamics of Aβ aggregation. In this work, de novo replica exchange molecular dynamics (REMD) simulations on the microseconds-per-replica timescale are used to characterize the structural ensembles of Aβ42, Aβ40, and M35-oxidized Aβ42, three physiologically relevant isoforms with substantially different aggregation properties. J-coupling data calculated from the REMD trajectories were compared to corresponding NMR-derived values acquired through two different pulse sequences, revealing that all simulations converge on the order of hundreds of nanoseconds-per-replica toward ensembles that yield good agreement with experiment. Though all three Aβ species adopt highly heterogeneous ensembles, these are considerably more structured compared to simulations on shorter timescales. Prominent in the C-terminus are antiparallel β-hairpins between L17-A21, A30-L36, and V39-I41, similar to oligomer and fibril intrapeptide models that expose these hydrophobic side chains to solvent and may serve as hotspots for self-association. Compared to reduced Aβ42, the absence of a second β-hairpin in Aβ40 and the sampling of alternate β topologies by M35-oxidized Aβ42 may explain the reduced aggregation rates of these forms. A persistent V24-K28 bend motif, observed in all three species, is stabilized by buried backbone to side-chain hydrogen bonds with D23 and a cross-region salt bridge between E22 and K28, highlighting the role of the familial AD-linked E22 and D23 residues in Aβ monomer folding. These characterizations help illustrate the conformational landscapes of Aβ monomers at atomic resolution and provide insight into the early stages of Aβ aggregation pathways.

Structural, thermodynamical, and dynamical properties of oligomers formed by the amyloid NNQQ peptide: insights from coarse-grained simulations

Characterizing the early formed oligomeric intermediates of amyloid peptides is of particular interest due to their links with neurodegenerative diseases. Here we study the NNQQ peptide, known to display parallel β-strands in amyloid fibrils by x-ray microcrystallography, and investigate the structural, thermodynamical, and dynamical properties of 20 NNQQ peptides using molecular dynamics and replica exchange molecular dynamics simulations coupled to a coarse-grained force field. All simulations are initiated from randomized and fully dispersed monomeric conformations. Our simulations reveal that the phase transition is characterized by a change in the oligomer and β-sheet size distributions and the percentage of mixed parallel/antiparallel β-strands when the sheets are formed. At all temperatures, however, the fraction of parallel β-strands remains low, though there are many association/fragmentation events. This work and a growing body of computational studies provide strong evidence that the critical nucleus goes beyond 20 chains and reordering of the β-strands occurs in larger oligomers.

C-terminal turn stability determines assembly differences between Aβ40 and Aβ42

Oligomerization of the amyloid β-protein (Aβ) is a seminal event in Alzheimer's disease. Aβ42, which is only two amino acids longer than Aβ40, is particularly pathogenic. Why this is so has not been elucidated fully. We report here results of computational and experimental studies revealing a C-terminal turn at Val36-Gly37 in Aβ42 that is not present in Aβ40. The dihedral angles of residues 36 and 37 in an Ile31-Ala42 peptide were consistent with β-turns, and a β-hairpin-like structure was indeed observed that was stabilized by hydrogen bonds and by hydrophobic interactions between residues 31-35 and residues 38-42. In contrast, Aβ(31-40) mainly existed as a statistical coil. To study the system experimentally, we chemically synthesized Aβ peptides containing amino acid substitutions designed to stabilize or destabilize the hairpin. The triple substitution Gly33Val-Val36Pro-Gly38Val ("VPV") facilitated Aβ42 hexamer and nonamer formation, while inhibiting formation of classical amyloid-type fibrils. These assemblies were as toxic as were assemblies from wild-type Aβ42. When substituted into Aβ40, the VPV substitution caused the peptide to oligomerize similarly to Aβ42. The modified Aβ40 was significantly more toxic than Aβ40. The double substitution d-Pro36-l-Pro37 abolished hexamer and dodecamer formation by Aβ42 and produced an oligomer size distribution similar to that of Aβ40. Our data suggest that the Val36-Gly37 turn could be the sine qua non of Aβ42. If true, this structure would be an exceptionally important therapeutic target.

A key role for lysine residues in amyloid β-protein folding, assembly, and toxicity

A combination of hydrophobic and electrostatic interactions is important in initiating the aberrant self-assembly process that leads to formation of toxic oligomers and aggregates by multiple disease-related proteins, including amyloid β-protein (Aβ), whose self-assembly is believed to initiate brain pathogenesis in Alzheimer's disease. Lys residues play key roles in this process and participate in both types of interaction. They also are the target of our recently reported molecular tweezer inhibitors. To obtain further insight into the role of the two Lys residues in Aβ assembly and toxicity, here we substituted each by Ala in both Aβ40 and Aβ42 and studied the impact of the substitution on Aβ oligomerization, aggregation, and toxicity. Our data show that each substitution has a major impact on Aβ assembly and toxicity, with significant differences depending on peptide length (40 versus 42 amino acids) and the position of the substitution. In particular, Lys16→Ala substitution dramatically reduces Aβ toxicity. The data support the use of compounds targeting Lys residues specifically as inhibitors of Aβ toxicity and suggest that exploring the role of Lys residues in other disease-related amyloidogenic proteins may help understanding the mechanisms of aggregation and toxicity of these proteins.

Dynamics of metastable β-hairpin structures in the folding nucleus of amyloid β-protein

The amyloid β-protein (Aβ), which is present predominately as a 40- or 42-residue peptide, is postulated to play a seminal role in the pathogenesis of Alzheimer's disease (AD). Folding of the Aβ(21-30) decapeptide region is a critical step in the aggregation of Aβ. We report results of constant temperature all-atom molecular dynamics simulations in explicit water of the dynamics of monomeric Aβ(21-30) and its Dutch [Glu22Gln], Arctic [Glu22Gly], and Iowa [Asp23Asn] isoforms that are associated with familial forms of cerebral amyloid angiopathy and AD. The simulations revealed a variety of loop conformers that exhibited a hydrogen bond network involving the Asp23 and Ser26 amino acids. A population of conformers, not part of the loop population, was found to form metastable β-hairpin structures with the highest probability in the Iowa mutant. At least three β-hairpin structures were found that differed in their hydrogen bonding register, average number of backbone hydrogen bonds, and lifetimes. Analysis revealed that the Dutch mutant had the longest β-hairpin lifetime (≥500 ns), closely followed by the Iowa mutant (≈500 ns). Aβ(21-30) and the Arctic mutant had significantly lower lifetimes (≈200 ns). Hydrophobic packing of side chains was responsible for enhanced β-hairpin lifetimes in the Dutch and Iowa mutants, whereas lifetimes in Aβ(21-30) and its Arctic mutant were influenced by the backbone hydrogen bonding. The data suggest that prolonged β-hairpin lifetimes may impact peptide pathogenicity in vivo.

Large loop conformation sampling using the activation relaxation technique, ART-nouveau method

We present an adaptation of the ART-nouveau energy surface sampling method to the problem of loop structure prediction. This method, previously used to study protein folding pathways and peptide aggregation, is well suited to the problem of sampling the conformation space of large loops by targeting probable folding pathways instead of sampling exhaustively that space. The number of sampled conformations needed by ART nouveau to find the global energy minimum for a loop was found to scale linearly with the sequence length of the loop for loops between 8 and about 20 amino acids. Considering the linear scaling dependence of the computation cost on the loop sequence length for sampling new conformations, we estimate the total computational cost of sampling larger loops to scale quadratically compared to the exponential scaling of exhaustive search methods.

Structural dynamics of the amyloid β-protein monomer folding nucleus

Alzheimer's disease (AD) is linked to the aberrant assembly of the amyloid β-protein (Aβ). The (21)AEDVGSNKGA(30) segment, Aβ(21-30), forms a turn that acts as a monomer folding nucleus. Amino acid substitutions within this nucleus cause familial forms of AD. To determine the biophysical characteristics of the folding nucleus, we studied the biologically relevant acetyl-Aβ(21-30)-amide peptide using experimental techniques (limited proteolysis, thermal denaturation, urea denaturation followed by pulse proteolysis, and electron microscopy) and computational methods (molecular dynamics). Our results reveal a highly stable foldon and suggest new strategies for therapeutic drug development.

Distinct dimerization for various alloforms of the amyloid-beta protein: Aβ(1-40), Aβ(1-42), and Aβ(1-40)(D23N)

The Amyloid-beta protein is related to Alzheimer's disease, and various experiments have shown that oligomers as small as the dimer are cytotoxic. Two alloforms are mainly produced: Aβ(1-40) and Aβ(1-42). They have very different oligomer distributions, and it was recently suggested, from experimental studies, that this variation may originate from structural differences in their dimer structures. Little structural information is available on the Aβ dimer, however, and to complement experimental observations, we simulated the folding of the wild-type Aβ(1-40) and Aβ(1-42) dimers as well as the mutated Aβ(1-40)(D23N) dimer using an accurate coarse-grained force field coupled to Hamiltonian-temperature replica exchange molecular dynamics. The D23N variant impedes the salt-bridge formation between D23 and K28 seen in the wild-type Aβ, leading to very different fibrillation properties and final amyloid fibrils. Our results show that the Aβ(1-42) dimer has a higher propensity than the Aβ(1-40) dimer to form β-strands at the central hydrophobic core (residues 17-21) and at the C-terminal (residues 30-42), which are two segments crucial to the oligomerization of Aβ. The free energy landscape of the Aβ(1-42) dimer is also broader and more complex than that of the Aβ(1-40) dimer. Interestingly, D23N also impacts the free energy landscape by increasing the population of configurations with higher β-strand propensities when compared against Aβ(40). In addition, while Aβ(1-40)(D23N) displays a higher β-strand propensity at the C-terminal, its solvent accessibility does not change with respect to the wild-type sequence. Overall, our results show the strong impact of the two amino acids Ile41-Ala42 and the salt-bridge D23-K28 on the folding of the Aβ dimer.

Amyloid-β peptide structure in aqueous solution varies with fragment size

Various fragment sizes of the amyloid-β (Aβ) peptide have been utilized to mimic the properties of the full-length Aβ peptide in solution. Among these smaller fragments, Aβ16 and Aβ28 have been investigated extensively. In this work, we report the structural and thermodynamic properties of the Aβ16, Aβ28, and Aβ42 peptides in an aqueous solution environment. We performed replica exchange molecular dynamics simulations along with thermodynamic calculations for investigating the conformational free energies, secondary and tertiary structures of the Aβ16, Aβ28, and Aβ42 peptides. The results show that the thermodynamic properties vary from each other for these peptides. Furthermore, the secondary structures in the Asp1-Lys16 and Asp1-Lys28 regions of Aβ42 cannot be completely captured by the Aβ16 and Aβ28 fragments. For example, the β-sheet structures in the N-terminal region of Aβ16 and Aβ28 are either not present or the abundance is significantly decreased in Aβ42. The α-helix and β-sheet abundances in Aβ28 and Aβ42 show trends--to some extent--with the potential of mean forces but no such trend could be obtained for Aβ16. Interestingly, Arg5 forms salt bridges with large abundances in all three peptides. The formation of a salt bridge between Asp23-Lys28 is more preferred over the Glu22-Lys28 salt bridge in Aβ28 but this trend is vice versa for Aβ42. This study shows that the Asp1-Lys16 and Asp1-Lys28 regions of the full length Aβ42 peptide cannot be completely mimicked by studying the Aβ16 and Aβ28 peptides.

Structures of Aβ17-42 trimers in isolation and with five small-molecule drugs using a hierarchical computational procedure

The amyloid-β protein (Aβ) oligomers are believed to be the main culprits in the cytoxicity of Alzheimer's disease (AD) and p3 peptides (Aβ17-42 fragments) are present in AD amyloid plaques. Many small-molecule or peptide-based inhibitors are known to slow down Aβ aggregation and reduce the toxicity in vitro, but their exact modes of action remain to be determined since there has been no atomic level of Aβ(p3)-drug oligomers. In this study, we have determined the structure of Aβ17-42 trimers both in aqueous solution and in the presence of five small-molecule inhibitors using a multiscale computational study. These inhibitors include 2002-H20, curcumin, EGCG, Nqtrp, and resveratrol. First, we used replica exchange molecular dynamics simulations coupled to the coarse-grained (CG) OPEP force field. These CG simulations reveal that the conformational ensemble of Aβ17-42 trimer can be described by 14 clusters with each peptide essentially adopting turn/random coil configurations, although the most populated cluster is characterized by one peptide with a β-hairpin at Phe19-Leu31. Second, these 14 dominant clusters and the less-frequent fibril-like state with parallel register of the peptides were subjected to atomistic Autodock simulations. Our analysis reveals that the drugs have multiple binding modes with different binding affinities for trimeric Aβ17-42 although they interact preferentially with the CHC region (residues 17-21). The compounds 2002-H20 and Nqtrp are found to be the worst and best binders, respectively, suggesting that the drugs may interfere at different stages of Aβ oligomerization. Finally, explicit solvent molecular dynamics of two predicted Nqtrp-Aβ17-42 conformations describe at atomic level some possible modes of action for Nqtrp.

Flexibility and binding affinity in protein-ligand, protein-protein and multi-component protein interactions: limitations of current computational approaches

The recognition process between a protein and a partner represents a significant theoretical challenge. In silico structure-based drug design carried out with nothing more than the three-dimensional structure of the protein has led to the introduction of many compounds into clinical trials and numerous drug approvals. Central to guiding the discovery process is to recognize active among non-active compounds. While large-scale computer simulations of compounds taken from a library (virtual screening) or designed de novo are highly desirable in the post-genomic area, many technical problems remain to be adequately addressed. This article presents an overview and discusses the limits of current computational methods for predicting the correct binding pose and accurate binding affinity. It also presents the performances of the most popular algorithms for exploring binary and multi-body protein interactions.

A multiscale approach to characterize the early aggregation steps of the amyloid-forming peptide GNNQQNY from the yeast prion sup-35

The self-organization of peptides into amyloidogenic oligomers is one of the key events for a wide range of molecular and degenerative diseases. Atomic-resolution characterization of the mechanisms responsible for the aggregation process and the resulting structures is thus a necessary step to improve our understanding of the determinants of these pathologies. To address this issue, we combine the accelerated sampling properties of replica exchange molecular dynamics simulations based on the OPEP coarse-grained potential with the atomic resolution description of interactions provided by all-atom MD simulations, and investigate the oligomerization process of the GNNQQNY for three system sizes: 3-mers, 12-mers and 20-mers. Results for our integrated simulations show a rich variety of structural arrangements for aggregates of all sizes. Elongated fibril-like structures can form transiently in the 20-mer case, but they are not stable and easily interconvert in more globular and disordered forms. Our extensive characterization of the intermediate structures and their physico-chemical determinants points to a high degree of polymorphism for the GNNQQNY sequence that can be reflected at the macroscopic scale. Detailed mechanisms and structures that underlie amyloid aggregation are also provided.

The formation of amyloid fibrils is associated with many neurodegenerative diseases such as Alzheimer's, Creutzfeld-Jakob, Parkinson's, the Prion disease and diabetes mellitus. In all cases, proteins misfold to form highly ordered insoluble aggregates called amyloid fibrils that deposit intra- and extracellularly and are resistant to proteases. All proteins are believed to have the instrinsic capability of forming amyloid fibrils that share common specific structural properties that have been observed by X-ray crystallography and by NMR. However, little is known about the aggregation dynamics of amyloid assemblies, and their toxicity mechanism is therefore poorly understood. It is believed that small amyloid oligomers, formed on the aggregation pathway of full amyloid fibrils, are the toxic species. A detailed atomic characterization of the oligomerization process is thus necessary to further our understanding of the amyloid oligomer's toxicity. Our approach here is to study the aggregation dynamics of a 7-residue amyloid peptide GNNQQNY through a combination of numerical techniques. Our results suggest that this amyloid sequence can form fibril-like structures and is polymorphic, which agrees with recent experimental observations. The ability to fully characterize and describe the aggregation pathway of amyloid sequences numerically is key to the development of future drugs to target amyloid oligomers.

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

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

Effects of G33A and G33I mutations on the structures of monomer and dimer of the amyloid-β fragment 29-42 by replica exchange molecular dynamics simulations

The early formed oligomers of amyloid-β proteins with 40 and 42 amino acids are believed to be the culprits of Alzheimer's disease. Aβ1-42 peptides with alanine and isoleucine mutations of glycine 33 are known to be much less toxic than the wild-type Aβ1-42 and promote the aggregation process in vitro. The fragment Aβ29-42 has also been shown to form fibrils, disrupt Aβ1-42 oligomerization, and inhibit Aβ1-42-induced neurotoxicity. As a first step toward understanding the impact of G33A and G33I mutations on the earliest steps along the Aβ1-42 aggregation pathway, we have studied the structures of the monomer and dimer of Aβ29-42 and its two G33 variants using coarse-grained replica exchange molecular dynamics simulations. These simulations, totaling 15 μs, indicate that both substitutions impact the conformational ensemble of Aβ29-42. For the monomer, the population of the β-hairpin is high for wild-type Aβ29-42, but marginal for Aβ29-42 G33I mutant. The three dimers are also stabilized by different patterns of interaction. The data are discussed in terms of the differences in the aggregation characteristics between wild-type Aβ1-42 and its two G33A and G33I variants.

Principles governing oligomer formation in amyloidogenic peptides

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

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

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

Structural diversity of the soluble trimers of the human amylin(20-29) peptide revealed by molecular dynamics simulations

The human islet amyloid polypeptide (hIAPP) or amylin is a 37-residue hormone found as amyloid deposits in pancreatic extracts of nearly all type 2 diabetes patients. The fragment 20-29 of sequence SNNFGAILSS (hIAPP20-29) has been shown to be responsible for the amyloidogenic propensities of the full length protein. Various polymorphic forms of hIAPP20-29 fibrils were described by using Fourier transform infrared (FTIR) and solid-state NMR experiments: unseeded hIAPP20-29 fibril with out-of-register antiparallel beta-strands, and two forms of seeded hIAPP20-29 fibril, with in-register antiparallel or in-register parallel beta-strands. As a first step toward understanding this polymorphism, we explore the equilibrium structures of the soluble hIAPP20-29 trimer, using multiple molecular dynamics (MD) simulations with the Optimized Potential for Efficient structure Prediction (OPEP) coarse-grained implicit solvent force field for a total length of 3.2 micros. Although, the trimer is found mainly random coil, consistent with the signal measured experimentally during the lag phase of hIAPP20-29 fibril formation, the central FGAIL residues have a relative high propensity to form interpeptide beta-sheets and antiparallel beta-strands are more probable than parallel beta-strands. One MD-predicted out-of-register antiparallel three-stranded beta-sheet matches exactly the FTIR-derived unseeded hIAPP20-29 fibril model. Our simulations, however, do not reveal any evidence of in-register parallel or in-register antiparallel beta-sheets as reported for seeded hIAPP20-29 fibrils. All these results indicate that fibril polymorphism is partially encoded in a trimer.

A coarse-grained protein force field for folding and structure prediction

We have revisited the protein coarse-grained optimized potential for efficient structure prediction (OPEP). The training and validation sets consist of 13 and 16 protein targets. Because optimization depends on details of how the ensemble of decoys is sampled, trial conformations are generated by molecular dynamics, threading, greedy, and Monte Carlo simulations, or taken from publicly available databases. The OPEP parameters are varied by a genetic algorithm using a scoring function which requires that the native structure has the lowest energy, and the native-like structures have energy higher than the native structure but lower than the remote conformations. Overall, we find that OPEP correctly identifies 24 native or native-like states for 29 targets and has very similar capability to the all-atom discrete optimized protein energy model (DOPE), found recently to outperform five currently used energy models.

Replica exchange molecular dynamics simulations of coarse-grained proteins in implicit solvent

Current approaches aimed at determining the free energy surface of all-atom medium-size proteins in explicit solvent are slow and are not sufficient to converge to equilibrium properties. To ensure a proper sampling of the configurational space, it is preferable to use reduced representations such as implicit solvent and/or coarse-grained protein models, which are much lighter computationally. Each model must be verified, however, to ensure that it can recover experimental structures and thermodynamics. Here we test the coarse-grained implicit solvent OPEP model with replica exchange molecular dynamics (REMD) on six peptides ranging in length from 10 to 28 residues: two alanine-based peptides, the second beta-hairpin from protein G, the Trp-cage and zinc-finger motif, and a dimer of a coiled coil peptide. We show that REMD-OPEP recovers the proper thermodynamics of the systems studied, with accurate structural description of the beta-hairpin and Trp-cage peptides (within 1-2 A from experiments). The light computational burden of REMD-OPEP, which enables us to generate many hundred nanoseconds at each temperature and fully assess convergence to equilibrium ensemble, opens the door to the determination of the free energy surface of larger proteins and assemblies.

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

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

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

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

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

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

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

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

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

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

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

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