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.

The role of α-sheet structure in amyloidogenesis: characterization and implications

Amyloid diseases are linked to protein misfolding whereby the amyloidogenic protein undergoes a conformational change, aggregates and eventually forms amyloid fibrils. While the amyloid fibrils and plaques are hallmarks of these diseases, they typically form late in the disease process and do not correlate with disease. Instead, there is growing evidence that smaller, soluble toxic oligomers form prior and appear to be early triggers of the molecular pathology underlying these diseases. Nearly 20 years ago, we proposed the α-sheet hypothesis after discovering that the early conformational changes observed during atomistic molecular dynamics simulations involve the formation of a non-standard protein structure, α-sheet. Furthermore, we proposed that toxic oligomers contain α-sheet structure and that preferentially targeting this structure could neutralize the toxicity, prevent further aggregation and serve as the basis for early detection of disease. Here, we present the origin of the α-sheet hypothesis and describe α-sheet structure and the corresponding mechanisms of conversion. We discuss experimental studies demonstrating that both mammalian and bacterial amyloid systems form α-sheet oligomers before converting to conventional β-sheet fibrils. Furthermore, we show that the process can be inhibited with de novo designed α-sheet peptides complementary to the structure in the toxic oligomers.

Novel Stilbene-Nitroxyl Hybrid Compounds Display Discrete Modulation of Amyloid Beta Toxicity and Structure

Several neurodegenerative diseases are driven by misfolded proteins that assemble into soluble aggregates. These "toxic oligomers" have been associated with a plethora of cellular dysfunction and dysregulation, however the structural features underlying their toxicity are poorly understood. A major impediment to answering this question relates to the heterogeneous nature of the oligomers, both in terms of structural disorder and oligomer size. This not only complicates elucidating the molecular etiology of these disorders, but also the druggability of these targets as well. We have synthesized a class of bifunctional stilbenes to modulate both the conformational toxicity within amyloid beta oligomers (AβO) and the oxidative stress elicited by AβO. Using a neuronal culture model, we demonstrate this bifunctional approach has the potential to counter the molecular pathogenesis of Alzheimer's disease in a powerful, synergistic manner. Examination of AβO structure by various biophysical tools shows that each stilbene candidate uniquely alters AβO conformation and toxicity, providing insight towards the future development of structural correctors for AβO. Correlations of AβO structural modulation and bioactivity displayed by each provides insights for future testing in vivo. The multi-target activity of these hybrid molecules represents a highly advantageous feature for disease modification in Alzheimer's, which displays a complex, multifactorial etiology. Importantly, these novel small molecules intervene with intraneuronal AβO, a necessary feature to counter the cycle of dysregulation, oxidative stress and inflammation triggered during the earliest stages of disease progression.

Determination of amino acids that favour the αL region using Ramachandran propensity plots. Implications for α-sheet as the possible amyloid intermediate

In amyloid diseases an insoluble amyloid fibril forms via a soluble oligomeric intermediate. It is this intermediate that mediates toxicity and it has been suggested, somewhat controversially, that it has the α-sheet structure. Nests and α-strands are similar peptide motifs in that alternate residues lie in the α_R and γ_L regions of the Ramachandran plot for nests, or α_R and α_L regions for α-strands. In nests a concavity is formed by the main chain NH atoms whereas in α-strands the main chain is almost straight. Using "Ramachandran propensity plots" to focus on the α_L/γ_L region, it is shown that glycine favours γ_L (82% of amino acids are glycine), but disfavours α_L (3% are glycine). Most charged and polar amino acids favour α_L with asparagine having by far the highest propensity. Thus, glycine favours nests but, contrary to common expectation, should not favour α-sheet. By contrast most charged or polar amino acids should favour α-sheet by their propensity for the α_L conformation, which is more discriminating amongst amino acids than the α_R conformation. Thus, these results suggest the composition of sequences that favour α-sheet formation and point towards effective prediction of α-sheet from sequence.

Tumorigenic p53 mutants undergo common structural disruptions including conversion to α-sheet structure

The p53 protein is a commonly studied cancer target because of its role in tumor suppression. Unfortunately, it is susceptible to mutation-associated loss of function; approximately 50% of cancers are associated with mutations to p53, the majority of which are located in the central DNA-binding domain. Here, we report molecular dynamics simulations of wild-type (WT) p53 and 20 different mutants, including a stabilized pseudo-WT mutant. Our findings indicate that p53 mutants tend to exacerbate latent structural-disruption tendencies, or vulnerabilities, already present in the WT protein, suggesting that it may be possible to develop cancer therapies by targeting a relatively small set of structural-disruption motifs rather than a multitude of effects specific to each mutant. In addition, α-sheet secondary structure formed in almost all of the proteins. α-Sheet has been hypothesized and recently demonstrated to play a role in amyloidogenesis, and its presence in the reported p53 simulations coincides with the recent re-consideration of cancer as an amyloid disease.

Local structural motifs in proteins: Detection and characterization of fragments inserted in helices

The global/local fold of protein structures is stabilized by a variety of specific interactions. A primary role in this context is played by hydrogen bonds. In order to identify novel motifs in proteins, we searched Protein Data Bank structures looking for backbone H-bonds formed by NH groups of two (or more) consecutive residues with consecutive CO groups of distant residues in the sequence. The present analysis unravels the occurrence of recurrent structural motifs that, to the best of our knowledge, had not been characterized in literature. Indeed, these H-bonding patterns are found (i) in a specific parallel β-sheet capping, (ii) in linking of β-hairpins to α-helices, and (iii) in α-helix insertions. Interestingly, structural analyses of these motifs indicate that Gly residues frequently occupy prominent positions. The formation of these motifs is likely favored by the limited propensity of Gly to be embodied in helices/sheets. Of particular interest is the motif corresponding to insertions in helices that was detected in 1% of analyzed structures. Inserted fragments may assume different structures and aminoacid compositions and usually display diversified evolutionary conservation. Since inserted regions are physically separated from the rest of the protein structure, they represent hot spots for ad-hoc protein functionalization.

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.

Insights from molecular dynamics simulations for computational protein design

A grand challenge in the field of structural biology is to design and engineer proteins that exhibit targeted functions. Although much success on this front has been achieved, design success rates remain low, an ever-present reminder of our limited understanding of the relationship between amino acid sequences and the structures they adopt. In addition to experimental techniques and rational design strategies, computational methods have been employed to aid in the design and engineering of proteins. Molecular dynamics (MD) is one such method that simulates the motions of proteins according to classical dynamics. Here, we review how insights into protein dynamics derived from MD simulations have influenced the design of proteins. One of the greatest strengths of MD is its capacity to reveal information beyond what is available in the static structures deposited in the Protein Data Bank. In this regard simulations can be used to directly guide protein design by providing atomistic details of the dynamic molecular interactions contributing to protein stability and function. MD simulations can also be used as a virtual screening tool to rank, select, identify, and assess potential designs. MD is uniquely poised to inform protein design efforts where the application requires realistic models of protein dynamics and atomic level descriptions of the relationship between dynamics and function. Here, we review cases where MD simulations was used to modulate protein stability and protein function by providing information regarding the conformation(s), conformational transitions, interactions, and dynamics that govern stability and function. In addition, we discuss cases where conformations from protein folding/unfolding simulations have been exploited for protein design, yielding novel outcomes that could not be obtained from static structures.

The Ramachandran Number: An Order Parameter for Protein Geometry

Three-dimensional protein structures usually contain regions of local order, called secondary structure, such as α-helices and β-sheets. Secondary structure is characterized by the local rotational state of the protein backbone, quantified by two dihedral angles called ϕ and ψ. Particular types of secondary structure can generally be described by a single (diffuse) location on a two-dimensional plot drawn in the space of the angles ϕ and ψ, called a Ramachandran plot. By contrast, a recently-discovered nanomaterial made from peptoids, structural isomers of peptides, displays a secondary-structure motif corresponding to two regions on the Ramachandran plot [Mannige et al., Nature 526, 415 (2015)]. In order to describe such ‘higher-order’ secondary structure in a compact way we introduce here a means of describing regions on the Ramachandran plot in terms of a single Ramachandran number, R, which is a structurally meaningful combination of ϕ and ψ. We show that the potential applications of R are numerous: it can be used to describe the geometric content of protein structures, and can be used to draw diagrams that reveal, at a glance, the frequency of occurrence of regular secondary structures and disordered regions in large protein datasets. We propose that R might be used as an order parameter for protein geometry for a wide range of applications.

Amyloid Properties of Asparagine and Glutamine in Prion-like Proteins

Sequences rich in glutamine (Q) and asparagine (N) are intrinsically disordered in monomeric form, but can aggregate into highly ordered amyloids, as seen in Q/N-rich prion domains (PrDs). Amyloids are fibrillar protein aggregates rich in β-sheet structures that can self-propagate through protein-conformational chain reactions. Here, we present a comprehensive theoretical study of N/Q-rich peptides, including sequences found in the yeast Sup35 PrD, in parallel and antiparallel β-sheet aggregates, and probe via fully atomistic molecular dynamics simulations all their possible steric-zipper interfaces in order to determine their protofibril structure and their relative stability. Our results show that polyglutamine aggregates are more stable than polyasparagine aggregates. Enthalpic contributions to the free energy favor the formation of polyQ protofibrils, while entropic contributions favor the formation of polyN protofibrils. The considerably larger phase space that disordered polyQ must sample on its way to aggregation probably is at the root of the associated slower kinetics observed experimentally. When other amino acids are present, such as in the Sup35 PrD, their shorter side chains favor steric-zipper formation for N but not Q, as they preclude the in-register association of the long Q side chains.

Solid-State Nuclear Magnetic Resonance on the Static and Dynamic Domains of Huntingtin Exon-1 Fibrils

Amyloid-like fibrils formed by huntingtin exon-1 (htt_ex1) are a hallmark of Huntington's disease (HD). The structure of these fibrils is unknown, and determining their structure is an important step toward understanding the misfolding processes that cause HD. In HD, a polyglutamine (polyQ) domain in htt_ex1 is expanded to a degree that it gains the ability to form aggregates comprising the core of the resulting fibrils. Despite the simplicity of this polyQ sequence, the structure of htt_ex1 fibrils has been difficult to determine. This study provides a detailed structural investigation of fibrils formed by htt_ex1 using solid-state nuclear magnetic resonance (NMR) spectroscopy. We show that the polyQ domain of htt_ex1 forms the static amyloid core similar to polyQ model peptides. The Gln residues of this domain exist in two distinct conformations that are found in separate domains or monomers but are relatively close in space. The rest of htt_ex1 is relatively dynamic on an NMR time scale, especially the proline-rich C-terminus, which we found to be in a polyproline II helical and random coil conformation. We observed a similar dynamic C-terminus in a soluble form of htt_ex1, indicating that the conformation of this part of htt_ex1 is not changed upon its aggregation into an amyloid fibril. From these data, we propose a bottlebrush model for the fibrils formed by htt_ex1. In this model, the polyQ domains form the center and the proline-rich domains the bristles of the bottlebrush.

Structural determinants of polyglutamine protofibrils and crystallites

Nine inherited neurodegenerative diseases are associated with the expansion of the CAG codon. Once the translated polyglutamine expansion becomes longer than ~36 residues, it triggers the formation of intraneural protein aggregates that often display the signature of cross-β amyloid fibrils. Here, we use fully atomistic molecular dynamics simulations to probe the structural stability and conformational dynamics of both previously proposed and new polyglutamine aggregate models. We test the relative stability of parallel and antiparallel β sheets, and characterize possible steric interfaces between neighboring sheets and the effects of different alignments of the side-chain carboxamide dipoles. Results indicate that (i) different initial oligomer structures converge to crystals consistent with available diffraction data, after undergoing cooperative side-chain rotational transitions and quarter-stagger displacements on a microsecond time scale, (ii) structures previously deemed stable on a hundred nanosecond time scale are unstable over the microsecond time scale, and (iii) conversely, structures previously deemed unstable did not account for the correct side-chain packing and once the correct symmetry is considered the structures become stable for over a microsecond, due to tightly interdigitated side chains, which lock into highly regular polar zippers with inter-side-chain and backbone-side-chain hydrogen bonds. With these insights, we built Q40 monomeric models with different combinations of arc and hairpin turns and tested them for stability. The stable monomers were further probed as a function of repeat length. Our results are consistent with the aggregation threshold. These results explain and reconcile previously reported experimental and model discrepancies about polyglutamine aggregate structures.

Are long-range structural correlations behind the aggregration phenomena of polyglutamine diseases?

We have characterized the conformational ensembles of polyglutamine peptides of various lengths (ranging from to ), both with and without the presence of a C-terminal polyproline hexapeptide. For this, we used state-of-the-art molecular dynamics simulations combined with a novel statistical analysis to characterize the various properties of the backbone dihedral angles and secondary structural motifs of the glutamine residues. For (i.e., just above the pathological length for Huntington's disease), the equilibrium conformations of the monomer consist primarily of disordered, compact structures with non-negligible -helical and turn content. We also observed a relatively small population of extended structures suitable for forming aggregates including - and -strands, and - and -hairpins. Most importantly, for we find that there exists a long-range correlation (ranging for at least residues) among the backbone dihedral angles of the Q residues. For polyglutamine peptides below the pathological length, the population of the extended strands and hairpins is considerably smaller, and the correlations are short-range (at most residues apart). Adding a C-terminal hexaproline to suppresses both the population of these rare motifs and the long-range correlation of the dihedral angles. We argue that the long-range correlation of the polyglutamine homopeptide, along with the presence of these rare motifs, could be responsible for its aggregation phenomena.

Nine neurodegenerative diseases are caused by polyglutamine (polyQ) expansions greater than a given threshold in proteins with little or no homology except for the polyQ regions. The diseases all share a common feature: the formation of polyQ aggregates and eventual neuronal death. Using molecular dynamics simulations, we have explored the conformations of polyQ peptides. Results indicate that for peptides (i.e., just above the pathological length for Hungtington's disease), the equilibrium conformations were found to consist primarily of disordered, compact structures with a non-negligible -helical and turn content. We also observed a small population of extended structures suitable for forming aggregates. For peptides below the pathological length, the population of these structures was found to be considerably lower. For longer peptides, we found evidence for long-range correlations among the dihedral angles. This correlation turns out to be short-range for the smaller polyQ peptides, and is suppressed (along with the extended structural motifs) when a C-terminal polyproline tail is added to the peptides. We believe that the existence of these long-range correlations in above-threshold polyQ peptides, along with the presence of rare motifs, could be responsible for the experimentally observed aggregation phenomena associated with polyQ diseases.

Energetic selection of topology in ferredoxins

Models of early protein evolution posit the existence of short peptides that bound metals and ions and served as transporters, membranes or catalysts. The Cys-X-X-Cys-X-X-Cys heptapeptide located within bacterial ferredoxins, enclosing an Fe₄S₄ metal center, is an attractive candidate for such an early peptide. Ferredoxins are ancient proteins and the simple α+β fold is found alone or as a domain in larger proteins throughout all three kingdoms of life. Previous analyses of the heptapeptide conformation in experimentally determined ferredoxin structures revealed a pervasive right-handed topology, despite the fact that the Fe₄S₄ cluster is achiral. Conformational enumeration of a model CGGCGGC heptapeptide bound to a cubane iron-sulfur cluster indicates both left-handed and right-handed folds could exist and have comparable stabilities. However, only the natural ferredoxin topology provides a significant network of backbone-to-cluster hydrogen bonds that would stabilize the metal-peptide complex. The optimal peptide configuration (alternating α_L,α_R) is that of an α-sheet, providing an additional mechanism where oligomerization could stabilize the peptide and facilitate iron-sulfur cluster binding.

The ferredoxin fold is one of the oldest structures capable of catalyzing electron transfer reactions. In nature, only a right-handed topology exists in the ferredoxin fold. To understand how a specific fold-handedness was selected, we analyzed the structural motif using the tools of de novo protein design, searching in an unbiased fashion for backbone geometries that can favorably interact with the tetrahedral iron-sulfur cluster. In silico, we found both left-handed and right-handed folds can be formed, however the right-handed folds provide up to six hydrogen bonds that can stabilize the reduced iron-sulfur cluster, whereas left-handed folds at most form three hydrogen bonds. The difference in electrostatic conformational energy may have influenced selection of topology early in the evolution of iron-sulfur cluster containing proteins. This observation led us to establish a fundamental protein design principle that only right-handed peptide folds can properly interact while maintain redox function. Our results provide guidance in the creation of artificial proteins capable of carrying out redox reactions.

PPII propensity of multiple-guest amino acids in a proline-rich environment

There has been considerable debate about the intrinsic PPII propensity of amino acid residues in denatured polypeptides. Experimentally, this scale is based on the behavior of guest amino acid residues placed in the middle of proline-based hosts. We have used classical molecular dynamics simulations combined with replica-exchange methods to carry out a comprehensive analysis of the conformational equilibria of proline-based host oligopeptides with multiple guest amino acids including alanine, glutamine, valine, and asparagine. The tracked structural characteristics include the secondary structural motifs based on the Ramachandran angles and the cis/trans isomerization of the prolyl bonds. In agreement with our recent study of single amino acid guests, we did not observe an intrinsic PPII propensity in any of the guest amino acids in a multiple-guest setting. Instead, the experimental results can be explained in terms of (i) the steric restrictions imposed on the C-terminal guest amino acid that is immediately followed by a proline residue and (ii) an increase in the trans content of the prolyl bonds due to the presence of guest residues. In terms of the latter, we found that the more guests added to the system, the larger the increase in the trans content of the prolyl bonds, which results in an effective increase in the PPII content of the peptide.