Tau R3-R4 Domain Dimer of the Wild Type and Phosphorylated Ser356 Sequences. I. In Solution by Atomistic Simulations

Influence of ALS-linked M337V mutation on the conformational ensembles of TDP-43321-340 peptide monomer and dimer

The transactive response (TAR) DNA/RNA-binding protein 43 (TDP-43) can self-assemble into both functional stress granules via liquid-liquid phase separation (LLPS) and pathogenic amyloid fibrillary aggregates that are closely linked to amyotrophic lateral sclerosis. Previous experimental studies reported that the low complexity domain (LCD) of TDP-43 plays an essential role in the LLPS and aggregation of the full-length protein, and it alone can also undergo LLPS to form liquid droplets mainly via intermolecular interactions in the 321-340 region. And the ALS-associated M337V mutation impairs LCD's LLPS and facilitates liquid-solid phase transition. However, the underlying atomistic mechanism is not well understood. Herein, as a first step to understand the M337V-caused LLPS disruption of TDP-43 LCD mediated by the 321-340 region and the fibrillization enhancement, we investigated the conformational properties of monomer/dimer of TDP-43321-340 peptide and its M337V mutant by performing extensive all-atom explicit-solvent replica exchange molecular dynamic simulations. Our simulations demonstrate that M337V mutation alters the residue regions with high helix/β-structure propensities and thus affects the conformational ensembles of both monomer and dimer. M337V mutation inhibits helix formation in the N-terminal Ala-rich region and the C-terminal mutation site region, while facilitating their long β-sheet formation, albeit with a minor impact on the average probability of both helix structure and β-structure. Further analysis of dimer system shows that M337V mutation disrupts inter-molecular helix-helix interactions and W334-W334 π-π stacking interactions which were reported to be important for the LLPS of TDP-43 LCD, whereas enhances the overall peptide residue-residue interactions and weakens peptide-water interactions, which is conducive to peptide fibrillization. This study provides mechanistic insights into the M337V-mutation-induced impairment of phase separation and facilitation of fibril formation of TDP-43 LCD.

A theoretical study on the activity and selectivity of IDO/TDO inhibitors

Indoleamine 2,3-dioxygenase 1 (IDO) is a tryptophan (Trp) metabolic enzyme along the kynurenine (NFK) pathway. Under pathological conditions, IDO overexpressed by tumor cells causes depletion of tryptophan and the accumulation of metabolic products, which inhibit the local immune response and form immune escape. Therefore, the suppression of IDO activity is one of the strategies for tumor immunotherapy, and drug design for this target has been the focus of research for more than two decades. Apart from IDO, tryptophan dioxygenase (TDO) of the same family can also catalyze the same biochemical reaction in the human body, but it has different tissue distribution and substrate selectivity from IDO. Based on the principle of drug design with high potency and low cross-reactivity to specific targets, in this subject, the activity and selectivity of IDO and TDO toward small molecular inhibitors were studied from the perspective of thermodynamics and kinetics. The aim was to elucidate the structural requirements for achieving favorable biological activity and selectivity of IDO and TDO inhibitors. Specifically, the interactions of inhibitors from eight families with IDO and TDO were initially investigated through molecular docking and molecular dynamics simulations, and the thermodynamic data for binding of inhibitors were predicted by the molecular mechanics/generalized Born surface area (MM/GBSA) method. Secondly, we explored the free energy landscape of JKloops, the kinetic control element of IDO/TDO, using temperature replica exchange molecular dynamics (T-REMD) simulations and elucidated the connection between the rules of IDO/TDO conformational changes and the inhibitor selectivity mechanism. Furthermore, the binding and dissociation processes of the C1 inhibitor (NLG919) were simulated by the adaptive steering molecular dynamics (ASMD) method, which not only addressed the possible stable, metastable, and transition states for C1 inhibitor-IDO/TDO interactions, but also accurately predicted kinetic data for C1 inhibitor binding and dissociation. In conclusion, we have constructed a complete process from enzyme (IDO/TDO) conformational activation to inhibitor binding/dissociation and used the thermodynamic and kinetic data of each link as clues to verify the control mechanism of IDO/TDO on inhibitor selectivity. This is of great significance for us to understand the design principles of tumor immunotherapy drugs and to avoid drug resistance of immunotherapy drugs.

Role of G326 in Determining the Aggregation Propensity of R3 Tau Repeat: Insights from Studies on R1R3 Tau Construct

The Microtubule-binding repeat region (MTBR) of Tau has been studied extensively due to its pathological implications in neurodegenerative diseases like Alzheimer's disease. The pathological property of MTBR is mainly due to the R3 repeat's high propensity for self-aggregation, highlighting the critical molecular grammar of the repeat. Utilizing the R1R3 construct (WT) and its G326E mutant (EE), we determine the distinct characteristics of various peptide segments that modulate the aggregation propensity of the R3 repeat using NMR spectroscopy. Through time-dependent experiments, we have identified 317KVTSKCGS324 in R3 repeat as the aggregation initiating motif (AIM) due to its role at the initial stages of aggregation. The G326E mutation induces changes in conformation and dynamics at the AIM, thereby effectively abrogating the aggregation propensity of the R1R3 construct. We further corroborate our findings through MD simulations and propose that AIM is a robust site of interest for tauopathy drug design.

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.

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.

Molecular Mechanism in the Disruption of Chronic Traumatic Encephalopathy-Related R3-R4 Tau Protofibril by Quercetin and Gallic Acid: Similarities and Differences

Chronic traumatic encephalopathy (CTE) is a unique progressive neurodegenerative tauopathy pathologically related to the aggregation of the tau protein to neurofibrillary tangles. Disrupting tau oligomers (protofibril) is a promising strategy to prevent CTE. Quercetin (QE) and gallic acid (GA), two polyphenol small molecules abundant in natural crops, were proved to inhibit recombinant tau and the R3 fragment of human full-length tau in vitro. However, their disruptive effect on CTE-related protofibril and the underlying molecular mechanism remain elusive. Cryo-electron microscopy resolution reveals that the R3-R4 fragment of tau forms the core of the CTE-related tau protofibril. In this study, we conducted extensive all-atom molecular dynamics simulations on CTE-related R3-R4 tau protofibril with and without QE/GA molecules. The results disclose that both QE and GA can disrupt the global structure of the protofibril, while GA shows a relatively strong effect. The binding sites, exact binding patterns, and disruptive modes for the two molecules show similarities and differences. Strikingly, both QE and GA can insert into the hydrophobic cavity of the protofibril, indicating they have the potential to compete for the space in the cavity with aggregation cofactors unique to CTE-related protofibril and thus impede the further aggregation of the tau protein. Due to relatively short time scale, our study captures the early disruptive mechanism of CTE-related R3-R4 tau protofibril by QE/GA. However, our research does provide valuable knowledge for the design of supplements or drugs to prevent or delay the development of CTE.

Phosphorylation at Ser289 Enhances the Oligomerization of Tau Repeat R2

In tauopathies such as Alzheimer's disease (AD), aberrant phosphorylation causes the dissociation of tau proteins from microtubules. The dissociated tau then aggregates into sequent forms from soluble oligomers to paired helical filaments and insoluble neurofibrillary tangles (NFTs). NFTs is a hallmark of AD, while oligomers are found to be the most toxic form of the tau aggregates. Therefore, understanding tau oligomerization with regard to abnormal phosphorylation is important for the therapeutic development of AD. In this study, we investigated the impact of phosphorylated Ser289, one of the 40 aberrant phosphorylation sites of full-length tau proteins, on monomeric and dimeric structures of tau repeat R2 peptides. We carried out intensive replica exchange molecular dynamics simulation with a total simulation time of up to 0.1 ms. Our result showed that the phosphorylation significantly affected the structures of both the monomer and the dimer. For the monomer, the phosphorylation enhanced ordered-disordered structural transition and intramolecular interaction, leading to more compactness of the phosphorylated R2 compared to the wild-type one. As to the dimer, the phosphorylation increased intermolecular interaction and β-sheet formation, which can accelerate the oligomerization of R2 peptides. This result suggests that the phosphorylation at Ser289 is likely to promote tau aggregation. We also observed a phosphorylated Ser289-Na+-phosphorylated Ser289 bridge in the phosphorylated R2 dimer, suggesting an important role of cation ions in tau aggregation. Our findings suggest that phosphorylation at Ser289 should be taken into account in the inhibitor screening of tau oligomerization.

Phosphorylation of Tau R2 Repeat Destabilizes Its Binding to Microtubules: A Molecular Dynamics Simulation Study

Phosphorylation, the most popular post-translational modification of tau protein, plays an important role in regulating tau physiological functions. However, aberrant phosphorylation attenuates the binding affinity of tau to a microtubule (MT), resulting in MT destabilization followed by accumulation of neurofibrillary tangles in the brain. There are in total 85 potential phosphorylation sites in a full-length tau protein, and about half of them are abnormal as they occur in tau of Alzheimer's disease (AD) brain only. In this work, we investigated the impact of abnormal Ser289, Ser293, and Ser289/Ser293 phosphorylation on tau R2-MT binding and the conformation of tau R2 using molecular dynamics simulation. We found that the phosphorylation significantly affected R2-MT interaction and reduced the binding affinity of tau R2 peptides to MTs. Free energy decomposition analysis suggested that the post-translational modified residues themselves made a significant contribution to destabilize tau repeat R2-MT binding. Therefore, the phosphorylation may attenuate the binding affinity of tau to MTs. Additionally, the phosphorylation also enhanced helix-coil transition of monomeric R2 peptides, which may result in the acceleration of tau aggregation. Since these phosphorylated sites have not been examined in previous experimental studies, our finding through all-atom molecular dynamics simulations and free energy analysis can inspire experimental scientists to investigate the impact of the phosphorylation on MT binding and aggregation of full-length tau and the pathological roles of the phosphorylation at those sites in AD development through in vitro/in vivo assays.

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.

Sequence-Dependent Conformational Properties of PGGG Motif in Tau Repeats: Insights from Molecular Dynamics Simulations of Narrow Pick Filament

Tauopathies are a class of neurodegenerative diseases correlated with the presence of pathological Tau fibrils as a diagnostic marker. The microtubule-binding repeat region of Tau protein, which includes R1, R2, R3, and R4 repeats, constitutes the core of these fibrils. Each repeat consists of a semiconserved C-terminal hexapeptide flanked by KxGS and PGGG motifs. Previous studies have shown the influence of these peptides on protein aggregation, yet their repeat-specific properties are less explored. Using molecular dynamics, we probed the sequence-specific influence of the C-terminal hexapeptide (²⁶⁴ENLKHQ²⁶⁹) in determining the compact local conformation of the R1 repeat of the narrow Pick filament (NPF) with a homologous E264G mutation. In addition to that, we also studied the influence of ²⁶²S phosphorylation on this conformation as the phosphorylation is proposed to alleviate the pathogenesis of Pick's disease. Interestingly, we determined that E264G mutation induces a conformational shift of ²⁷⁰PGGG²⁷³ from a turn to a random coil. This conformational dependence is experimentally verified with the R1R3-E264G mutant construct, which displayed accelerated aggregation compared with the R1R3 wild-type construct. A significant delay in aggregation of the R1R3-G326E mutant further demonstrates the importance of ³²⁶G in determining the conformation of the R3 repeat. Thus, we conclude that the conformational properties of the PGGG motif in Tau repeats are strongly dependent on the repeat-specific sequence of the C-terminal hexapeptide.

Interactions of Polyphenolic Gallotannins with Amyloidogenic Polypeptides Associated with Alzheimer's Disease: From Molecular Insights to Physiological Significance

Polyphenols are natural compounds abundantly found in plants. They are known for their numerous benefits to human health, including antioxidant properties and anti-inflammatory activities. Interestingly, many studies have revealed that polyphenols can also modulate the formation of amyloid fibrils associated with disease states and can prevent the formation of cytotoxic oligomer species. In this review, we underline the numerous effects of four hydrolysable gallotannins (HGTs) with high conformational flexibility, low toxicity, and multi-targeticity, e.g., tannic acid, pentagalloyl glucose, corilagin, and 1,3,6-tri-O-galloyl-β-D-glucose, on the aggregation of amyloidogenic proteins associated with the Alzheimer's Disease (AD). These HGTs have demonstrated interesting abilities to reduce, at different levels, the formation of amyloid fibrils involved in AD, including those assembled from the amyloid β-peptide, the tubulin-associated unit, and the islet amyloid polypeptide. HGTs were also shown to disassemble pre-formed fibrils and to diminish cognitive decline in mice. Finally, this manuscript highlights the importance of further investigating these naturally occurring HGTs as promising scaffolds to design molecules that can interfere with the formation of proteotoxic oligomers and aggregates associated with AD pathogenesis.

Mechanistic insight into the disruption of Tau R3-R4 protofibrils by curcumin and epinephrine: an all-atom molecular dynamics study

The accumulation of Tau protein aggregates is a pathological hallmark of tauopathy, including chronic traumatic encephalopathy (CTE). Inhibiting Tau aggregation or disrupting preformed Tau fibrils is considered one of the rational therapeutic strategies to combat tauopathy. Previous studies reported that curcumin (Cur, a molecule of a labile natural product) and epinephrine (EP, an important neurotransmitter) could effectively inhibit the formation of Tau fibrillar aggregates and disassociate preformed fibrils. However, the underlying molecular mechanisms remain elusive. In this study, we performed multiple molecular dynamics simulations for 17.5 μs in total to investigate the influence of Cur and EP on the C-shaped Tau protofibril associated with CTE. Our simulations show that the protofibrillar pentamer is the smallest stable Tau R3-R4 protofibril. Taking the pentamer as a protofibril model, we found that both Cur and EP molecules could affect the shape of the Tau pentamer by changing the β2-β3 and β7-β8 angles, leading to a more extended structure. Cur and EP display a disruptive effect on the local β-sheets and the formation of hydrogen bonds, and thus destabilize the global protofibril structure. The contact number analysis shows that Cur has a higher binding affinity with the Tau pentamer than EP, especially in the nucleating segment PHF6. Hydrophobic, π-π and cation-π interactions together facilitate the binding of Cur and EP with the Tau pentamer. Cur exhibits stronger hydrophobic and π-π interactions with Tau than EP, and EP displays a stronger cation-π interaction. Our findings provide molecular insights into the disruptive mechanisms of the Tau R3-R4 protofibrils by curcumin and epinephrine, which may be useful for the design of effective drug candidates for the treatment of CTE.

Molecular Dynamics Simulations of the Tau Amyloid Fibril Core Dimer at the Surface of a Lipid Bilayer Model: I. In Alzheimer's Disease

A tau R3-R4 domain spanning residues 306-378 was shown to form an amyloid fibril core of a full-length tau in the brain of patients with Alzheimer's disease. Recently, we studied the dynamics of a tau R3-R4 monomer at the surface of a lipid bilayer model and revealed deep insertion of the amino acids spanning the PHF6 motif (residues 306-311) and its flanking residues. Here, we explore the membrane-associated conformational ensemble of a tau R3-R4 dimer by means of atomistic molecular dynamics. Similar to the monomer simulation, the R3-R4 dimer has the propensity to form β-hairpin-like conformation. Unlike the monomer, the dimer shows insertion of the C-terminal R4 region and transient adsorption of the PHF6 motif. Taken together, these results reveal the multiplicity of adsorption and insertion modes of tau into membranes depending on its oligomer size.

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.

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.

Dynamics of Amyloid Formation from Simplified Representation to Atomistic Simulations

Amyloid fibril formation is an intrinsic property of short peptides, non-disease proteins, and proteins associated with neurodegenerative diseases. Aggregates of the Aβ and tau proteins, the α-synuclein protein, and the prion protein are observed in the brain of Alzheimer's, Parkinson's, and prion disease patients, respectively. Due to the transient short-range and long-range interactions of all species and their high aggregation propensities, the conformational ensemble of these devastating proteins, the exception being for the monomeric prion protein, remains elusive by standard structural biology methods in bulk solution and in lipid membranes. To overcome these limitations, an increasing number of simulations using different sampling methods and protein models have been performed. In this chapter, we first review our main contributions to the field of amyloid protein simulations aimed at understanding the early aggregation steps of short linear amyloid peptides, the conformational ensemble of the Aβ40/42 dimers in bulk solution, and the stability of Aβ aggregates in lipid membrane models. Then we focus on our studies on the interactions of amyloid peptides/inhibitors to prevent aggregation, and long amyloid sequences, including new results on a monomeric tau construct.

Unraveling the Influence of K280 Acetylation on the Conformational Features of Tau Core Fragment: A Molecular Dynamics Simulation Study

Abnormal aggregation of the microtubule-associated protein Tau is closely associated with tauopathies, including Alzheimer's disease and chronic traumatic encephalopathy. The hexapeptide 275VQIINK280 (PHF6*), a fibril-nucleating core motif of Tau, has been shown to play a vital role in the aggregation of Tau. Mounting experiment evidence demonstrated the acetylation of a single-lysine residue K280 in the PHF6* was a critical event for the formation of pathological Tau amyloid deposits. However, the underlying mechanisms by which K280 acetylation affects Tau aggregation at the atomic level remain elusive. In this work, we performed replica exchange molecular dynamics simulations to investigate the influence of acetylation of K280 on the aggregation of PHF6*. Our simulations show that acetylation of K280 not only enhances the self-assembly capability of PHF6* peptides but also increases the β-sheet structure propensity of the PHF6*. The inter-molecular interactions among PHF6* peptides are strengthened by the acetylation of K280, resulting in an increased ordered β-sheet-rich conformations of the PHF6* assemblies along with a decrease of the structural diversity. The residue-pairwise contact frequency analysis shows that K280 acetylation increases the interactions among the hydrophobic chemical groups from PHF6* peptides, which promotes the aggregation of PHF6*. This study offers mechanistic insights into the effects of acetylation on the aggregation of PHF6*, which will be helpful for an in-depth understanding of the relationship between acetylation and Tau aggregation at the molecular level.

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.

Mechanisms of melatonin binding and destabilizing the protofilament and filament of tau R3-R4 domains revealed by molecular dynamics simulation

The accumulation of β-amyloid (Aβ) and tau protein is considered to be an important pathological characteristic of Alzheimer's disease (AD). Failure of medicine targeting Aβ has drawn more attention to the influence of tau protein and its fibrillization on neurodegeneration. Increasing evidence shows that melatonin (Mel) can effectively inhibit the formation of tau fibrils and disassemble preformed tau fibrils. However, the underlying mechanism is poorly understood. In this work, we investigated the kinetics of melatonin binding and destabilizing the tetrameric protofilament and octameric filament of tau R3-R4 domains by performing microsecond all-atom molecular dynamics simulations. Our results show that Mel is able to disrupt the C-shaped structure of the tau protofilament and filament, and destabilizes the association between N- and C-termini. Mel predominantly binds to β1 and β6-β8 regions and favors contact with the elongation surface, which is dominantly driven by hydrogen bonding interactions and facilitated by other interactions. The strong π-π stacking interaction of Mel with Y310 impedes the intramolecular CH-π interaction between I308 and Y310, and the cation-π interaction of Mel with R379 interferes with the formation of the D348-R379 salt bridge. Moreover, Mel occupies the protofilament surface in the tetrameric protofilament and prevents the formation of intermolecular hydrogen bonds between residues K331 and Q336 in the octameric filament. Our work provides molecular insights into Mel hindering tau fibrillization or destabilizing the protofilament and filament, and the revealed inhibitory mechanisms provide useful clues for the design of efficient anti-amyloid agents.

Protein Aggregation Landscape in Neurodegenerative Diseases: Clinical Relevance and Future Applications

Intrinsic disorder is a natural feature of polypeptide chains, resulting in the lack of a defined three-dimensional structure. Conformational changes in intrinsically disordered regions of a protein lead to unstable β-sheet enriched intermediates, which are stabilized by intermolecular interactions with other β-sheet enriched molecules, producing stable proteinaceous aggregates. Upon misfolding, several pathways may be undertaken depending on the composition of the amino acidic string and the surrounding environment, leading to different structures. Accumulating evidence is suggesting that the conformational state of a protein may initiate signalling pathways involved both in pathology and physiology. In this review, we will summarize the heterogeneity of structures that are produced from intrinsically disordered protein domains and highlight the routes that lead to the formation of physiological liquid droplets as well as pathogenic aggregates. The most common proteins found in aggregates in neurodegenerative diseases and their structural variability will be addressed. We will further evaluate the clinical relevance and future applications of the study of the structural heterogeneity of protein aggregates, which may aid the understanding of the phenotypic diversity observed in neurodegenerative disorders.

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.

Phosphorylation-Induced Structural Reorganization in Tau-Paired Helical Filaments

Taupathies involve the deposition of abnormal tau protein into neurofibrillary tangles (NFTs) in the human brain. The abnormally hyperphosphorylated tau dissociates from microtubules and forms insoluble aggregates known as paired helical filaments (PHFs), highlighting the importance of post-translational modifications in taupathies. The present study examines the factors responsible for the structural stability of PHFs in native as well as in phosphorylated and O-GlcNAcylated tau. We carried out molecular dynamics simulations on the R3-R4 repeat domains of the human tau protein to gain atomic insights into the key noncovalent interactions responsible for their unique dimeric C-shaped structure. The structural effects upon post-translational modification were found to be prominent for phosphorylation when compared with O-GlcNAcylation. O-GlcNAcylated tau was found to retain the "C conformation" observed in the native tau PHF, whereas upon phosphorylation, we observed a conformational transition to a more opened "H conformation". We found that this conformational transition is brought about by the loss of a key salt bridge between Lys353 and Asp358 due to the phosphorylation at Ser356 that results in the reorganization of the dimeric interface.

Insight Into Seeded Tau Fibril Growth From Molecular Dynamics Simulation of the Alzheimer's Disease Protofibril Core

Aggregates of the microtubule associated tau protein are a major constituent of neurofibrillary lesions that define Alzheimer’s disease (AD) pathology. Increasing experimental evidence suggests that the spread of tau neurofibrillary tangles results from a prion-like seeding mechanism in which small oligomeric tau fibrils template the conversion of native, intrinsically disordered, tau proteins into their pathological form. By using atomistic molecular dynamics (MD) simulations, we investigate the stability and dissociation thermodynamics of high-resolution cryo-electron microscopy (cryo-EM) structures of both the AD paired-helical filament (PHF) and straight filament (SF). Non-equilibrium steered MD (SMD) center-of-mass pulling simulations are used to probe the stability of the protofibril structure and identify intermolecular contacts that must be broken before a single tau peptide can dissociate from the protofibril end. Using a combination of exploratory metadynamics and umbrella sampling, we investigate the complete dissociation pathway and compute a free energy profile for the dissociation of a single tau peptide from the fibril end. Different features of the free energy surface between the PHF and SF protofibril result from a different mechanism of tau unfolding. Comparison of wild-type tau PHF and post-translationally modified pSer356 tau shows that phosphorylation at this site changes the dissociation free energy surface of the terminal peptide. These results demonstrate how different protofibril morphologies template the folding of endogenous tau in distinct ways, and how post-translational modification can perturb the folding mechanism.

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.

Molecular dynamics simulations reveal the destabilization mechanism of Alzheimer's disease-related tau R3-R4 Protofilament by norepinephrine

Aggregation of Tau protein into neurofibrillary tangles is associated with the pathogenesis of Alzheimer's disease (AD) which has no cure yet. Clearing neurofibrillary tangles is one of major therapeutic strategies. Experimental studies reported that norepinephrine (NE) has the ability to disrupt Tau filament and cause Tau degradation. However, the underlying mechanism remains elusive. Herein, we performed molecular dynamic simulations to investigate the influence of NE on the C-shaped Tau R3-R4 protofilament. Our simulations show that NE compound destabilizes Tau protofilament by mostly disrupting β6/β8 and altering the β2-β3 and β6-β7 angles. NE binds mainly with aromatic residues Y310/P312/H374/F378 through ππ stacking and charged residues E338/E342/D348/D358/E372 via hydrogen-bonding interactions. Our results, together with the findings that exercise can markedly increase NE level, suggest that exercise might be a potent therapy against AD. This study reveals the disruptive mechanism of Tau protofilament by NE molecules, which may provide new clues for AD drug candidate design.

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.

Structures of the intrinsically disordered Aβ, tau and α-synuclein proteins in aqueous solution from computer simulations

Intrinsically disordered proteins (IDPs) play many biological roles in the human proteome ranging from vesicular transport, signal transduction to neurodegenerative diseases. The Aβ and tau proteins, and the α-synuclein protein, key players in Alzheimer's and Parkinson's diseases, respectively are fully disordered at the monomer level. The structural heterogeneity of the monomeric and oligomeric states and the high self-assembly propensity of these three IDPs have precluded experimental structural determination. Simulations have been used to determine the atomic structures of these IDPs. In this article, we review recent computer models to capture the equilibrium ensemble of Aβ, tau and α-synuclein proteins at different association steps in aqueous solution and present new results of the PEP-FOLD framework on α-synuclein monomer.