Terminal Capping of an Amyloidogenic Tau Fragment Modulates Its Fibrillation Propensity

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.

Effect of Cosolutes on the Aggregation of a Tau Fragment: A Combined Experimental and Simulation Approach

The intrinsically disordered protein Tau represents the main component of neurofibrillary tangles that are a hallmark of Alzheimer's disease. A small fragment of Tau, known as paired helical filament 6 (PHF6), is considered to be important for the formation of the β-structure core of the fibrils. Here we study the aggregation of this fragment in the presence of different cosolutes, including urea, TMAO, sucrose and 2-hydroxypropyl-β-cyclodextrin (2-HPβCD), using both experiments and molecular dynamics simulations. A novel implicit solvation approach (MIST - Model with Implicit Solvation Thermodynamics) is used, where an energetic contribution based on the concept of transfer free energies describes the effect of the cosolutes. The simulation predictions are compared to thioflavin-T and atomic force microscopy results, and the good agreement observed confirms the predictive ability of the computational approach herein proposed. Both simulations and experiments indicate that PHF6 aggregation is inhibited in the presence of urea and 2-HPβCD, while TMAO and sucrose stabilize associated conformations. The remarkable ability of HPβCD to inhibit aggregation represents an extremely promising result for future applications, especially considering the widespread use of this molecule as a drug carrier to the brain and as a solubilizer/excipient in pharmaceutical formulations.

Differential Binding and Conformational Dynamics of Tau Microtubule-Binding Repeats with a Preformed Amyloid-β Fibril Seed

Both senile plaques formed by amyloid-β (Aβ) and neurofibrillary tangles (NFTs) comprised of tau are pathological hallmarks of Alzheimer's disease (AD). The accumulation of NFTs better correlates with the loss of cognitive function than senile plaques, but NFTs are rarely observed without the presence of senile plaques. Hence, cross-seeding of tau by preformed Aβ amyloid fibril seeds has been proposed to drive the aggregation of tau and exacerbate AD progression, but the molecular mechanism remains unknown. Here, we first identified cross-interaction hotspots between Aβ and tau using atomistic discrete molecular dynamics simulations (DMD) and confirmed the critical role of the four microtubule-binding repeats of tau (R1-R4) in the cross-interaction with Aβ. We further investigated the binding structure and dynamics of each tau repeat with a preformed Aβ fibril seed. Specifically, R1 and R3 preferred to bind the Aβ fibril lateral surface instead of the elongation end. In contrast, R2 and R4 had higher binding propensities to the fibril elongation end than the lateral surface, enhancing β-sheet content by forming hydrogen bonds with the exposed hydrogen bond donors and acceptors. Together, our results suggest that the four repeats play distinct roles in driving the binding of tau to different surfaces of an Aβ fibril seed. Binding of tau to the lateral surface of Aβ fibril can increase the local concentration, while the binding to the elongation surface promotes β-sheet formation, both of which reduce the free energy barrier for tau aggregation nucleation and subsequent fibrillization.

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.

On the Tracks of the Aggregation Mechanism of the PHF6 Peptide from Tau Protein: Molecular Dynamics, Energy, and Interaction Network Investigations

The formation of neurofibrillary tangles (NFTs), composed of tau protein aggregates, is a hallmark of some neurodegenerative diseases called tauopathies. NFTs are composed of paired helical filaments (PHFs) of tau protein with a dominant β-sheet secondary structuration. The NFT formation mechanism is not known yet. This study focuses on PHF6, a crucial hexapeptide responsible for tau aggregation. A 2 μs molecular dynamics simulation was launched to determine the keys of the PHF6 aggregation mechanism. Hydrogen bonding, van der Waals, and other non-covalent interactions as π-stacking were investigated. Parallel aggregation was slightly preferred due to its adaptability, but antiparallel aggregation remained widely present during the PHF6 aggregation. The analysis highlighted the leading role of hydrogen bonds identified at the atomic level for each aggregation process. The aggregation study emphasized the importance of Tyr310 during the β-sheets' complexation through π-stacking.

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.

Molecular Dynamics Simulations of the Tau R3-R4 Domain Monomer in the Bulk Solution and at the Surface of a Lipid Bilayer Model

The aggregation of the tau protein plays a significant role in Alzheimer's disease, and the tau R3-R4 domain spanning residues 306-378 was shown to form the amyloid fibril core of a full-length tau. The conformations of the tau R3-R4 monomer in the bulk solution and at the surface of membranes are unknown. In this study, we address these questions by means of atomistic molecular dynamics. The simulations in the bulk solution show a very heterogeneous ensemble of conformations with low β and helical contents. The tau R3-R4 monomer has the propensity to form transient β-hairpins within the R3 repeat and between the R3 and R4 repeats and parallel β-sheets spanning the R3 and R4 repeats. The simulations also show that the surface of the membrane does not induce β-sheet insertion and leads to an ensemble of structures very different from those in the bulk solution. They also reveal the dynamical properties of the membrane-bound state of the tau R3-R4 monomer, enabling insertion of the residues 306-318 and 376-378.

Structurally Distinct Polymorphs of Tau Aggregates Revealed by Nanoscale Infrared Spectroscopy

Aggregation of the tau protein plays a central role in several neurodegenerative diseases collectively known as tauopathies, including Alzheimer's and Parkinson's disease. Tau misfolds into fibrillar β sheet structures that constitute the paired helical filaments found in neurofibrillary tangles. It is known that there can be significant structural heterogeneities in tau aggregates associated with different diseases. However, while structures of mature fibrils have been studied, the structural distributions in early-stage tau aggregates is not well-understood. In the present study, we use atomic force microscopy-IR to investigate nanoscale spectra of individual tau fibrils at different stages of aggregation and demonstrate the presence of multiple fibrillar polymorphs that exhibit different secondary structures. We further show that mature fibrils contain significant amounts of antiparallel β sheets. Our results are the very first application of nanoscale infrared spectroscopy to tau aggregates and underscore the promise of spatially resolved infrared spectroscopy for investigating protein aggregation.

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.

Heparin remodels the microtubule-binding repeat R3 of Tau protein towards fibril-prone conformations

Abnormal aggregation of proteins into pathological amyloid fibrils is implicated in a wide range of devastating human neurodegenerative diseases. Intracellular fibrillary inclusions formed by Tau protein are characterized as the hallmark of tauopathies, including Alzheimer's disease and frontotemporal dementia. Heparin has been often used to trigger Tau aggregation in in vitro studies. However, the conformational changes induced by heparin and the underlying mechanism of promotion of Tau aggregation by heparin are not well understood. Structural characterization of Tau oligomers in the early stage of fibrillation is of great importance but remains challenging due to their dynamic and heterogeneous nature. R3, the third microtubule-binding repeat of Tau, contains the fibril-nucleating core (PHF6) and is crucial for Tau aggregation. In this study, utilizing extensive all-atom replica-exchange molecular dynamic simulations, we explored the conformational ensembles of R3 monomer/dimer in the absence and presence of heparin. Our results show that without heparin, both monomeric and dimeric R3 preferentially adopt collapsed β-sheet-containing conformations and PHF6 plays an important role in the formation of interchain β-sheet structures, while in the presence of heparin, R3 can populate relatively extended disordered states where chain dimension is similar to that of R3 in Tau filaments. Through electrostatic, hydrogen-bonding and hydrophobic interactions, heparin has a preference for interacting with residues V306/Q307/K317/K321/H329/H330/K331 which distribute throughout the entire sequence of R3, in turn acting as a template to extend R3 conformations. More importantly, heparin alters intramolecular/intermolecular interaction patterns of R3 and increases the intermolecular contact regions. Our results suggest that heparin remodels the conformations of R3 towards fibril-prone structures by increasing chain dimension and intermolecular contact regions, which may shed light on the atomic mechanism of heparin-induced amyloid fibrillization of Tau protein.

Liquid-liquid phase separation of Tau by self and complex coacervation

The liquid-liquid phase separation (LLPS) of Tau has been postulated to play a role in modulating the aggregation property of Tau, a process known to be critically associated with the pathology of a broad range of neurodegenerative diseases including Alzheimer's Disease. Tau can undergo LLPS by homotypic interaction through self-coacervation (SC) or by heterotypic association through complex-coacervation (CC) between Tau and binding partners such as RNA. What is unclear is in what way the formation mechanisms for self and complex coacervation of Tau are similar or different, and the addition of a binding partner to Tau alters the properties of LLPS and Tau. A combination of in vitro experimental and computational study reveals that the primary driving force for both Tau CC and SC is electrostatic interactions between Tau-RNA or Tau-Tau macromolecules. The liquid condensates formed by the complex coacervation of Tau and RNA have distinctly higher micro-viscosity and greater thermal stability than that formed by the SC of Tau. Our study shows that subtle changes in solution conditions, including molecular crowding and the presence of binding partners, can lead to the formation of different types of Tau condensates with distinct micro-viscosity that can coexist as persistent and immiscible entities in solution. We speculate that the formation, rheological properties and stability of Tau droplets can be readily tuned by cellular factors, and that liquid condensation of Tau can alter the conformational equilibrium of Tau.

Misfolding and Self-Assembly Dynamics of Microtubule-Binding Repeats of the Alzheimer-Related Protein Tau

Pathological aggregation of intrinsically disordered tau protein, driven by the interactions between microtubule-binding (MTB) domains, is associated with Alzheimer's disease. The MTB domain contains either three or four repeats with sequence similarities. Compared to amyloid β, many aspects of the misfolding and aggregation mechanisms of tau are largely unknown. In this study, we systematically investigated the dynamics of monomer misfolding and dimerization of each MTB repeat using atomistic discrete molecular dynamic simulations. Our results revealed that all the four repeat monomers (R1-R4) were very dynamic, featuring frequent conformational conversion and lacking stable conformations. While R1, R2, and R4 monomers occasionally adopted partially helical conformations, R3 monomers frequently formed β-sheets. In dimerization simulations, R3 displayed the strongest aggregation propensity with high β-sheet contents, while R1 was the least prone to aggregation. The R2 and R4 dimers contained both helix and β-sheet structures. The β-sheets in R4 assemblies were dominant with β-hairpin conformation. In R2 and R3 dimers, intermolecular β-sheets were mainly driven by residues around the paired helical filament (PHF) regions. Residues around the PHF6* in R2 and PHF6 in R3 had significantly higher intermolecular contacts than other regions, suggesting that these residues play a key role in the amyloid aggregation of tau. Our results on the structural ensembles and early aggregation dynamics of each tau MTB repeat will help understand the nucleation and fibrillization of tau.

Liquid-Liquid Phase Separation of Tau Protein Is Encoded at the Monomeric Level

Liquid-liquid phase separation (LLPS) is involved in both physiological and pathological processes. The intrinsically disordered protein Tau and its K18 construct can undergo LLPS in a distinct temperature-dependent manner, and the LLPS of Tau protein can initiate Tau aggregation. However, the underlying mechanism driving Tau LLPS remains largely elusive. To understand the temperature-dependent LLPS behavior of Tau at the monomeric level, we explored the conformational ensemble of Tau at different temperatures by performing all-atom replica-exchange molecular dynamic simulation on K18 monomer with an accumulated simulation time of 26.4 μs. Our simulation demonstrates that the compactness, β-structure propensity, and intramolecular interaction of K18 monomer exhibit nonlinear temperature-dependent behavior. ²⁹⁵DNIKHV³⁰⁰/³²⁶GNIHHK³³¹/³³⁷VEVKSE³⁴² make significant contributions to the temperature dependence of the β propensity of K18 monomer, while the two fibril-nucleating cores display relatively high β propensity at all temperatures. At a specific temperature, K18 monomer adopts the most collapsed state with exposed sites for both persistent and transient interactions. Given that more collapsed polypeptide chains were reported to be more prone to phase separate, our results suggest that K18 monomer inherently possesses conformational characteristics favoring LLPS. Our simulation predicts the importance of ²⁹⁵DNIKHV³⁰⁰/³²⁶GNIHHK³³¹/³³⁷VEVKSE³⁴² to the temperature-dependent conformational properties of K18, which is corroborated by CD spectra, turbidity assays, and DIC microscopy. Taken together, we offer a computational and experimental approach to comprehend the structural basis for LLPS by amyloidal building blocks.