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.     

Binding Modes of Phthalocyanines to Amyloid β Peptide and Their Effects on Amyloid Fibril Formation

The inherent tendency of proteins to convert from their native states into amyloid aggregates is associated with a range of human disorders, including Alzheimer’s and Parkinson’s diseases. In that sense, the use of small molecules as probes for the structural and toxic mechanism related to amyloid aggregation has become an active area of research. Compared with other compounds, the structural and molecular basis behind the inhibitory interaction of phthalocyanine tetrasulfonate (PcTS) with proteins such as αS and tau has been well established, contributing to a better understanding of the amyloid aggregation process in these proteins. We present here the structural characterization of the binding of PcTS and its Cu(II) and Zn(II)-loaded forms to the amyloid β-peptide (Aβ) and the impact of these interactions on the peptide amyloid fibril assembly. Elucidation of the PcTS binding modes to Aβ₄₀ revealed the involvement of specific aromatic and hydrophobic interactions in the formation of the Aβ₄₀-PcTS complex, ascribed to a binding mode in which the planarity and hydrophobicity of the aromatic ring system in the phthalocyanine act as main structural determinants for the interaction. Our results demonstrated that formation of the Aβ₄₀-PcTS complex does not interfere with the progression of the peptide toward the formation of amyloid fibrils. On the other hand, conjugation of Zn(II) but not Cu(II) at the center of the PcTS macrocyclic ring modified substantially the binding profile of this phthalocyanine to Aβ₄₀ and became crucial to reverse the effects of metal-free PcTS on the fibril assembly of the peptide. Overall, our results provide a firm basis to understand the structural rules directing phthalocyanine-protein interactions and their implications on the amyloid fibril assembly of the target proteins; in particular, our results contradict the hypothesis that PcTS might have similar mechanisms of action in slowing the formation of a variety of pathological aggregates.     

Simulations of cross-amyloid aggregation of amyloid-β and islet amyloid polypeptide fragments

Amyloid-β (Aβ) and islet amyloid polypeptide (IAPP) are small peptides, classified as amyloids, that have the potential to self-assemble and form cytotoxic species, such as small soluble oligomers and large insoluble fibrils. The formation of Aβ aggregates facilitates the progression of Alzheimer’s disease (AD), while IAPP aggregates induce pancreatic β-cell apoptosis, leading to exacerbation of type 2 diabetes (T2D). Cross-amyloid interactions between Aβ and IAPP have been described both in vivo and in vitro, implying the role of Aβ or IAPP as modulators of cytotoxic self-aggregation of each species, and suggesting that Aβ-IAPP interactions are a potential molecular link between AD and T2D. Using molecular dynamics (MD) simulations, “hotspot” regions of the two peptides were studied to understand the formation of hexamers in a heterogeneous and homogeneous peptide-containing environment. Systems of only Aβ_(16–22) peptides formed antiparallel, β-barrel-like structures, while systems of only IAPP_(20–29) peptides formed stacked, parallel β-sheets and had relatively unstable aggregation structures after 2 μs of simulation time. Systems containing both Aβ and IAPP (1:1 ratio) hexamers showed antiparallel, β-barrel-like structures, with an interdigitated arrangement of Aβ_(16–22) and IAPP_(20–29). These β-barrel structures have features of cytotoxic amyloid species identified in previous literature. Ultimately, this work seeks to provide atomistic insight into both the mechanism behind cross-amyloid interactions and structural morphologies of these toxic amyloid species.     

Structural Determination of Aβ25-35 Micelles by Molecular Dynamics Simulations

Amyloid-β (Aβ) peptides and other amyloidogenic proteins can form a wide range of soluble oligomers of varied morphologies at the early aggregation stage, and some of these oligomers are biologically relevant to the pathogenesis of Alzheimer's disease. Spherical micelle-like oligomers have been often observed for many different types of amyloids. Here, we report a hybrid computational approach to systematically construct, search, optimize, and rank soluble micelle-like Aβ_25-35 structures with different side-chain packings at the atomic level. Simulations reveal for the first time, to our knowledge, that two Aβ micelles with antiparallel peptide organization and distinct surface hydrophobicity display high structural stability. Stable micelles experience a slow secondary structural transition from turn to α-helix. Energetic analysis coupled with computational mutagenesis reveals that van der Waals and solvation energies play a more pronounced role in stabilizing the micelles, whereas the electrostatic energies present a stable but minor energetic contribution to peptide assemblies. Modeled Aβ micelles with shapes and dimensions similar to those of experimentally derived spherical structures also provide detailed information about the roles of structural dynamics and transition in the formation of amyloid fibrils. The strong binding affinity of our micelles to antibodies implies that micelles may be a biologically relevant species.     

Induced β-Barrel Formation of the Alzheimer's Aβ25-35 Oligomers on Carbon Nanotube Surfaces: Implication for Amyloid Fibril Inhibition

Recent experimental studies show that carbon nanotubes impact the aggregation process of proteins associated with neurodegenerative diseases. However, the details of molecular interactions between proteins and carbon nanotubes are still not well understood. In this study, we investigate the initial adsorption features and dynamics of the Alzheimer's amyloid-β peptide spanning residues 25-35 (Aβ25-35) on a single-walled carbon nanotube (SWNT) surface using fully atomic molecular dynamics simulations (MD) in explicit solvent. The initial configurations of the Aβ25-35 peptides consist of two preformed bilayer β-sheets, each with four or five β-strands in parallel or mixed antiparallel-parallel orientations. Our simulations show, for what we believe is the first time, that two disjointed Aβ25-35 β-sheets with mixed antiparallel-parallel strands can assemble into β-barrels wrapping the SWNT. In contrast, both simulations of Aβ25-35 without SWNT, and simulations of SWNT−Aβ25-35 with purely parallel β-strands, lead to disordered aggregates. We find that Aβ25-35 β-barrel formation involves at least two steps: i), curving of the Aβ25-35 β-sheets as a result of strong hydrophobic interactions with carbon nanotube concomitantly with dehydration of the SWNT-peptide interface; and ii), intersheet backbone hydrogen bond formation with fluctuating intrasheet hydrogen bonds. Detailed analysis of the conversion shows that β-barrel formation on SWNT surface results from the interplay of dehydration and peptide-SWNT/peptide-peptide interactions. Implications of our results on amyloid fibril inhibition are discussed.     

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.     

Intrinsic Determinants of Neurotoxic Aggregate Formation by the Amyloid β Peptide

The extent to which proteins aggregate into distinct structures ranging from prefibrillar oligomers to amyloid fibrils is key to the pathogenesis of many age-related degenerative diseases. We describe here for the Alzheimer's disease-related amyloid β peptide (Aβ) an investigation of the sequence-based determinants of the balance between the formation of prefibrillar aggregates and amyloid fibrils. We show that by introducing single-point mutations, it is possible to convert the normally harmless Aβ₄₀ peptide into a pathogenic species by increasing its relative propensity to form prefibrillar but not fibrillar aggregates, and, conversely, to abolish the pathogenicity of the highly neurotoxic E22G Aβ₄₂ peptide by reducing its relative propensity to form prefibrillar species rather than mature fibrillar ones. This observation can be rationalized by the demonstration that whereas regions of the sequence of high aggregation propensity dominate the overall tendency to aggregate, regions with low intrinsic aggregation propensities exert significant control over the balance of the prefibrillar and fibrillar species formed, and therefore play a major role in determining the neurotoxicity of the Aβ peptide.     

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.     

Temperature-induced dissociation of Abeta monomers from amyloid fibril

Using all-atom molecular dynamics, we study the temperature-induced dissociation of Aβ monomers from the fibril protofilament. To accelerate conformational sampling, simulations are performed at elevated temperatures and peptide concentrations. By computing free energy disconnectivity graphs we mapped the free energy landscape of monomers on the surface of Aβ fibril. We found that Aβ monomers sample diverse sets of low free energy states with different degrees of association with the fibril. Some of these states have residual amounts of fibril interactions, whereas others lack fibril-like content. Generally, Aβ monomers with partially formed fibril-like interactions have the lowest free energies, but their backbone conformations may differ considerably from those in the fibril interior. Overall, Aβ amyloid protofilaments seem to be highly resistant to thermal dissociation. Monomer dissociation from the fibril edge proceeds via multiple stages and pathways. Our simulation findings are discussed in the context of recent experimental results.     

Highly Disordered Amyloid-β Monomer Probed by Single-Molecule FRET and MD Simulation

Monomers of amyloid-β (Aβ) protein are known to be disordered, but there is considerable controversy over the existence of residual or transient conformations that can potentially promote oligomerization and fibril formation. We employed single-molecule Förster resonance energy transfer (FRET) spectroscopy with site-specific dye labeling using an unnatural amino acid and molecular dynamics simulations to investigate conformations and dynamics of Aβ isoforms with 40 (Aβ40) and 42 residues (Aβ42). The FRET efficiency distributions of both proteins measured in phosphate-buffered saline at room temperature show a single peak with very similar FRET efficiencies, indicating there is apparently only one state. 2D FRET efficiency-donor lifetime analysis reveals, however, that there is a broad distribution of rapidly interconverting conformations. Using nanosecond fluorescence correlation spectroscopy, we measured the timescale of the fluctuations between these conformations to be ～35 ns, similar to that of disordered proteins. These results suggest that both Aβ40 and Aβ42 populate an ensemble of rapidly reconfiguring unfolded states, with no long-lived conformational state distinguishable from that of the disordered ensemble. To gain molecular-level insights into these observations, we performed molecular dynamics simulations with a force field optimized to describe disordered proteins. We find, as in experiments, that both peptides populate configurations consistent with random polymer chains, with the vast majority of conformations lacking significant secondary structure, giving rise to very similar ensemble-averaged FRET efficiencies.     

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.     

Scavenging amyloid oligomers from neurons with silica nanobowls: Implications for amyloid diseases

Amyloid-β (Aβ) oligomers are toxic species implicated in Alzheimer’s disease (AD). The prevailing hypothesis implicates a major role of membrane-associated amyloid oligomers in AD pathology. Our silica nanobowls (NB) coated with lipid-polymer have submicromolar affinity for Aβ binding. We demonstrate that NB scavenges distinct fractions of Aβs in a time-resolved manner from amyloid precursor protein-null neuronal cells after incubation with Aβ. At short incubation times in cell culture, NB-Aβ seeds have aggregation kinetics resembling that of extracellular fraction of Aβ, whereas at longer incubation times, NB-Aβ seeds scavenge membrane-associated Aβ. Aβ aggregates can be eluted from NB surfaces by mechanical agitation and appear to retain their aggregation driving domains as seen in seeding aggregation experiments. These results demonstrate that the NB system can be used for time-resolved separation of toxic Aβ species from biological samples for characterization and in diagnostics. Scavenging membrane-associated amyloids using lipid-functionalized NB without chemical manipulation has wide applications in the diagnosis and therapy of AD and other neurodegenerative diseases, cancer, and cardiovascular conditions.     

The amyloid precursor protein forms plasmalemmal clusters via its pathogenic amyloid-β domain

The amyloid precursor protein (APP) is a large, ubiquitous integral membrane protein with a small amyloid-β (Aβ) domain. In the human brain, endosomal processing of APP produces neurotoxic Aβ-peptides, which are involved in Alzheimer's disease. Here, we show that the Aβ sequence exerts a physiological function when still present in the unprocessed APP molecule. From the extracellular site, Aβ concentrates APP molecules into plasmalemmal membrane protein clusters. Moreover, Aβ stabilization of clusters is a prerequisite for their targeting to endocytic clathrin structures. Therefore, we conclude that the Aβ domain directly mediates a central step in APP trafficking, driving its own conversion into neurotoxic peptides.     

Atomic-Scale Simulations Confirm that Soluble β-Sheet-Rich Peptide Self-Assemblies Provide Amyloid Mimics Presenting Similar Conformational Properties

The peptide self-assembly mimic (PSAM) from the outer surface protein A (OspA) can form highly stable but soluble β-rich self-assembly-like structures similar to those formed by native amyloid-forming peptides. However, unlike amyloids that predominantly form insoluble aggregates, PSAMs are highly water-soluble. Here, we characterize the conformations of these soluble β-sheet-rich assemblies. We simulate PSAMs with different-sized β-sheets in the presence and absence of end-capping proteins using all-atom explicit-solvent molecular dynamics, comparing the structural stability, conformational dynamics, and association force. Structural and free-energy comparisons among β-sheets with different numbers of layers and sequences indicate that in similarity to amyloids, the intersheet side chain-side chain interactions and hydrogen bonds combined with intrasheet salt bridges are the major driving forces in stabilizing the overall structural organization. A detailed structural analysis shows that in similarity to amyloid fibrils, all wild-type and mutated PSAM structures display twisted and bent β-sheets to some extent, implying that a twisted and bent β-sheet is a general motif of β-rich assemblies. Thus, our studies indicate that soluble β-sheet-rich peptide self-assemblies can provide good amyloid mimics, and as such confirm on the atomic scale that they are excellent systems for amyloid studies. These results provide further insight into the usefulness of such mimics for nanostructure design.     

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.     

An Unusual Hydrophobic Core Confers Extreme Flexibility to HEAT Repeat Proteins

Alpha-solenoid proteins are suggested to constitute highly flexible macromolecules, whose structural variability and large surface area is instrumental in many important protein-protein binding processes. By equilibrium and nonequilibrium molecular dynamics simulations, we show that importin-β, an archetypical α-solenoid, displays unprecedentedly large and fully reversible elasticity. Our stretching molecular dynamics simulations reveal full elasticity over up to twofold end-to-end extensions compared to its bound state. Despite the absence of any long-range intramolecular contacts, the protein can return to its equilibrium structure to within 3 Å backbone RMSD after the release of mechanical stress. We find that this extreme degree of flexibility is based on an unusually flexible hydrophobic core that differs substantially from that of structurally similar but more rigid globular proteins. In that respect, the core of importin-β resembles molten globules. The elastic behavior is dominated by nonpolar interactions between HEAT repeats, combined with conformational entropic effects. Our results suggest that α-solenoid structures such as importin-β may bridge the molecular gap between completely structured and intrinsically disordered proteins.     

Structural and Dynamic Properties of the Human Prion Protein

Prion diseases involve the conformational conversion of the cellular prion protein (PrP^C) to its misfolded pathogenic form (PrP^Sc). To better understand the structural mechanism of this conversion, we performed extensive all-atom, explicit-solvent molecular-dynamics simulations for three structures of the wild-type human PrP (huPrP) at different pH values and temperatures. Residue 129 is polymorphic, being either Met or Val. Two of the three structures have Met in position 129 and the other has Val. Lowering the pH or raising the temperature induced large conformational changes of the C-terminal globular domain and increased exposure of its hydrophobic core. In some simulations, HA and its preceding S1-HA loop underwent large displacements. The C-terminus of HB was unstable and sometimes partially unfolded. Two hydrophobic residues, Phe-198 and Met-134, frequently became exposed to solvent. These conformational changes became more dramatic at lower pH or higher temperature. Furthermore, Tyr-169 and the S2-HB loop, or the X-loop, were different in the starting structures but converged to common conformations in the simulations for the Met-129, but not the Val-129, protein. α-Strands and β-strands formed in the initially unstructured N-terminus. α-Strand propensity in the N-terminus was different between the Met-129 and Val129 proteins, but β-strand propensity was similar. This study reveals detailed structural and dynamic properties of huPrP, providing insight into the mechanism of the conversion of PrP^C to PrP^Sc.     

Alzheimer's Disease Drug Candidates Stabilize A-β Protein Native Structure by Interacting with the Hydrophobic Core

Deposition of amyloid fibrils, consisting primarily of Aβ₄₀ and Aβ₄₂ peptides, in the extracellular space in the brain is a major characteristic of Alzheimer's disease (AD). We recently developed new (to our knowledge) drug candidates for AD that inhibit the fibril formation of Aβ peptides and eliminate their neurotoxicity. We performed all-atom molecular-dynamics simulations on the Aβ₄₂ monomer at its α-helical conformation and a pentamer fibril fragment of Aβ₄₂ peptide with or without LRL and fluorene series compounds to investigate the mechanism of inhibition. The results show that the active drug candidates, LRL22 (EC₅₀ = 0.734 μM) and K162 (EC₅₀ = 0.080 μM), stabilize hydrophobic core I of Aβ₄₂ peptide (residues 17–21) to its α-helical conformation by interacting specifically in this region. The nonactive drug candidates, LRL27 (EC₅₀ > 10 μM) and K182 (EC₅₀ > 5 μM), have little to no similar effect. This explains the different behavior of the drug candidates in experiments. Of more importance, this phenomenon indicates that hydrophobic core I of the Aβ₄₂ peptide plays a major mechanistic role in the formation of amyloid fibrils, and paves the way for the development of new drugs against AD.     

Binding and Cleavage of E. coli HUβ by the E. coli Lon Protease

The Escherichia coli Lon protease degrades the E. coli DNA-binding protein HUβ, but not the related protein HUα. Here we show that the Lon protease binds to both HUβ and HUα, but selectively degrades only HUβ in the presence of ATP. Mass spectrometry of HUβ peptide fragments revealed that region K18-G22 is the preferred cleavage site, followed in preference by L36-K37. The preferred cleavage site was further refined to A20-A21 by constructing and testing mutant proteins; Lon degraded HUβ-A20Q and HUβ-A20D more slowly than HUβ. We used optical tweezers to measure the rupture force between HU proteins and Lon; HUα, HUβ, and HUβ-A20D can bind to Lon, and in the presence of ATP, the rupture force between each of these proteins and Lon became weaker. Our results support a mechanism of Lon protease cleavage of HU proteins in at least three stages: binding of Lon with the HU protein (HUβ, HUα, or HUβ-A20D); hydrolysis of ATP by Lon to provide energy to loosen the binding to the HU protein and to allow an induced-fit conformational change; and specific cleavage of only HUβ.