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

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.

Transient interdomain interactions modulate the monomeric structural ensemble and oligomerization landscape of Huntingtin Exon 1

Polyglutamine expansion (≥ 36 residues) within the N-terminal exon-1 of Huntingtin (Httex1) leads to Huntington's disease, a neurogenerative condition marked by the presence of intranuclear Htt inclusions. Notably, the polyglutamine tract in Httex1 is flanked by an N-terminal coiled-coil domain - N17 (17 amino acids), which undergoes self-association to promote the formation of soluble Httex1 oligomers and brings the aggregation-prone polyQ tracts in close spatial proximity. However, the mechanisms underlying the subsequent conversion of soluble oligomers into insoluble β-rich aggregates with increasing polyQ length, remain unclear. Current knowledge suggests that expansion of the polyQ tract increases its helicity, and this favors its oligomerization and aggregation. In addition, studies utilizing conformation-specific antibodies and a stable coiled-coil heterotetrametric system fused to polyQ indicate that domain "cross-talk" (i.e., interdomain interactions) may be necessary to efficiently promote the emergence of toxic conformations (in monomers and oligomers) and fibrillar aggregation. Here, we performed extensive atomistic molecular dynamics (MD) simulations (aggregate time ∼ 0.7 ms) of N17-polyQ fragments to uncover the interplay between structural transformation and domain "cross-talk" on the monomeric structural ensemble and oligomerization landscape of Httex1. Our simulation ensembles of N17-polyQ monomers validated against 13 C NMR chemical shifts indicated that in addition to elevated α-helicity, polyQ expansion also favors transient, interdomain (N17-polyQ) interactions which result in the emergence of β-conformations. Further, interdomain interactions decreased the overall stability of N17-mediated dimers by counteracting the stabilizing effect of increased α-helicity and promoted a heterogenous oligomerization landscape on the sub-microsecond timescale. Overall, our study uncovers the significance of domain "cross-talk" in modulating the monomeric conformational ensemble and oligomerization landscape of Httex1 to favor the formation of amyloid aggregates.

Neuropathogenesis-on-chips for neurodegenerative diseases

Developing diagnostics and treatments for neurodegenerative diseases (NDs) is challenging due to multifactorial pathogenesis that progresses gradually. Advanced in vitro systems that recapitulate patient-like pathophysiology are emerging as alternatives to conventional animal-based models. In this review, we explore the interconnected pathogenic features of different types of ND, discuss the general strategy to modelling NDs using a microfluidic chip, and introduce the organoid-on-a-chip as the next advanced relevant model. Lastly, we overview how these models are being applied in academic and industrial drug development. The integration of microfluidic chips, stem cells, and biotechnological devices promises to provide valuable insights for biomedical research and developing diagnostic and therapeutic solutions for NDs.

Repurposing Amphotericin B: anti-microbial, molecular docking and molecular dynamics simulation studies suggest inhibition potential of Amphotericin B against MRSA

Amphotericin B (AMPH) is an anti-fungal drug and this study, for the first time as best of our knowledge, reports the repurposing of the Amphotericin B. The drug was found to show significant antibacterial potential revealed by antimicrobial screening, molecular docking, and mode of action analysis targeting Penicillin Binding Protein 2a (PBP 2a protein) which is target of β-lactam drugs and is involved in cell wall synthesis. Mode of action analysis showed the drug to have hydrophobic and hydrophilic interactions with both C-terminal, trans-peptidase and non-penicillin binding domain of the protein. Additionally, to evaluate the impact of ligand binding on the protein's conformational dynamics, molecular dynamics (MD) simulations were used. Comparative Dynamical flexibility (RMSF) and Dynamics Cross Correlation (DCCM) followed by MD simulations revealed the complex formation significantly effecting structural dynamics of the enzyme significantly in the non-penicillin binding domain (327-668) and slightly in trans peptidase domain. Radius of gyration assessment further showed ligand binding also decreasing over all compactness of protein. Secondary structure analysis indicated the complex formation changing the conformational integrity in non-penicillin binding domain. Hydrogen bond analysis and MMPBSA, free energy of calculations followed by MD simulations, also complemented the antimicrobial and molecular docking revelations suggesting Amphotericin B to have substantial antibacterial potential.

The structure of pathogenic huntingtin exon 1 defines the bases of its aggregation propensity

Huntington's disease is a neurodegenerative disorder caused by a CAG expansion in the first exon of the HTT gene, resulting in an extended polyglutamine (poly-Q) tract in huntingtin (httex1). The structural changes occurring to the poly-Q when increasing its length remain poorly understood due to its intrinsic flexibility and the strong compositional bias. The systematic application of site-specific isotopic labeling has enabled residue-specific NMR investigations of the poly-Q tract of pathogenic httex1 variants with 46 and 66 consecutive glutamines. Integrative data analysis reveals that the poly-Q tract adopts long α-helical conformations propagated and stabilized by glutamine side chain to backbone hydrogen bonds. We show that α-helical stability is a stronger signature in defining aggregation kinetics and the structure of the resulting fibrils than the number of glutamines. Our observations provide a structural perspective of the pathogenicity of expanded httex1 and pave the way to a deeper understanding of poly-Q-related diseases.

Structural systems biology approach delineate the functional implications of SNPs in exon junction complex interaction network

In eukaryotes, transcripts that carry premature termination codons (PTC) leading to truncated proteins are degraded by the Nonsense Mediated Decay (NMD) machinery. Missense and nonsense Single Nucleotide Polymorphisms (SNPs) in proteins belonging to Exon junction complex (EJC) and up-frameshift protein (UPF) will compromise NMD leading to the accumulation of truncated proteins in various diseases. The EJC and UPF which are involved in NMD is a good model system to study the effect of SNPs at a system level. Despite the availability of crystal structures, computational tools, and data on mutational and deletion studies, with functional implications, an integrated effort to understand the impact of SNPs at the systems level is lacking. To study the functional consequences of missense SNPs, sequence-based techniques like SIFT and PolyPhen which classify SNPs as deleterious or non-deleterious and structure-based methods like FoldX which calculate the Delta Delta G, (ddGs, ∆∆G) are used. Using FoldX, the ddG for mutations with experimentally validated functional effects is calculated and compared with those calculated for SNPs in the same protein-protein interaction interface. Further, a model is conceived to explain the functional implications of SNPs based on the effects observed for known mutants. The results are visualized in a network format. The effects of nonsense mutations are discerned by comparing with deletion mutation studies and loss of interaction in the crystal structure. The present work not only integrates genomics, proteomics, and classical genetics with 'Structural Biology' but also helps to integrate it into a 'systems-level functional network'.Communicated by Ramaswamy H. Sarma.

Comparative molecular dynamics simulations of pathogenic and non-pathogenic huntingtin protein monomers and dimers

Polyglutamine expansion at the N-terminus of the huntingtin protein exon 1 (Htt-ex1) is closely associated with a number of neurodegenerative diseases, which result from the aggregation of the increased polyQ repeat. However, the underlying structures and aggregation mechanism are still poorly understood. We performed microsecond-long all-atom molecular dynamics simulations to study the folding and dimerization of Htt-ex1 (about 100 residues) with non-pathogenic and pathogenic polyQ lengths, and uncovered substantial differences. The non-pathogenic monomer adopts a long α-helix that includes most of the polyQ residues, which forms the interaction interface for dimerization, and a PPII-turn-PPII motif in the proline-rich region. In the pathogenic monomer, the polyQ region is disordered, leading to compact structures with many intra-protein interactions and the formation of short β-sheets. Dimerization can proceed via different modes, where those involving the N-terminal headpiece bury more hydrophobic residues and are thus more stable. Moreover, in the pathogenic Htt-ex1 dimers the proline-rich region interacts with the polyQ region, which slows the formation of β-sheets.

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.

Possible Co-Evolution of Polyglutamine and Polyproline in Huntingtin Protein: Proline-Rich Domain as Transient Folding Chaperone

Huntington's disease is an inherited neurodegenerative disorder caused by the overduplication of CAG repeats in the Huntingtin gene. Recent findings revealed that among the orthologs, the expansion of CAG repeats (polyQ) in the Huntingtin gene occurs in tandem with the duplication of CCG repeats (polyP). However, the molecular mechanism of this possible co-evolution remains unknown. We examined the structures of Huntingtin exon 1 (HttEx1) from six species along with five designed mutants. We found that the polyP segments "chaperone" the rest of the HttEx1 by forming ad hoc polyP binding grooves. Such a process elongates the otherwise poorly solvated polyQ domain, while modulating its secondary structure propensity from β-strands to α-helices. This chaperoning effect is achieved mostly through transient hydrogen bond interactions between polyP and the rest of HttEx1, resulting in a striking golden ratio of ∼2:1 between the chain lengths of polyQ and polyP.

Synthesis and Evaluation of a Fluorine-18 Radioligand for Imaging Huntingtin Aggregates by Positron Emission Tomographic Imaging

Mutations in the huntingtin gene (HTT) triggers aggregation of huntingtin protein (mHTT), which is the hallmark pathology of neurodegenerative Huntington's disease (HD). Development of a high affinity 18F radiotracer would enable the study of Huntington's disease pathology using a non-invasive imaging modality, positron emission tomography (PET) imaging. Herein, we report the first synthesis of fluorine-18 imaging agent, 6-(5-((5-(2,2-difluoro-2-(fluoro-18F)ethoxy)pyridin-2-yl)methoxy)benzo[d]oxazol-2-yl)-2-methylpyridazin-3(2H)-one ([18F]1), a radioligand for HD and its preclinical evaluation in vitro (autoradiography of post-mortem HD brains) and in vivo (rodent and non-human primate brain PET). [18F]1 was synthesized in a 4.1% RCY (decay corrected) and in an average molar activity of 16.5 ± 12.5 GBq/μmol (445 ± 339 Ci/mmol). [18F]1 penetrated the blood-brain barrier of both rodents and primates, and specific saturable binding in post-mortem brain slices was observed that correlated to mHTT aggregates identified by immunohistochemistry.

miR-302 Attenuates Mutant Huntingtin-Induced Cytotoxicity through Restoration of Autophagy and Insulin Sensitivity

Huntington’s disease (HD) is an autosomal-dominant brain disorder caused by mutant huntingtin (mHtt). Although the detailed mechanisms remain unclear, the mutational expansion of polyglutamine in mHtt is proposed to induce protein aggregates and neuronal toxicity. Previous studies have shown that the decreased insulin sensitivity is closely related to mHtt-associated impairments in HD patients. However, how mHtt interferes with insulin signaling in neurons is still unknown. In the present study, we used a HD cell model to demonstrate that the miR-302 cluster, an embryonic stem cell-specific polycistronic miRNA, is significantly downregulated in mHtt-Q74-overexpressing neuronal cells. On the contrary, restoration of miR-302 cluster was shown to attenuate mHtt-induced cytotoxicity by improving insulin sensitivity, leading to a reduction of mHtt aggregates through the enhancement of autophagy. In addition, miR-302 also promoted mitophagy and stimulated Sirt1/AMPK-PGC1α pathway thereby preserving mitochondrial function. Taken together, these results highlight the potential role of miR-302 cluster in neuronal cells, and provide a novel mechanism for mHtt-impaired insulin signaling in the pathogenesis of HD.

Specific Key-Point Mutations along the Helical Conformation of Huntingtin-Exon 1 Protein Might Have an Antagonistic Effect on the Toxic Helical Content's Formation

The polyglutamine tract length represents a key regulator for the Huntington's disease toxicity level and its aggregation rates, often being related to helical structural conformations. In this study, we performed all-atom MD simulations on mutant Huntingtin-Exon1 protein with additional mutation spots, aiming to observe the corresponding structural and dynamical changes at the level of the helix. The simulated structures consist of three sets of Q residue mutations into P residues (4P, 7P, and 9P), with each set including different spots of mutations: random along the mutant sequence (R models), at the edges of the helix (E models), as well as at the edges and in the middle of the helix (EM models). At the helical level, our results predict less compactness profiles for a higher number of P mutations (7P and 9P models) with particular mutation spots at the edges and at the edges-middle of the helix. Moreover, the C-alpha atom distances decreased for 7P and 9P models in comparison to 4P models, and the RMSF values show the highest fluctuation rates for 9P models with point mutations at the edges and in the middle of the helix. The secondary structure analysis suggests greater structural transitions from α-helices to bends, turns, and random coils for 7P and 9P models, particularly for point mutations considered at the edges and in the middle of the helical content. The obtained results support our hypothesis that specific key-point mutations along the helical conformation might have an antagonistic effect on the toxic helical content's formation.