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

Recent Computational Advances Regarding Amyloid-β and Tau Membrane Interactions in Alzheimer's Disease

The interactions of amyloid proteins with membranes have been subject to many experimental and computational studies, as these interactions contribute in part to neurodegenerative diseases. In this review, we report on recent simulations that have focused on the adsorption and insertion modes of amyloid-β and tau proteins in membranes. The atomistic-resolution characterization of the conformational changes of these amyloid proteins upon lipid cell membrane and free lipid interactions is of interest to rationally design drugs targeting transient oligomers in Alzheimer's disease.

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.

Interaction of the tau fibrils with the neuronal membrane

Tau proteins are recently gaining a lot of interest due to their active role in causing a range of tauopathies. Molecular mechanisms underlying the tau interaction with the neuronal membrane are hitherto unknown and difficult to characterize using experimental methods. Using the cryo-EM structure of the tau-fibrils we have used atomistic molecular dynamics simulation to model the tau fibril and neuronal membrane interaction using explicit solvation. The dynamics and structural characteristics of the tau fibril with the neuronal membrane are compared to the tau fibril in the aqueous phase to corroborate the effect of the neuronal membrane in the tau structure. Tau fibrils have been modelled using CHARMM-36m force field and the six component neuronal membrane composition is taken from the earlier simulation results. The timescale conceivable in our molecular dynamics simulations is of the order of microseconds which captures the onset of the interaction of the tau fibrils with the neuronal membrane. This interaction is found to impact the tau pathogenesis that finally causes neuronal toxicity. Our study initiates the understanding of tau conformational ensemble in the presence of neuronal membrane and sheds the light on the significant tau-membrane interactions.

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.

Recent Advances in Molecular Dynamics Simulations of Tau Fibrils and Oligomers

The study of tau protein aggregation and interactions with other molecules or solvents using molecular dynamics simulations (MDs) is of interest to many researchers to propose new mechanism-based therapeutics for neurodegenerative diseases such as Alzheimer's disease, Pick's disease, chronic traumatic encephalopathy, and other tauopathies. In this review, we present recent MD simulation studies of tau oligomers and fibrils such as tau-NPK, tau-PHF, tau-K18, and tau-R3-R4 monomers and dimers. All-atom simulations by replica exchange MDs and coarse-grained MDs in lipid bilayers and in solution were used. The simulations revealed different mechanisms in the binding of tau in bilayers and in solutions, depending on the peptide size. Phosphorylation is also an important factor in MD simulations. The use of steered MDs was also included to simulate the dissociation of tau fibrils. The exponential improvement in the computing power of computers has led to an increasing number of scientists and engineers using a cost-effective, high-performance computing platform to study how the tau protein interacts and the effects of changing its structure, such as the phosphorylation of tau fibrils.

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.

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.

Self-Assembly of Amyloid-Beta (Aβ) Peptides from Solution to Near In Vivo Conditions

Understanding the atomistic resolution changes during the self-assembly of amyloid peptides or proteins is important to develop compounds or conditions to alter the aggregation pathways and suppress the toxicity and potentially aid in the development of drugs. However, the complexity of protein aggregation and the transient order/disorder of oligomers along the pathways to fibril are very challenging. In this Perspective, we discuss computational studies of amyloid polypeptides carried out under various conditions, including conditions closely mimicking in vivo and point out the challenges in obtaining physiologically relevant results, focusing mainly on the amyloid-beta Aβ peptides.

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.