Differences in β-strand populations of monomeric Aβ40 and Aβ42

Effect of a Paramagnetic Spin Label on the Intrinsically Disordered Peptide Ensemble of Amyloid-β

Paramagnetic relaxation enhancement is an NMR technique that has yielded important insight into the structure of folded proteins, although the perturbation introduced by the large spin probe might be thought to diminish its usefulness when applied to characterizing the structural ensembles of intrinsically disordered proteins (IDPs). We compare the computationally generated structural ensembles of the IDP amyloid-β42 (Aβ42) to an alternative sequence in which a nitroxide spin label attached to cysteine has been introduced at its N-terminus. Based on this internally consistent computational comparison, we find that the spin label does not perturb the signature population of the β-hairpin formed by residues 16-21 and 29-36 that is dominant in the Aβ42 reference ensemble. However, the presence of the tag induces a strong population shift in a subset of the original Aβ42 structural sub-populations, including a sevenfold enhancement of the β-hairpin formed by residues 27-31 and 33-38. Through back-calculation of NMR observables from the computational structural ensembles, we show that the structural differences between the labeled and unlabeled peptide would be evident in local residual dipolar couplings, and possibly differences in homonuclear ¹H-¹H nuclear Overhauser effects (NOEs) and heteronuclear ¹H-¹⁵N NOEs if the paramagnetic contribution to the longitudinal relaxation does not suppress the NOE intensities in the real experiment. This work shows that molecular simulation provides a complementary approach to resolving the potential structural perturbations introduced by reporter tags that can aid in the interpretation of paramagnetic relaxation enhancement, double electron-electron resonance, and fluorescence resonance energy transfer experiments applied to IDPs.

Comparing the Aggregation Free Energy Landscapes of Amyloid Beta(1-42) and Amyloid Beta(1-40)

Using a predictive coarse-grained protein force field, we compute and compare the free energy landscapes and relative stabilities of amyloid-β protein (1-42) and amyloid-β protein (1-40) in their monomeric and oligomeric forms up to the octamer. At the same concentration, the aggregation free energy profile of Aβ42 is more downhill, with a computed solubility that is about 10 times smaller than that of Aβ40. At a concentration of 40 μM, the clear free energy barrier between the pre-fibrillar tetramer form and the fibrillar pentamer in the Aβ40 aggregation landscape disappears for Aβ42, suggesting that the Aβ42 tetramer has a more diverse structural range. To further compare the landscapes, we develop a cluster analysis based on the structural similarity between configurations and use it to construct an oligomerization map that captures the paths of easy interconversion between different but structurally similar states of oligomers for both species. A taxonomy of the oligomer species based on β-sheet stacking topologies is proposed. The comparison of the two oligomerization maps highlights several key differences in the landscapes that can be attributed to the two additional C-terminal residues that Aβ40 lacks. In general, the two terminal residues strongly stabilize the oligomeric structures for Aβ42 relative to Aβ40, and greatly facilitate the conversion from pre-fibrillar trimers to fibrillar tetramers.

Self-Assembly Pathways of β-Sheet-Rich Amyloid-β(1-40) Dimers: Markov State Model Analysis on Millisecond Hybrid-Resolution Simulations

Early oligomerization during amyloid-β (Aβ) aggregation is essential for Aβ neurotoxicity. Understanding how unstructured Aβs assemble into oligomers, especially those rich in β-sheets, is essential but remains challenging as the assembly process is too transient for experimental characterization and too slow for molecular dynamics simulations. So far, atomic simulations are limited only to studies of either oligomer structures or assembly pathways for short Aβ segments. To overcome the computational challenge, we combine in this study a hybrid-resolution model and adaptive sampling techniques to perform over 2.7 ms of simulations of formation of full-length Aβ40 dimers that are the earliest toxic oligomeric species. The Markov state model is further employed to characterize the transition pathways and associated kinetics. Our results show that for two major forms of β-sheet-rich structures reported experimentally, the corresponding assembly mechanisms are markedly different. Hairpin-containing structures are formed by direct binding of soluble Aβ in β-hairpin-like conformations. Formation of parallel, in-register structures resembling fibrils occurs ∼100-fold more slowly and involves a rapid encounter of Aβ in arbitrary conformations followed by a slow structural conversion. The structural conversion proceeds via diverse pathways but always requires transient unfolding of encounter complexes. We find that the transition kinetics could be affected differently by intra-/intermolecular interactions involving individual residues in a conformation-dependent manner. In particular, the interactions involving Aβ's N-terminal part promote the assembly into hairpin-containing structures but delay the formation of fibril-like structures, thus explaining puzzling observations reported previously regarding the roles of this region in the early assembly process.

Bexarotene Does Not Clear Amyloid Beta Plaques but Delays Fibril Growth: Molecular Mechanisms

In 2012, it was reported that anticancer drug bexarotene reduced amyloid plaque and improved mental functioning in a small sample of mice engineered to exhibit Alzheimer's like symptoms. It has been suggested that bexarotene stimulates expression of apolipoprotein E (ApoE) leading to intracellular clearance of amyloid beta (Aβ). However, the effect of bexarotene on clearance of plaques has not been seen in some mouse models. Two interesting questions include whether bexarotene can destroy Aβ fibrils via direct interaction with them and how this compound impacts the lag phase in the fibril growth process. By the Thioflavin T fluorescence assay and atomic force microscopy, we have shown that bexarotene prolongs the lag phase, but it does not degrade Aβ fibrils. The impotence of bexarotene in destroying fibrils means that this compound is weakly bound to Aβ. On the other hand, the weak binding would prevent bexarotene from prolonging the lag phase. Thus, our two main in vitro observations seem to contradict each other. In order to settle this problem at the atomic level, we have performed all-atom molecular dynamics simulations in explicit water. We have demonstrated that bexarotene is not capable to reduce amyloid deposits due to weak binding to Aβ fibrils. However, it delays the self-assembly through reduction of the β-content of Aβ monomers at high enough ligand concentrations. Bexarotene is the first compound which displays such an unusual behavior. We have also shown that bexarotene has a low binding propensity to Aβ monomer and dimer.

Protein folding, misfolding and aggregation: The importance of two-electron stabilizing interactions

Proteins associated with neurodegenerative diseases are highly pleiomorphic and may adopt an all-α-helical fold in one environment, assemble into all-β-sheet or collapse into a coil in another, and rapidly polymerize in yet another one via divergent aggregation pathways that yield broad diversity of aggregates’ morphology. A thorough understanding of this behaviour may be necessary to develop a treatment for Alzheimer’s and related disorders. Unfortunately, our present comprehension of folding and misfolding is limited for want of a physicochemical theory of protein secondary and tertiary structure. Here we demonstrate that electronic configuration and hyperconjugation of the peptide amide bonds ought to be taken into account to advance such a theory. To capture the effect of polarization of peptide linkages on conformational and H-bonding propensity of the polypeptide backbone, we introduce a function of shielding tensors of the C^α atoms. Carrying no information about side chain-side chain interactions, this function nonetheless identifies basic features of the secondary and tertiary structure, establishes sequence correlates of the metamorphic and pH-driven equilibria, relates binding affinities and folding rate constants to secondary structure preferences, and manifests common patterns of backbone density distribution in amyloidogenic regions of Alzheimer’s amyloid β and tau, Parkinson’s α-synuclein and prions. Based on those findings, a split-intein like mechanism of molecular recognition is proposed to underlie dimerization of Aβ, tau, αS and PrP^C, and divergent pathways for subsequent association of dimers are outlined; a related mechanism is proposed to underlie formation of PrP^Sc fibrils. The model does account for: (i) structural features of paranuclei, off-pathway oligomers, non-fibrillar aggregates and fibrils; (ii) effects of incubation conditions, point mutations, isoform lengths, small-molecule assembly modulators and chirality of solid-liquid interface on the rate and morphology of aggregation; (iii) fibril-surface catalysis of secondary nucleation; and (iv) self-propagation of infectious strains of mammalian prions.

Emergence of Alternative Structures in Amyloid Beta 1-42 Monomeric Landscape by N-terminal Hexapeptide Amyloid Inhibitors

Alzheimer’s disease (AD) is characterized by deposition of amyloid beta (Aβ) peptides into senile plaques in the brain. While most familial mutations are associated with early-onset AD, recent studies report the AD-protective nature of two genetic human Aβ variants, i.e. A2T and A2V, in the heterozygous state. The mixture of A2V Aβ1-6 (Aβ₆) hexapeptide and WT Aβ1–42 (Αβ₄₂) is also found neuroprotective. Motivated by these findings, in this study we investigate the effects of WT, A2V, and A2T Aβ₆ hexapeptide binding on the monomeric WT Aβ₄₂ landscape. For this purpose, we have performed extensive atomistic Replica Exchange Molecular Dynamics simulations, elucidating preferential binding of Aβ₄₂ with the A2V and A2T hexapeptides compared to WT Aβ₆. A notable reorganization of the Aβ₄₂ landscape is revealed due to hexapeptide association, as manifested by lowering of transient interactions between the central and C-terminal hydrophobic patches. Concurrently, Aβ₆-bound Aβ₄₂ monomer exhibits alternative structural features that are strongly dependent on the hexapeptide sequence. For example, a central helix is more frequently populated within the A2T-bound monomer, while A2V-bound Aβ₄₂ is often enhanced in overall disorder. Taken together, the present simulations offer novel molecular insights onto the effect of the N-terminal hexapeptide binding on the Aβ₄₂ monomer structure, which might help in explaining their reported amyloid inhibition properties.

Study of structural stability and damaging effect on membrane for four Aβ42 dimers

Increasing evidence shows that Aβ oligomers are key pathogenic molecules in Alzheimer’s disease. Among Aβ oligomers, dimer is the smallest aggregate and toxic unit. Therefore, understanding its structural and dynamic properties is quite useful to prevent the formation and toxicity of the Aβ oligomers. In this study, we performed molecular dynamic simulations on four Aβ42 dimers, 2NCb, CNNC, NCNC and NCCN, within the hydrated DPPC membrane. Four Aβ42 dimers differ in the arrangements of two Aβ42 peptides. This study aims to investigate the impact of aggregation pattern of two Aβ peptides on the structural stability of the Aβ42 dimer and its disruption to the biological membrane. The MD results demonstrate that the NCCN, CNNC and NCNC have the larger structural fluctuation at the N-terminus of Aβ42 peptide, where the β-strand structure converts into the coil structure. The loss of the N-terminal β-strand further impairs the aggregate ability of Aβ42 dimer. In addition, inserting Aβ42 dimer into the membrane can considerably decrease the average APL of DPPC membrane. Moreover this decrease effect is largely dependent on the distance to the location of Aβ42 dimer and its secondary structure forms. Based on the results, the 2NCb is considered as a stable dimeric unit for aggregating the larger Aβ42 oligomer, and has a potent ability to disrupt the membrane.

Amyloidogenicity at a Distance: How Distal Protein Regions Modulate Aggregation in Disease

The misfolding of proteins to form amyloid is a key pathological feature of several progressive, and currently incurable, diseases. A mechanistic understanding of the pathway from soluble, native protein to insoluble amyloid is crucial for therapeutic design, and recent efforts have helped to elucidate the key molecular events that trigger protein misfolding. Generally, either global or local structural perturbations occur early in amyloidogenesis to expose aggregation-prone regions of the protein that can then self-associate to form toxic oligomers. Surprisingly, these initiating structural changes are often caused or influenced by protein regions distal to the classically amyloidogenic sequences. Understanding the importance of these distal regions in the pathogenic process has highlighted many remaining knowledge gaps regarding the precise molecular events that occur in classic aggregation pathways. In this review, we discuss how these distal regions can influence aggregation in disease and the recent technical and conceptual advances that have allowed this insight.

Alzheimer's Protective Cross-Interaction between Wild-Type and A2T Variants Alters Aβ42 Dimer Structure

Whole genome sequencing has recently revealed the protective effect of a single A2T mutation in heterozygous carriers against Alzheimer's disease (AD) and age-related cognitive decline. The impact of the protective cross-interaction between the wild-type (WT) and A2T variants on the dimer structure is therefore of high interest, as the Aβ dimers are the smallest known neurotoxic species. Toward this goal, extensive atomistic replica exchange molecular dynamics simulations of the solvated WT homo- and A2T hetero- Aβ_1-42 dimers have been performed, resulting into a total of 51 μs of sampling for each system. Weakening of a set of transient, intrachain contacts formed between the central and C-terminal hydrophobic residues is observed in the heterodimeric system. The majority of the heterodimers with reduced interaction between central and C-terminal regions lack any significant secondary structure and display a weak interchain interface. Interestingly, the A2T N-terminus, particularly residue F4, is frequently engaged in tertiary and quaternary interactions with central and C-terminal hydrophobic residues in those distinct structures, leading to hydrophobic burial. This atypical involvement of the N-terminus within A2T heterodimer revealed in our simulations implies possible interference on Aβ₄₂ aggregation and toxic oligomer formation, which is consistent with experiments. In conclusion, the present study provides detailed structural insights onto A2T Aβ₄₂ heterodimer, which might provide molecular insights onto the AD protective effect of the A2T mutation in the heterozygous state.

Fibrillation-prone conformations of the amyloid-β-42 peptide at the gold/water interface

Proteins in the proximity of inorganic surfaces and nanoparticles may undergo profound adjustments that trigger biomedically relevant processes, such as protein fibrillation. The mechanisms that govern protein-surface interactions at the molecular level are still poorly understood. In this work, we investigate the adsorption onto a gold surface, in water, of an amyloid-β (Aβ) peptide, which is the amyloidogenic peptide involved in Alzheimer's disease. The entire adsorption process, from the peptide in bulk water to its conformational relaxation on the surface, is explored by large-scale atomistic molecular dynamics (MD) simulations. We start by providing a description of the conformational ensemble of Aβ in solution by a 22 μs temperature replica exchange MD simulation, which is consistent with previous results. Then, we obtain a statistical description of how the peptide approaches the gold surface by multiple MD simulations, identifying the preferential gold-binding sites and giving a kinetic picture of the association process. Finally, relaxation of the Aβ conformations at the gold/water interface is performed by a 19 μs Hamiltonian-temperature replica exchange MD simulation. We find that the conformational ensemble of Aβ is strongly perturbed by the presence of the surface. In particular, at the gold/water interface the population of the conformers akin to amyloid fibrils is significantly enriched, suggesting that this extended contact geometry may promote fibrillation.

Computer Simulations of Intrinsically Disordered Proteins

The investigation of intrinsically disordered proteins (IDPs) is a new frontier in structural and molecular biology that requires a new paradigm to connect structural disorder to function. Molecular dynamics simulations and statistical thermodynamics potentially offer ideal tools for atomic-level characterizations and thermodynamic descriptions of this fascinating class of proteins that will complement experimental studies. However, IDPs display sensitivity to inaccuracies in the underlying molecular mechanics force fields. Thus, achieving an accurate structural characterization of IDPs via simulations is a challenge. It is also daunting to perform a configuration-space integration over heterogeneous structural ensembles sampled by IDPs to extract, in particular, protein configurational entropy. In this review, we summarize recent efforts devoted to the development of force fields and the critical evaluations of their performance when applied to IDPs. We also survey recent advances in computational methods for protein configurational entropy that aim to provide a thermodynamic link between structural disorder and protein activity.

Comparison of force fields for Alzheimer's A β42: A case study for intrinsically disordered proteins

Intrinsically disordered proteins are essential for biological processes such as cell signalling, but are also associated to devastating diseases including Alzheimer's disease, Parkinson's disease or type II diabetes. Because of their lack of a stable three-dimensional structure, molecular dynamics simulations are often used to obtain atomistic details that cannot be observed experimentally. The applicability of molecular dynamics simulations depends on the accuracy of the force field chosen to represent the underlying free energy surface of the system. Here, we use replica exchange molecular dynamics simulations to test five modern force fields, OPLS, AMBER99SB, AMBER99SB*ILDN, AMBER99SBILDN-NMR and CHARMM22*, in their ability to model Aβ42 , an intrinsically disordered peptide associated with Alzheimer's disease, and compare our results to nuclear magnetic resonance (NMR) experimental data. We observe that all force fields except AMBER99SBILDN-NMR successfully reproduce local NMR observables, with CHARMM22* being slightly better than the other force fields.

Is the Conformational Ensemble of Alzheimer's Aβ10-40 Peptide Force Field Dependent?

By applying REMD simulations we have performed comparative analysis of the conformational ensembles of amino-truncated Aβ10-40 peptide produced with five force fields, which combine four protein parameterizations (CHARMM36, CHARMM22*, CHARMM22/cmap, and OPLS-AA) and two water models (standard and modified TIP3P). Aβ10-40 conformations were analyzed by computing secondary structure, backbone fluctuations, tertiary interactions, and radius of gyration. We have also calculated Aβ10-40 ³J_HNHα-coupling and RDC constants and compared them with their experimental counterparts obtained for the full-length Aβ1-40 peptide. Our study led us to several conclusions. First, all force fields predict that Aβ adopts unfolded structure dominated by turn and random coil conformations. Second, specific TIP3P water model does not dramatically affect secondary or tertiary Aβ10-40 structure, albeit standard TIP3P model favors slightly more compact states. Third, although the secondary structures observed in CHARMM36 and CHARMM22/cmap simulations are qualitatively similar, their tertiary interactions show little consistency. Fourth, two force fields, OPLS-AA and CHARMM22* have unique features setting them apart from CHARMM36 or CHARMM22/cmap. OPLS-AA reveals moderate β-structure propensity coupled with extensive, but weak long-range tertiary interactions leading to Aβ collapsed conformations. CHARMM22* exhibits moderate helix propensity and generates multiple exceptionally stable long- and short-range interactions. Our investigation suggests that among all force fields CHARMM22* differs the most from CHARMM36. Fifth, the analysis of ³J_HNHα-coupling and RDC constants based on CHARMM36 force field with standard TIP3P model led us to an unexpected finding that in silico Aβ10-40 and experimental Aβ1-40 constants are generally in better agreement than these quantities computed and measured for identical peptides, such as Aβ1-40 or Aβ1-42. This observation suggests that the differences in the conformational ensembles of Aβ10-40 and Aβ1-40 are small and the former can be used as proxy of the full-length peptide. Based on this argument, we concluded that CHARMM36 force field with standard TIP3P model produces the most accurate representation of Aβ10-40 conformational ensemble.

Dependence of protein conformational ensembles on force field parameterizations limits the predictive power of molecular dynamics simulations. To address this problem, we evaluated five all-atom force fields for their consistency in reproducing the conformational ensemble of Alzheimer’s Aβ10-40 peptide. To generate conformational ensembles, we have used replica exchange molecular dynamics and computed Aβ10-40 secondary and tertiary structures. We found that, although all force fields predict Aβ10-40 unfolded structure, they strongly disagree on helix and β propensities and tertiary structure distributions. We have also calculated Aβ10-40 J-coupling and residual dipolar coupling constants and compared them with the experimental data for the full-length Aβ1-40 peptide. Unexpectedly, we determined that in silico Aβ10-40 and experimental Aβ1-40 constants are in better agreement than these quantities computed and measured previously for identical peptides, such as Aβ1-40 or Aβ1-42. We then concluded that the conformational ensembles of Aβ10-40 and Aβ1-40 are similar and on this basis argue that CHARMM36 force field with standard TIP3P water model provides the most accurate description of Aβ10-40. Although our objective was not to evaluate the biomolecular force fields in general, our study is expected to facilitate their proper selection for the simulations of Alzheimer’s peptides.

Molecular Dynamics Simulation of Tau Peptides for the Investigation of Conformational Changes Induced by Specific Phosphorylation Patterns

The Tau protein plays an important role due to its biomolecular interactions in neurodegenerative diseases. The lack of stable structure and various posttranslational modifications such as phosphorylation at various sites in the Tau protein pose a challenge for many experimental methods that are traditionally used to study protein folding and aggregation. Atomistic molecular dynamics (MD) simulations can help around deciphering relationship between phosphorylation and various intermediate and stable conformations of the Tau protein which occur on longer timescales. This chapter outlines protocols for the preparation, execution, and analysis of all-atom MD simulations of a 21-amino acid-long phosphorylated Tau peptide with the aim of generating biologically relevant structural and dynamic information. The simulations are done in explicit solvent and starting from nearly extended configurations of the peptide. The scaled MD method implemented in AMBER14 was chosen to achieve enhanced conformational sampling in addition to a conventional MD approach, thereby allowing the characterization of folding for such an intrinsically disordered peptide at 293 K. Emphasis is placed on the analysis of the simulation trajectories to establish correlations with NMR data (i.e., chemical shifts and NOEs). Finally, in-depth discussions are provided for commonly encountered problems.

Mechanism of Nucleation and Growth of Aβ40 Fibrils from All-Atom and Coarse-Grained Simulations

In this work, we characterize the nucleation and elongation mechanisms of the "diseased" polymorph of the amyloid-β 40 (Aβ40) fibril using an off-lattice coarse-grained (CG) protein model. After determining the nucleation size and subsequent stable protofibrillar structure from the CG model, validated with all-atom simulations, we consider the "lock and dock" and "activated monomer" fibril elongation mechanisms for the protofibril by statistical additions of a monomer drawn from four different ensembles of the free Aβ40 peptide to grow the fibril. Our CG model shows that the dominant mechanism for fibril elongation is the lock and dock mechanism across all monomer ensembles, even when the monomer is in the activated form. Although our CG model finds no thermodynamic difference between the two fibril elongation mechanisms, the activated monomer is found to be kinetically faster by a factor of 2 for the "locking" step compared with all other structured or unstructured monomer ensembles.

Comparing generalized ensemble methods for sampling of systems with many degrees of freedom

We compare two standard replica exchange methods using temperature and dielectric constant as the scaling variables for independent replicas against two new corresponding enhanced sampling methods based on non-equilibrium statistical cooling (temperature) or descreening (dielectric). We test the four methods on a rough 1D potential as well as for alanine dipeptide in water, for which their relatively small phase space allows for the ability to define quantitative convergence metrics. We show that both dielectric methods are inferior to the temperature enhanced sampling methods, and in turn show that temperature cool walking (TCW) systematically outperforms the standard temperature replica exchange (TREx) method. We extend our comparisons of the TCW and TREx methods to the 5 residue met-enkephalin peptide, in which we evaluate the Kullback-Leibler divergence metric to show that the rate of convergence between two independent trajectories is faster for TCW compared to TREx. Finally we apply the temperature methods to the 42 residue amyloid-β peptide in which we find non-negligible differences in the disordered ensemble using TCW compared to the standard TREx. All four methods have been made available as software through the OpenMM Omnia software consortium (http://www.omnia.md/).

Preparation of pure populations of covalently stabilized amyloid β-protein oligomers of specific sizes

Evidence suggests that amyloid β-protein (Aβ) oligomers may be seminal pathogenic agents in Alzheimer's disease (AD). If so, developing oligomer-targeted therapeutics requires an understanding of oligomer structure. This has been difficult due to the instability of these non-covalently associated Aβ assemblies. We previously used rapid, zero-length, in situ chemical cross-linking to stabilize oligomers of Aβ40. These enabled us to isolate pure, stable populations of dimers, trimers, and tetramers and to determine their structure-activity relationships. However, equivalent methods applied to Aβ42 did not produce stable oligomers. We report here that the use of an Aβ42 homologue, [F10, Y42]Aβ42, coupled with sequential denaturation/dissociation and gel electrophoresis procedures, provides the means to produce highly pure, stable populations of oligomers of sizes ranging from dimer through dodecamer that are suitable for structure-activity relationship determination.

Protein Aggregation and Molecular Crowding: Perspectives From Multiscale Simulations

Cells are extremely crowded environments, thus the use of diluted salted aqueous solutions containing a single protein is too simplistic to mimic the real situation. Macromolecular crowding might affect protein structure, folding, shape, conformational stability, binding of small molecules, enzymatic activity, interactions with cognate biomolecules, and pathological aggregation. The latter phenomenon typically leads to the formation of amyloid fibrils that are linked to several lethal neurodegenerative diseases, but that can also play a functional role in certain organisms. The majority of molecular simulations performed before the last few years were conducted in diluted solutions and were restricted both in the timescales and in the system dimensions by the available computational resources. In recent years, several computational solutions were developed to get close to physiological conditions. In this review we summarize the main computational techniques used to tackle the issue of protein aggregation both in a diluted and in a crowded environment.

Impact of Cu(II) Binding on Structures and Dynamics of Aβ42 Monomer and Dimer: Molecular Dynamics Study

The classical force field, which is compatible with the Amber force field 99SB, has been obtained for the interaction of Cu(II) with monomer and dimers of amyloid-β peptides using the coordination where Cu(II) is bound to His6, His13 (or His14), and Asp1 with distorted planar geometry. The newly developed force field and molecular dynamics simulation were employed to study the impact of Cu(II) binding on structures and dynamics of Aβ₄₂ monomer and dimers. It was shown that in the presence of Cu(II) the β content of monomer is reduced substantially compared with the wild-type Aβ₄₂ suggesting that, in accord with experiments, metal ions facilitate formation of amorphous aggregates rather than amyloid fibrils with cross-β structures. In addition, one possible mechanism for amorphous assembly is that the Asp23-Lys28 salt bridge, which plays a crucial role in β sheet formation, becomes more flexible upon copper ion binding to the Aβ N-terminus. The simulation of dimers was conducted with the Cu(II)/Aβ stoichiometric ratios of 1:1 and 1:2. For the 1:1 ratio Cu(II) delays the Aβ dimerization process as observed in a number of experiments. The mechanism underlying this phenomenon is associated with slow formation of interchain salt bridges in dimer as well as with decreased hydrophobicity of monomer upon Cu-binding.

Assembly of Peptides Derived from β-Sheet Regions of β-Amyloid

In Alzheimer’s disease, aggregation of the β-amyloid peptide (Aβ) results in the formation of oligomers and fibrils that are associated with neurodegeneration. Aggregation of Aβ occurs through interactions between different regions of the peptide. This paper and the accompanying paper constitute a two-part investigation of two key regions of Aβ: the central region and the C-terminal region. These two regions promote aggregation and adopt β-sheet structure in the fibrils, and may also do so in the oligomers. In this paper, we study the assembly of macrocyclic β-sheet peptides that contain residues 17–23 (LVFFAED) from the central region and residues 30–36 (AIIGLMV) from the C-terminal region. These peptides assemble to form tetramers. Each tetramer consists of two hydrogen-bonded dimers that pack through hydrophobic interactions in a sandwich-like fashion. Incorporation of a single ¹⁵N isotopic label into each peptide provides a spectroscopic probe with which to elucidate the β-sheet assembly and interaction: ¹H,¹⁵N HSQC studies facilitate the identification of the monomers and tetramers; ¹⁵N-edited NOESY studies corroborate the pairing of the dimers within the tetramers. In the following paper, J. Am. Chem. Soc.2016, DOI: 10.1021/jacs.6b06001, we will extend these studies to elucidate the coassembly of the peptides to form heterotetramers.

Exploring the aggregation free energy landscape of the amyloid-β protein (1-40)

A predictive coarse-grained protein force field [associative memory, water-mediated, structure, and energy model for molecular dynamics (AWSEM)-MD] is used to study the energy landscapes and relative stabilities of amyloid-β protein (1-40) in the monomer and all of its oligomeric forms up to an octamer. We find that an isolated monomer is mainly disordered with a short α-helix formed at the central hydrophobic core region (L17-D23). A less stable hairpin structure, however, becomes increasingly more stable in oligomers, where hydrogen bonds can form between neighboring monomers. We explore the structure and stability of both prefibrillar oligomers that consist of mainly antiparallel β-sheets and fibrillar oligomers with only parallel β-sheets. Prefibrillar oligomers are polymorphic but typically take on a cylindrin-like shape composed of mostly antiparallel β-strands. At the concentration of the simulation, the aggregation free energy landscape is nearly downhill. We use umbrella sampling along a structural progress coordinate for interconversion between prefibrillar and fibrillar forms to identify a conversion pathway between these forms. The fibrillar oligomer only becomes favored over its prefibrillar counterpart in the pentamer where an interconversion bottleneck appears. The structural characterization of the pathway along with statistical mechanical perturbation theory allow us to evaluate the effects of concentration on the free energy landscape of aggregation as well as the effects of the Dutch and Arctic mutations associated with early onset of Alzheimer's disease.

Monomer Dynamics of Alzheimer Peptides and Kinetic Control of Early Aggregation in Alzheimer's Disease

The rate of reconfiguration-or intramolecular diffusion-of monomeric Alzheimer (Aβ) peptides is measured and, under conditions that aggregation is more likely, peptide diffusion slows down significantly, which allows bimolecular associations to be initiated. By using the method of Trp-Cys contact quenching, the rate of reconfiguration is observed to be about five times faster for Aβ₄₀ , which aggregates slowly, than that for Aβ₄₂ , which aggregates quickly. Furthermore, the rate of reconfiguration for Aβ₄₂ speeds up at higher pH, which slows aggregation, and in the presence of the aggregation inhibitor curcumin. The measured reconfiguration rates are able to predict the early aggregation behavior of the Aβ peptide and provide a kinetic basis for why Aβ₄₂ is more prone to aggregation than Aβ₄₀ , despite a difference of only two amino acids.

Finding Our Way in the Dark Proteome

The traditional structure-function paradigm has provided significant insights for well-folded proteins in which structures can be easily and rapidly revealed by X-ray crystallography beamlines. However, approximately one-third of the human proteome is comprised of intrinsically disordered proteins and regions (IDPs/IDRs) that do not adopt a dominant well-folded structure, and therefore remain "unseen" by traditional structural biology methods. This Perspective considers the challenges raised by the "Dark Proteome", in which determining the diverse conformational substates of IDPs in their free states, in encounter complexes of bound states, and in complexes retaining significant disorder requires an unprecedented level of integration of multiple and complementary solution-based experiments that are analyzed with state-of-the art molecular simulation, Bayesian probabilistic models, and high-throughput computation. We envision how these diverse experimental and computational tools can work together through formation of a "computational beamline" that will allow key functional features to be identified in IDP structural ensembles.

Molecular architecture of Aβ fibrils grown in cerebrospinal fluid solution and in a cell culture model of Aβ plaque formation

OBJECTIVES

The detailed structure of brain-derived Aβ amyloid fibrils is unknown. To approach this issue, we investigate the molecular architecture of Aβ(1-40) fibrils grown in either human cerebrospinal fluid solution, in chemically simple phosphate buffer in vitro or extracted from a cell culture model of Aβ amyloid plaque formation.

METHODS

We have used hydrogen-deuterium exchange (HX) combined with nuclear magnetic resonance, transmission electron microscopy, seeding experiments both in vitro and in cell culture as well as several other spectroscopic measurements to compare the morphology and residue-specific conformation of these different Aβ fibrils.

RESULTS AND CONCLUSIONS

Our data reveal that, despite considerable variations in morphology, the spectroscopic properties and the pattern of slowly exchanging backbone amides are closely similar in the fibrils investigated. This finding implies that a fundamentally conserved molecular architecture of Aβ peptide fold is common to Aβ fibrils.

Slide-and-exchange mechanism for rapid and selective transport through the nuclear pore complex

Nucleocytoplasmic transport is mediated by the interaction of transport factors (TFs) with disordered phenylalanine-glycine (FG) repeats that fill the central channel of the nuclear pore complex (NPC). However, the mechanism by which TFs rapidly diffuse through multiple FG repeats without compromising NPC selectivity is not yet fully understood. In this study, we build on our recent NMR investigations showing that FG repeats are highly dynamic, flexible, and rapidly exchanging among TF interaction sites. We use unbiased long timescale all-atom simulations on the Anton supercomputer, combined with extensive enhanced sampling simulations and NMR experiments, to characterize the thermodynamic and kinetic properties of FG repeats and their interaction with a model transport factor. Both the simulations and experimental data indicate that FG repeats are highly dynamic random coils, lack intrachain interactions, and exhibit significant entropically driven resistance to spatial confinement. We show that the FG motifs reversibly slide in and out of multiple TF interaction sites, transitioning rapidly between a strongly interacting state and a weakly interacting state, rather than undergoing a much slower transition between strongly interacting and completely noninteracting (unbound) states. In the weakly interacting state, FG motifs can be more easily displaced by other competing FG motifs, providing a simple mechanism for rapid exchange of TF/FG motif contacts during transport. This slide-and-exchange mechanism highlights the direct role of the disorder within FG repeats in nucleocytoplasmic transport, and resolves the apparent conflict between the selectivity and speed of transport.

Experimental Inferential Structure Determination of Ensembles for Intrinsically Disordered Proteins

We develop a Bayesian approach to determine the most probable structural ensemble model from candidate structures for intrinsically disordered proteins (IDPs) that takes full advantage of NMR chemical shifts and J-coupling data, their known errors and variances, and the quality of the theoretical back-calculation from structure to experimental observables. Our approach differs from previous formulations in the optimization of experimental and back-calculation nuisance parameters that are treated as random variables with known distributions, as opposed to structural or ensemble weight optimization or use of a reference ensemble. The resulting experimental inferential structure determination (EISD) method is size extensive with O(N) scaling, with N = number of structures, that allows for the rapid ranking of large ensemble data comprising tens of thousands of conformations. We apply the EISD approach on singular folded proteins and a corresponding set of ∼25 000 misfolded states to illustrate the problems that can arise using Boltzmann weighted priors. We then apply the EISD method to rank IDP ensembles most consistent with the NMR data and show that the primary error for ranking or creating good IDP ensembles resides in the poor back-calculation from structure to simulated experimental observable. We show that a reduction by a factor of 3 in the uncertainty of the back-calculation error can improve the discrimination among qualitatively different IDP ensembles for the amyloid-beta peptide.

Advances in the Simulation of Protein Aggregation at the Atomistic Scale

Protein aggregation into highly structured amyloid fibrils is associated with various diseases including Alzheimer's disease, Parkinson's disease, and type II diabetes. Amyloids can also have normal biological functions and, in the future, could be used as the basis for novel nanoscale materials. However, a full understanding of the physicochemical forces that drive protein aggregation is still lacking. Such understanding is crucial for the development of drugs that can effectively inhibit aberrant amyloid aggregation and for the directed design of functional amyloids. Atomistic simulations can help understand protein aggregation. In particular, atomistic simulations can be used to study the initial formation of toxic oligomers which are hard to characterize experimentally and to understand the difference in aggregation behavior between different amyloidogenic peptides. Here, we review the latest atomistic simulations of protein aggregation, concentrating on amyloidogenic protein fragments, and provide an outlook for the future in this field.

Monomeric Aβ(1-40) and Aβ(1-42) Peptides in Solution Adopt Very Similar Ramachandran Map Distributions That Closely Resemble Random Coil

The pathogenesis of Alzheimer’s disease is characterized by the aggregation and fibrillation of amyloid peptides Aβ^1–40 and Aβ^1–42 into amyloid plaques. Despite strong potential therapeutic interest, the structural pathways associated with the conversion of monomeric Aβ peptides into oligomeric species remain largely unknown. In particular, the higher aggregation propensity and associated toxicity of Aβ^1–42 compared to that of Aβ^1–40 are poorly understood. To explore in detail the structural propensity of the monomeric Aβ^1–40 and Aβ^1–42 peptides in solution, we recorded a large set of nuclear magnetic resonance (NMR) parameters, including chemical shifts, nuclear Overhauser effects (NOEs), and J couplings. Systematic comparisons show that at neutral pH the Aβ^1–40 and Aβ^1–42 peptides populate almost indistinguishable coil-like conformations. Nuclear Overhauser effect spectra collected at very high resolution remove assignment ambiguities and show no long-range NOE contacts. Six sets of backbone J couplings (³J_HNHα, ³J_C′C′, ³J_C′Hα, ¹J_HαCα, ²J_NCα, and ¹J_NCα) recorded for Aβ^1–40 were used as input for the recently developed MERA Ramachandran map analysis, yielding residue-specific backbone ϕ/ψ torsion angle distributions that closely resemble random coil distributions, the absence of a significantly elevated propensity for β-conformations in the C-terminal region of the peptide, and a small but distinct propensity for α_L at K28. Our results suggest that the self-association of Aβ peptides into toxic oligomers is not driven by elevated propensities of the monomeric species to adopt β-strand-like conformations. Instead, the accelerated disappearance of Aβ NMR signals in D₂O over H₂O, particularly pronounced for Aβ^1–42, suggests that intermolecular interactions between the hydrophobic regions of the peptide dominate the aggregation process.

Accurate measurement of (3)J(HNHα) couplings in small or disordered proteins from WATERGATE-optimized TROSY spectra

Provided that care is taken in adjusting the WATERGATE element of a (1)H-(15)N TROSY-HSQC experiment, such that neither the water magnetization nor the (1)H(α) protons are inverted by its final 180° pulse, (3)JHNHα couplings can be measured directly from splittings in the (1)H dimension of the spectrum. With band-selective (1)H decoupling, very high (15)N resolution can be achieved. A complete set of (3)JHNHα values, ranging from 3.4 to 10.1 Hz was measured for the 56-residue third domain of IgG-binding protein G (GB3). Using the H-N-C(α)-H(α) dihedral angles extracted from a RDC-refined structure of GB3, (3)JHNHα values predicted by a previously parameterized Karplus equation agree to within a root-mean-square deviation (rmsd) of 0.37 Hz with the experimental data. Values measured for the Alzheimer's implicated Aβ(1-40) peptide fit to within an rmsd of 0.45 Hz to random coil (3)JHNHα values.

Transient β-hairpin formation in α-synuclein monomer revealed by coarse-grained molecular dynamics simulation

Parkinson's disease, originating from the intrinsically disordered peptide α-synuclein, is a common neurodegenerative disorder that affects more than 5% of the population above age 85. It remains unclear how α-synuclein monomers undergo conformational changes leading to aggregation and formation of fibrils characteristic for the disease. In the present study, we perform molecular dynamics simulations (over 180 μs in aggregated time) using a hybrid-resolution model, Proteins with Atomic details in Coarse-grained Environment (PACE), to characterize in atomic detail structural ensembles of wild type and mutant monomeric α-synuclein in aqueous solution. The simulations reproduce structural properties of α-synuclein characterized in experiments, such as secondary structure content, long-range contacts, chemical shifts, and (3)J(HNHCα )-coupling constants. Most notably, the simulations reveal that a short fragment encompassing region 38-53, adjacent to the non-amyloid-β component region, exhibits a high probability of forming a β-hairpin; this fragment, when isolated from the remainder of α-synuclein, fluctuates frequently into its β-hairpin conformation. Two disease-prone mutations, namely, A30P and A53T, significantly accelerate the formation of a β-hairpin in the stated fragment. We conclude that the formation of a β-hairpin in region 38-53 is a key event during α-synuclein aggregation. We predict further that the G47V mutation impedes the formation of a turn in the β-hairpin and slows down β-hairpin formation, thereby retarding α-synuclein aggregation.

Characterization of Aβ Monomers through the Convergence of Ensemble Properties among Simulations with Multiple Force Fields

Amyloid β (Aβ) monomers represent a base state in the pathways of aggregation that result in the fibrils and oligomers implicated in the pathogenesis of Alzheimer's disease (AD). The structural properties of these intrinsically disordered peptides remain unclear despite extensive experimental and computational investigations. Further, there are mutations within Aβ that change the way the peptide aggregates and are known to cause familial AD (FAD). Here, we analyze the ensembles of different isoforms (Aβ42 and Aβ40) and mutants (E22Δ, D23N, E22K, E22G, and A2T in Aβ40) of Aβ generated with all-atom replica exchange molecular dynamics (REMD) simulations on the μs/replica time scale. These were run using three different force field/water model combinations: OPLS-AA/L and TIP3P ("OPLS"), AMBER99sb-ILDN and TIP4P-Ew ("ILDN"), as well as CHARMM22* and TIP3SP ("CHARMM"). Despite fundamental changes in simulation parameters, we find that the resulting ensembles demonstrate a strong convergence in structural properties. In particular, antiparallel contacts between L17-A21 and A30-L34 are prevalent in ensembles of Aβ40, directly forming β sheets in the OPLS and ILDN combinations. A21-A30 commonly forms an interceding region that rarely interacts with the rest of the peptide. Further, Aβ42 contributes new β hairpin motifs involving V40-I41 in both OPLS and ILDN. However, the structural flexibility of the central region and the electrostatic interactions that characterize it are notably different between the different conditions. Further, for OPLS, each of the FAD mutations disrupts central bend character and increases the polymorphism of antiparallel contacts across the central region. However, the studied mutations in the ILDN set primarily encourage more global contacts involving the N-terminus and the central region, and promote the formation of new β topologies that may seed different aggregates involved in disease phenotypes. These differences aside, the large degree of agreement between simulation sets across multiple force fields provides a generalizable characterization of Aβ that is also consistent with experimental data and models.

Investigating protein conformational energy landscapes and atomic resolution dynamics from NMR dipolar couplings: a review

Nuclear magnetic resonance spectroscopy is exquisitely sensitive to protein dynamics. In particular inter-nuclear dipolar couplings, that become measurable in solution when the protein is dissolved in a dilute liquid crystalline solution, report on all conformations sampled up to millisecond timescales. As such they provide the opportunity to describe the Boltzmann distribution present in solution at atomic resolution, and thereby to map the conformational energy landscape in unprecedented detail. The development of analytical methods and approaches based on numerical simulation and their application to numerous biologically important systems is presented.

A facile method for expression and purification of (15)N isotope-labeled human Alzheimer's β-amyloid peptides from E. coli for NMR-based structural analysis

Alzheimer's disease (AD) is a progressive neurodegenerative disease affecting millions of people worldwide. AD is characterized by the presence of extracellular plaques composed of aggregated/oligomerized β-amyloid peptides with Aβ42 peptide representing a major isoform in the senile plaques. Given the pathological significance of Aβ42 in the progression of AD, there is considerable interest in understanding the structural ensembles for soluble monomer and oligomeric forms of Aβ42. This report describes an efficient method to express and purify high quality (15)N isotope-labeled Aβ42 for structural studies by NMR. The protocol involves utilization of an auto induction system with (15)N isotope labeled medium, for high-level expression of Aβ42 as a fusion with IFABP. After the over-expression of the (15)N isotope-labeled IFABP-Aβ42 fusion protein in the inclusion bodies, pure (15)N isotope-labeled Aβ42 peptide is obtained following a purification method that is streamlined and improved from the method originally developed for the isolation of unlabeled Aβ42 peptide (Garai et al., 2009). We obtain a final yield of ∼ 6 mg/L culture for (15)N isotope-labeled Aβ42 peptide. Mass spectrometry and (1)H-(15)N HSQC spectra of monomeric Aβ42 peptide validate the uniform incorporation of the isotopic label. The method described here is equally applicable for the uniform isotope labeling with (15)N and (13)C in Aβ42 peptide as well as its other variants including any Aβ42 peptide mutants.

Structural Ensembles of Intrinsically Disordered Proteins Depend Strongly on Force Field: A Comparison to Experiment

Intrinsically disordered proteins (IDPs) are notoriously challenging to study both experimentally and computationally. The structure of IDPs cannot be described by a single conformation but must instead be described as an ensemble of interconverting conformations. Atomistic simulations are increasingly used to obtain such IDP conformational ensembles. Here, we have compared the IDP ensembles generated by eight all-atom empirical force fields against primary small-angle X-ray scattering (SAXS) and NMR data. Ensembles obtained with different force fields exhibit marked differences in chain dimensions, hydrogen bonding, and secondary structure content. These differences are unexpectedly large: changing the force field is found to have a stronger effect on secondary structure content than changing the entire peptide sequence. The CHARMM 22* ensemble performs best in this force field comparison: it has the lowest error in chemical shifts and J-couplings and agrees well with the SAXS data. A high population of left-handed α-helix is present in the CHARMM 36 ensemble, which is inconsistent with measured scalar couplings. To eliminate inadequate sampling as a reason for differences between force fields, extensive simulations were carried out (0.964 ms in total); the remaining small sampling uncertainty is shown to be much smaller than the observed differences. Our findings highlight how IDPs, with their rugged energy landscapes, are highly sensitive test systems that are capable of revealing force field deficiencies and, therefore, contributing to force field development.

Structures of the Alzheimer's Wild-Type Aβ1-40 Dimer from Atomistic Simulations

We have studied the dimer of amyloid beta peptide Aβ of 40 residues by means of all-atom replica exchange molecular dynamics. The Aβ-dimers have been found to be the smallest toxic species in Alzheimer's disease, but their inherent flexibilities have precluded structural characterization by experimental methods. Though the 24-μs-scale simulation reveals a mean secondary structure of 18% β-strand and 10% α helix, we find transient configurations with an unstructured N-terminus and multiple β-hairpins spanning residues 17-21 and 30-36, but the antiparallel and perpendicular peptide orientations are preferred over the parallel organization. Short-lived conformational states also consist of all α topologies, and one compact peptide with β-sheet structure stabilized by a rather extended peptide with α-helical content. Overall, this first all-atom study provides insights into the equilibrium structure of the Aβ1-40 dimer in aqueous solution, opening a new avenue for a comprehensive understanding of the impact of pathogenic and protective mutations in early-stage Alzheimer's disease on a molecular level.

High-resolution NMR characterization of low abundance oligomers of amyloid-β without purification

Alzheimer's disease is characterized by the misfolding and self-assembly of the amyloidogenic protein amyloid-β (Aβ). The aggregation of Aβ leads to diverse oligomeric states, each of which may be potential targets for intervention. Obtaining insight into Aβ oligomers at the atomic level has been a major challenge to most techniques. Here, we use magic angle spinning recoupling (1)H-(1)H NMR experiments to overcome many of these limitations. Using (1)H-(1)H dipolar couplings as a NMR spectral filter to remove both high and low molecular weight species, we provide atomic-level characterization of a non-fibrillar aggregation product of the Aβ1-40 peptide using non-frozen samples without isotopic labeling. Importantly, this spectral filter allows the detection of the specific oligomer signal without a separate purification procedure. In comparison to other solid-state NMR techniques, the experiment is extraordinarily selective and sensitive. A resolved 2D spectra could be acquired of a small population of oligomers (6 micrograms, 7% of the total) amongst a much larger population of monomers and fibers (93% of the total). By coupling real-time (1)H-(1)H NMR experiments with other biophysical measurements, we show that a stable, primarily disordered Aβ1-40 oligomer 5-15 nm in diameter can form and coexist in parallel with the well-known cross-β-sheet fibrils.

Exploring the Alzheimer amyloid-β peptide conformational ensemble: A review of molecular dynamics approaches

Alzheimer's disease is one of the most common dementia among elderly worldwide. There is no therapeutic drugs until now to treat effectively this disease. One main reason is due to the poorly understood mechanism of Aβ peptide aggregation, which plays a crucial role in the development of Alzheimer's disease. It remains challenging to experimentally or theoretically characterize the secondary and tertiary structures of the Aβ monomer because of its high flexibility and aggregation propensity, and its conformations that lead to the aggregation are not fully identified. In this review, we highlight various structural ensembles of Aβ peptide revealed and characterized by computational approaches in order to find converging structures of Aβ monomer. Understanding how Aβ peptide forms transiently stable structures prior to aggregation will contribute to the design of new therapeutic molecules against the Alzheimer's disease.

Coupling of Zinc-Binding and Secondary Structure in Nonfibrillar Aβ40 Peptide Oligomerization

Nonfibrillar neurotoxic amyloid β (Aβ) oligomer structures are typically rich in β-sheets, which could be promoted by metal ions like Zn(2+). Here, using molecular dynamics (MD) simulations, we systematically examined combinations of Aβ40 peptide conformations and Zn(2+) binding modes to probe the effects of secondary structure on Aβ dimerization energies and kinetics. We found that random conformations do not contribute to dimerization either thermodynamically or kinetically. Zn(2+) couples with preformed secondary structures (α-helix and β-hairpin) to speed dimerization and stabilize the resulting dimer. Partial α-helices increase the dimerization speed, and dimers with α-helix rich conformations have the lowest energy. When Zn(2+) coordinates with residues D1, H6, H13, and H14, Aβ40 β-hairpin monomers have the fastest dimerization speed. Dimers with experimentally observed zinc coordination (E11, H6, H13, and H14) form with slower rate but have lower energy. Zn(2+) cannot stabilize fibril-like β-arch dimers. However, Zn(2+)-bound β-arch tetramers have the lowest energy. Collectively, zinc-stabilized β-hairpin oligomers could be important in the nucleation-polymerization of cross-β structures. Our results are consistent with experimental findings that α-helix to β-structural transition should accompany Aβ aggregation in the presence of zinc ions and that Zn(2+) stabilizes nonfibrillar Aβ oligomers and, thus, inhibits formation of less toxic Aβ fibrils.

Varied Probability of Staying Collapsed/Extended at the Conformational Equilibrium of Monomeric Aβ40 and Aβ42

In present study, we set out to investigate the conformation dynamics of Aβ40 and Aβ42 through exploring the impact of intra-molecular interactions on conformation dynamics using equilibrium molecular dynamics simulations. Our 40 microsecond-scale simulations reveal heterogeneous conformation ensembles of Aβ40 and Aβ42 that encompass ~35% β-strand and ~60% unstructured coils. Two conformational states were identified in both alloforms: a collapsed state (CS) that resembles the structural motif of face-to-face hydrophobic clustering in amyloid fibrils, and an extended state (ES) that features the structural characteristics of anti-parallel β-sheets in amyloid oligomers. In Aβ40, the C-terminus remains unstructured and rarely interacts with other parts, thereof the hydrophobic clustering is in loose contact and the peptide assumes ES with high probability. In contrast, the C-terminus of Aβ42 adopts a β-strand structure that strongly interacts with segments E3-R5 and V18-A21. The active association leads to a more compact hydrophobic collapse and refrain the alloform from ES. Based on the structural characterization, we propose that the fibril and oligomer assembly pathways could respectively take off from CS and ES, and their aggregation propensity may be governed by the probability of visiting the corresponding conformational states at the equilibrium.

Effects of the Pathogenic Mutation A117V and the Protective Mutation H111S on the Folding and Aggregation of PrP106-126: Insights from Replica Exchange Molecular Dynamics Simulations

The fragment 106-126 of prion protein exhibits similar properties to full-length prion. Experiments have shown that the A117V mutation enhances the aggregation of PrP106-126, while the H111S mutation abolishes the assembly. However, the mechanism of the change in the aggregation behavior of PrP106-126 upon the two mutations is not fully understood. In this study, replica exchange molecular dynamics simulations were performed to investigate the conformational ensemble of the WT PrP106-126 and its two mutants A117V and H111S. The obtained results indicate that the three species are all intrinsically disordered but they have distinct morphological differences. The A117V mutant has a higher propensity to form β-hairpin structures than the WT, while the H111S mutant has a higher population of helical structures. Furthermore, the A117V mutation increases the hydrophobic solvent accessible surface areas of PrP106-126 and the H111S mutation reduces the exposure of hydrophobic residues. It can be concluded that the difference in populations of β-hairpin structures and the change of hydrophobic solvent accessible areas may induce the different aggregation behaviors of the A117V and the H111S mutated PrP106-126. Understanding why the two mutations have contrary effects on the aggregation of PrP106-126 is very meaningful for further elucidation of the mechanism underlying aggregation and design of inhibitor against aggregation process.

Amyloid misfolding, aggregation, and the early onset of protein deposition diseases: insights from AFM experiments and computational analyses

The development of Alzheimer's disease is believed to be caused by the assembly of amyloid β proteins into aggregates and the formation of extracellular senile plaques. Similar models suggest that structural misfolding and aggregation of proteins are associated with the early onset of diseases such as Parkinson's, Huntington's, and other protein deposition diseases. Initially, the aggregates were structurally characterized by traditional techniques such as x-ray crystallography, NMR, electron microscopy, and AFM. However, data regarding the structures formed during the early stages of the aggregation process were unknown. Experimental models of protein deposition diseases have demonstrated that the small oligomeric species have significant neurotoxicity. This highlights the urgent need to discover the properties of these species, to enable the development of efficient diagnostic and therapeutic strategies. The oligomers exist transiently, making it impossible to use traditional structural techniques to study their characteristics. The recent implementation of single-molecule imaging and probing techniques that are capable of probing transient states have enabled the properties of these oligomers to be characterized. Additionally, powerful computational techniques capable of structurally analyzing oligomers at the atomic level advanced our understanding of the amyloid aggregation problem. This review outlines the progress in AFM experimental studies and computational analyses with a primary focus on understanding the very first stage of the aggregation process. Experimental approaches can aid in the development of novel sensitive diagnostic and preventive strategies for protein deposition diseases, and several examples of these approaches will be discussed.