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.

Effect of the R126C mutation on the structure and function of the glucose transporter GLUT1: A molecular dynamics simulation study

Glucose transporter 1 (GLUT1) is responsible for basal glucose uptake and is expressed in most tissues under normal conditions. GLUT1 mutations can cause early-onset absence epilepsy and myoclonus dystonia syndrome (MDS), with MDS potentially lethal. In this study, the effect of the R126C mutation, which is associated with MDS, on structural stability and substrate transport of GLUT1 was investigated. Various bioinformatics tools were used to predict the stability of GLUT1, revealing that the R126C mutation reduces the structural stability of GLUT1. Molecular dynamics (MD) simulations were used to further characterize the effect of the R126C mutation on GLUT1 structural stability. Based on the MD simulations, specific conformational changes and dominant motions of the GLUT1 mutant were characterized by Principal component analysis (PCA). The mutation disrupts hydrogen bonds between substrate-binding residues and glucose, thus likely reducing substrate affinity. The R126C mutation reduces the conformational stability of the protein, and fewer intramolecular hydrogen bonds were present in the mutated GLUT1 when compared with that of wild-type GLUT1. The mutation increased the free energy of glucose transport through GLUT1 significantly, especially at the mutation site, indicating that passage of glucose through the channel is hindered, and this mutant may even release cytoplasmic glucose. This study provides a detailed atomic-level explanation for the reduced structural stability and substrate transport capacity of a GLUT1 mutant. The results aid our understanding of the structure of GLUT1 and provide a framework for developing drugs to treat GLUT1-related diseases, such as MDS.

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 ¹⁷¹IPAKTPPAPK¹⁸⁰, 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.

Unraveling the Histidine Tautomerism Effect on the Initial Stages of Prion Misfolding: New Insights from a Computational Perspective

The aggregation and structural conversion of normal prion peptide (PrP^C) into the pathogenic scrapie form (PrP^Sc), which can act as a seed to enhance prion amyloid fiber formation, is believed to be a crucial event in prionopathies. Previous research suggests that the prion monomer may play an important role in oligomer generation during disease pathogenesis. In the present study, extensive replica-exchange molecular dynamics (REMD) simulations were conducted to explore the conformational characteristics of the huPrP (125-160) monomer under the histidine tautomerism effect. Investigating the structural characteristics and fibrilization process is challenging because two histidine tautomers [N_ε2-H (ε) and N_δ1-H (δ)] can occur in the open neutral state. Molecular dynamics (MD) simulation outcomes have shown that the toxic εδ and δδ isomer (containing several and broader local minima) had the highest α-helix structures, with contents of 21.11% and 21.01%, respectively, and may have a strong influence on the organizational behavior of a monomeric prion. The amino acids aspartate 20 (D20)-asparagine 29 (N29) and isoleucine 15 (I15)-histidine 16 (H16), D20-arginine 27 (R27) as well as N29 formed α-helix with the highest probabilities in the δδ and εδ isomer, accordingly. On the basis of our findings, we propose the histidine tautomerization hypothesis as a new prion accumulation mechanism, which may exist to induce the formation of prion accumulates. Overall, our tautomerism hypothesis constitutes a promising perspective for enhancing understanding of prion disease pathobiology and may help in the design of a good inhibitor.

Dimeric allostery mechanism of the plant circadian clock photoreceptor ZEITLUPE

In Arabidopsis thaliana, the Light-Oxygen-Voltage (LOV) domain containing protein ZEITLUPE (ZTL) integrates light quality, intensity, and duration into regulation of the circadian clock. Recent structural and biochemical studies of ZTL indicate that the protein diverges from other members of the LOV superfamily in its allosteric mechanism, and that the divergent allosteric mechanism hinges upon conservation of two signaling residues G46 and V48 that alter dynamic motions of a Gln residue implicated in signal transduction in all LOV proteins. Here, we delineate the allosteric mechanism of ZTL via an integrated computational approach that employs atomistic simulations of wild type and allosteric variants of ZTL in the functional dark and light states, together with Markov state and supervised machine learning classification models. This approach has unveiled key factors of the ZTL allosteric mechanisms, and identified specific interactions and residues implicated in functional allosteric changes. The final results reveal atomic level insights into allosteric mechanisms of ZTL function that operate via a non-trivial combination of population-shift and dynamics-driven allosteric pathways.

UMAP as a Dimensionality Reduction Tool for Molecular Dynamics Simulations of Biomacromolecules: A Comparison Study

Proteins are the molecular machines of life. The multitude of possible conformations that proteins can adopt determines their free-energy landscapes. However, the inherently high dimensionality of a protein free-energy landscape poses a challenge to deciphering how proteins perform their functions. For this reason, dimensionality reduction is an active field of research for molecular biologists. The uniform manifold approximation and projection (UMAP) is a dimensionality reduction method based on a fuzzy topological analysis of data. In the present study, the performance of UMAP is compared with that of other popular dimensionality reduction methods such as t-distributed stochastic neighbor embedding (t-SNE), principal component analysis (PCA), and time-structure independent components analysis (tICA) in the context of analyzing molecular dynamics simulations of the circadian clock protein VIVID. A good dimensionality reduction method should accurately represent the data structure on the projected components. The comparison of the raw high-dimensional data with the projections obtained using different dimensionality reduction methods based on various metrics showed that UMAP has superior performance when compared with linear reduction methods (PCA and tICA) and has competitive performance and scalable computational cost.

Tunneling nanotubes: A novel pharmacological target for neurodegenerative diseases?

Diversiform ways of intercellular communication are vital links in maintaining homeostasis and disseminating physiological states. Among intercellular bridges, tunneling nanotubes (TNTs) discovered in 2004 were recognized as potential pharmacology targets related to the pathogenesis of common or infrequent neurodegenerative disorders. The neurotoxic aggregates in neurodegenerative diseases including scrapie prion protein (PrPSc), mutant tau protein, amyloid-beta (Aβ) protein, alpha-synuclein (α-syn) as well as mutant Huntington (mHTT) protein could promote TNT formation via certain physiological mechanisms, in turn, mediating the intercellular transmission of neurotoxicity. In this review, we described in detail the skeleton, the formation, the physicochemical properties, and the functions of TNTs, while paying particular attention to the key role of TNTs in the transport of pathological proteins during neurodegeneration.

The metastable states of proteins

The intriguing process of protein folding comprises discrete steps that stabilize the protein molecules in different conformations. The metastable state of protein is represented by specific conformational characteristics, which place the protein in a local free energy minimum state of the energy landscape. The native-to-metastable structural transitions are governed by transient or long-lived thermodynamic and kinetic fluctuations of the intrinsic interactions of the protein molecules. Depiction of the structural and functional properties of metastable proteins is not only required to understand the complexity of folding patterns but also to comprehend the mechanisms of anomalous aggregation of different proteins. In this article, we review the properties of metastable proteins in context of their stability and capability of undergoing atypical aggregation in physiological conditions.

How do mutations affect the structural characteristics and substrate binding of CYP21A2? An investigation by molecular dynamics simulations

CYP21A2 mutations affect the activity of the protein leading to CAH disease.

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 ³⁰⁶VQIVYK³¹¹ (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.

Improving the productivity of monodisperse polyhedral cages by the rational design of kinetic self-assembly pathways

Hollow polyhedral cages hold great potential for application in nanotechnological and biomedical fields. Understanding the formation mechanism of these self-assembled structures could provide guidance for the rational design of the desired polyhedral cages. Here, by constructing kinetic network models from extensive coarse-grained molecular dynamics simulations, we elucidated the formation mechanism of the dodecahedral cage, which is formed by the self-assembly of patchy particles. We found that the dodecahedral cage is formed through increasing the aggregate size followed by structure rearrangement. Based on this mechanistic understanding, we improved the productivity of the dodecahedral cage through the rational design of the patch arrangement of patchy particles, which promotes the structural rearrangement process. Our results demonstrate that it should be a feasible strategy to achieve the rational design of the desired nanostructures via the kinetic analysis. We anticipate that this methodology could be extended to other self-assembly systems for the fabrication of functional nanomaterials.

Experimental and Theoretical Insights into the Inhibition Mechanism of Prion Fibrillation by Resveratrol and its Derivatives

Resveratrol and its derivatives have been shown to display beneficial effects to neurodegenerative diseases. However, the molecular mechanism of resveratrol and its derivatives on prion conformational conversion is poorly understood. In this work, the interaction mechanism between prion and resveratrol as well as its derivatives was investigated using steady-state fluorescence quenching, Thioflavin T binding assay, Western blotting, and molecular dynamics simulation. Protein fluorescence quenching method and Thioflavin T assay revealed that resveratrol and its derivatives could interact with prion and interrupt prion fibril formation. Molecular dynamics simulation results indicated that resveratrol can stabilize the PrP_127-147 peptide mainly through π-π stacking interactions between resveratrol and Tyr128. The hydrogen bonds interactions between resveratrol and the PrP_127-147 peptide could further reduce the flexibility and the propensity to aggregate. The results of this study not only can provide useful information about the interaction mechanism between resveratrol and prion, but also can provide useful clues for further design of new inhibitors inhibiting prion aggregation.