Mechanical Resistance in Unstructured Proteins

Single-molecule pulling experiments on unstructured proteins linked to neurodegenerative diseases have measured rupture forces comparable to those for stable folded proteins. To investigate the structural mechanisms of this unexpected force resistance, we perform pulling simulations of the amyloid β-peptide (Aβ) and α-synuclein (αS), starting from simulated conformational ensembles for the free monomers. For both proteins, the simulations yield a set of rupture events that agree well with the experimental data. By analyzing the conformations occurring shortly before rupture in each event, we find that the mechanically resistant structures share a common architecture, with similarities to the folds adopted by Aβ and αS in amyloid fibrils. The disease-linked Arctic mutation of Aβ is found to increase the occurrence of highly force-resistant structures. Our study suggests that the high rupture forces observed in Aβ and αS pulling experiments are caused by structures that might have a key role in amyloid formation.     

The Structure of Misfolded Amyloidogenic Dimers: Computational Analysis of Force Spectroscopy Data

Progress in understanding the molecular mechanism of self-assembly of amyloidogenic proteins and peptides requires knowledge about their structure in misfolded states. Structural studies of amyloid aggregates formed during the early aggregation stage are very limited. Atomic force microscopy (AFM) spectroscopy is widely used to analyze misfolded proteins and peptides, but the structural characterization of transiently formed misfolded dimers is limited by the lack of computational approaches that allow direct comparison with AFM experiments. Steered molecular dynamics (SMD) simulation is capable of modeling force spectroscopy experiments, but the modeling requires pulling rates 10⁷ times higher than those used in AFM experiments. In this study, we describe a computational all-atom Monte Carlo pulling (MCP) approach that enables us to model results at pulling rates comparable to those used in AFM pulling experiments. We tested the approach by modeling pulling experimental data for I91 from titin I-band (PDB ID: 1TIT) and ubiquitin (PDB ID: 1UBQ). We then used MCP to analyze AFM spectroscopy experiments that probed the interaction of the peptides [Q6C] Sup35 (6–13) and [H13C] Aβ (13–23). A comparison of experimental results with the computational data for the Sup35 dimer with out-of-register and in-register arrangements of β-sheets suggests that Sup35 monomers adopt an out-of-register arrangement in the dimer. A similar analysis performed for Aβ peptide demonstrates that the out-of-register antiparallel β-sheet arrangement of monomers also occurs in this peptide. Although the rupture of hydrogen bonds is the major contributor to dimer dissociation, the aromatic-aromatic interaction also contributes to the dimer rupture process.     

Effects of Hydrophobic Macromolecular Crowders on Amyloid β (16–22) Aggregation

In Alzheimer’s disease (AD), the amyloid β (Aβ) peptide aggregates in the brain to form progressively larger oligomers, fibrils, and plaques. The aggregation process is strongly influenced by the presence of other macromolecular species, called crowders, that can exert forces on the proteins. One very common attribute of macromolecular crowders is their hydrophobicity. We examined the effect of hydrophobic crowders on protein aggregation by using discontinuous molecular dynamics (DMD) simulations in combination with an intermediate resolution protein model, PRIME20. The systems considered contained 48 Aβ (16–22) peptides and crowders with diameters of 5 Å, 20 Å, and 40 Å, represented by hard spheres or spheres with square-well/square-shoulder interactions, at a crowder volume fraction of ϕ = 0.10. Results show that low levels of crowder hydrophobicity are capable of increasing the fibrillation lag time and high levels of crowder hydrophobicity can fully prevent the formation of fibrils. The types of structures that remain during the final stages of the simulations are summarized in a global phase diagram that shows fibril, disordered oligomer, or β-sheet phases in the space spanned by crowder size and crowder hydrophobicity. In particular, at high levels of hydrophobicity, simulations with 5 Å crowders result in only disordered oligomers and simulations with 40 Å crowders result in only β-sheets. The presence of hydrophobic crowders reduces the antiparallel β-sheet content of fibrils, whereas hard sphere crowders increase it. Finally, strong hydrophobic crowders alter the secondary structure of the Aβ (16–22) monomers, bending them into a shape that is incapable of forming ordered β-sheets or fibrils. These results qualitatively agree with previous theoretical and experimental work.     

In Silico Analysis of the Apolipoprotein E and the Amyloid β Peptide Interaction: Misfolding Induced by Frustration of the Salt Bridge Network

The relationship between Apolipoprotein E (ApoE) and the aggregation processes of the amyloid β (Aβ) peptide has been shown to be crucial for Alzheimer''s disease (AD). The presence of the ApoE4 isoform is considered to be a contributing risk factor for AD. However, the detailed molecular properties of ApoE4 interacting with the Aβ peptide are unknown, although various mechanisms have been proposed to explain the physiological and pathological role of this relationship. Here, computer simulations have been used to investigate the process of Aβ interaction with the N-terminal domain of the human ApoE isoforms (ApoE2, ApoE3 and ApoE4). Molecular docking combined with molecular dynamics simulations have been undertaken to determine the Aβ peptide binding sites and the relative stability of binding to each of the ApoE isoforms. Our results show that from the several ApoE isoforms investigated, only ApoE4 presents a misfolded intermediate when bound to Aβ. Moreover, the initial α-helix used as the Aβ peptide model structure also becomes unstructured due to the interaction with ApoE4. These structural changes appear to be related to a rearrangement of the salt bridge network in ApoE4, for which we propose a model. It seems plausible that ApoE4 in its partially unfolded state is incapable of performing the clearance of Aβ, thereby promoting amyloid forming processes. Hence, the proposed model can be used to identify potential drug binding sites in the ApoE4-Aβ complex, where the interaction between the two molecules can be inhibited.     

Increasing the Amphiphilicity of an Amyloidogenic Peptide Changes the β-Sheet Structure in the Fibrils from Antiparallel to Parallel

Solid-state NMR measurements have been reported for four peptides derived from β-amyloid peptide Aβ(1–42): Aβ(1–40), Aβ(10–35), Aβ(16–22), and Aβ(34–42). Of these, the first two are predicted to be amphiphilic and were reported to form parallel β-sheets, whereas the latter two peptides appear nonamphiphilic and adopt an antiparallel β-sheet organization. These results suggest that amphiphilicity may be significant in determining fibril structure. Here, we demonstrate that acylation of Aβ(16–22) with octanoic acid increases its amphiphilicity and changes the organization of fibrillar β-sheet from antiparallel to parallel. Electron microscopy, Congo Red binding, and one-dimensional ¹³C NMR measurements demonstrate that octanoyl-Aβ(16–22) forms typical amyloid fibrils. Based on the stability of monolayers at the air-water interface, octanoyl-Aβ(16–22) is more amphiphilic than Aβ(16–22). Measurements of ¹³C-¹³C and ¹⁵N-¹³C nuclear magnetic dipole-dipole couplings in isotopically labeled fibril samples, using the constant-time finite-pulse radiofrequency-driven recoupling (fpRFDR-CT) and rotational echo double resonance (REDOR) solid-state NMR techniques, demonstrate that octanoyl-Aβ(16–22) fibrils are composed of parallel β-sheets, whereas Aβ(16–22) fibrils are composed of antiparallel β-sheets. These data demonstrate that amphiphilicity is critical in determining the structural organization of β-sheets in the amyloid fibril. This work also shows that all amyloid fibrils do not share a common supramolecular structure, and suggests a method for controlling the structure of amyloid fibrils.     

Differences in β-strand Populations of Monomeric Aβ40 and Aβ42

Using homonuclear ¹H NOESY spectra, with chemical shifts, ³J_H^N_H^α scalar couplings, residual dipolar couplings, and ¹H-¹⁵N NOEs, we have optimized and validated the conformational ensembles of the amyloid-β 1–40 (Aβ40) and amyloid-β 1–42 (Aβ42) peptides generated by molecular dynamics simulations. We find that both peptides have a diverse set of secondary structure elements including turns, helices, and antiparallel and parallel β-strands. The most significant difference in the structural ensembles of the two peptides is the type of β-hairpins and β-strands they populate. We find that Aβ42 forms a major antiparallel β-hairpin involving the central hydrophobic cluster residues (16–21) with residues 29–36, compatible with known amyloid fibril forming regions, whereas Aβ40 forms an alternative but less populated antiparallel β-hairpin between the central hydrophobic cluster and residues 9–13, that sometimes forms a β-sheet by association with residues 35–37. Furthermore, we show that the two additional C-terminal residues of Aβ42, in particular Ile-41, directly control the differences in the β-strand content found between the Aβ40 and Aβ42 structural ensembles. Integrating the experimental and theoretical evidence accumulated over the last decade, it is now possible to present monomeric structural ensembles of Aβ40 and Aβ42 consistent with available information that produce a plausible molecular basis for why Aβ42 exhibits greater fibrillization rates than Aβ40.     

Amyloids of Alpha-Synuclein Affect the Structure and Dynamics of Supported Lipid Bilayers

Interactions of monomeric alpha-synuclein (αS) with lipid membranes have been suggested to play an important role in initiating aggregation of αS. We have systematically analyzed the distribution and self-assembly of monomeric αS on supported lipid bilayers. We observe that at protein/lipid ratios higher than 1:10, αS forms micrometer-sized clusters, leading to observable membrane defects and decrease in lateral diffusion of both lipids and proteins. An αS deletion mutant lacking amino-acid residues 71–82 binds to membranes, but does not observably affect membrane integrity. Although this deletion mutant cannot form amyloid, significant amyloid formation is observed in the wild-type αS clusters. These results suggest that the process of amyloid formation, rather than binding of αS on membranes, is crucial in compromising membrane integrity.     

F-actin asymmetry and the endoplasmic reticulum–associated TCC-1 protein contribute to stereotypic spindle movements in the Caenorhabditis elegans embryo

The microtubule spindle apparatus dictates the plane of cell cleavage in animal cells. During development, dividing cells control the position of the spindle to determine the size, location, and fate of daughter cells. Spindle positioning depends on pulling forces that act between the cell periphery and astral microtubules. This involves dynein recruitment to the cell cortex by a heterotrimeric G-protein α subunit in complex with a TPR-GoLoco motif protein (GPR-1/2, Pins, LGN) and coiled-coil protein (LIN-5, Mud, NuMA). In this study, we searched for additional factors that contribute to spindle positioning in the one-cell Caenorhabditis elegans embryo. We show that cortical actin is not needed for Gα–GPR–LIN-5 localization and pulling force generation. Instead, actin accumulation in the anterior actually reduces pulling forces, possibly by increasing cortical rigidity. Examining membrane-associated proteins that copurified with GOA-1 Gα, we found that the transmembrane and coiled-coil domain protein 1 (TCC-1) contributes to proper spindle movements. TCC-1 localizes to the endoplasmic reticulum membrane and interacts with UNC-116 kinesin-1 heavy chain in yeast two-hybrid assays. RNA interference of tcc-1 and unc-116 causes similar defects in meiotic spindle positioning, supporting the concept of TCC-1 acting with kinesin-1 in vivo. These results emphasize the contribution of membrane-associated and cortical proteins other than Gα–GPR–LIN-5 in balancing the pulling forces that position the spindle during asymmetric cell division.     

Cross Dimerization of Amyloid-β and αSynuclein Proteins in Aqueous Environment: A Molecular Dynamics Simulations Study

Self-assembly of the intrinsically unstructured proteins, amyloid beta (Aβ) and alpha synclein (αSyn), are associated with Alzheimer’s Disease, and Parkinson’s and Lewy Body Diseases, respectively. Importantly, pathological overlaps between these neurodegenerative diseases, and the possibilities of interactions between Aβ and αSyn in biological milieu emerge from several recent clinical reports and in vitro studies. Nevertheless, there are very few molecular level studies that have probed the nature of spontaneous interactions between these two sequentially dissimilar proteins and key characteristics of the resulting cross complexes. In this study, we have used atomistic molecular dynamics simulations to probe the possibility of cross dimerization between αSyn_1–95 and Aβ _1–42, and thereby gain insights into their plausible early assembly pathways in aqueous environment. Our analyses indicate a strong probability of association between the two sequences, with inter-protein attractive electrostatic interactions playing dominant roles. Principal component analysis revealed significant heterogeneity in the strength and nature of the associations in the key interaction modes. In most, the interactions of repeating Lys residues, mainly in the imperfect repeats ‘KTKEGV’ present in αSyn_1–95 were found to be essential for cross interactions and formation of inter-protein salt bridges. Additionally, a hydrophobicity driven interaction mode devoid of salt bridges, where the non-amyloid component (NAC) region of αSyn_1–95 came in contact with the hydrophobic core of Aβ _1–42 was observed. The existence of such hetero complexes, and therefore hetero assembly pathways may lead to polymorphic aggregates with variations in pathological attributes. Our results provide a perspective on development of therapeutic strategies for preventing pathogenic interactions between these proteins.     

Interaction of Amyloid Inhibitor Proteins with Amyloid Beta Peptides: Insight from Molecular Dynamics Simulations

Knowledge of the detailed mechanism by which proteins such as human αB- crystallin and human lysozyme inhibit amyloid beta (Aβ) peptide aggregation is crucial for designing treatment for Alzheimer''s disease. Thus, unconstrained, atomistic molecular dynamics simulations in explicit solvent have been performed to characterize the Aβ_17–42 assembly in presence of the αB-crystallin core domain and of lysozyme. Simulations reveal that both inhibitor proteins compete with inter-peptide interaction by binding to the peptides during the early stage of aggregation, which is consistent with their inhibitory action reported in experiments. However, the Aβ binding dynamics appear different for each inhibitor. The binding between crystallin and the peptide monomer, dominated by electrostatics, is relatively weak and transient due to the heterogeneous amino acid distribution of the inhibitor surface. The crystallin-bound Aβ oligomers are relatively long-lived, as they form more extensive contact surface with the inhibitor protein. In contrast, a high local density of arginines from lysozyme allows strong binding with Aβ peptide monomers, resulting in stable complexes. Our findings not only illustrate, in atomic detail, how the amyloid inhibitory mechanism of human αB-crystallin, a natural chaperone, is different from that of human lysozyme, but also may aid de novo design of amyloid inhibitors.     

Single-Molecule Studies on PolySUMO Proteins Reveal Their Mechanical Flexibility

Proteins with β-sandwich and β-grasp topologies are resistant to mechanical unfolding as shown by single-molecule force spectroscopy studies. Their high mechanical stability has generally been associated with the mechanical clamp geometry present at the termini. However, there is also evidence for the importance of interactions other than the mechanical clamp in providing mechanical stability, which needs to be tested thoroughly. Here, we report the mechanical unfolding properties of ubiquitin-like proteins (SUMO1 and SUMO2) and their comparison with those of ubiquitin. Although ubiquitin and SUMOs have similar size and structural topology, they differ in their sequences and structural contacts, making them ideal candidates to understand the variations in the mechanical stability of a given protein topology. We observe a two-state unfolding pathway for SUMO1 and SUMO2, similar to that of ubiquitin. Nevertheless, the unfolding forces of SUMO1 (∼130 pN) and SUMO2 (∼120 pN) are lower than that of ubiquitin (∼190 pN) at a pulling speed of 400 nm/s, indicating their lower mechanical stability. The mechanical stabilities of SUMO proteins and ubiquitin are well correlated with the number of interresidue contacts present in their structures. From pulling speed-dependent mechanical unfolding experiments and Monte Carlo simulations, we find that the unfolding potential widths of SUMO1 (∼0.51 nm) and SUMO2 (∼0.33 nm) are much larger than that of ubiquitin (∼0.19 nm), indicating that SUMO1 is six times and SUMO2 is three times mechanically more flexible than ubiquitin. These findings might also be important in understanding the functional differences between ubiquitin and SUMOs.     

Thermodynamic Selection of Steric Zipper Patterns in the Amyloid Cross-β Spine

At the core of amyloid fibrils is the cross-β spine, a long tape of β-sheets formed by the constituent proteins. Recent high-resolution x-ray studies show that the unit of this filamentous structure is a β-sheet bilayer with side chains within the bilayer forming a tightly interdigitating “steric zipper” interface. However, for a given peptide, different bilayer patterns are possible, and no quantitative explanation exists regarding which pattern is selected or under what condition there can be more than one pattern observed, exhibiting molecular polymorphism. We address the structural selection mechanism by performing molecular dynamics simulations to calculate the free energy of incorporating a peptide monomer into a β-sheet bilayer. We test filaments formed by several types of peptides including GNNQQNY, NNQQ, VEALYL, KLVFFAE and STVIIE, and find that the patterns with the lowest binding free energy correspond to available atomistic structures with high accuracy. Molecular polymorphism, as exhibited by NNQQ, is likely because there are more than one most stable structures whose binding free energies differ by less than the thermal energy. Detailed analysis of individual energy terms reveals that these short peptides are not strained nor do they lose much conformational entropy upon incorporating into a β-sheet bilayer. The selection of a bilayer pattern is determined mainly by the van der Waals and hydrophobic forces as a quantitative measure of shape complementarity among side chains between the β-sheets. The requirement for self-complementary steric zipper formation supports that amyloid fibrils form more easily among similar or same sequences, and it also makes parallel β-sheets generally preferred over anti-parallel ones. But the presence of charged side chains appears to kinetically drive anti-parallel β-sheets to form at early stages of assembly, after which the bilayer formation is likely driven by energetics.     

Effects of Familial Mutations on the Monomer Structure of Aβ42

Amyloid beta (Aβ) peptide plays an important role in Alzheimer’s disease. A number of mutations in the Aβ sequence lead to familial Alzheimer’s disease, congophilic amyloid angiopathy, or hereditary cerebral hemorrhage with amyloid. Using molecular dynamics simulations of ∼200 μs for each system, we characterize and contrast the consequences of four pathogenic mutations (Italian, Dutch, Arctic, and Iowa) for the structural ensemble of the Aβ monomer. The four familial mutations are found to have distinct consequences for the monomer structure.Amyloid beta (Aβ) peptides have long been thought to play a central role in Alzheimer’s disease (AD). Usually 40 or 42 residues in length, Aβ peptides are proteolytic products of the Aβ precursor protein and they aggregate to form the fibrillar plaques in AD patients’ brains. Besides fibrillar plaques, Aβ oligomers are also neurotoxic. The significance and nature of Aβ oligomerization has recently become a focus of intensive research studies and debates (1,2). Notably, numerous pathogenic mutations have been identified in the Aβ precursor protein sequence and in the enzymes involved in Aβ processing (3). These mutations generally lead to early onset of AD or cerebral amyloid angiopathy. Understanding how the pathogenic mutations alter Aβ oligomerization/aggregation is essential to our understanding of the disease mechanism.Four of these pathogenic mutations (Italian E22K, Dutch E22Q, Arctic E22G, and Iowa D23N) cluster in the region of E22 and D23 in the Aβ sequence (distal from proteolytic cleavage sites) and they have higher neurotoxicity compared to wild-type (WT) Aβ (4). These mutations are thought to modify the physicochemistry of the peptide. For example, kinetic studies (4) show that the E22K and E22Q mutations lead to faster peptide aggregation, whereas the E22G and D23N mutations result in slightly slower aggregation than WT Aβ₄₂ (although the E22G mutation shows increased protofibril formation (5)). Recent solid-state NMR studies also suggest that rather than the in-register β-sheet conformation adopted by WT Aβ, the Iowa D23N mutant forms amyloid fibrils with antiparallel β-sheet structure (6).To understand how the mutations modify the peptide oligomerization/aggregation it is critical to characterize the starting point of the process, the monomers. Unfortunately, investigating the early phase of the oligomerization process experimentally is a challenging task due to the high aggregation propensity of Aβ and its intrinsic disorder. Therefore, a number of computational approaches have been adopted to investigate the consequences of mutations for the monomer structure (7–16). However, due to the high computational demands of explicit-solvent molecular dynamics (MD) simulations to simulate full-length Aβ peptides, most of these computational studies are either on Aβ fragments (to decrease the system size) using explicit-solvent simulations (8–12) or on full-length Aβ using implicit-solvent simulations (which are less computationally demanding and enable longer simulation times, but lack explicit water molecules in the simulations to fully describe water-peptide interactions) (13–15). In a very recent report, explicit-solvent simulations were used to study the effects of the E22Q mutation on full-length Aβ; however, rather limited data (<10 μs) were collected (16). Thus, characterizing full-length Aβ monomers remains quite a daunting task even with simulations.To characterize the effects of mutations on full-length Aβ monomer using explicit-solvent MD simulations, we employed distributed computing (17) to simulate the WT Aβ₄₂, Aβ₄₂-E22K, Aβ₄₂-E22Q, Aβ₄₂-E22G, and Aβ₄₂-D23N monomers. MD simulations of >200 μs were performed for each system and AMBER ff99sb (18) and the tip3p water model (19) were used for force field parameters. Peptide configurations in the MD trajectories were clustered with the root mean-square deviation metric to identify representative conformations (i.e., states) and transitions between these states were counted. Markov state model analysis was then performed where the master equations were solved and the equilibrium population of each state deduced (20). Details of the MD simulation procedures and Markov state model analysis can be found in the Supporting Material.Each of the five Aβ monomer systems exhibits great structural diversity and can only be characterized in an ensemble fashion (rather than described by a handful of representative configurations). This is in accord with the notion that full-length Aβ peptides are intrinsically disordered (21,22). Using the Dictionary of Secondary Structure of Proteins program (23) to assign secondary structure, it is clear that the five Aβ monomer systems are found overall not well structured, although small β-hairpins and α-helices are observed. In Fig. 1 we plot the residue-dependent extended β propensity and α-helix propensity, in the top and bottom panels, respectively, for each Aβ monomer system. Although we are reasonably confident of the convergence behavior of the α-helix propensity, we note that the convergence of the extended β-propensity might be more challenging and demand a much longer sampling time than the current aggregate simulation time of ∼200 μs (24).Open in a separate windowFigure 1Ensemble-averaged %population of β-strand (top) and α-helix (bottom) propensity for all five monomer systems. The sequence of the WT Aβ₄₂ is given on the x axis.We observe in Fig. 1 that all five Aβ monomer systems share a rather similar residue-dependent tendency to form an extended β-structure, although minor differences are present. On the other hand, these pathogenic mutations alter the α-helix propensity quite significantly. The E22K and E22Q mutations increase the α-helix propensity in the region of residues 20–23. All four mutations (E22K, E22Q, E22G, and D23N) decrease the α-helix propensity in the region of residues 33–36.Notably, we find that in all five systems only short stretches of α-helices are formed. That is, when a residue is involved in α-helix formation, it participates in forming mostly short helical segments (consisting of only four helical residues). To provide more insight into the changes of α-helix propensity due to the mutations, in Fig. S1 we plot the tendency of forming short α-helices along the sequence for all five systems. Each data point in Fig. S1 represents the propensity to form an α-helix of four residues in length, ending at the specific residue. For example, in the structural ensemble adopted by the WT peptide, ∼5.5% of the conformations have a short α-helix of size four, involving residues 15–18. We see from Fig. S1 that the E22K and E22Q mutations induce the formation of two short helices in residues 19–22 and 20–23. The higher α-helix propensity in this region for the E22K mutant compared to the WT was previously attributed to the elimination of the electrostatic repulsion between E22 and D23 in the WT by the mutation and the longer aliphatic chain of K22 in the mutant compared to E22 in the WT (9,22). This is consistent with the observation that the E22Q mutation also induces helix formation in this region (by eliminating the electrostatic repulsion between E22 and D23 in the WT) but to a lesser extent, possibly due to the shorter aliphatic chain of Q22 compared to K22.In the E22G mutant, although the mutation eliminates the electrostatic repulsion between E22 and D23 in the WT peptide, glycine is known to be a helix breaker (25), leading to diminished α-helix propensity in the region around residue G22 seen in Fig. S1.In the D23N mutant, although the mutation eliminates the electrostatic repulsion between E22 and D23 in the WT peptide, it does not induce (or rather even slightly decreases) helix formation around residue 23. This may be due to the short aliphatic chain of N23 but it is possible that the mutation induces some nonlocal effects on the peptide structure, disfavoring helix formation in this region.It is worth noting that all four mutations (E22K, E22Q, E22G, and D23N) virtually eliminate the α-helix propensity in the region of residues 33–36. This region is rather far away from the mutation sites in sequence but its α-helix propensity is nonetheless affected. The origin of such a nonlocal effect is less straightforward to explain and further analysis will aid untangling this behavior. Nonetheless, the diminished α-helix propensity in the region of residues 33–36 appears to be a consistent feature across all four mutants.The four mutations studied here (E22K, E22Q, E22G, and D23N) have been thought to modify the physicochemistry of the peptide and alter the oligomerization/aggregation process, leading to higher neurotoxicity. In predicting intrinsic aggregation propensities using peptide sequences, all four mutants are suggested to be more aggregation prone (26). On the other hand, kinetic studies show that only the E22K and E22Q mutants aggregate more quickly, whereas the E22G and D23N mutations result in slightly slower aggregation than WT Aβ₄₂ (4). Our simulation results suggest these pathogenic mutations have complicated effects on the monomer structure—all four mutations decrease helix propensity in residues 33–36, whereas only the E22K and E22Q mutations increase helix propensity in residues 20–23. It is interesting to note that α-helix propensity is generally thought to anticorrelate with aggregation propensity; however, recent studies have suggested an important role of α-helical intermediates in amyloid oligomerization (27–29). Our studies suggest that it would be of great value to investigate how the distinct patterns of α-helix propensity in these five systems may propagate to give rise to different oligomerization kinetics or even mechanisms. The pathogenic mutations studied here have complex effects on the oligomerization of the peptide. The characterization of the monomer structural ensembles reported here should aid understanding of such an important and complicated process.     

Effects of Hydroxylated Carbon Nanotubes on the Aggregation of Aβ16–22 Peptides: A Combined Simulation and Experimental Study

The pathogenesis of Alzheimer’s disease (AD) is associated with the aggregation of amyloid-β (Aβ) peptides into toxic aggregates with β-sheet character. In a previous computational study, we showed that pristine single-walled carbon nanotubes (SWCNTs) can inhibit the formation of β-sheet-rich oligomers in the central hydrophobic core fragment of Aβ (Aβ_16–22). However, the poor solubility of SWCNTs in water hinders their use in biomedical applications and nanomedicine. Here, we investigate the influence of hydroxylated SWCNT, a water-soluble SWCNT derivative, on the aggregation of Aβ_16–22 peptides using all-atom explicit-water replica exchange molecular dynamics simulations. Our results show that hydroxylated SWCNTs can significantly inhibit β-sheet formation and shift the conformations of Aβ_16–22 oligomers from ordered β-sheet-rich structures toward disordered coil aggregates. Detailed analyses of the SWCNT-Aβ interaction reveal that the inhibition of β-sheet formation by hydroxylated SWCNTs mainly results from strong electrostatic interactions between the hydroxyl groups of SWCNTs and the positively charged residue K16 of Aβ_16–22 and hydrophobic and aromatic stacking interactions between SWCNTs and F19 and F20. In addition, our atomic force microscopy and thioflavin T fluorescence experiments confirm the inhibitory effect of both pristine and hydroxylated SWCNTs on Aβ_16–22 fibrillization, in support of our previous and present replica exchange molecular dynamics simulation results. These results demonstrate that hydroxylated SWCNTs efficiently inhibit the aggregation of Aβ_16–22; in addition, they offer molecular insight into the inhibition mechanism, thus providing new clues for the design of therapeutic drugs against amyloidosis.     

A Flexible Nanoarray Approach for the Assembly and Probing of Molecular Complexes

Immobilization is a key step involved in probing molecular interactions using single-molecule force spectroscopy methods, including atomic force microscopy (AFM). To our knowledge, we describe a novel approach termed flexible nanoarray (FNA) in which the interaction between the two internally immobilized amyloid β peptides is measured by pulling of the tether. The FNA tether was synthesized with nonnucleotide phosphoramidite monomers using the DNA synthesis chemistry. The two anchoring points for immobilization of the peptides inside the tether were incorporated at defined distances between them and from the ends of the polymer. Decamers of amyloid β peptide capable of dimer formation were selected as a test system. The formation of the peptide dimers was verified by AFM force spectroscopy by pulling the tether at the ends. In these experiments, the thiolated end of the FNA tether was covalently immobilized on the AFM substrate functionalized with maleimide. The other end of the FNA tether was functionalized with biotin to form a noncovalent link with the streptavidin functionalized AFM tip during the approach stage. The dimers’ rupture fingerprint was unambiguously identified on the force curves by its position and the force value. The FNA design allowed reversible experiments in which the monomers were allowed to associate after the rupture of the dimers by performing the approach stage before the rupture of the biotin-streptavidin link. This suggests that the FNA technique is capable of analyzing multiple intermolecular interactions in the same molecular complex. The computational analysis showed that the tethered peptides assemble into the same dimer structure as that formed by nontethered peptides, suggesting that the FNA tether has the necessary flexibility to enable assembly of the dimer even during the course of the force spectroscopy experiment.     

Stability of Transmembrane Amyloid β-Peptide and Membrane Integrity Tested by Molecular Modeling of Site-Specific Aβ 42 Mutations

Interactions of the amyloid β-protein (Aβ) with neuronal cell membranes, leading to the disruption of membrane integrity, are considered to play a key role in the development of Alzheimer’s disease. Natural mutations in Aβ ₄₂, such as the Arctic mutation (E22G) have been shown to increase Aβ ₄₂ aggregation and neurotoxicity, leading to the early-onset of Alzheimer’s disease. A correlation between the propensity of Aβ ₄₂ to form protofibrils and its effect on neuronal dysfunction and degeneration has been established. Using rational mutagenesis of the Aβ ₄₂ peptide it was further revealed that the aggregation of different Aβ ₄₂ mutants in lipid membranes results in a variety of polymorphic aggregates in a mutation dependent manner. The mutant peptides also have a variable ability to disrupt bilayer integrity. To further test the connection between Aβ ₄₂ mutation and peptide–membrane interactions, we perform molecular dynamics simulations of membrane-inserted Aβ ₄₂ variants (wild-type and E22G, D23G, E22G/D23G, K16M/K28M and K16M/E22G/D23G/K28M mutants) as β-sheet monomers and tetramers. The effects of charged residues on transmembrane Aβ ₄₂ stability and membrane integrity are analyzed at atomistic level. We observe an increased stability for the E22G Aβ ₄₂ peptide and a decreased stability for D23G compared to wild-type Aβ ₄₂, while D23G has the largest membrane-disruptive effect. These results support the experimental observation that the altered toxicity arising from mutations in Aβ is not only a result of the altered aggregation propensity, but also originates from modified Aβ interactions with neuronal membranes.     

Alzheimer’s Protective A2T Mutation Changes the Conformational Landscape of the Aβ1–42 Monomer Differently Than Does the A2V Mutation

The aggregation of amyloid-β (Aβ) peptides plays a crucial role in the etiology of Alzheimer’s disease (AD). Recently, it has been reported that an A2T mutation in Aβ can protect against AD. Interestingly, a nonpolar A2V mutation also has been found to offer protection against AD in the heterozygous state, although it causes early-onset AD in homozygous carriers. Since the conformational landscape of the Aβ monomer is known to directly contribute to the early-stage aggregation mechanism, it is important to characterize the effects of the A2T and A2V mutations on Aβ_1–42 monomer structure. Here, we have performed extensive atomistic replica-exchange molecular dynamics simulations of the solvated wild-type (WT), A2V, and A2T Aβ_1–42 monomers. Our simulations reveal that although all three variants remain as collapsed coils in solution, there exist significant structural differences among them at shorter timescales. A2V exhibits an enhanced double-hairpin population in comparison to the WT, similar to those reported in toxic WT Aβ_1–42 oligomers. Such double-hairpin formation is caused by hydrophobic clustering between the N-terminus and the central and C-terminal hydrophobic patches. In contrast, the A2T mutation causes the N-terminus to engage in unusual electrostatic interactions with distant residues, such as K16 and E22, resulting in a unique population comprising only the C-terminal hairpin. These findings imply that a single A2X (where X = V or T) mutation in the primarily disordered N-terminus of the Aβ_1–42 monomer can dramatically alter the β-hairpin population and switch the equilibrium toward alternative structures. The atomistically detailed, comparative view of the structural landscapes of A2V and A2T variant monomers obtained in this study can enhance our understanding of the mechanistic differences in their early-stage aggregation.     

Effect of Proline Mutations on the Monomer Conformations of Amylin

The formation of human islet amyloid polypeptide (hIAPP) is implicated in the loss of pancreatic β-cells in type II diabetes. Rat amylin, which differs from human amylin at six residues, does not lead to formation of amyloid fibrils. Pramlintide is a synthetic analog of human amylin that shares three proline substitutions with rat amylin. Pramlintide has a much smaller propensity to form amyloid aggregates and has been widely prescribed in amylin replacement treatment. It is known that the three prolines attenuate β-sheet formation. However, the detailed effects of these proline substitutions on full-length hIAPP remain poorly understood. In this work, we use molecular simulations and bias-exchange metadynamics to investigate the effect of proline substitutions on the conformation of the hIAPP monomer. Our results demonstrate that hIAPP can adopt various β-sheet conformations, some of which have been reported in experiments. The proline substitutions perturb the formation of long β-sheets and reduce their stability. More importantly, we find that all three proline substitutions of pramlintide are required to inhibit β conformations and stabilize the α-helical conformation. Fewer substitutions do not have a significant inhibiting effect.     

Force Measurements of the α5β1 Integrin–Fibronectin Interaction

The interaction of the α₅β₁ integrin and its ligand, fibronectin (FN), plays a crucial role in the adhesion of cells to the extracellular matrix. An important intrinsic property of the α₅β₁/FN interaction is the dynamic response of the complex to a pulling force. We have carried out atomic force microscopy measurements of the interaction between α₅β₁ and a fibronectin fragment derived from the seventh through tenth type III repeats of FN (i.e., FN7-10) containing both the arg-gly-asp (RGD) sequence and the synergy site. Direct force measurements obtained from an experimental system consisting of an α₅β₁ expressing K562 cell attached to the atomic force microscopy cantilever and FN7-10 adsorbed on a substrate were used to determine the dynamic response of the α₅β₁/FN7-10 complex to a pulling force. The experiments were carried out over a three-orders-of-magnitude change in loading rate and under conditions that allowed for detection of individual α₅β₁/FN7-10 interactions. The dynamic rupture force of the α₅β₁/FN7-10 complex revealed two regimes of loading: a fast loading regime (>10,000 pN/s) and a slow loading regime (<10,000 pN/s) that characterize the inner and outer activation barriers of the complex, respectively. Activation by TS2/16 antibody increased both the frequency of adhesion and elevated the rupture force of the α₅β₁/wild type FN7-10 complex to higher values in the slow loading regime. In experiments carried out with a FN7-10 RGD deleted mutant, the force measurements revealed that both inner and outer activation barriers were suppressed by the mutation. Mutations to the synergy site of FN, however, suppressed only the outer barrier activation of the complex. For both the RGD and synergy deletions, the frequency of adhesion was less than that of the wild type FN7-10, but was increased by integrin activation. The rupture force of these mutants was only slightly less than that of the wild type, and was not increased by activation. These results suggest that integrin activation involved a cooperative interaction with both the RGD and synergy sites.     

β-Hairpin-Mediated Formation of Structurally Distinct Multimers of Neurotoxic Prion Peptides

Protein misfolding disorders are associated with conformational changes in specific proteins, leading to the formation of potentially neurotoxic amyloid fibrils. During pathogenesis of prion disease, the prion protein misfolds into β-sheet rich, protease-resistant isoforms. A key, hydrophobic domain within the prion protein, comprising residues 109–122, recapitulates many properties of the full protein, such as helix-to-sheet structural transition, formation of fibrils and cytotoxicity of the misfolded isoform. Using all-atom, molecular simulations, it is demonstrated that the monomeric 109–122 peptide has a preference for α-helical conformations, but that this peptide can also form β-hairpin structures resulting from turns around specific glycine residues of the peptide. Altering a single amino acid within the 109–122 peptide (A₁₁₇V, associated with familial prion disease) increases the prevalence of β-hairpin formation and these observations are replicated in a longer peptide, comprising residues 106–126. Multi-molecule simulations of aggregation yield different assemblies of peptide molecules composed of conformationally-distinct monomer units. Small molecular assemblies, consistent with oligomers, comprise peptide monomers in a β-hairpin-like conformation and in many simulations appear to exist only transiently. Conversely, larger assemblies are comprised of extended peptides in predominately antiparallel β-sheets and are stable relative to the length of the simulations. These larger assemblies are consistent with amyloid fibrils, show cross-β structure and can form through elongation of monomer units within pre-existing oligomers. In some simulations, assemblies containing both β-hairpin and linear peptides are evident. Thus, in this work oligomers are on pathway to fibril formation and a preference for β-hairpin structure should enhance oligomer formation whilst inhibiting maturation into fibrils. These simulations provide an important new atomic-level model for the formation of oligomers and fibrils of the prion protein and suggest that stabilization of β-hairpin structure may enhance cellular toxicity by altering the balance between oligomeric and fibrillar protein assemblies.