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

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

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.

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.

pH-Induced Misfolding Mechanism of Prion Protein: Insights from Microsecond-Accelerated Molecular Dynamics Simulations.

The conformational transition of prion protein (PrP) from a native form PrPC to a pathological isoform PrPSc is the main cause of a number of prion diseases in human and animals. Thus, understanding the molecular basis of conformational transition of PrP will be valuable for unveiling the etiology of PrP-related diseases. Here, to explore the potential misfolding mechanism of PrP under the acidic condition, which is known to promote PrP misfolding and trigger its aggregation, the conventional and accelerated molecular dynamics (MD) simulations combined with the Markov state model (MSM) analysis were performed. The conventional MD simulations reveal that, at an acidic pH, the globular domain of PrP is partially unfolded, particularly for the α2 C-terminus. Structural analysis of the key macrostates obtained by MSM indicates that the α2 C-terminus and the ß2-α2 loop may serve as important sites for the pH-induced PrP misfolding. Meanwhile, the α1 may also participate in the pH-induced structural conversion by moving away from the α2-α3 subdomain. Notably, dynamical network analysis of the key metastable states indicates that the protonated H187 weakens the interactions between the α2 C-terminus, α1-ß2 loop, and α2-α3 loop, leading these domains, especially the α2 C-terminus, to become unstable and to begin to misfold. Therefore, the α2 C-terminus plays a key role in the PrP misfolding process and serves as a potential site for drug targeting. Overall, our findings can deepen the understanding of the pathogenesis related to PrP and provide useful guidance for the future drug discovery.

The folding mechanism and key metastable state identification of the PrP127-147 monomer studied by molecular dynamics simulations and Markov state model analysis.

The structural transition of prion proteins from a native α-helix (PrPC) to a misfolded ß-sheet-rich conformation (PrPSc) is believed to be the main cause of a number of prion diseases in humans and animals. Understanding the molecular basis of misfolding and aggregation of prion proteins will be valuable for unveiling the etiology of prion diseases. However, due to the limitation of conventional experimental techniques and the heterogeneous property of oligomers, little is known about the molecular architecture of misfolded PrPSc and the mechanism of structural transition from PrPC to PrPSc. The prion fragment 127-147 (PrP127-147) has been reported to be a critical region for PrPSc formation in Gerstmann-Straussler-Scheinker (GSS) syndrome and thus has been used as a model for the study of prion aggregation. In the present study, we employ molecular dynamics (MD) simulation techniques to study the conformational change of this fragment that could be relevant to the PrPC-PrPSc transition. Employing extensive replica exchange molecular dynamics (REMD) and conventional MD simulations, we sample a huge number of conformations of PrP127-147. Using the Markov state model (MSM), we identify the metastable conformational states of this fragment and the kinetic network of transitions between the states. The resulting MSM reveals that disordered random-coiled conformations are the dominant structures. A key metastable folded state with typical extended ß-sheet structures is identified with Pro137 being located in a turn region, consistent with a previous experimental report. Conformational analysis reveals that intrapeptide hydrophobic interaction and two key residue interactions, including Arg136-His140 and Pro137-His140, contribute a lot to the formation of ordered extended ß-sheet states. However, network pathway analysis from the most populated disordered state indicates that the formation of extended ß-sheet states is quite slow (at the millisecond level), as large structural rearrangement is needed from disordered states. We speculate that the formation process of the extended ß-sheet folded states may represent an important event during the early formation of prion oligomers and the results of our study provide insights into the molecular details of the early stage of prion aggregation.

Ligand induced change of ß2 adrenergic receptor from active to inactive conformation and its implication for the closed/open state of the water channel: insight from molecular dynamics simulation, free energy calculation and Markov state model analysis.

The reported crystal structures of ß2 adrenergic receptor (ß2AR) reveal that the open and closed states of the water channel are correlated with the inactive and active conformations of ß2AR. However, more details about the process by which the water channel states are affected by the active to inactive conformational change of ß2AR remain illusive. In this work, molecular dynamics simulations are performed to study the dynamical inactive and active conformational change of ß2AR induced by inverse agonist ICI 118,551. Markov state model analysis and free energy calculation are employed to explore the open and close states of the water channel. The simulation results show that inverse agonist ICI 118,551 can induce water channel opening during the conformational transition of ß2AR. Markov state model (MSM) analysis proves that the energy contour can be divided into seven states. States S1, S2 and S5, which represent the active conformation of ß2AR, show that the water channel is in the closed state, while states S4 and S6, which correspond to the intermediate state conformation of ß2AR, indicate the water channel opens gradually. State S7, which represents the inactive structure of ß2AR, corresponds to the full open state of the water channel. The opening mechanism of the water channel is involved in the ligand-induced conformational change of ß2AR. These results can provide useful information for understanding the opening mechanism of the water channel and will be useful for the rational design of potent inverse agonists of ß2AR.