The Possible Structural Models for Polyglutamine Aggregation: A Molecular Dynamics Simulations Study

Reviewing the Structure-Function Paradigm in Polyglutamine Disorders: A Synergistic Perspective on Theoretical and Experimental Approaches

Polyglutamine (polyQ) disorders are a group of neurodegenerative diseases characterized by the excessive expansion of CAG (cytosine, adenine, guanine) repeats within host proteins. The quest to unravel the complex diseases mechanism has led researchers to adopt both theoretical and experimental methods, each offering unique insights into the underlying pathogenesis. This review emphasizes the significance of combining multiple approaches in the study of polyQ disorders, focusing on the structure-function correlations and the relevance of polyQ-related protein dynamics in neurodegeneration. By integrating computational/theoretical predictions with experimental observations, one can establish robust structure-function correlations, aiding in the identification of key molecular targets for therapeutic interventions. PolyQ proteins' dynamics, influenced by their length and interactions with other molecular partners, play a pivotal role in the polyQ-related pathogenic cascade. Moreover, conformational dynamics of polyQ proteins can trigger aggregation, leading to toxic assembles that hinder proper cellular homeostasis. Understanding these intricacies offers new avenues for therapeutic strategies by fine-tuning polyQ kinetics, in order to prevent and control disease progression. Last but not least, this review highlights the importance of integrating multidisciplinary efforts to advancing research in this field, bringing us closer to the ultimate goal of finding effective treatments against polyQ disorders.

Advances in the understanding of protein misfolding and aggregation through molecular dynamics simulation

Aberrant protein folding known as protein misfolding is counted as one of the striking factors of neurodegenerative diseases. The extensive range of pathologies caused by protein misfolding, aggregation and subsequent accumulation are mainly classified into either gain of function diseases or loss of function diseases. In order to seek for novel strategies for treatment and diagnosis of neurodegenerative diseases, insights into the mechanism of misfolding and aggregation is essential. A comprehensive knowledge on the factors influencing misfolding and aggregation is required as well. An extensive experimental study on protein aggregation is somewhat challenging due to the insoluble and noncrystalline nature of amyloid fibrils. Thus there has been a growing use of computational approaches including Monte Carlo simulation, docking simulation, molecular dynamics simulation in the study of protein misfolding and aggregation. The review presents a discussion on molecular dynamics simulation alone as to how it has emerged as a promising tool in the understanding of protein misfolding and aggregation in general, detailing upon three different aspects considering four misfold prone proteins in particular. It is noticeable that all four proteins considered in this review i.e prion, superoxide dismutase1, huntingtin and amyloid β are linked to chronic neurodegenerative diseases with debilitating effects. Initially the review elaborates on the factors influencing the misfolding and aggregation. Next, it addresses our current understanding of the amyloid structures and the associated aggregation mechanisms, finally, summarizing the contribution of this computational tool in the search for therapeutic strategies against the respective protein-deposition diseases.

Molecular Dynamics Simulations Applied to Structural and Dynamical Transitions of the Huntingtin Protein: A Review

Over the recent years, Huntington's disease (HD) has become widely discussed in the scientific literature especially because at the mutant level there are several contradictions regarding the aggregation mechanism. The specific role of the physiological huntingtin protein remains unknown, due to the lack of characterization of its entire crystallographic structure, making the experimental and theoretical research even harder when taking into consideration its involvement in multiple biological functions and its high affinity for different interacting partners. Different types of models, containing fewer (not more than 35 Qs) polyglutamine residues for the WT structure and above 35 Qs for the mutants, were subjected to classical or advanced MD simulations to establish the proteins' structural stability by evaluating their conformational changes. Outside the polyQ tract, there are two other regions of interest (the N17 domain and the polyP rich domain) considered to be essential for the aggregation kinetics at the mutant level. The polymerization process is considered to be dependent on the polyQ length. As the polyQ tract's dimension increases, the structures present more β-sheet conformations. Contrarily, it is also considered that the aggregation stability is not necessarily dependent on the number of Qs, while the initial stage of the aggregation seed might play the decisive role. A general assumption regarding the polyP domain is that it might preserve the polyQ structures soluble by acting as an antagonist for β-sheet formation.

Molecular dynamics analysis of the aggregation propensity of polyglutamine segments

Protein misfolding and aggregation is a pathogenic feature shared among at least ten polyglutamine (polyQ) neurodegenerative diseases. While solvent-solution interaction is a key factor driving protein folding and aggregation, the solvation properties of expanded polyQ tracts are not well understood. By using GPU-enabled all-atom molecular dynamics simulations of polyQ monomers in an explicit solvent environment, this study shows that solvent-polyQ interaction propensity decreases as the lengths of polyQ tract increases. This study finds a predominance in long-distance interactions between residues far apart in polyQ sequences with longer polyQ segments, that leads to significant conformational differences. This study also indicates that large loops, comprised of parallel β-structures, appear in long polyQ tracts and present new aggregation building blocks with aggregation driven by long-distance intra-polyQ interactions. Finally, consistent with previous observations using coarse-grain simulations, this study demonstrates that there is a gain in the aggregation propensity with increased polyQ length, and that this gain is correlated with decreasing ability of solvent-polyQ interaction. These results suggest the modulation of solvent-polyQ interactions as a possible therapeutic strategy for treating polyQ diseases.

Homology modeling and docking studies of FabH (β-ketoacyl-ACP synthase III) enzyme involved in type II fatty acid biosynthesis of Chlorella variabilis: a potential algal feedstock for biofuel production

The concept of using microalgae as an alternative renewable source of biofuel has gained much importance in recent years. However, its commercial feasibility is still an area of concern for researchers. Unraveling the fatty acid metabolic pathway and understanding structural features of various key enzymes regulating the process will provide valuable insights to target microalgae for augmented oil content. FabH (β-ketoacyl-acyl carrier protein synthase; KAS III) is a condensing enzyme catalyzing the initial elongation step of type II fatty acid biosynthetic process and acyl carrier protein (ACP) facilitates the shuttling of the fatty acyl intermediates to the active site of the respective enzymes in the pathway. In the present study, a reliable three-dimensional structure of FabH from Chlorella variabilis, an oleaginous green microalga was modeled and subsequently the key residues involved in substrate binding were determined by employing protein-protein docking and molecular dynamics (MD) simulation protocols. The FabH-ACP complex having the lowest docking energy score showed the binding of ACP to the electropositive FabH surface with strong hydrogen bond interactions. The MD simulation results indicated that the substrate-complexed FabH adopted a more stable conformation than the free enzyme. Further, the FabH structure retained its stability throughout the simulation although noticeable displacements were observed in the loop regions. Molecular simulation studies suggested the importance of crucial hydrogen bonding of the conserved Arg(91) of FabH with Glu(53) and Asp(56) of ACP for exhibiting high affinity between the enzyme and substrate. The molecular modeling results are consistent with available experimental results on the flexibility of FabH and the present study provides first in silico insights into the structural and dynamical aspect of catalytic mechanism of FabH, which could be used for further site-specific mutagenic experiments to develop engineered high oil-yielding microalgal strains for biofuel production.

Longitudinal changes of outcome measures in spinal and bulbar muscular atrophy

Spinal and bulbar muscular atrophy is an adult-onset, hereditary motor neuron disease caused by the expansion of a trinucleotide CAG repeat within the gene encoding the androgen receptor. To date, several agents have been shown to prevent or slow disease progression in animal models of this disease. For the translational research of these agents, it is necessary to perform the detailed analysis of natural history with quantitative outcome measures and to establish sensitive and validated disease-specific endpoints in the clinical trials. To this end, we performed a prospective observation of disease progression over 3 years in 34 genetically confirmed Japanese patients with spinal and bulbar muscular atrophy by using quantitative outcome measures, including functional and blood parameters. The baseline evaluation revealed that CAG repeat length in the androgen receptor gene correlated not only with the age of onset but also with the timing of substantial changes in activity of daily living. Multiple regression analyses indicated that the serum level of creatinine is the most useful blood parameter that reflects the severity of motor dysfunction in spinal and bulbar muscular atrophy. In 3-year prospective analyses, a slow but steady progression was affirmed in most of the outcome measures we examined. In the analyses using random coefficient models that summarize the individual data into a representative line, disease progression was not affected by CAG repeat length or onset age. These models showed large interindividual variation, which was also independent of the differences of CAG repeat size. Analyses using these models also demonstrated that the subtle neurological deficits at an early or preclinical stage were more likely to be detected by objective motor functional tests such as the 6-min walk test and grip power or serum creatinine levels than by functional rating scales, such as the revised amyotrophic lateral sclerosis functional rating scale or modified Norris scale. Categorization of the clinical phenotypes using factor analysis showed that upper limb function is closely related to bulbar function, but not to lower limb function at baseline, whereas the site of onset had no substantial effects on disease progression. These results suggest that patients with spinal and bulbar muscular atrophy show a slow but steady progression of motor dysfunction over time that is independent of CAG repeat length or clinical phenotype, and that objective outcome measures may be used to evaluate disease severity at an early stage of this disease.

The effect of fucoidan on tyrosinase: computational molecular dynamics integrating inhibition kinetics

Fucoidan is a complex sulfated polysaccharide extracted from brown seaweed and has a wide variety of biological activities. In this study, we investigated the inhibitory effect of fucoidan on tyrosinase via a combination of inhibition kinetics and computational simulations. Fucoidan reversibly inhibited tyrosinase in a mixed-type manner. Time-interval kinetics showed that the inhibition was processed as first order with biphasic processes. For further insight, we simulated dockings with various sizes of molecular models (monomer to decamer) of fucoidan and showed that the best binding energy change results were obtained from the pentamer (-1.89 kcal/mol) and the hexamer (-1.97 kcal/mol) models of AutoDock Vina. The molecular dynamics simulation confirmed the binding mechanisms between tyrosinase and fucoidan and suggested that fucoidan mostly interacts with several residues including copper ions located in the active site. Our study suggests that fucoidan might be a potential natural antipigment agent.

Sequence-dependent stability test of a left-handed β-helix motif

The left-handed β-helix (LHBH) is an intriguing, rare structural pattern in polypeptides that has been implicated in the formation of amyloid aggregates. We used accurate all-atom replica-exchange molecular dynamics (REMD) simulations to study the relative stability of diverse sequences in the LHBH conformation. Ensemble-average coordinates from REMD served as a scoring criterion to identify sequences and threadings optimally suited to the LHBH, as in a fold recognition paradigm. We examined the repeatability of our REMD simulations, finding that single simulations can be reliable to a quantifiable extent. We find expected behavior for the positive and negative control cases of a native LHBH and intrinsically disordered sequences, respectively. Polyglutamine and a designed hexapeptide repeat show remarkable affinity for the LHBH motif. A structural model for misfolded murine prion protein was also considered, and showed intermediate stability under the given conditions. Our technique is found to be an effective probe of LHBH stability, and promises to be scalable to broader studies of this and potentially other novel or rare motifs. The superstable character of the designed hexapeptide repeat suggests theoretical and experimental follow-ups.

Physical chemistry of polyglutamine: intriguing tales of a monotonous sequence

Polyglutamine (polyQ) sequences of unknown normal function are present in a significant number of proteins, and their repeat expansion is associated with a number of genetic neurodegenerative diseases. PolyQ solution structure and properties are important not only because of the normal and abnormal biology associated with these sequences but also because they represent an interesting case of a biologically relevant homopolymer. As the common thread in expanded polyQ repeat diseases, it is important to understand the structure and properties of simple polyQ sequences. At the same time, experience has shown that sequences attached to polyQ, whether in artificial constructs or in disease proteins, can influence structure and properties. The two major contenders for the molecular source of the neurotoxicity implicit in polyQ expansion within disease proteins are a populated toxic conformation in the monomer ensemble and a toxic aggregated species. This review summarizes experimental and computational studies on the solution structure and aggregation properties of both simple and complex polyQ sequences, and their repeat-length dependence. As a representative of complex polyQ proteins, the behavior of huntingtin N-terminal fragments, such as exon-1, receives special attention.