Deciphering the Effect of Lysine Acetylation on the Misfolding and Aggregation of Human Tau Fragment 171IPAKTPPAPK180 Using Molecular Dynamic Simulation and the Markov State Model.

The formation of neurofibrillary tangles (NFT) with ß-sheet-rich structure caused by abnormal aggregation of misfolded microtubule-associated protein Tau is a hallmark of tauopathies, including Alzheimer's Disease. It has been reported that acetylation, especially K174 located in the proline-rich region, can largely promote Tau aggregation. So far, the mechanism of the abnormal acetylation of Tau that affects its misfolding and aggregation is still unclear. Therefore, revealing the effect of acetylation on Tau aggregation could help elucidate the pathogenic mechanism of tauopathies. In this study, molecular dynamics simulation combined with multiple computational analytical methods were performed to reveal the effect of K174 acetylation on the spontaneous aggregation of Tau peptide 171IPAKTPPAPK180, and the dimerization mechanism as an early stage of the spontaneous aggregation was further specifically analyzed by Markov state model (MSM) analysis. The results showed that both the actual acetylation and the mutation mimicking the acetylated state at K174 induced the aggregation of the studied Tau fragment; however, the effect of actual acetylation on the aggregation was more pronounced. In addition, acetylated K174 plays a major contributing role in forming and stabilizing the antiparallel ß-sheet dimer by forming several hydrogen bonds and side chain van der Waals interactions with residues I171, P172, A173 and T175 of the corresponding chain. In brief, this study uncovered the underlying mechanism of Tau peptide aggregation in response to the lysine K174 acetylation, which can deepen our understanding on the pathogenesis of tauopathies.

Assessment of Amyloid Forming Tendency of Peptide Sequences from Amyloid Beta and Tau Proteins Using Force-Field, Semi-Empirical, and Density Functional Theory Calculations.

A wide variety of neurodegenerative diseases are characterized by the accumulation of protein aggregates in intraneuronal or extraneuronal brain regions. In Alzheimer's disease (AD), the extracellular aggregates originate from amyloid-ß proteins, while the intracellular aggregates are formed from microtubule-binding tau proteins. The amyloid forming peptide sequences in the amyloid-ß peptides and tau proteins are responsible for aggregate formation. Experimental studies have until the date reported many of such amyloid forming peptide sequences in different proteins, however, there is still limited molecular level understanding about their tendency to form aggregates. In this study, we employed umbrella sampling simulations and subsequent electronic structure theory calculations in order to estimate the energy profiles for interconversion of the helix to ß-sheet like secondary structures of sequences from amyloid-ß protein (KLVFFA) and tau protein (QVEVKSEKLD and VQIVYKPVD). The study also included a poly-alanine sequence as a reference system. The calculated force-field based free energy profiles predicted a flat minimum for monomers of sequences from amyloid and tau proteins corresponding to an α-helix like secondary structure. For the parallel and anti-parallel dimer of KLVFFA, double well potentials were obtained with the minima corresponding to α-helix and ß-sheet like secondary structures. A similar double well-like potential has been found for dimeric forms for the sequences from tau fibril. Complementary semi-empirical and density functional theory calculations displayed similar trends, validating the force-field based free energy profiles obtained for these systems.

The misfolding mechanism of the key fragment R3 of tau protein: a combined molecular dynamics simulation and Markov state model study.

The formation of neurofibrillary tangles (NFT) by abnormal aggregation of misfolded microtubule-associated protein tau is a hallmark of tauopathies, including Alzheimer's disease. However, it remains unclear how tau monomers undergo conformational changes and further lead to the abnormal aggregation. In this work, molecular dynamics simulation combined with the Markov state model (MSM) analysis was used to uncover the misfolding progress and structural characteristics of the key R3 fragment of tau protein at the atomic level. The simulation results show that R3 exists in disordered structures mainly, which is consistent with the experimental results. The MSM analysis identified multiple ß-sheet conformations of R3. The residues involved in the ß-sheet structure formation are mainly located in three regions: PHF6 at the N-terminal, S324 to N327 at the middle of R3, and K331 to G334 at the C-terminal. In addition, the path analysis of the formation of the ß-sheet structure by transition path theory (TPT) revealed that there are multiple paths to form ß-sheet structures from the disordered state, and the timescales are at the millisecond level, indicating that a large number of structural rearrangements occur during the formation of ß-sheet structures. It is interesting to note that S19 is a critical intermediate state for the formation of two target ß-sheet structures, S23 and S4. In S19, three regions of V306 to K311, C322 to G326, and K331 to G334 form a turn structure, the regions that form the ß-sheet structure in target states S23 and S4, indicating that the formation of a turn structure is necessary to form a ß-sheet structure and then the turn structure will eventually transform into the ß-sheet structure through key hydrogen bonding interactions. These findings can provide insights into the kinetics of tau protein misfolding.

Disclosing the Mechanism of Spontaneous Aggregation and Template-Induced Misfolding of the Key Hexapeptide (PHF6) of Tau Protein Based on Molecular Dynamics Simulation.

The microtubule-associated protein tau is critical for the development and maintenance of the nervous system. Tau dysfunction is associated with a variety of neurodegenerative diseases called tauopathies, which are characterized by neurofibrillary tangles formed by abnormally aggregated tau protein. Studying the aggregation mechanism of tau protein is of great significance for elucidating the etiology of tauopathies. The hexapeptide 306VQIVYK311 (PHF6) of R3 has been shown to play a vital role in promoting tau aggregation. In this study, long-term all-atom molecular dynamics simulations in explicit solvent were performed to investigate the mechanisms of spontaneous aggregation and template-induced misfolding of PHF6, and the dimerization at the early stage of nucleation was further specifically analyzed by the Markov state model (MSM). Our results show that PHF6 can spontaneously aggregate to form multimers enriched with ß-sheet structure and the ß-sheets in multimers prefer to exist in a parallel way. It is observed that PHF6 monomer can be induced to form a ß-sheet structure on either side of the template but in a different way. In detail, the ß-sheet structure is easier to form on the left side but does not extend well, but on the right side, the monomer can form the extended ß-sheet structure. Furthermore, MSM analysis shows that the formation of dimer mainly occurs in three steps. First, the separated monomers collide with each other at random orientations, and then a dimer with short ß-sheet structure at the N-terminal forms; finally, ß-sheets elongate to form an extended parallel ß-sheet dimer. During these processes, multiple intermediate states are identified and multiple paths can form a parallel ß-sheet dimer from the disordered coil structure. Moreover, the residues I308, V309, and Y310 play an essential role in the dimerization. In a word, our results uncover the aggregation and misfolding mechanism of PHF6 from the atomic level, which can provide useful theoretical guidance for rational design of effective therapeutic drugs against tauopathies.

Assessment of the efficacy of different procedures that remove and disassemble alpha-synuclein, tau and A-beta fibrils from laboratory material and surfaces.

α-synuclein fibrillar polymorphs, Tau and Aß 1-42 fibrillar assemblies have been shown to propagate, amplify and trigger the formation of protein deposits reminiscent of those present within the central nervous system of patients developing synucleinopathies, tauopathies and amyloid plaques after injection intracerebrally, intramuscularly, intraperitoneally or within the blood stream of model animals. They are thus hazardous and there is need for decontamination and inactivation procedures for laboratory surfaces and non-disposable material. We assessed the effectiveness of different reagents to clean and disassemble potentially pathogenic assemblies adsorbed on non-disposable materials in laboratories. We show that commercial detergents and SDS are way more suited to detach α-synuclein fibrillar polymorphs, Tau and Aß 1-42 fibrillar assemblies from contaminated surfaces and disassemble the fibrils than methods designed to decrease PrP prion infectivity. Our observations reveal that the choice of the most adapted cleaning procedure for one given protein assembly or fibrillar polymorph should integrate detergent's cleaning efficiency, material compatibility and capacity to dismantle assemblies. We provide an integrated representation where desorption and neutralization efficacy and surface compatibility are combined to facilitate the choice of the most adapted decontamination procedure. This representation, together with good laboratory practices, contributes to reducing potential health hazards associated to manipulating protein assemblies with prion-like properties.

The conformational ensembles of α-synuclein and tau: combining single-molecule FRET and simulations.

Intrinsically disordered proteins (IDPs) are increasingly recognized for their important roles in a range of biological contexts, both in normal physiological function and in a variety of devastating human diseases. However, their structural characterization by traditional biophysical methods, for the purposes of understanding their function and dysfunction, has proved challenging. Here, we investigate the model IDPs α-Synuclein (αS) and tau, that are involved in major neurodegenerative conditions including Parkinson's and Alzheimer's diseases, using excluded volume Monte Carlo simulations constrained by pairwise distance distributions from single-molecule fluorescence measurements. Using this, to our knowledge, novel approach we find that a relatively small number of intermolecular distance constraints are sufficient to accurately determine the dimensions and polymer conformational statistics of αS and tau in solution. Moreover, this method can detect local changes in αS and tau conformations that correlate with enhanced aggregation. Constrained Monte Carlo simulations produce ensembles that are in excellent agreement both with experimental measurements on αS and tau and with all-atom, explicit solvent molecular dynamics simulations of αS, with much lower configurational sampling requirements and computational expense.

Monte carlo simulations of tau proteins: effect of phosphorylation.

We perform Monte Carlo simulations of tau proteins bound to a cylinder that mimics a microtubule (MT), and then study them in solution. Tau protein binds to a highly anionic MT surface to stabilize the cylindrical structure of MT. The negatively charged tail domain floats away from the anionic MT surface while positively charged tau segments localize near the MT surface. Monte Carlo simulations demonstrate that, in 3RS tau isoform (which has three imperfect repeats (R) short (S) isoform), amino acids are more condensed near a highly charged interface compared to 4RL isoform (which has four imperfect repeats (R) long (L) isoform). In 4RL isoform, amino acids in tail domain stay mostly apart from the MT surface. In the bulk solution, dephosphorylated taus are separated due to Coulomb repulsion between similarly charged isoforms. Moderate phosphorylation of 3RS isoform decreases average intermolecular distance between dephosphorylated and phosphorylated taus and lead to their overlap. Further phosphorylation does not change noticeably the intermolecular distances.

A thorough assessment of benign genetic variability in GRN and MAPT.

Mutations in APP, PSEN1, MAPTand GRNare the most common genetic causes of dementia. The previous miss-assignment of pathogenicity to benign variants in these genes stresses the importance of discerning between disease causing mutations and benign variants with no pathogenic effect on the function of the respective protein. In this study we sequenced GRNand MAPTin 282 samples from the Centre d'Etude du Polymorphisme Humain - Human Genome Diversity Cell Line Panel, in order to identify benign variants that could otherwise be mistaken for pathogenic mutations. We found sixteen different non-synonymous changes, eleven of which are novel variants.

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

Small oligomers formed early in the process of amyloid fibril formation may be the major toxic species in Alzheimer's disease. We investigate the early stages of amyloid aggregation for the tau fragment AcPHF6 (Ac-VQIVYK-NH2) using an implicit solvent all-atom model and extensive Monte Carlo simulations of 12, 24, and 36 chains. A variety of small metastable aggregates form and dissolve until an aggregate of a critical size and conformation arises. However, the stable oligomers, which are beta-sheet-rich and feature many hydrophobic contacts, are not always growth-ready. The simulations indicate instead that these supercritical oligomers spend a lengthy period in equilibrium in which considerable reorganization takes place accompanied by exchange of chains with the solution. Growth competence of the stable oligomers correlates with the alignment of the strands in the beta-sheets. The larger aggregates seen in our simulations are all composed of two twisted beta-sheets, packed against each other with hydrophobic side chains at the sheet-sheet interface. These beta-sandwiches show similarities with the proposed steric zipper structure for PHF6 fibrils but have a mixed parallel/antiparallel beta-strand organization as opposed to the parallel organization found in experiments on fibrils. Interestingly, we find that the fraction of parallel beta-sheet structure increases with aggregate size. We speculate that the reorganization of the beta-sheets into parallel ones is an important rate-limiting step in the formation of PHF6 fibrils.