Structural transitions and oligomerization along polyalanine fibril formation pathways from computer simulations

Advanced computational approaches to understand protein aggregation

Protein aggregation is a widespread phenomenon implicated in debilitating diseases like Alzheimer's, Parkinson's, and cataracts, presenting complex hurdles for the field of molecular biology. In this review, we explore the evolving realm of computational methods and bioinformatics tools that have revolutionized our comprehension of protein aggregation. Beginning with a discussion of the multifaceted challenges associated with understanding this process and emphasizing the critical need for precise predictive tools, we highlight how computational techniques have become indispensable for understanding protein aggregation. We focus on molecular simulations, notably molecular dynamics (MD) simulations, spanning from atomistic to coarse-grained levels, which have emerged as pivotal tools in unraveling the complex dynamics governing protein aggregation in diseases such as cataracts, Alzheimer's, and Parkinson's. MD simulations provide microscopic insights into protein interactions and the subtleties of aggregation pathways, with advanced techniques like replica exchange molecular dynamics, Metadynamics (MetaD), and umbrella sampling enhancing our understanding by probing intricate energy landscapes and transition states. We delve into specific applications of MD simulations, elucidating the chaperone mechanism underlying cataract formation using Markov state modeling and the intricate pathways and interactions driving the toxic aggregate formation in Alzheimer's and Parkinson's disease. Transitioning we highlight how computational techniques, including bioinformatics, sequence analysis, structural data, machine learning algorithms, and artificial intelligence have become indispensable for predicting protein aggregation propensity and locating aggregation-prone regions within protein sequences. Throughout our exploration, we underscore the symbiotic relationship between computational approaches and empirical data, which has paved the way for potential therapeutic strategies against protein aggregation-related diseases. In conclusion, this review offers a comprehensive overview of advanced computational methodologies and bioinformatics tools that have catalyzed breakthroughs in unraveling the molecular basis of protein aggregation, with significant implications for clinical interventions, standing at the intersection of computational biology and experimental research.

Theoretical and computational advances in protein misfolding

Misfolded proteins escape the cellular quality control mechanism and fail to fold properly or remain correctly folded leading to a loss in their functional specificity. Thus misfolding of proteins cause a large number of very different diseases ranging from errors in metabolism to various types of complex neurodegenerative diseases. A theoretical and computational perspective of protein misfolding is presented with a special emphasis on its salient features, mechanism and consequences. These insights quantitatively analyze different determinants of misfolding, that may be applied to design disease specific molecular targets.

Aggregation of Chameleon Peptides: Implications of α-Helicity in Fibril Formation

We investigate the relationship between the inherent secondary structure and aggregation propensity of peptides containing chameleon sequences (i.e., sequences that can adopt either α or β structure depending on context) using a combination of replica exchange molecular dynamics simulations, ion-mobility mass spectrometry, circular dichroism, and transmission electron microscopy. We focus on an eight-residue long chameleon sequence that can adopt an α-helical structure in the context of the iron-binding protein from Bacillus anthracis (PDB id 1JIG ) and a β-strand in the context of the baculovirus P35 protein (PDB id 1P35 ). We show that the isolated chameleon sequence is intrinsically disordered, interconverting between α-helical and β-rich conformations. The inherent conformational plasticity of the sequence can be constrained by addition of flanking residues with a given secondary structure propensity. Intriguingly, we show that the chameleon sequence with helical flanking residues aggregates rapidly into fibrils, whereas the chameleon sequence with flanking residues that favor β-conformations has weak aggregation propensity. This work sheds new insights into the possible role of α-helical intermediates in fibril formation.

Conformational behavior of polyalanine peptides with and without protecting groups of varying chain lengths: population of PP-II structure!

Oculopharyngeal muscular dystrophy (OPMD), a polyalanine myopathy, occurs due to expansion of homo-polyalanine stretch in normal polyadenylating binding protein nuclear 1 (PABPN1) protein from Ala10 to Ala11-17. Therefore, the conformational behavior of polyalanine peptides with n = 10-17, with and without terminal protecting groups, have been investigated with different starting geometries in water by molecular dynamics simulation studies. Alanine peptides are shown to give rise to unordered structure irrespective of starting geometry and not more than two residues at a stretch have the same/similar set of φ, ψ values. However, the final structure with terminal protecting groups look like β-strand. Unprotected poly-Ala peptides adopt twisted β-hairpin/multi hairpin like structure with increasing chain length. The number of residues having φ, ψ values in collagen region is found to be less in peptides with unprotected termini as compared to peptides with protected termini of same chain length. The results have been supported by recent synchrotron radiation circular dichroism spectroscopy of polyproline II and unordered secondary structures. Opening of the helical structure in poly-Ala peptides with protecting groups has been shown to take place from C-terminal and in peptides without protecting groups opening of helix starts from both terminals. Further, opening of helix takes more time in poly-Ala peptides without terminal protecting groups. The deviations in amide bond planarity have been discussed and compared with available experimental and computational results.

The effect of glycosaminoglycans (GAGs) on amyloid aggregation and toxicity

Amyloidosis is a protein folding disorder in which normally soluble proteins are deposited extracellularly as insoluble fibrils, impairing tissue structure and function. Charged polyelectrolytes such as glycosaminoglycans (GAGs) are frequently found associated with the proteinaceous deposits in tissues of patients affected by amyloid diseases. Experimental evidence indicate that they can play an active role in favoring amyloid fibril formation and stabilization. Binding of GAGs to amyloid fibrils occurs mainly through electrostatic interactions involving the negative polyelectrolyte charges and positively charged side chains residues of aggregating protein. Similarly to catalyst for reactions, GAGs favor aggregation, nucleation and amyloid fibril formation functioning as a structural templates for the self-assembly of highly cytotoxic oligomeric precursors, rich in β-sheets, into harmless amyloid fibrils. Moreover, the GAGs amyloid promoting activity can be facilitated through specific interactions via consensus binding sites between amyloid polypeptide and GAGs molecules. We review the effect of GAGs on amyloid deposition as well as proteins not strictly related to diseases. In addition, we consider the potential of the GAGs therapy in amyloidosis.

Competition between surface adsorption and folding of fibril-forming polypeptides

Self-assembly of polypeptides into fibrillar structures can be initiated by planar surfaces that interact favorably with certain residues. Using a coarse-grained model, we systematically studied the folding and adsorption behavior of a β-roll forming polypeptide. We find that there are two different folding pathways depending on the temperature: (i) at low temperature, the polypeptide folds in solution into a β-roll before adsorbing onto the attractive surface; (ii) at higher temperature, the polypeptide first adsorbs in a disordered state and folds while on the surface. The folding temperature increases with increasing attraction as the folded β-roll is stabilized by the surface. Surprisingly, further increasing the attraction lowers the folding temperature again, as strong attraction also stabilizes the adsorbed disordered state, which competes with folding of the polypeptide. Our results suggest that to enhance the folding, one should use a weakly attractive surface. They also explain the recent experimental observation of the nonmonotonic effect of charge on the fibril formation on an oppositely charged surface [C. Charbonneau et al., ACS Nano 8, 2328 (2014)].

Involvement of non-polyalanine (polyA) residues in aggregation of polyA proteins: Clue for inhibition of aggregation

Presence of polyalanine (polyA) stretches in some proteins is found to be associated with their aggregation, which causes disorders in various developmental processes. In this work, inherent propensities towards aggregation of some residues, which are not part of the polyA stretches, have been identified by using the primary sequences of seven polyA proteins with the help of Betascan, PASTA and Tango programs and explored unambiguously. This provides a basis for proposing molecular mechanism of this type of aggregation. Reported suppression of aggregation of polyA proteins by chaperones like HSP40 and HSP70 is substantiated through molecular docking. The hydrophobic residues of identified aggregating region are found to be interacting with hydrophobic surface of chaperones. This suggests a crucial clue for possible way to inhibit the aggregation of such proteins.

A coarse-grained model for peptide aggregation on a membrane surface

The aggregation of peptides on a lipid bilayer is studied using coarse-grained molecular dynamics in implicit solvent. Peptides bind to and self-assemble on the membrane surface into β-rich fibrillar aggregates, even under conditions where only disordered oligomers form in bulk solution. Relative to a solid surface, the membrane surface facilitates peptide mobility and a more complex network of morphology transitions as aggregation proceeds. Additionally, final aggregate structures realized on the membrane surface are distinct from those observed on a comparable solid surface. The aggregated fibrils alter the local structure and material properties of the lipid bilayer in their immediate vicinity but have only a modest effect on the overall bending rigidity of the bilayer.

Conformational transitions of the cross-linking domains of elastin during self-assembly

Elastin is the intrinsically disordered polymeric protein imparting the exceptional properties of extension and elastic recoil to the extracellular matrix of most vertebrates. The monomeric precursor of elastin, tropoelastin, as well as polypeptides containing smaller subsets of the tropoelastin sequence, can self-assemble through a colloidal phase separation process called coacervation. Present understanding suggests that self-assembly is promoted by association of hydrophobic domains contained within the tropoelastin sequence, whereas polymerization is achieved by covalent joining of lysine side chains within distinct alanine-rich, α-helical cross-linking domains. In this study, model elastin polypeptides were used to determine the structure of cross-linking domains during the assembly process and the effect of sequence alterations in these domains on assembly and structure. CD temperature melts indicated that partial α-helical structure in cross-linking domains at lower temperatures was absent at physiological temperature. Solid-state NMR demonstrated that β-strand structure of the cross-linking domains dominated in the coacervate state, although α-helix was predominant after subsequent cross-linking of lysine side chains with genipin. Mutation of lysine residues to hydrophobic amino acids, tyrosine or alanine, leads to increased propensity for β-structure and the formation of amyloid-like fibrils, characterized by thioflavin-T binding and transmission electron microscopy. These findings indicate that cross-linking domains are structurally labile during assembly, adapting to changes in their environment and aggregated state. Furthermore, the sequence of cross-linking domains has a dramatic effect on self-assembly properties of elastin-like polypeptides, and the presence of lysine residues in these domains may serve to prevent inappropriate ordered aggregation.

Interplay between folding and assembly of fibril-forming polypeptides

Polypeptides can self-assemble into hierarchically organized fibrils consisting of a stack of individually folded polypeptides driven together by hydrophobic interaction. Using a coarse-grained model, we systematically studied this self-assembly as a function of temperature and hydrophobicity of the residues on the outside of the building block. We find the self-assembly can occur via two different pathways-a random aggregation-folding route and a templated-folding process-thus indicating a strong coupling between folding and assembly. The simulation results can explain experimental evidence that assembly through stacking of folded building blocks is rarely observed, at the experimental concentrations. The model thus provides a generic picture of hierarchical fibril formation.