Formation and growth of oligomers: a Monte Carlo study of an amyloid tau fragment

Structures of Oligomeric States of Tau Protein, Amyloid-β, α-Synuclein and Prion Protein Implicated in Alzheimer's Disease, Parkinson's Disease and Prionopathies

In this review, we focus on the biophysical and structural aspects of the oligomeric states of physiologically intrinsically disordered proteins and peptides tau, amyloid-β and α-synuclein and partly disordered prion protein and their isolations from animal models and human brains. These protein states may be the most toxic agents in the pathogenesis of Alzheimer's and Parkinson's disease. It was shown that oligomers are important players in the aggregation cascade of these proteins. The structural information about these structural states has been provided by methods such as solution and solid-state NMR, cryo-EM, crosslinking mass spectrometry, AFM, TEM, etc., as well as from hybrid structural biology approaches combining experiments with computational modelling and simulations. The reliable structural models of these protein states may provide valuable information for future drug design and therapies.

Energy gap of conformational transition related with temperature for the NACore of α-synuclein

Pathological aggregation of α-synuclein (α-syn) into amyloid fibrils is a major feature of Parkinson's disease (PD). The self-assembly of α-syn is mainly governed by a non-amyloid-β component core (NACore). However, the effects of concentrations and temperatures on their conformational transition remain unclear. To answer this question, we investigated the aggregation kinetics of NACore oligomers in silico by performing several independent all-atom molecular dynamics simulations. The simulation results show that tetramers are more prone to form β-sheets at 300 K than dimers and octamers. We also found that the NACore oligomers had higher β-sheet and β-barrel contents at 310 K. The inter-chain hydrophobic interactions, the backbone hydrogen bonding, the residue-residue interactions between V70-V77 as well as V77-V77 play important roles in the aggregation tendency of NACore octamers at 310 K. Interestingly, the energy gap analysis revealed that the conformational transition of NACore oligomers from intermediate states (β-barrel conformation) to stable structures (β-sheet layers) was dependent on the temperatures. In short, our study provides insight into the kinetic and thermodynamic mechanisms of the conformational transition of NACore at different concentrations and temperatures, contributing to a better understanding of the aggregation process of α-syn in Parkinson's disease.

Molecular Dynamics Simulation Study of the Self-Assembly of Tau-Derived PHF6 and Its Inhibition by Oleuropein Aglycone from Extra Virgin Olive Oil

Alzheimer's disease (AD) and other taupathies are neurodegenerative disorders associated with the amyloid deposition of the Tau protein in the brain. This amyloid formation may be inhibited by small molecules, which is recognized as one of the best therapeutic strategies to stop the progression of the disease. This work focuses on the small nucleating segment, hexapeptide-paired helical filament 6 (PHF6), responsible for Tau aggregation. Using computational modeling and classical molecular dynamics simulations, we show that PHF6 monomers collapse in water to form β-sheet rich structures, and the main olive oil polyphenol oleuropein aglycone (OleA) prevents peptide aggregation significantly. We gradually increase the ratio of the PHF6-OleA from 1:1 to 1:3 and find that for the 1:1 ratio, the peptide monomers are prone to form aggregated structures, while for the 1:2 ratio, the formation of the extended β-sheet structure is significantly less. For a 1:3 ratio of protein/OleA, the peptide residues are sufficiently crowded by OleA molecules through hydrogen bonding, hydrophobic interactions, and π-π stacking; hence, the peptide chains prefer to exist in a monomeric random coil conformation.

Structures and Dynamics of β-Rich Oligomers of ATTR (105-115) Assembly

Transthyretin (TTR) is a tetrameric homologous protein that can dissociate into monomers. Misfolding and aggregation of TTR can lead to amyloid transthyretin amyloidosis (ATTR), which can cause many diseases (e.g., senile systemic amyloidosis, familial amyloid cardiomyopathy, and familial amyloid polyneuropathy). Despite growing evidence indicating that small oligomers play a critical role in regulating cytotoxicity, the structures of these oligomeric intermediates and their conformational transformations are still unclear, impeding our understanding of neurodegenerative mechanisms and the development of therapeutics targeting early aggregation species. The TTR monomer protein consists of various fragments prone to self-aggregation, including the residue 105-115 sequence. Therefore, our study investigated the assembly progress of ATTR (105-115) peptides using all-atom molecular dynamics simulations. The findings indicate that the probability of β-sheet content increases with increasing numbers of peptides. Additionally, interactions between hydrophobic residues L110 and L111 are crucial for the formation of a β-rich oligomer formation. These β-rich oligomers may adopt β-barrel conformations, potentially toxic oligomer species. Free-energy analysis reveals that β-barrel conformations serve as intermediates for these β-rich oligomers. Our insights into the structural ensemble dynamics of ATTR (105-115) contribute to understanding the physical mechanisms underlying the β-barrel oligomers of ATTR. These findings may shed light on the pathological role of ATTR in neurodegenerative diseases and offer potential therapeutic targets.

Exploring the Aggregation Propensity of PHF6 Peptide Segments of the Tau Protein Using Ion Mobility Mass Spectrometry Techniques

Peptide and protein aggregation involves the formation of oligomeric species, but the complex interplay between oligomers of different conformations and sizes complicates their structural elucidation. Using ion mobility mass spectrometry (IM-MS), we aim to reveal these early steps of aggregation for the Ac-PHF6-NH2 peptide segment from tau protein, thereby distinguishing between different oligomeric species and gaining an understanding of the aggregation pathway. An important factor that is often neglected, but which can alter the aggregation propensity of peptides, is the terminal capping groups. Here, we demonstrate the use of IM-MS to probe the early stages of aggregate formation of Ac-PHF6-NH2, Ac-PHF6, PHF6-NH2, and uncapped PHF6 peptide segments. The aggregation propensity of the four PHF6 segments is confirmed using thioflavin T fluorescence assays and transmission electron microscopy. A novel approach based on post-IM fragmentation and quadrupole selection on the TIMS-Qq-ToF (trapped ion mobility) spectrometer was developed to enhance oligomer assignment, especially for the higher-order aggregates. This approach pushes the limits of IM identification of isobaric species, whose signatures appear closer to each other with increasing oligomer size, and provides new insights into the interpretation of IM-MS data. In addition, TIMS collision cross section values are compared with traveling wave ion mobility (TWIMS) data to evaluate potential instrumental bias in the trapped ion mobility results. The two IM-MS instrumental platforms are based on different ion mobility principles and have different configurations, thereby providing us with valuable insight into the preservation of weakly bound biomolecular complexes such as peptide aggregates.

Mapping the configurational landscape and aggregation phase behavior of the tau protein fragment PHF6

The PHF6 (Val-Gln-Ile-Val-Tyr-Lys) motif, found in all isoforms of the microtubule-associated protein tau, forms an integral part of ordered cores of amyloid fibrils formed in tauopathies and is thought to play a fundamental role in tau aggregation. Because PHF6 as an isolated hexapeptide assembles into ordered fibrils on its own, it is investigated as a minimal model for insight into the initial stages of aggregation of larger tau fragments. Even for this small peptide, however, the large length and time scales associated with fibrillization pose challenges for simulation studies of its dynamic assembly, equilibrium configurational landscape, and phase behavior. Here, we develop an accurate, bottom-up coarse-grained model of PHF6 for large-scale simulations of its aggregation, which we use to uncover molecular interactions and thermodynamic driving forces governing its assembly. The model, not trained on any explicit information about fibrillar structure, predicts coexistence of formed fibrils with monomers in solution, and we calculate a putative equilibrium phase diagram in concentration-temperature space. We also characterize the configurational and free energetic landscape of PHF6 oligomers. Importantly, we demonstrate with a model of heparin that this widely studied cofactor enhances the aggregation propensity of PHF6 by ordering monomers during nucleation and remaining associated with growing fibrils, consistent with experimentally characterized heparin-tau interactions. Overall, this effort provides detailed molecular insight into PHF6 aggregation thermodynamics and pathways and, furthermore, demonstrates the potential of modern multiscale modeling techniques to produce predictive models of amyloidogenic peptides simultaneously capturing sequence-specific effects and emergent aggregate structures.

Identification of Aggregation Mechanism of Acetylated PHF6* and PHF6 Tau Peptides Based on Molecular Dynamics Simulations and Markov State Modeling

The microtubule-associated protein tau (MAPT) has a critical role in the development and preservation of the nervous system. However, tau's dysfunction and accumulation in the human brain can lead to several neurodegenerative diseases, such as Alzheimer's disease, Down's syndrome, and frontotemporal dementia. The microtubule binding (MTB) domain plays a significant, important role in determining the tau's pathophysiology, as the core of paired helical filaments PHF6* (275VQIINK280) and PHF6 (306VQIVYK311) of R2 and R3 repeat units, respectively, are formed in this region, which promotes tau aggregation. Post-translational modifications, and in particular lysine acetylation at K280 of PHF6* and K311 of PHF6, have been previously established to promote tau misfolding and aggregation. However, the exact aggregation mechanism is not known. In this study, we established an atomic-level nucleation-extension mechanism of the separated aggregation of acetylated PHF6* and PHF6 hexapeptides, respectively, of tau. We show that the acetylation of the lysine residues promotes the formation of β-sheet enriched high-ordered oligomers. The Markov state model analysis of ac-PHF6* and ac-PHF6 aggregation revealed the formation of an antiparallel dimer nucleus which could be extended from both sides in a parallel manner to form mixed-oriented and high-ordered oligomers. Our study describes the detailed mechanism for acetylation-driven tau aggregation, which provides valuable insights into the effect of post-translation modification in altering the pathophysiology of tau hexapeptides.

Dose-dependent binding behavior of anthraquinone derivative purpurin interacting with tau-derived peptide protofibril

Alzheimer's disease is hallmarked by microtubule-associated protein tau tangles and amyloid-β plaques. The β-structure propensity of tau inclusions is closely related to the hexapeptide motif VQIVYK (termed PHF6), and disruption of this motif prevents tau aggregation. Small-molecule inhibitors are considered a promising therapeutic strategy, but the molecular mechanisms underlying the correlation between dose and inhibitory effects are still unclear. In this work, we investigated the dose-induced influence of purpurin, an anthraquinone derivative, on the structural stability of the PHF6 fibrillar nucleus by performing microsecond all-atom molecular dynamics simulations in explicit water. The stability of PHF6 protofibrils of different sizes was first examined, and it was found that the structural stability of fibrillar oligomers increases with oligomer size, and that the octamer is the minimal stable nucleus for fibril formation. When purpurin molecules were added to the protofibril octamer at a low purpurin/peptide ratio, they bound to the octamer with different coupling states, and the different states may transition to each of the other states through an uncoupling state or directly through a short-time transition. With increasing purpurin/peptide ratio, purpurins tend to self-aggregate rather than bind to the protein surface. Interestingly, the contacts between individual purpurins and the octamer as a function of the purpurin number show a power-law behavior, which may serve as a useful indicator to reflect the binding efficiency of ligands to proteins in drug screening. The interaction analysis reveals that purpurin prefers to bind to the hydrophilic and aromatic Tyr and has the lowest probability with the hydrophobic Val located in the middle of PHF6. Aromatic stacking plays a key role in the octamer-purpurin interaction, in which the three aromatic rings of purpurin have different contributions. In addition, purpurin shows a remarkable disruptive effect on the protofibril octamer when the molar ratio of purpurin to peptide is 1 : 2; above this ratio, the binding mode and disruption effect of purpurin do not change significantly. Our work provides a detailed picture of the dynamics and interactions of purpurin binding to the PHF6 protofibril and expands the understanding of the dose-induced inhibitory mechanism.

Multistep molecular mechanisms of Aβ16-22 fibril formation revealed by lattice Monte Carlo simulations

As a model of self-assembly from disordered monomers to fibrils, the amyloid-β fragment Aβ16-22 was subject to past numerous experimental and computational studies. Because dynamics information between milliseconds and seconds cannot be assessed by both studies, we lack a full understanding of its oligomerization. Lattice simulations are particularly well suited to capture pathways to fibrils. In this study, we explored the aggregation of 10 Aβ16-22 peptides using 65 lattice Monte Carlo simulations, each simulation consisting of 3 × 109 steps. Based on a total of 24 and 41 simulations that converge and do not converge to the fibril state, respectively, we are able to reveal the diversity of the pathways leading to fibril structure and the conformational traps slowing down the fibril formation.

Aggregation of Amyloidogenic Peptide Uperin-Molecular Dynamics Simulations

Uperin 3.5 is a remarkable natural peptide obtained from the skin of toadlets comprised of 17 amino acids which exhibits both antimicrobial and amyloidogenic properties. Molecular dynamics simulations were performed to study the β-aggregation process of uperin 3.5 as well as two of its mutants, in which the positively charged residues Arg7 and Lys8 have been replaced by alanine. All three peptides rapidly underwent spontaneous aggregation and conformational transition from random coils to beta-rich structures. The simulations reveal that the initial and essential step of the aggregation process involves peptide dimerization and the formation of small beta-sheets. A decrease in positive charge and an increase in the number of hydrophobic residues in the mutant peptides lead to an increase in the rate of their aggregation.

Dissecting how ALS-associated D290V mutation enhances pathogenic aggregation of hnRNPA2286-291 peptides: Dynamics and conformational ensembles

The aggregation of RNA binding proteins, including hnRNPA1/2, TDP-43 and FUS, is heavily implicated in causing or increasing disease risk for a series of neurodegenerative diseases such as amyotrophic lateral sclerosis (ALS). A recent experimental study demonstrated that an ALS-related D290V mutation in the low complexity domain (LCD) of hnRNPA2 can enhance the aggregation propensity of wild type (WT) hnRNPA2_286-291 peptide. However, the underlying molecular mechanisms remain elusive. Herein, we investigated effects of D290V mutation on aggregation dynamics of hnRNPA2_286-291 peptide and the conformational ensemble of hnRNPA2_286-291 oligomers by performing all-atom molecular dynamic and replica-exchange molecular dynamic simulations. Our simulations demonstrate that D290V mutation greatly reduces the dynamics of hnRNPA2_286-291 peptide and that D290V oligomers possess higher compactness and β-sheet content than WT, indicative of mutation-enhanced aggregation capability. Specifically, D290V mutation strengthens inter-peptide hydrophobic, main-chain hydrogen bonding and side-chain aromatic stacking interactions. Those interactions collectively lead to the enhancement of aggregation capability of hnRNPA2_286-291 peptides. Overall, our study provides insights into the dynamics and thermodynamic mechanisms underlying D290V-induced disease-causing aggregation of hnRNPA2_286-291, which could contribute to better understanding of the transitions from reversible condensates to irreversible pathogenic aggregates of hnRNPA2 LCD in ALS-related diseases.

Investigation of the Structure of Full-Length Tau Proteins with Coarse-Grained and All-Atom Molecular Dynamics Simulations

Tau proteins not only have many important biological functions but also are associated with several neurodegenerative diseases, such as Parkinson's disease and Alzheimer's disease (AD). However, it is still a challenge to identify the atomic structure of full-length tau proteins due to their lengthy and disordered characteristics and the factor that there are no crystal structures of full-length tau proteins available. We performed multi- and large-scale molecular dynamics simulations of the full-length tau monomer (the 2N4R isoform and 441 residues) in aqueous solution under biological conditions with coarse-grained and all-atom force fields. The obtained atomic structures produced radii of gyration and chemical shifts that are in excellent agreement with those of experiment. The generated monomer structure ensemble would be very useful for further studying the oligomerization mechanism and discovering tau oligomerization inhibitors, which are important events in AD drug development.

Atomistic Insights into the Inhibitory Mechanism of Tyrosine Phosphorylation against the Aggregation of Human Tau Fragment PHF6

Abnormal aggregation of the microtubule-associated protein tau into intracellular fibrillary inclusions is characterized as the hallmark of tauopathies, including Alzheimer's disease and chronic traumatic encephalopathy. The hexapeptide ³⁰⁶VQIVYK³¹¹ (PHF6) of R3 plays an important role in the aggregation of tau. Recent experimental studies reported that phosphorylation of residue tyrosine 310 (Y310) could decrease the propensity of PHF6 to form fibrils and inhibit tau aggregation. However, the underlying inhibitory mechanism is not well understood. In this work, we systematically investigated the influences of phosphorylation on the conformational ensembles and oligomerization dynamics of PHF6 by performing extensive all-atom molecular dynamics (MD) simulations. Our replica exchange MD simulations demonstrate that Y310 phosphorylation could effectively suppress the formation of β-structure and shift PHF6 oligomers toward coil-rich aggregates. The interaction analyses show that hydrogen bonding and hydrophobic interactions among PHF6 peptides, as well as Y310-Y310 π-π stacking and I308-Y310 CH-π interactions, are weakened by phosphorylation. Additional microsecond MD simulations show that Y310 phosphorylation could inhibit the oligomerization of PHF6 by preventing the formation of large β-sheet oligomers and multi-layer β-sheet aggregates. This study provides mechanistic insights into the phosphorylation-inhibited tau aggregation, which may be helpful for the in-depth understanding of the pathogenesis of tauopathies.

Design of short peptides and peptide amphiphiles as collagen mimics and an investigation of their interactions with collagen using molecular dynamics simulations and docking studies

Short peptide sequences and bolaamphiphiles derived from natural proteins are gaining importance due to their ability to form unique nanoscale architectures for a variety of biological applications. In this work, we have designed six short peptides (triplet or monomeric forms) and two peptide bolaamphiphiles that either incorporate the bioactive collagen motif (Gly-X-Y) or sequences where Gly, Pro, or hydroxyproline (Hyp) are replaced by Ala or His. For the bolaamphiphiles, a malate moiety was used as the aliphatic linker for connecting His with Hyp to create collagen mimics. Stability of the assemblies was assessed through molecular dynamics simulations and results indicated that (Pro-Ala-His)₃ and (Ala-His-Hyp)₃ formed the most stable structures, while the amphiphiles and the monomers showed some disintegration over the course of the 200 ns simulation, though most regained structural integrity and formed fibrillar structures, and micelles by the end of the simulation, likely due to the formation of more thermodynamically stable conformations. Multiple replica simulations (REMD) were also conducted where the sequences were simulated at different temperatures. Our results showed excellent convergence in most cases compared to constant temperature molecular dynamics simulation. Furthermore, molecular docking and MD simulations of the sequences bound to collagen triple helix structure revealed that several of the sequences had a high binding affinity and formed stable complexes, particularly (Pro-Ala-His)₃ and (Ala-His-Hyp)₃. Thus, we have designed new hybrid-peptide-based sequences which may be developed for potential applications as biomaterials for tissue engineering or drug delivery.

Deciphering the Role of ATP on PHF6 Aggregation

The aggregation of Tau protein, which are involved in Alzheimer's disease, are associated with the self-assembly of the hexapeptide sequence, paired helical filament 6 (PHF6) from repeat 3 of Tau. In order to treat Alzheimer's disease and other such tauopathies, one of the therapeutic strategies is to inhibit aggregation of Tau and its nucleating segments. Therefore, we have studied the effect of adenosine triphosphate (ATP) on the aggregation of PHF6. ATP has, interestingly, demonstrated its ability to inhibit and dissolve protein aggregates. Using classical molecular dynamics simulations, we observed that the hydrophobic core of PHF6 segment displays extended β-sheet conformation, which stabilizes PHF6 aggregates. However, the distribution of ATP around the vicinity of the peptides enables PHF6 to remain discrete and attain random coil conformers. The interpeptide interactions are substituted by PHF6-ATP interactions through hydrogen bonding and hydrophobic interactions (including π-π stacking). Furthermore, the adenosine moiety of ATP contributes more than the triphosphate chain toward PHF6-ATP interaction. Ultimately, this work establishes the inhibitory activity of ATP against Tau aggregation; hence, the therapeutic effect of ATP should be explored further in regard to the effective treatment of Alzheimer's disease.

Development of a Hybrid-Resolution Force Field for Peptide Self-Assembly Simulations: Optimizing Peptide-Peptide and Peptide-Solvent Interactions

Atomic descriptions of peptide self-assembly are crucial to an understanding of disease-related peptide aggregation and the design of peptide-assembled materials. Obtaining these descriptions through computer simulation is challenging because current force fields, which were not designed for this process and are often unable to describe correctly peptide self-assembly behavior and the sequence dependence. Here, we developed a framework using dipeptide aggregation as a model system to improve force fields for simulations of self-assembly. Aggregation-related structural properties were designed and used to guide the optimization of peptide-peptide and peptide-solvent interactions. With this framework, we developed a self-assembly force field, termed PACE-ASM, by reoptimizing a hybrid-resolution force field that was originally developed for folding simulation. With its applicability in folding simulations, the new PACE was used to simulate the self-assembly of two disease-related short peptides, Aβ_16-21 and PHF6, into β-sheet-rich cross-β amyloids. These simulations reproduced the crystal structures of Aβ_16-21 and PHF6 amyloids at near-atomic resolution and captured the difference in packing orientations between the two sequences, a task which is challenging even with all-atom force fields. Apart from cross-β amyloids, the self-assembly of emerging helix-rich cross-α amyloids by another peptide PSMα3 can also be correctly described with the new PACE, manifesting the versatility of the force field. We demonstrated that the ability of the PACE-ASM to model peptide self-assembly is based largely on its improved description of peptide-peptide and peptide-solvent interactions. This was achieved with our optimization framework that can readily identify and address the deficiency in describing these interactions.

Tau aggregation

Here we revisit tau protein aggregation at primary, secondary, tertiary and quaternary structures. In addition, the presence of non-aggregated tau protein, which has been recently discovered, is also commented on.

Computer Simulations Aimed at Exploring Protein Aggregation and Dissociation

Protein aggregation can lead to well-defined structures that are functional, but is also the cause of the death of neuron cells in many neurodegenerative diseases. The complexity of the molecular events involved in the aggregation kinetics of amyloid proteins and the transient and heterogeneous characters of all oligomers prevent high-resolution structural experiments. As a result, computer simulations have been used to determine the atomic structures of amyloid proteins at different association stages as well as to understand fibril dissociation. In this chapter, we first review the current computer simulation methods used for aggregation with some atomistic and coarse-grained results aimed at better characterizing the early formed oligomers and amyloid fibril formation. Then we present the applications of non-equilibrium molecular dynamics simulations to comprehend the dissociation of protein assemblies.

Effects of All-Atom Molecular Mechanics Force Fields on Amyloid Peptide Assembly: The Case of PHF6 Peptide of Tau Protein

Molecular dynamics (MD) simulations play a vital role in revealing the mechanism of amyloid aggregation that is crucial to the therapeutic agent development for Alzheimer's Disease. However, the accuracy of MD simulation results strongly depends on the force field employed. In our previous benchmark for 17 all-atom force fields on modeling of amyloid aggregation using the Aβ_16-22 dimer, we showed that AMBER14SB and CHARMM36m are suitable force fields for amyloid aggregation simulation, while GROMOS54a7 and OPLSAA are not good for the task. In this work, we continue assessing the applicability of atomistic force fields on amyloid aggregation using the VQIVYK (PHF6) peptide which is essential for tau-protein aggregation. Although, both Aβ_16-22 and PHF6 peptides formed fibrils in vitro, the PHF6 fibrils are parallel β-sheets, while the Aβ_16-22 fibrils are antiparallel β-sheets. We performed an all-atom large-scale MD simulation in explicit water on the PHF6 dimer and octa-peptides systems using five mainstream force fields, including AMBER99SB-disp, AMBER14SB, CHARMM36m, GROMOS54a7, and OPLSAA. The accumulated simulation time is 0.2 ms. Our result showed that the β-sheet structures of PHF6 peptides sampled by AMBER99SB-disp, AMBER14SB, GROMOS54a7, and OPLSAA are in favor of the antiparallel β-sheets, while the dominant type of β-sheet structures is parallel β-sheet by using CHARMM36m. Among the five force fields, CHARMM36m provides the strongest CH-π interaction that was observed in an NMR study. The comparison between our results and experimental observation indicates that CHARMM36m achieved the best performance on modeling the aggregation of PHF6 peptides. In summary, CHARMM36m is currently the most suitable force field for studying the aggregation of both amyloid-β and Tau through MD simulations.

Oxidative Stress Conditions Result in Trapping of PHF-Core Tau (297-391) Intermediates

The self-assembly of tau into paired helical filaments (PHFs) in neurofibrillary tangles (NFTs) is a significant event in Alzheimer's disease (AD) pathogenesis. Numerous post-translational modifications enhance or inhibit tau assembly into NFTs. Oxidative stress, which accompanies AD, induces multiple post-translational modifications in proteins, including the formation of dityrosine (DiY) cross-links. Previous studies have revealed that metal-catalysed oxidation (MCO) using Cu²⁺ and H₂O₂ leads to the formation of DiY cross-links in two misfolding proteins, Aβ and α-synuclein, associated with AD and Parkinson's disease respectively. The effect of MCO on tau remains unknown. Here, we examined the effect of MCO and ultra-violet oxidation to study the influence of DiY cross-linking on the self-assembly of the PHF-core tau fragment. We report that DiY cross-linking facilitates tau assembly into tau oligomers that fail to bind thioflavin S, lack β-sheet structure and prevents their elongation into filaments. At a higher concentration, Cu²⁺ (without H₂O₂) also facilitates the formation of these tau oligomers. The DiY cross-linked tau oligomers do not cause cell death. Our findings suggest that DiY cross-linking of pre-assembled tau promotes the formation of soluble tau oligomers that show no acute impact on cell viability.

Amyloid Oligomers: A Joint Experimental/Computational Perspective on Alzheimer's Disease, Parkinson's Disease, Type II Diabetes, and Amyotrophic Lateral Sclerosis

Protein misfolding and aggregation is observed in many amyloidogenic diseases affecting either the central nervous system or a variety of peripheral tissues. Structural and dynamic characterization of all species along the pathways from monomers to fibrils is challenging by experimental and computational means because they involve intrinsically disordered proteins in most diseases. Yet understanding how amyloid species become toxic is the challenge in developing a treatment for these diseases. Here we review what computer, in vitro, in vivo, and pharmacological experiments tell us about the accumulation and deposition of the oligomers of the (Aβ, tau), α-synuclein, IAPP, and superoxide dismutase 1 proteins, which have been the mainstream concept underlying Alzheimer's disease (AD), Parkinson's disease (PD), type II diabetes (T2D), and amyotrophic lateral sclerosis (ALS) research, respectively, for many years.

A theoretical study of polymorphism in VQIVYK fibrils

The VQIVYK fragment from the Tau protein, also known as PHF6, is essential for aggregation of Tau into neurofibrillary lesions associated with neurodegenerative diseases. VQIVYK itself forms amyloid fibrils composed of paired β-sheets. Therefore, the full Tau protein and VQIVYK fibrils have been intensively investigated. A central issue in these studies is polymorphism, the ability of a protein to fold into more than one structure. Using all-atom molecular simulations, we generate five stable polymorphs of VQIVYK fibrils, establish their relative free energy with umbrella sampling methods, and identify the side chain interactions that provide stability. The two most stable polymorphs, which have nearly equal free energy, are formed by interdigitation of the mostly hydrophobic VIY "face" sides of the β-sheets. Another stable polymorph is formed by interdigitation of the QVK "back" sides. When we turn to examine structures from cryo-electron microscopy experiments on Tau filaments taken from diseased patients or generated in vitro, we find that the pattern of side chain interactions found in the two most stable face-to-face as well as the back-to-back polymorphs are recapitulated in amyloid structures of the full protein. Thus, our studies suggest that the interactions stabilizing PHF6 fibrils explain the amyloidogenicity of the VQIVYK motif within the full Tau protein and provide justification for the use of VQIVYK fibrils as a test bed for the design of molecules that identify or inhibit amyloid structures.

From beta amyloid to altered proteostasis in Alzheimer's disease

Alzheimer's disease (AD) is an age related neurodegenerative disorder causing severe disability and important socio-economic burden, but with no cure available to date. To disentangle this puzzling disease genetic studies represented an important way for the comprehension of pathogenic mechanisms. Abnormal processing and accumulation of amyloid-β peptide (Aβ) has been considered the main cause and trigger factor of the disease. The amyloid cascade theory has fallen into crisis because the failure of several anti-amyloid drugs trials and because of the simple equation AD = abnormal Aβ deposition is not always the case. We now know that multiple neurodegenerative diseases share common pathogenic mechanisms leading to accumulation of misfolded protein species. Genome Wide Association studies (GWAS) led to the identification of large numbers of DNA common variants (SNPs) distributed on different chromosomes and modulating the Alzheimer's risk. GWAS genes fall into several common pathways such as immune system and neuroinflammation, lipid metabolism, synaptic dysfunction and endocytosis, all of them addressing to novel routes for different pathogenic mechanisms. Other hints could be derived from epidemiological and experimental studies showing some lifestyles may have a major role in the pathogenesis of many age-associated diseases by modifying cell metabolism, proteostasis and microglia mediated neuroinflammation.

Terminal Capping of an Amyloidogenic Tau Fragment Modulates Its Fibrillation Propensity

Aberrant protein folding leading to the formation of characteristic cross-β-sheet-rich amyloid structures is well known for its association with a variety of debilitating human diseases. Often, depending upon amino acid composition, only a small segment of a large protein participates in amyloid formation and is in fact capable of self-assembling into amyloid, independent of the rest of the protein. Therefore, such peptide fragments serve as useful model systems for understanding the process of amyloid formation. An important factor that has often been overlooked while using peptides to mimic full-length protein is the charge on the termini of these peptides. Here, we show the influence of terminal charges on the aggregation of an amyloidogenic peptide from microtubule-associated protein Tau, implicated in Alzheimer's disease and tauopathies. We found that modification of terminal charges by capping the peptide at one or both of the termini drastically modulates the fibrillation of the hexapeptide sequence paired helical filament 6 (PHF6) from repeat 3 of Tau, both with and without heparin. Without heparin, the PHF6 peptide capped at both termini and PHF6 capped only at the N-terminus self-assembled to form amyloid fibrils. With heparin, all capping variants of PHF6, except for PHF6 with both termini free, formed typical amyloid fibrils. However, the rate and extent of aggregation both with and without heparin as well as the morphology of aggregates were found to be highly dependent on the terminal charges. Our molecular dynamics simulations on PHF6 capping variants corroborated our experiments and provided critical insights into the mechanism of PHF6 self-assembly. Overall, our results emphasize the importance of terminal modifications in fibrillation of small peptide fragments and provide significant insights into the aggregation of a small Tau fragment, which is considered essential for Tau filament assembly.

A two-step biopolymer nucleation model shows a nonequilibrium critical point

Biopolymer self-assembly pathways are complicated by the ability of their monomeric subunits to adopt different conformational states. This means nucleation often involves a two-step mechanism where the monomers first condense to form a metastable intermediate, which then converts to a stable polymer by conformational rearrangement of constituent monomers. Nucleation intermediates play a causative role in amyloid diseases such as Alzheimer's and Parkinson's. While existing mathematical models neglect the conversion dynamics, experiments show that conversion events frequently occur on comparable timescales to the condensation of intermediates and growth of mature polymers and thus cannot be ignored. We present a model that explicitly accounts for simultaneous assembly and conversion. To describe conversion, we propose an experimentally motivated initiation-propagation mechanism in which the stable phase arises locally within the intermediate and then spreads by nearest-neighbor interactions, in a manner analogous to one-dimensional Glauber dynamics. Our analysis shows that the competing timescales of assembly and conversion result in a nonequilibrium critical point, separating a regime where intermediates are kinetically unstable from one where conformationally mixed intermediates accumulate. This strongly affects the accumulation rate of the stable biopolymer phase. Our model is uniquely able to explain experimental phenomena such as the formation of mixed intermediates and abrupt changes in the scaling exponent γ, which relates the total monomer concentration to the accumulation rate of the stable phase. This provides a first step toward a general model of two-step biopolymer nucleation, which can quantitatively predict the concentration and composition of biologically crucial intermediates.

Purpurin modulates Tau-derived VQIVYK fibrillization and ameliorates Alzheimer's disease-like symptoms in animal model

Neurofibrillary tangles of the Tau protein and plaques of the amyloid β peptide are hallmarks of Alzheimer's disease (AD), which is characterized by the conversion of monomeric proteins/peptides into misfolded β-sheet rich fibrils. Halting the fibrillation process and disrupting the existing aggregates are key challenges for AD drug development. Previously, we performed in vitro high-throughput screening for the identification of potent inhibitors of Tau aggregation using a proxy model, a highly aggregation-prone hexapeptide fragment 306VQIVYK311 (termed PHF6) derived from Tau. Here we have characterized a hit molecule from that screen as a modulator of Tau aggregation using in vitro, in silico, and in vivo techniques. This molecule, an anthraquinone derivative named Purpurin, inhibited ~ 50% of PHF6 fibrillization in vitro at equimolar concentration and disassembled pre-formed PHF6 fibrils. In silico studies showed that Purpurin interacted with key residues of PHF6, which are responsible for maintaining its β-sheets conformation. Isothermal titration calorimetry and surface plasmon resonance experiments with PHF6 and full-length Tau (FL-Tau), respectively, indicated that Purpurin interacted with PHF6 predominantly via hydrophobic contacts and displayed a dose-dependent complexation with FL-Tau. Purpurin was non-toxic when fed to Drosophila and it significantly ameliorated the AD-related neurotoxic symptoms of transgenic flies expressing WT-FL human Tau (hTau) plausibly by inhibiting Tau accumulation and reducing Tau phosphorylation. Purpurin also reduced hTau accumulation in cell culture overexpressing hTau. Importantly, Purpurin efficiently crossed an in vitro human blood-brain barrier model. Our findings suggest that Purpurin could be a potential lead molecule for AD therapeutics.

A computational approach to get insights into multiple faces of additives in modulation of protein aggregation pathways

An enormous population worldwide is presently confronted with debilitating neurodegenerative diseases. The etiology of the disease is connected to protein aggregation and the events involved therein. Thus, a complete understanding of an inhibitor at different stages in the process is imperative for the formulation of a drug molecule. This review presents a detailed summary of the current status of different cosolvents. It further develops how the complex aggregation pathway can be simplified into three steps common to all proteins and the way computer simulations can be exploited to gain insights into the ways by which known inhibitors can affect all these stages. Computation of theoretical parameters in this regard and their correlation with experimental techniques is accentuated. In addition to providing an outline of the scope of different additives, this review showcases the way by which the problem of analyzing an effect of an additive can be addressed effectively via MD simulations.

Probing the Basis of α-Synuclein Aggregation by Comparing Simulations to Single-Molecule Experiments

Intrinsically disordered proteins often play an important role in protein aggregation. However, it is challenging to determine the structures and interactions that drive the early stages of aggregation because they are transient and obscured in a heterogeneous mixture of disordered states. Even computational methods are limited because the lack of ordered structure makes it difficult to ensure that the relevant conformations are sampled. We address these challenges by integrating atomistic simulations with high-resolution single-molecule measurements reported previously, using the measurements to help discern which parts of the disordered ensemble of structures in the simulations are most probable while using the simulations to identify residues and interactions that are important for oligomer stability. This approach was applied to α-synuclein, an intrinsically disordered protein that aggregates in the context of Parkinson's disease. We simulated single-molecule pulling experiments on dimers, the minimal oligomer, and compared them to force spectroscopy measurements. Force-extension curves were simulated starting from a set of 66 structures with substantial structured content selected from the ensemble of dimer structures generated at zero force via Monte Carlo simulations. The pattern of contour length changes as the structures unfolded through intermediate states was compared to the results from optical trapping measurements on the same dimer to discern likely structures occurring in the measurements. Simulated pulling curves were generally consistent with experimental data but with a larger number of transient intermediates. We identified an ensemble of β-rich dimer structures consistent with the experimental data from which dimer interfaces could be deduced. These results suggest specific druggable targets in the structural motifs of α-synuclein that may help prevent the earliest steps of oligomerization.

Amyloid Oligomers, Protofibrils and Fibrils

Amyloid diseases are of major concern all over the world due to a number of factors including: (i) aging population, (ii) increasing life span and (iii) lack of effective pharmacotherapy options. The past decade has seen intense research in discovering disease-modifying multi-targeting small molecules as therapeutic options. In recent years, targeting the amyloid cascade has emerged as an attractive strategy to discover novel neurotherapeutics. Formation of amyloid species, with different degrees of solubility and neurotoxicity is associated with the gradual decline in cognition leading to dementia/cell dysfunction. Here, in this chapter, we have described the recent scenario of amyloid diseases with a great deal of information about the structural features of oligomers, protofibrils and fibrils. Also, comprehensive details have been provided to differentiate the degree of toxicity associated with prefibrillar aggregates. Moreover, a review of the technologies that aid characterisation of oligomer, protofibrils and fibrils as well as various inhibition strategies to overcome protein fibrillation are also discussed.

Secondary nucleation in amyloid formation

Nucleation of new peptide and protein aggregates on the surfaces of amyloid fibrils of the same peptide or protein has emerged in the past two decades as a major pathway for both the generation of molecular species responsible for cellular toxicity and for the autocatalytic proliferation of peptide and protein aggregates. A key question in current research is the molecular mechanism and driving forces governing such processes, known as secondary nucleation. In this context, the analogies with other self-assembling systems for which monomer-dependent secondary nucleation has been studied for more than a century provide a valuable source of inspiration. Here, we present a short overview of this background and then review recent results regarding secondary nucleation of amyloid-forming peptides and proteins, focusing in particular on the amyloid β peptide (Aβ) from Alzheimer's disease, with some examples regarding α-synuclein from Parkinson's disease. Monomer-dependent secondary nucleation of Aβ was discovered using a combination of kinetic experiments, global analysis, seeding experiments and selective isotope-enrichment, which pinpoint the monomer as the origin of new aggregates in a fibril-catalyzed reaction. Insights into driving forces are gained from variations of solution conditions, temperature and peptide sequence. Selective inhibition of secondary nucleation is explored as an effective means to limit oligomer production and toxicity. We also review experiments aimed at finding interaction partners of oligomers generated by secondary nucleation in an ongoing aggregation process. At the end of this feature article we bring forward outstanding questions and testable mechanistic hypotheses regarding monomer-dependent secondary nucleation in amyloid formation.

Carbamylation promotes amyloidogenesis and induces structural changes in Tau-core hexapeptide fibrils

BACKGROUND

Carbamylation is a non-enzymatic post-translational modification (PTM), which involves the covalent modification of N-terminus of protein or ε-amino group of Lys. The role of carbamylation in several age-related disorders is well documented, however, the relationship between carbamylation and neurodegenerative disorders including Alzheimer's disease remains uncharted.

METHODS

In the present study, using aggregation-prone tau-core hexapeptide fragments ³⁰⁶VQIVYK³¹¹ (PHF6) and ²⁷⁵VQIINK²⁸⁰ (PHF6*) as models, we have elucidated the effect of carbamylation on aggregation kinetics and the changes occurring in the 3-dimensional architecture of fibrils using biophysical assays and molecular dynamics simulations.

RESULTS

We found that carbamylation aids in amyloid formation and can convert the unstructured off-pathway aggregates into robust amyloids, which were toxic to cells. Electron microscopy images and molecular dynamics simulations of PHF6 fibrils showed that carbamylated peptides can form excess hydrogen bonds and modulate the pitch length and twist of peptides fibrils. We have also compared N-terminal carbamylation to acetylation and further extended our finding to full length tau that exhibits aggregation upon carbamylation even in the absence of any external inducer.

CONCLUSION

Our in vitro and in silico results together suggest that carbamylation can modulate the aggregation pathway of the amyloidegenic sequences and cause structural changes in fibril assemblies.

GENERAL SIGNIFICANCE

Carbamylation acts as a switch, which triggers the aggregation in short amyloidogenic peptide fragments and modulate the structural changes in resulting amyloid fibrils.

Integrating in vitro and in silico approaches to evaluate the "dual functionality" of palmatine chloride in inhibiting and disassembling Tau-derived VQIVYK peptide fibrils

BACKGROUND

Alzheimer's disease (AD) is the most common neurodegenerative disorder which is characterized by the deposits of intra-cellular tau protein and extra-cellular amyloid-β (Aβ) peptides in the human brain. Understanding the mechanism of protein aggregation and finding compounds that are capable of inhibiting its aggregation is considered to be highly important for disease therapy.

METHODS

We used an in vitro High-Throughput Screening for the identification of potent inhibitors of tau aggregation using a proxy model; a highly aggregation-prone hexapeptide fragment ³⁰⁶VQIVYK³¹¹ derived from tau. Using ThS fluorescence assay we screened a library of 2401 FDA approved, bio-active and natural compounds in attempt to find molecules which can efficiently modulate tau aggregation.

RESULTS

Among the screened compounds, palmatine chloride (PC) alkaloid was able to dramatically reduce the aggregation propensity of PHF6 at sub-molar concentrations. PC was also able to disassemble preformed aggregates of PHF6 and reduce the amyloid content in a dose-dependent manner. Insights obtained from MD simulation showed that PC interacted with the key residues of PHF6 responsible for β-sheet formation, which could likely be the mechanism of inhibition and disassembly. Furthermore, PC could effectively inhibit the aggregation of full-length tau and disassemble preformed aggregates.

CONCLUSIONS

We found that PC possesses "dual functionality" towards PHF6 and full-length tau, i.e. inhibit their aggregation and disassemble pre-formed fibrils.

GENERAL SIGNIFICANCE

The "dual functionality" of PC is valuable as a disease modifying strategy for AD, and other tauopathies, by inhibiting their progress and reducing the effect of fibrils already present in the brain.

Inhibition of GNNQQNY prion peptide aggregation by trehalose: a mechanistic view

Deposition of amyloid fibrils is the seminal event in the pathogenesis of numerous neurodegenerative diseases. The formation of this amyloid assembly is the manifestation of a cascade of structural transitions including toxic oligomer formation in the early stages of aggregation. Thus a viable therapeutic strategy involves the use of small molecular ligands to interfere with this assembly. In this perspective, we have explored the kinetics of aggregate formation of the fibril forming GNNQQNY peptide fragment from the yeast prion protein SUP35 using multiple all atom MD simulations with explicit solvent and provided mechanistic insights into the way trehalose, an experimentally known aggregation inhibitor, modulates the aggregation pathway. The results suggest that the assimilation process is impeded by different barriers at smaller and larger oligomeric sizes: the initial one being easily surpassed at higher temperatures and peptide concentrations. The kinetic profile demonstrates that trehalose delays the aggregation process by increasing both these activation barriers, specifically the latter one. It increases the sampling of small-sized aggregates that lack the beta sheet conformation. Analysis reveals that the barrier in the growth of larger stable oligomers causes the formation of multiple stable small oligomers which then fuse together bimolecularly. The PCA of 26 properties was carried out to deconvolute the events within the temporary lag phases, which suggested dynamism in lags involving an increase in interchain contacts and burial of SASA. The predominant growth route is monomer addition, which changes to condensation on account of a large number of depolymerisation events in the presence of trehalose. The favourable interaction of trehalose specifically with the sidechain of the peptide promotes crowding of trehalose molecules in its vicinity - the combination of both these factors imparts the observed behaviour. Furthermore, increasing trehalose concentration leads to faster expulsion of water molecules than interpeptide interactions. These expelled water molecules have larger translational movement, suggesting an entropy factor to favor the assembly process. Different conformations observed under this condition suggest the role of water molecules in guiding the morphology of the aggregates as well. A similar scenario exists on increasing peptide concentration.

Lattice model for amyloid peptides: OPEP force field parametrization and applications to the nucleus size of Alzheimer's peptides

Coarse-grained protein lattice models approximate atomistic details and keep the essential interactions. They are, therefore, suitable for capturing generic features of protein folding and amyloid formation at low computational cost. As our aim is to study the critical nucleus sizes of two experimentally well-characterized peptide fragments Aβ16-22 and Aβ37-42 of the full length Aβ1-42 Alzheimer's peptide, it is important that simulations with the lattice model reproduce all-atom simulations. In this study, we present a comprehensive force field parameterization based on the OPEP (Optimized Potential for Efficient protein structure Prediction) force field for an on-lattice protein model, which incorporates explicitly the formation of hydrogen bonds and directions of side-chains. Our bottom-up approach starts with the determination of the best lattice force parameters for the Aβ16-22 dimer by fitting its equilibrium parallel and anti-parallel β-sheet populations to all-atom simulation results. Surprisingly, the calibrated force field is transferable to the trimer of Aβ16-22 and the dimer and trimer of Aβ37-42. Encouraged by this finding, we characterized the free energy landscapes of the two decamers. The dominant structure of the Aβ16-22 decamer matches the microcrystal structure. Pushing the simulations for aggregates between 4-mer and 12-mer suggests a nucleus size for fibril formation of 10 chains. In contrast, the Aβ37-42 decamer is largely disordered with mixed by parallel and antiparallel chains, suggesting that the nucleus size is >10 peptides. Our refined force field coupled to this on-lattice model should provide useful insights into the critical nucleation number associated with neurodegenerative diseases.

Characterizing Structural Stability of Amyloid Motif Fibrils Mediated by Water Molecules

In biological systems, structural confinements of amyloid fibrils can be mediated by the role of water molecules. However, the underlying effect of the dynamic behavior of water molecules on structural stabilities of amyloid fibrils is still unclear. By performing molecular dynamics simulations, we investigate the dynamic features and the effect of interior water molecules on conformations and mechanical characteristics of various amyloid fibrils. We find that a specific mechanism induced by the dynamic properties of interior water molecules can affect diffusion of water molecules inside amyloid fibrils, inducing their different structural stabilities. The conformation of amyloid fibrils induced by interior water molecules show the fibrils' different mechanical features. We elucidate the role of confined and movable interior water molecules in structural stabilities of various amyloid fibrils. Our results offer insights not only in further understanding of mechanical features of amyloids as mediated by water molecules, but also in the fine-tuning of the functional abilities of amyloid fibrils for applications.

Primary Fibril Nucleation of Aggregation Prone Tau Fragments PHF6 and PHF6

We performed replica exchange molecular dynamics and forward flux sampling simulations of hexapeptide VQIINK and VQIVYK systems, also known as, respectively, fragments PHF6* and PHF6 from the tau protein. Being a part of the microtubule binding region, these fragments are known to be aggregation prone, and at least one of them is a prerequisite for fibril formation of the tau protein. Using a coarse-grained force field, we establish the phase behavior of both fragments, and investigate the nucleation kinetics for the conversion into a β-sheet fibril. As the conversion is, in principle, a reversible process, we predict the rate constants for both the fibril formation and melting, and examine the corresponding mechanisms. Our simulations indicate that, while both fragments form disordered aggregates, only PHF6 is able to form β-sheet fibrils. This observation provides a possible explanation for the lack of available steric zipper crystal structures for PHF6*.

An Atomistic View of Amyloidogenic Self-assembly: Structure and Dynamics of Heterogeneous Conformational States in the Pre-nucleation Phase

The formation of well-defined filamentous amyloid structures involves a polydisperse collection of oligomeric states for which relatively little is known in terms of structural organization. Here we use extensive, unbiased explicit solvent molecular dynamics (MD) simulations to investigate the structural and dynamical features of oligomeric aggregates formed by a number of highly amyloidogenic peptides at atomistic resolution on the μs time scale. A consensus approach has been adopted to analyse the simulations in multiple force fields, yielding an in-depth characterization of pre-fibrillar oligomers and their global and local structure properties. A collision cross section analysis revealed structurally heterogeneous aggregate ensembles for the individual oligomeric states that lack a single defined quaternary structure during the pre-nucleation phase. To gain insight into the conformational space sampled in early aggregates, we probed their substructure and found emerging β-sheet subunit layers and a multitude of ordered intermolecular β-structure motifs with growing aggregate size. Among those, anti-parallel out-of-register β-strands compatible with toxic β-barrel oligomers were particularly prevalent already in smaller aggregates and formed prior to ordered fibrillar structure elements. Notably, also distinct fibril-like conformations emerged in the oligomeric state and underscore the notion that pre-nucleated oligomers serve as a critical intermediate step on-pathway to fibrils.

Thermodynamics of amyloid formation and the role of intersheet interactions

The self-assembly of proteins into β-sheet-rich amyloid fibrils has been observed to occur with sigmoidal kinetics, indicating that the system initially is trapped in a metastable state. Here, we use a minimal lattice-based model to explore the thermodynamic forces driving amyloid formation in a finite canonical (NVT) system. By means of generalized-ensemble Monte Carlo techniques and a semi-analytical method, the thermodynamic properties of this model are investigated for different sets of intersheet interaction parameters. When the interactions support lateral growth into multi-layered fibrillar structures, an evaporation/condensation transition is observed, between a supersaturated solution state and a thermodynamically distinct state where small and large fibril-like species exist in equilibrium. Intermediate-size aggregates are statistically suppressed. These properties do not hold if aggregate growth is one-dimensional.

Peptide folding in the presence of interacting protein crowders

Using Monte Carlo methods, we explore and compare the effects of two protein crowders, BPTI and GB1, on the folding thermodynamics of two peptides, the compact helical trp-cage and the β-hairpin-forming GB1m3. The thermally highly stable crowder proteins are modeled using a fixed backbone and rotatable side-chains, whereas the peptides are free to fold and unfold. In the simulations, the crowder proteins tend to distort the trp-cage fold, while having a stabilizing effect on GB1m3. The extent of the effects on a given peptide depends on the crowder type. Due to a sticky patch on its surface, BPTI causes larger changes than GB1 in the melting properties of the peptides. The observed effects on the peptides stem largely from attractive and specific interactions with the crowder surfaces, and differ from those seen in reference simulations with purely steric crowder particles.

Tau assembly: the dominant role of PHF6 (VQIVYK) in microtubule binding region repeat R3

Self-aggregation of the microtubule-binding protein Tau reduces its functionality and is tightly associated with Tau-related diseases, termed tauopathies. Tau aggregation is also strongly associated with two nucleating six-residue segments, namely PHF6 (VQIVYK) and PHF6* (VQIINK). In this paper, using experiments and computational modeling, we study the self-assembly of individual and binary mixtures of Tau fragments containing PHF6* (R2/wt; (273)GKVQIINKKLDL(284)) and PHF6 (R3/wt; (306)VQIVYKPVDLSK(317)) and a mutant R2/ΔK280 associated with a neurodegenerative tauopathy. The initial stage of aggregation is probed by ion-mobility mass spectrometry, the kinetics of aggregation monitored with Thioflavin T assays, and the morphology of aggregates visualized by transmission electron microscopy. Insights into the structure of early aggregates and the factors stabilizing the aggregates are obtained from replica exchange molecular dynamics simulations. Our data suggest that R3/wt has a much stronger aggregation propensity than either R2/wt or R2/ΔK280. Heterodimers containing R3/wt are less stable than R3/wt homodimers but much more stable than homodimers of R2/wt and R2/ΔK280, suggesting a possible role of PHF6*-PHF6 interactions in initiating the aggregation of full-length Tau. Lastly, R2/ΔK280 binds more strongly to R3/wt than R2/wt, suggesting a possible mechanism for a pathological loss of normal Tau function.

Molecular basis for structural heterogeneity of an intrinsically disordered protein bound to a partner by combined ESI-IM-MS and modeling

Intrinsically disordered proteins (IDPs) form biologically active complexes that can retain a high degree of conformational disorder, escaping structural characterization by conventional approaches. An example is offered by the complex between the intrinsically disordered N(TAIL) domain and the phosphoprotein X domain (P(XD)) from measles virus (MeV). Here, distinct conformers of the complex are detected by electrospray ionization-mass spectrometry (ESI-MS) and ion mobility (IM) techniques yielding estimates for the solvent-accessible surface area (SASA) in solution and the average collision cross-section (CCS) in the gas phase. Computational modeling of the complex in solution, based on experimental constraints, provides atomic-resolution structural models featuring different levels of compactness. The resulting models indicate high structural heterogeneity. The intermolecular interactions are predominantly hydrophobic, not only in the ordered core of the complex, but also in the dynamic, disordered regions. Electrostatic interactions become involved in the more compact states. This system represents an illustrative example of a hydrophobic complex that could be directly detected in the gas phase by native mass spectrometry. This work represents the first attempt to modeling the entire N(TAIL) domain bound to P(XD) at atomic resolution.

Kinetic Monte Carlo simulations of the assembly of filamentous biomacromolecules by dimer addition mechanism

In cells, several important biomacromolecules form oligomers through a dimer addition mechanism. Rate equations based on mean field approximations are usually employed to describe the assembly process. However, such equations often require multiple assumptions that mask some detailed changes of the biomolecular configurations during aggregations. Here, we present a Kinetic Monte Carlo simulation scheme to account for the diffusion and rotation of dimers on two-dimensional hexagonal lattices while naturally including the stochastic features. We investigate the effects of the interaction energy between dimers, the diffusion coefficient and the concentration of dimers on the aggregation by dimer addition mechanism. Our simulations identified unusual double-S shape evolutions of aggregation kinetics, which are probably associated with the formation of metastable clusters.