Length-dependent aggregation of uninterrupted polyalanine peptides

Dynamics and Machine Learning Reveal the Link between Tripeptide Sequences and Evaporation-Driven Material Properties

Previous research showed that a peptide composed of three tyrosines (YYY) can turn into organic glass and cause strong adhesion between substrates via evaporation. However, the mechanisms of these processes remain unclear, and the exploration of applications of other peptide sequences is necessary. In this study, an optimized evaporation method was employed in molecular dynamics. It was found that YYY evaporation products possess abundant internal hydrogen bonds, which may facilitate the amorphous glass state formation. Moderate hydrophilicity of a peptide enhances molecular mobility and self-healing ability, while excessive hydrophilicity causes a water plasticizing effect. Stronger hydrophilicity also brings a larger curvature of evaporation products on polydimethylsiloxane (PDMS) substrate. A machine learning model was developed to predict the evaporation contact angle of peptide evaporation products and agrees well with the experiment. This research aims to improve understanding of peptide structure-function relationships and aid in designing custom organic optical devices based on peptide sequences.

Liquid-liquid crystalline phase separation of spider silk proteins

Liquid-liquid phase separation (LLPS) of proteins can be considered an intermediate solubility regime between disperse solutions and solid fibers. While LLPS has been described for several pathogenic amyloids, recent evidence suggests that it is similarly relevant for functional amyloids. Here, we review the evidence that links spider silk proteins (spidroins) and LLPS and its role in the spinning process. Major ampullate spidroins undergo LLPS mediated by stickers and spacers in their repeat regions. During spinning, the spidroins droplets shift from liquid to crystalline states. Shear force, altered ion composition, and pH changes cause micelle-like spidroin assemblies to form an increasingly ordered liquid-crystalline phase. Interactions between polyalanine regions in the repeat regions ultimately yield the characteristic β-crystalline structure of mature dragline silk fibers. Based on these findings, we hypothesize that liquid-liquid crystalline phase separation (LLCPS) can describe the molecular and macroscopic features of the phase transitions of major ampullate spidroins during spinning and speculate whether other silk types may use a similar mechanism to convert from liquid dope to solid fiber.

AlphaFold with conformational sampling reveals the structural landscape of homorepeats

Homorepeats are motifs with reiterations of the same amino acid. They are prevalent in proteins associated with diverse physiological functions but also linked to several pathologies. Structural characterization of homorepeats has remained largely elusive, primarily because they generally occur in the disordered regions or proteins. Here, we address this subject by combining structures derived from machine learning with conformational sampling through physics-based simulations. We find that hydrophobic homorepeats have a tendency to fold into structured secondary conformations, while hydrophilic ones predominantly exist in unstructured states. Our data show that the flexibility rendered by disorder is a critical component besides the chemical feature that drives homorepeats composition toward hydrophilicity. The formation of regular secondary structures also influences their solubility, as pathologically relevant homorepeats display a direct correlation between repeat expansion, induction of helicity, and self-assembly. Our study provides critical insights into the conformational landscape of protein homorepeats and their structure-activity relationship.

An Undergraduate Laboratory Module Integrating Organic Chemistry and Polymer Science

Polymer science is receiving wider acceptance in the organic chemistry community; thus, it is imperative to include it in the undergraduate organic chemistry curriculum. Despite the ever-increasing popularity of the topic of polymer chemistry in undergraduate curricula, a comprehensive laboratory experiment module describing a polypeptide synthesis by ring-opening polymerization of N-carboxyanhydride (NCA ROP) and a homopolymer synthesis by activators-regenerated by electron-transfer for atom transfer radical polymerization (ARGET ATRP) has yet to be proposed. Herein, we report a semester-long, ten week undergraduate laboratory module focusing on the synthesis and analytical characterization of polyalanine and polystyrene for an advanced organic chemistry class. Students received hands-on-experiences in synthesizing polymers followed by their characterization via proton nuclear magnetic resonance (1H NMR) spectroscopy, electrospray ionization-mass spectrometry (ESI-MS), gel permeation chromatography (GPC), thermogravimetry (TGA), differential scanning calorimetry (DSC), and transmission electron microscopy (TEM), which are not well-presented in the typical organic chemistry curricula. These engaging hands-on lessons in the newly designed laboratory module not only increase students' interests in an interdisciplinary environment of organic chemistry and polymer science but also cultivate their research interests and communication skills and promote critical thinking.

Backbone Hydration of α-Helical Peptides: Hydrogen-Bonding and Surface Hydrophobicity/Hydrophilicity

The stability of proteins and small peptides depends on the way they interact with the surrounding water molecules. For small peptides, such as α-helical polyalanine (polyALA), water molecules can weaken the intramolecular hydrogen-bonds (HB) formed between the peptide backbone O and NH groups which are responsible for the α-helix structure. Here, we perform molecular dynamics simulations to study the hydration of polyALA, polyserine (polySER), and other homopolymer peptide α-helices at different temperatures and pressures. We find that water molecules form HB with most polyALA carbonyl O atoms, despite ALA hydrophobic CH3 side chain. Similar water-peptide backbone HB are found in other (hydrophobic and hydrophilic) homopolymer α-helices with large side chains, including polyvaline, polyleucine, and polyphenyalanine. A novel hydration mechanism is observed in polyserine (polySER): the backbone peptide rarely forms HB with water and, instead, the carbonyl O atoms tend to form HB with polySER side chain OH groups. We also quantify the hydrophobicity/hydrophilicity of polyALA and polySER by calculating the contact angle θ c of a water droplet pierced by a long polyALA/polySER α-helix. Unexpectedly, even when polyALA α-helix is supposed to be hydrophobic (θ c > 90°), we find that θ c ≈ 79°. For polySER, θ c ≈ 70°, consistent with α-helical polySER being hydrophilic.

Alternative low-populated conformations prompt phase transitions in polyalanine repeat expansions

Abnormal trinucleotide repeat expansions alter protein conformation causing malfunction and contribute to a significant number of incurable human diseases. Scarce structural insights available on disease-related homorepeat expansions hinder the design of effective therapeutics. Here, we present the dynamic structure of human PHOX2B C-terminal fragment, which contains the longest polyalanine segment known in mammals. The major α-helical conformation of the polyalanine tract is solely extended by polyalanine expansions in PHOX2B, which are responsible for most congenital central hypoventilation syndrome cases. However, polyalanine expansions in PHOX2B additionally promote nascent homorepeat conformations that trigger length-dependent phase transitions into solid condensates that capture wild-type PHOX2B. Remarkably, HSP70 and HSP90 chaperones specifically seize PHOX2B alternative conformations preventing phase transitions. The precise observation of emerging polymorphs in expanded PHOX2B postulates unbalanced phase transitions as distinct pathophysiological mechanisms in homorepeat expansion diseases, paving the way towards the search of therapeutics modulating biomolecular condensates in central hypoventilation syndrome.

Dissociation process of polyalanine aggregates by free electron laser irradiation

Polyalanine (polyA) disease-causative proteins with an expansion of alanine repeats can be aggregated. Although curative treatments for polyA diseases have not been explored, the dissociation of polyA aggregates likely reduces the cytotoxicity of polyA. Mid-infrared free electron laser (FEL) successfully dissociated multiple aggregates. However, whether the FEL dissociates polyA aggregates like other aggregates has not been tested. Here, we show that FEL at 6.1 μm experimentally weakened the extent of aggregation of a peptide with 13 alanine repeats (13A), and the irradiated 13A exerted lesser cytotoxicity to neuron-like cells than non-irradiated 13A. Then, we applied molecular dynamics (MD) simulation to follow the dissociation process by FEL. We successfully observed how the intermolecular β-sheet of polyA aggregates was dissociated and separated into monomers with helix structures upon FEL irradiation. After the dissociation by FEL, water molecules inhibited the reformation of polyA aggregates. We recently verified the same dissociation process using FEL-treated amyloid-β aggregates. Thus, a common mechanism underlies the dissociation of different protein aggregates that cause different diseases, polyA disease and Alzheimer's disease. However, MD simulation indicated that polyA aggregates are less easily dissociated than amyloid-β aggregates and require longer laser irradiation due to hydrophobic alanine repeats.

Probing protein misfolding and dissociation with an infrared free-electron laser

Misfolding is observed in the mutant proteins that are causative for neurodegenerative disorders such as polyglutamine diseases. These proteins are prone to aggregate in the cytoplasm and nucleus of cells. To reproduce cells with the aggregated proteins, gene expression system is usually applied, in which the expression construct having the mutated DNA sequence of the interest is transfected into cells. The transfected DNA is finally converted into the mutant protein, which is gradually aggregated in the cells. In addition, a simple method to prepare the cells having aggregates inside has been recently applied. Peptides were first aggregated by incubating them in water. The aggregates are spontaneously taken up by cells because aggregated proteins generally transfer between cells. Peptides with different degrees of aggregation can be made by changing the incubation times and temperatures, which enables to examine contribution of aggregation to the toxicity to the recipient cells. Moreover, such cells can be used for therapeutic researches of diseases in which aggregates are involved. In this chapter, we show methods to induce aggregation of peptides. The functional analyses of the cells with aggregates are also described. Then, experimental dissociation of the aggregates produced using this method by mid infrared free electron laser irradiation and its theoretical support by molecular dynamics simulation are introduced as the therapeutic research for neurodegenerative disorders.

The sequence context in poly-alanine regions: structure, function and conservation

MOTIVATION

Poly-alanine (polyA) regions are protein stretches mostly composed of alanines. Despite their abundance in eukaryotic proteomes and their association to nine inherited human diseases, the structural and functional roles exerted by polyA stretches remain poorly understood. In this work we study how the amino acid context in which polyA regions are settled in proteins influences their structure and function.

RESULTS

We identified glycine and proline as the most abundant amino acids within polyA and in the flanking regions of polyA tracts, in human proteins as well as in 17 additional eukaryotic species. Our analyses indicate that the non-structuring nature of these two amino acids influences the α-helical conformations predicted for polyA, suggesting a relevant role in reducing the inherent aggregation propensity of long polyA. Then, we show how polyA position in protein N-termini relates with their function as transit peptides. PolyA placed just after the initial methionine is often predicted as part of mitochondrial transit peptides, whereas when placed in downstream positions, polyA are part of signal peptides. A few examples from known structures suggest that short polyA can emerge by alanine substitutions in α-helices; but evolution by insertion is observed for longer polyA. Our results showcase the importance of studying the sequence context of homorepeats as a mechanism to shape their structure-function relationships.

AVAILABILITY AND IMPLEMENTATION

The datasets used and/or analyzed during the current study are available from the corresponding author onreasonable request.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A chemical probe for BAG1 targets androgen receptor-positive prostate cancer through oxidative stress signaling pathway

BAG1 is a family of polypeptides with a conserved C-terminal BAG domain that functions as a nucleotide exchange factor for the molecular chaperone HSP70. BAG1 proteins also control several signaling processes including proteostasis, apoptosis, and transcription. The largest isoform, BAG1L, controls the activity of the androgen receptor (AR) and is upregulated in prostate cancer. Here, we show that BAG1L regulates AR dynamics in the nucleus and its ablation attenuates AR target gene expression especially those involved in oxidative stress and metabolism. We show that a small molecule, A4B17, that targets the BAG domain downregulates AR target genes similar to a complete BAG1L knockout and upregulates the expression of oxidative stress-induced genes involved in cell death. Furthermore, A4B17 outperformed the clinically approved antagonist enzalutamide in inhibiting cell proliferation and prostate tumor development in a mouse xenograft model. BAG1 inhibitors therefore offer unique opportunities for antagonizing AR action and prostate cancer growth.

•

BAG1L interacts with a sequence overlapping a polyalanine tract in the AR NTD

•

Knockdown of BAG1L increase AR dynamics in the nucleus

•

BAG1L uses ROS pathway to regulate AR⁺ prostate cancer cell proliferation

•

A small molecule BAG1 inhibitor inhibits prostate tumor growth in mouse xenografts

Exogenous polyserine and polyleucine are toxic to recipient cells

Repeat-associated non-AUG (RAN) translation of mRNAs/transcripts responsible for polyglutamine (polyQ) diseases may generate peptides containing different mono amino acid tracts such as polyserine (polyS) and polyleucine (polyL). The propagation of aggregated polyQ from one cell to another is also an intriguing feature of polyQ proteins. However, whether the RAN translation-related polyS and polyL have the ability to propagate remains unclear, and if they do, whether the exogenous polyS and polyL exert toxicity on the recipient cells is also not known yet. In the present study, we found that aggregated polyS and polyL peptides spontaneously enter neuron-like cells and astrocytes in vitro. Aggregated polyS led to the degeneration of the differentiated neuron-like cultured cells. Likewise, the two types of aggregates taken up by astrocytes induced aberrant differentiation and cell death in vitro. Furthermore, injection of each of the two types of aggregates into the ventricles of adult mice resulted in their behavioral changes. The polyS-injected mice showed extensive vacuolar degeneration in the brain. Thus, the RAN translation-related proteins containing polyS and polyL have the potential to propagate and the proteins generated by all polyQ diseases might exert universal toxicity in the recipient cells.

Toxicity of internalized polyalanine to cells depends on aggregation

In polyalanine (PA) diseases, the disease-causing transcription factors contain an expansion of alanine repeats. While aggregated proteins that are responsible for the pathogenesis of neurodegenerative disorders show cell-to-cell propagation and thereby exert toxic effects on the recipient cells, whether this is also the case with expanded PA has not been studied. It is also not known whether the internalized PA is toxic to recipient cells based on the degree of aggregation. In this study, we therefore prepared different degrees of aggregation of a peptide having 13 alanine repeats without flanking sequences of PA disease-causative proteins (13A). The aggregated 13A was spontaneously taken up by neuron-like cultured cells. Functionally, strong aggregates but not weak aggregates displayed a deficit in neuron-like differentiation in vitro. Moreover, the injection of strong but not weak 13A aggregates into the ventricle of mice during the neonatal stage led to enhanced spontaneous motor activity later in life. Thus, PA in the extracellular space has the potential to enter adjacent cells, and may exert toxicity depending on the degree of aggregation.

Recombinant Silk Proteins with Additional Polyalanine Have Excellent Mechanical Properties

This paper explores the structures of exogenous protein molecules that can effectively improve the mechanical properties of silkworm silk. Several transgenic vectors fused with the silkworm fibroin light chain and type 3 repeats in different multiples of the ampullate dragline silk protein 1 (MaSp1) from black widow spider with different lengths of the polyalanine motifs were constructed for this study. Transgenic silkworms were successfully obtained by piggyBac-mediated microinjection. Molecular detection showed that foreign proteins were successfully secreted and contained within the cocoon shells. According to the prediction of PONDR^® VSL2 and PONDR^® VL-XT, the type 3 repeats and the polyalanine motif of the MaSp1 protein were amorphous. The results of FTIR analysis showed that the content of β-sheets in the silk of transgenic silkworms engineered with transgenic vectors with additional polyalanine was significantly higher than that of wild-type silkworm silk. Additionally, silk with a higher β-sheet content had better fracture strength and Young’s modulus. The mechanical properties of silk with longer chains of exogenous proteins were improved. In general, our results provide theoretical guidance and technical support for the large-scale production of excellent bionic silk.

Characterizing aggregate growth and morphology of alanine-rich polypeptides as a function of sequence chemistry and solution temperature from scattering, spectroscopy, and microscopy

The aggregation behavior and stability of a series of alanine-rich peptides, which are included as components of peptide-polymer conjugates, were characterized using a combination of biophysical techniques. Light scattering techniques were used to monitor changes in peptide morphology and size distributions as a function of time and temperature. The results show large particles immediately upon dissolution in buffer. At room temperature, these particles relaxed to reach a mostly monomeric peptide state, while at higher temperatures, they grew to form aggregates. Circular dichroism spectroscopy (CD) was used to monitor temperature- and time-dependent conformational changes as a function of peptide sequence and incubation time. CD measurements reveal that all of the sequences are helical at low temperatures with transitions to non-helical conformation with increased temperature. Samples incubated at room temperature were able to recover their original helicity. At increased temperature, the shorter and longer peptide sequences showed notable changes in conformation, and were not able to recover their original helicity after 72 h. After incubation for up to one week, β-sheet conformations were observed in these two cases, while only α-helical conformation loss was observed for the peptide of intermediate molecular weight. Transmission electron microscopy measurements reveal the formation of fibrils after 72 h of incubation at 60 °C for all samples, in agreement with the scattering measurements. Additional quenching experiments show that peptide aggregation can be stalled when solutions are cooled to room temperature.

Molecular insights into the role of the polyalanine region in mediating PHOX2B aggregation

About 90% of congenital central hypoventilation syndrome (CCHS) patients show polyalanine triplet expansions in the coding region of transcription factor PHOX2B, which renders this protein an intriguing target to understand the insurgence of this syndrome and for the design of a novel therapeutical approach. Consistently with the role of PHOX2B as a transcriptional regulator, it is reasonable that a general transcriptional dysregulation caused by the polyalanine expansion might represent an important mechanism underlying CCHS pathogenesis. Therefore, this study focused on the biochemical characterization of different PHOX2B variants, such as a variant containing the correct C-terminal (20 alanines) stretch, one of the most frequent polyalanine expansions (+7 alanines), and a variant lacking the complete alanine stretch (0 alanines). Comparison of the different variants by a multidisciplinary approach based on different methodologies (including circular dichroism, spectrofluorimetry, light scattering, and Atomic Force Microscopy studies) highlighted the propensity to aggregate for the PHOX2B variant containing the polyalanine expansion (+7-alanines), especially in the presence of DNA, while the 0-alanines variant resembled the protein with the correct polyalanine length. Moreover, and unexpectedly, the formation of fibrils was revealed only for the pathological variant, suggesting a plausible role of such fibrils in the insurgence of CCHS.

Is there codon usage bias for poly-Q stretches in the human proteome?

We have analyzed codon usage for poly-Q stretches of different lengths for the human proteome. First, we have obtained that all long poly-Q stretches in Protein Data Bank (PDB) belong to the disordered regions. Second, we have found the bias for codon usage for glutamine homo-repeats in the human proteome. In the cases when the same codon is used for poly-Q stretches only CAG triplets are found. Similar results are obtained for human proteins with glutamine homo-repeats associated with diseases. Moreover, for proteins associated with diseases (from the HraDis database), the fraction of proteins for which the same codon is used for glutamine homo-repeats is less (22%) than for proteins from the human proteome (26%). We have demonstrated for poly-Q stretches in the human proteome that in some cases (28) the splicing sites correspond to the homo-repeats and in 11 cases, these sites appear at the C -terminal part of the homo-repeats with statistical significance 10 -8 .

Injectable tissue integrating networks from recombinant polypeptides with tunable order

Emergent properties of natural biomaterials result from the collective effects of nanoscale interactions among ordered and disordered domains. Here, using recombinant sequence design, we have created a set of partially ordered polypeptides to study emergent hierarchical structures by precisely encoding nanoscale order-disorder interactions. These materials, which combine the stimuli-responsiveness of disordered elastin-like polypeptides and the structural stability of polyalanine helices, are thermally responsive with tunable thermal hysteresis and the ability to reversibly form porous, viscoelastic networks above threshold temperatures. Through coarse-grain simulations, we show that hysteresis arises from physical crosslinking due to mesoscale phase separation of ordered and disordered domains. On injection of partially ordered polypeptides designed to transition at body temperature, they form stable, porous scaffolds that rapidly integrate into surrounding tissue with minimal inflammation and a high degree of vascularization. Sequence-level modulation of structural order and disorder is an untapped principle for the design of functional protein-based biomaterials.

Glutamine Codon Usage and polyQ Evolution in Primates Depend on the Q Stretch Length

Amino acid usage in a proteome depends mostly on its taxonomy, as it does the codon usage in transcriptomes. Here, we explore the level of variation in the codon usage of a specific amino acid, glutamine, in relation to the number of consecutive glutamine residues. We show that CAG triplets are consistently more abundant in short glutamine homorepeats (polyQ, four to eight residues) than in shorter glutamine stretches (one to three residues), leading to the evolutionary growth of the repeat region in a CAG-dependent manner. The length of orthologous polyQ regions is mostly stable in primates, particularly the short ones. Interestingly, given a short polyQ the CAG usage is higher in unstable-in-length orthologous polyQ regions. This indicates that CAG triplets produce the necessary instability for a glutamine stretch to grow. Proteins related to polyQ-associated diseases behave in a more extreme way, with longer glutamine stretches in human and evolutionarily closer nonhuman primates, and an overall higher CAG usage. In the light of our results, we suggest an evolutionary model to explain the glutamine codon usage in polyQ regions.

Intrinsic Disorder in Proteins with Pathogenic Repeat Expansions

Intrinsically disordered proteins and proteins with intrinsically disordered regions have been shown to be highly prevalent in disease. Furthermore, disease-causing expansions of the regions containing tandem amino acid repeats often push repetitive proteins towards formation of irreversible aggregates. In fact, in disease-relevant proteins, the increased repeat length often positively correlates with the increased aggregation efficiency and the increased disease severity and penetrance, being negatively correlated with the age of disease onset. The major categories of repeat extensions involved in disease include poly-glutamine and poly-alanine homorepeats, which are often times located in the intrinsically disordered regions, as well as repeats in non-coding regions of genes typically encoding proteins with ordered structures. Repeats in such non-coding regions of genes can be expressed at the mRNA level. Although they can affect the expression levels of encoded proteins, they are not translated as parts of an affected protein and have no effect on its structure. However, in some cases, the repetitive mRNAs can be translated in a non-canonical manner, generating highly repetitive peptides of different length and amino acid composition. The repeat extension-caused aggregation of a repetitive protein may represent a pivotal step for its transformation into a proteotoxic entity that can lead to pathology. The goals of this article are to systematically analyze molecular mechanisms of the proteinopathies caused by the poly-glutamine and poly-alanine homorepeat expansion, as well as by the polypeptides generated as a result of the microsatellite expansions in non-coding gene regions and to examine the related proteins. We also present results of the analysis of the prevalence and functional roles of intrinsic disorder in proteins associated with pathological repeat expansions.

Structure Formation in Langmuir Peptide Films As Revealed from Coarse-Grained Molecular Dynamics Simulations

Molecular dynamics simulations in conjunction with the Martini coarse-grained model have been used to investigate the (nonequilibrium) behavior of helical 22-residue poly(γ-benzyl-l-glutamate) (PBLG) peptides at the water/vapor interface. Preformed PBLG mono- or bilayers homogeneously covering the water surface laterally collapse in tens of nanoseconds, exposing significant proportions of empty water surface. This behavior was also observed in recent AFM experiments at similar areas per monomer, where a complete coverage had been assumed in earlier work. In the simulations, depending on the area per monomer, either elongated clusters or fibrils form, whose heights (together with the portion of empty water surface) increase over time. Peptides tend to align with respect to the fiber axis or with the major principal axis of the cluster, respectively. The aspect ratio of the cluster observed is 1.7 and, hence, comparable to though somewhat smaller than the aspect ratio of the peptides in α-helical conformation, which is 2.2. The heights of the fibrils is 3 nm after 20 ns and increases to 4.5 nm if the relaxation time is increased by 2 orders of magnitude, in agreement with the experiment. Aggregates with heights of about 3 or 4.5 nm are found to correspond to local bi- or trilayer structures, respectively.

Buckling transition in long α-helices

The treatment of bending and buckling of stiff biopolymer filaments by the popular worm-like chain model does not provide adequate understanding of these processes at the microscopic level. Thus, we have used the atomistic molecular-dynamic simulations and the Amber03 force field to examine the compression buckling of α-helix (AH) filaments at room temperature. It was found that the buckling instability occurs in AHs at the critical force f(c) in the range of tens of pN depending on the AH length. The decrease of the force f(c) with the contour length follows the prediction of the classic thin rod theory. At the force f(c) the helical filament undergoes the swift and irreversible transition from the smoothly bent structure to the buckled one. A sharp kink in the AH contour arises at the transition, accompanied by the disruption of the hydrogen bonds in its vicinity. The kink defect brings in an effective softening of the AH molecule at buckling. Nonbonded interactions between helical branches drive the rearrangement of a kinked AH into the ultimate buckled structure of a compact helical hairpin described earlier in the literature.

Automated synthesis of backbone protected peptides

The automated introduction of removable substitution along a peptide backbone prevents chain-association and synthesis failure.

The synthesis of peptides rich in aggregation prone sequences can be improved with backbone protection. We report the automated introduction of backbone protection to a peptide. This new method was applied in a fully-automated synthesis, giving improved handling, quality and yield of several challenging target sequences.

Role of side-chain interactions on the formation of α-helices in model peptides

The role played by side-chain interactions on the formation of α-helices is studied using extensive all-atom molecular dynamics simulations of polyalanine-like peptides in explicit TIP4P water. The peptide is described by the OPLS-AA force field except for the Lennard-Jones interaction between Cβ-Cβ atoms, which is modified systematically. We identify values of the Lennard-Jones parameter that promote α-helix formation. To rationalize these results, potentials of mean force (PMF) between methane-like molecules that mimic side chains in our polyalanine-like peptides are computed. These PMF exhibit a complex distance dependence where global and local minima are separated by an energy barrier. We show that α-helix propensity correlates with values of these PMF at distances corresponding to Cβ-Cβ of i-i+3 and other nearest neighbors in the α-helix. In particular, the set of Lennard-Jones parameters that promote α-helices is characterized by PMF that exhibit a global minimum at distances corresponding to i-i+3 neighbors in α-helices. Implications of these results are discussed.

Stimuli responsive fibrous hydrogels from hierarchical self-assembly of a triblock copolypeptide

In this work, the self-assembly behavior and pH responsiveness of a triblock copolypeptide in aqueous media are demonstrated. The copolypeptide was composed of a central pH responsive poly(l-glutamic acid) (PGA), flanked by two hydrophobic poly(l-alanine) blocks (PAla) (PAla5-PGA11-PAla5). This system showed a pH-responsive transition from short tapes to spherical aggregates by increasing the pH, as a result of deprotonation of the PGA block and a conformational change from α-helix to random coil. Increasing the ionic strength to physiological conditions (0.15 M) has triggered fibrillar self-assembly through intermolecular hydrogen bonding of PAla end-blocks that form β-sheet nanostructures, in conjunction with charge screening of the central random coil PGA segments. At elevated concentrations a thermo-responsive free supporting hydrogel was obtained, consisting of rigid β-sheet based twisted superfibers, resulting from hierarchical self-assembly of the copolypeptide. Yet, morphological transformation of this nanostructure was observed upon switching the pH from physiological conditions to pH 4. An unexpected morphology constituted of α-helix-based giant nanobelts was observed as a consequence of the secondary peptide transitions.

Effects of hydrophobic interaction strength on the self-assembled structures of model peptides

Stable and ordered self-assembled peptide nanostructures are formed as a result of cooperative effects of various relatively weak intermolecular interactions. We systematically studied the influence of hydrophobic interaction strength and temperature on the self-assembly of peptides with a coarse-grained model by Monte Carlo simulations. The simulation results show a rich phase behavior of peptide self-assembly, indicating that the formation and morphology of peptide assemblies may be tuned by varying the temperature and the strength of hydrophobic interactions. There exist optimal combinations of temperature and hydrophobic interaction strength where ordered fibrillar nanostructures are readily formed. Our simulation results not only facilitate the understanding of the self-assembly behavior of peptides at the molecular level, but also provide useful insights into the development of fabrication strategies for high-quality peptide fibrils.

Association of polyalanine and polyglutamine coiled coils mediates expansion disease-related protein aggregation and dysfunction

The expansion of homopolymeric glutamine (polyQ) or alanine (polyA) repeats in certain proteins owing to genetic mutations induces protein aggregation and toxicity, causing at least 18 human diseases. PolyQ and polyA repeats can also associate in the same proteins, but the general extent of their association in proteomes is unknown. Furthermore, the structural mechanisms by which their expansion causes disease are not well understood, and these repeats are generally thought to misfold upon expansion into aggregation-prone β-sheet structures like amyloids. However, recent evidence indicates a critical role for coiled-coil (CC) structures in triggering aggregation and toxicity of polyQ-expanded proteins, raising the possibility that polyA repeats may as well form these structures, by themselves or in association with polyQ. We found through bioinformatics screenings that polyA, polyQ and polyQA repeats have a phylogenetically graded association in human and non-human proteomes and associate/overlap with CC domains. Circular dichroism and cross-linking experiments revealed that polyA repeats can form—alone or with polyQ and polyQA—CC structures that increase in stability with polyA length, forming higher-order multimers and polymers in vitro. Using structure-guided mutagenesis, we studied the relevance of polyA CCs to the in vivo aggregation and toxicity of RUNX2—a polyQ/polyA protein associated with cleidocranial dysplasia upon polyA expansion—and found that the stability of its polyQ/polyA CC controls its aggregation, localization and toxicity. These findings indicate that, like polyQ, polyA repeats form CC structures that can trigger protein aggregation and toxicity upon expansion in human genetic diseases.

Length-dependent β-sheet growth mechanisms of polyalanine peptides in water and on hydrophobic surfaces

Fibrillar assemblies by peptides are becoming one of the most promising nanomaterials due to their exceptional properties. The self-assembly of peptides into β sheets is a critical step in the fibrillization pathway. We investigated the length-dependent β-sheet growth mechanisms of polyalanine [poly(A)] peptides consisting of 6 to 24 alanines (A6 to A24) in water and on the hydrophobic surface, respectively, by molecular dynamics simulations. β-sheet growth behavior in water fits negative exponential growth model, showing that β-sheet growth rate decays exponentially with time. Meanwhile, increasing chain length leads to an accelerated decay of the β-sheet growth rate. By contrast, β-sheet growth on the surface from A6 to A18 occurs in two consecutive stages, both of which fit linear growth models. β-sheet growth rate in the first stage increases as chain length is increased, while the intermediate length peptide A12 has the highest β-sheet growth rate in the second stage. β-sheet growth behavior of A24 on the surface still fits negative exponential model. Overall, the hydrophobic surface accelerates β-sheet growth by enhancing local concentration and reducing conformational entropy of poly(A) peptide, and the β-sheet growth of the intermediate length peptide A12 is the fastest on the surface. Our simulation results shed light on understanding the accelerated peptide fibrillization on the hydrophobic surface.

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.

Misfolded polyglutamine, polyalanine, and superoxide dismutase 1 aggregate via distinct pathways in the cell

Protein aggregation into intracellular inclusions is a key feature of many neurodegenerative disorders. A common theme has emerged that inappropriate self-aggregation of misfolded or mutant polypeptide sequences is detrimental to cell health. Yet protein quality control mechanisms may also deliberately cluster them together into distinct inclusion subtypes, including the insoluble protein deposit (IPOD) and the juxtanuclear quality control (JUNQ). Here we investigated how the intrinsic oligomeric state of three model systems of disease-relevant mutant protein and peptide sequences relates to the IPOD and JUNQ patterns of aggregation using sedimentation velocity analysis. Two of the models (polyalanine (37A) and superoxide dismutase 1 (SOD1) mutants A4V and G85R) accumulated into the same JUNQ-like inclusion whereas the other, polyglutamine (72Q), formed spatially distinct IPOD-like inclusions. Using flow cytometry pulse shape analysis (PulSA) to separate cells with inclusions from those without revealed the SOD1 mutants and 37A to have abruptly altered oligomeric states with respect to the nonaggregating forms, regardless of whether cells had inclusions or not, whereas 72Q was almost exclusively monomeric until inclusions formed. We propose that mutations leading to JUNQ inclusions induce a constitutively "misfolded" state exposing hydrophobic side chains that attract and ultimately overextend protein quality capacity, which leads to aggregation into JUNQ inclusions. Poly(Q) is not misfolded in this same sense due to universal polar side chains, but is highly prone to forming amyloid fibrils that we propose invoke a different engagement mechanism with quality control.

Factors that drive peptide assembly and fibril formation: experimental and theoretical analysis of Sup35 NNQQNY mutants

Residue mutations have substantial effects on aggregation kinetics and propensities of amyloid peptides and their aggregate morphologies. Such effects are attributed to conformational transitions accessed by various types of oligomers such as steric zipper or single β-sheet. We have studied the aggregation propensities of six NNQQNY mutants: NVVVVY, NNVVNV, NNVVNY, VIQVVY, NVVQIY, and NVQVVY in water using a combination of ion-mobility mass spectrometry, transmission electron microscopy, atomic force microscopy, and all-atom molecular dynamics simulations. Our data show a strong correlation between the tendency to form early β-sheet oligomers and the subsequent aggregation propensity. Our molecular dynamics simulations indicate that the stability of a steric zipper structure can enhance the propensity for fibril formation. Such stability can be attained by either hydrophobic interactions in the mutant peptide or polar side-chain interdigitations in the wild-type peptide. The overall results display only modest agreement with the aggregation propensity prediction methods such as PASTA, Zyggregator, and RosettaProfile, suggesting the need for better parametrization and model peptides for these algorithms.

Non-Arrhenius protein aggregation

Protein aggregation presents one of the key challenges in the development of protein biotherapeutics. It affects not only product quality but also potentially impacts safety, as protein aggregates have been shown to be linked with cytotoxicity and patient immunogenicity. Therefore, investigations of protein aggregation remain a major focus in pharmaceutical companies and academic institutions. Due to the complexity of the aggregation process and temperature-dependent conformational stability, temperature-induced protein aggregation is often non-Arrhenius over even relatively small temperature windows relevant for product development, and this makes low-temperature extrapolation difficult based simply on accelerated stability studies at high temperatures. This review discusses the non-Arrhenius nature of the temperature dependence of protein aggregation, explores possible causes, and considers inherent hurdles for accurately extrapolating aggregation rates from conventional industrial approaches for selecting accelerated conditions and from conventional or more advanced methods of analyzing the resulting rate data.

Thermodynamics of ideal proteinogenic homopolymer chains as a function of the energy spectrum E, helical propensity ω and enthalpic energy barrier

A reformulation and generalization of the Zwanzig model (ZW model) for ideal homopolymer chains poly-X, where X represents any of the twenty naturally occurring proteinogenic amino acid residues is presented. This reformulation and generalization provides a direct connection between coarse-grained parameters originally proposed in the ZW model with variables from the Lifson-Roig (LR) theory, such as the helical propensity per residue ω, and new variables introduced here, such as the energy gap Δ between unfolded and folded structures, as well as the ratio f of the energy scales involved. This enables us to discover the relevance of the energy spectrum E to the onset of configurational phase transitions. From the configurational partition function Q, thermodynamic properties such as the configurational entropy S, specific heat v and average energy <E> are calculated in terms of the number of residues K, temperature T, helical propensity ω and energy barrier ΔH for different poly-X chains in vacuo. Results obtained here provide substantial evidence that configurational phase transitions for ideal poly-X chains correspond to first-order phase transitions. An anomalous behavior of the thermodynamic functions <E>, Cv, S with respect to the number K of residues is also highlighted. On-going methods of solution are outlined.

Reduced atomic pair-interaction design (RAPID) model for simulations of proteins

Increasingly, theoretical studies of proteins focus on large systems. This trend demands the development of computational models that are fast, to overcome the growing complexity, and accurate, to capture the physically relevant features. To address this demand, we introduce a protein model that uses all-atom architecture to ensure the highest level of chemical detail while employing effective pair potentials to represent the effect of solvent to achieve the maximum speed. The effective potentials are derived for amino acid residues based on the condition that the solvent-free model matches the relevant pair-distribution functions observed in explicit solvent simulations. As a test, the model is applied to alanine polypeptides. For the chain with 10 amino acid residues, the model is found to reproduce properly the native state and its population. Small discrepancies are observed for other folding properties and can be attributed to the approximations inherent in the model. The transferability of the generated effective potentials is investigated in simulations of a longer peptide with 25 residues. A minimal set of potentials is identified that leads to qualitatively correct results in comparison with the explicit solvent simulations. Further tests, conducted for multiple peptide chains, show that the transferable model correctly reproduces the experimentally observed tendency of polyalanines to aggregate into β-sheets more strongly with the growing length of the peptide chain. Taken together, the reported results suggest that the proposed model could be used to succesfully simulate folding and aggregation of small peptides in atomic detail. Further tests are needed to assess the strengths and limitations of the model more thoroughly.

Kinetic analysis of aggregation data

Aggregation of repeat-containing proteins is associated with neurodegenerative disorders; a specific example is the established link between expansion of the polyglutamine domain in huntingtin and the appearance of nuclear inclusions in Huntington's disease. This connection between aggregation and pathology has motivated numerous investigations into the kinetics of aggregation. Quantitative analysis of kinetic data is needed both for comparative purposes (e.g., to compare the effect of different compounds on aggregation kinetics) and for mechanistic insight. Here we describe some analytical equations that can be used to model aggregation data and demonstrate appropriate and simple methods for extracting valid model parameters by fitting equations to kinetic data.

Conformational properties of polyglutamine sequences in guanidine hydrochloride solutions

Two sets of iso-1-cytochrome c variants have been prepared with N-terminal insertions of pure polyglutamine, i.e., PolyQ variants, or polyglutamine interrupted with lysine every sixth residue, i.e., Gln-rich variants. The polymer properties of these pure polyGln or Gln-rich sequences have been evaluated using equilibrium and kinetic His-heme loop formation methods for loop sizes ranging from 22 to 46 in 1.5, 3.0, and 6.0 M guanidine hydrochloride (GdnHCl). In 6.0 M GdnHCl, the scaling exponent, ν(3), for the pure polyGln sequences, is ~1.7--significantly less than ν(3) ≈ 2.15 for the Gln-rich sequences. The stability of the His-heme loops becomes progressively greater for the pure polyGln sequences relative to the Gln-rich sequences as GdnHCl concentration decreases from 6.0 to 1.5 M. Thus, the context of the sequence effects the polymer properties of Gln repeats even in denaturing concentrations of GdnHCl. Comparison of data for the Gln-rich variants with previous results for Gly-rich and Ala-rich variants shows that ν(3) ~ 2.2 for the Gln-rich, Gly-rich, and Ala-rich sequences in 6.0 M GdnHCl, whereas ν(3) remains unchanged at 3.0 M GdnHCl concentration for the Gln-rich and Ala-rich sequences but decreases to ~1.7 for the Gly-rich sequences. Thus, the polymer properties of Gln-rich and Ala-rich sequences are less sensitive to solvent quality in denaturing solutions of GdnHCl than Gly-rich sequences. Evaluation of Flory's characteristic ratio, C(n), for the Gln-rich and Ala-rich sequences relative to the Gly-rich sequences shows that Gln-rich sequences are stiffer than Ala-rich sequences at both 3.0 and 6.0 M GdnHCl.

Understanding protein aggregation from the view of monomer dynamics

Much work in recent years has been devoted to understanding the complex process of protein aggregation. This review looks at the earliest stages of aggregation, long before the formation of fibrils that are the hallmark of many aggregation-based diseases, and proposes that the first steps are controlled by the reconfiguration dynamics of the monomer. When reconfiguration is much faster or much slower than bimolecular diffusion, then aggregation is slow, but when they are similar, aggregation is fast. The experimental evidence for this model is reviewed and the prospects for small molecule aggregation inhibitors to prevent disease are discussed.

Explaining the length threshold of polyglutamine aggregation

The existence of a length threshold, of about 35 residues, above which polyglutamine repeats can give rise to aggregation and to pathologies, is one of the hallmarks of polyglutamine neurodegenerative diseases such as Huntington's disease. The reason why such a minimal length exists at all has remained one of the main open issues in research on the molecular origins of such classes of diseases. Following the seminal proposals of Perutz, most research has focused on the hunt for a special structure, attainable only above the minimal length, able to trigger aggregation. Such a structure has remained elusive and there is growing evidence that it might not exist at all. Here we review some basic polymer and statistical physics facts and show that the existence of a threshold is compatible with the modulation that the repeat length imposes on the association and dissociation rates of polyglutamine polypeptides to and from oligomers. In particular, their dramatically different functional dependence on the length rationalizes the very presence of a threshold and hints at the cellular processes that might be at play, in vivo, to prevent aggregation and the consequent onset of the disease.

Dynamic imaging by fluorescence correlation spectroscopy identifies diverse populations of polyglutamine oligomers formed in vivo

Protein misfolding and aggregation are exacerbated by aging and diseases of protein conformation including neurodegeneration, metabolic diseases, and cancer. In the cellular environment, aggregates can exist as discrete entities, or heterogeneous complexes of diverse solubility and conformational state. In this study, we have examined the in vivo dynamics of aggregation using imaging methods including fluorescence microscopy, fluorescence recovery after photobleaching (FRAP), and fluorescence correlation spectroscopy (FCS), to monitor the diverse biophysical states of expanded polyglutamine (polyQ) proteins expressed in Caenorhabditis elegans. We show that monomers, oligomers and aggregates co-exist at different concentrations in young and aged animals expressing different polyQ-lengths. During aging, when aggregation and toxicity are exacerbated, FCS-based burst analysis and purified single molecule FCS detected a populational shift toward an increase in the frequency of brighter and larger oligomeric species. Regardless of age or polyQ-length, oligomers were maintained in a heterogeneous distribution that spans multiple orders of magnitude in brightness. We employed genetic suppressors that prevent polyQ aggregation and observed a reduction in visible immobile species with the persistence of heterogeneous oligomers, yet our analysis did not detect the appearance of any discrete oligomeric states associated with toxicity. These studies reveal that the reversible transition from monomers to immobile aggregates is not represented by discrete oligomeric states, but rather suggests that the process of aggregation involves a more complex pattern of molecular interactions of diverse intermediate species that can appear in vivo and contribute to aggregate formation and toxicity.

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.