Structure-activity relationship of amyloid fibrils

Protein aggregation is a process in which proteins self-associate into imperfectly ordered macroscopic entities. Such aggregates are generally classified as either amorphous or highly ordered, the most common form of the latter being amyloid fibrils. Amyloid fibrils composed of cross-β-sheet structure are the pathological hallmarks of several diseases including Alzheimer’s disease, but are also associated with functional states such as the fungal HET-s prion. This review aims to summarize the recent high-resolution structural studies of amyloid fibrils in light of their (potential) activities. We propose that the repetitive nature of the cross-β-sheet structure of amyloids is key for their multiple properties: the repeating motifs can translate a rather non-specific interaction into a specific one through cooperativity.     

What can AlphaFold do for antimicrobial amyloids?

Amyloids, protein, and peptide assemblies in various organisms are crucial in physiological and pathological processes. Their intricate structures, however, present significant challenges, limiting our understanding of their functions, regulatory mechanisms, and potential applications in biomedicine and technology. This study evaluated the AlphaFold2 ColabFold method's structure predictions for antimicrobial amyloids, using eight antimicrobial peptides (AMPs), including those with experimentally determined structures and AMPs known for their distinct amyloidogenic morphological features. Additionally, two well-known human amyloids, amyloid-β and islet amyloid polypeptide, were included in the analysis due to their disease relevance, short sequences, and antimicrobial properties. Amyloids typically exhibit tightly mated β-strand sheets forming a cross-β configuration. However, certain amphipathic α-helical subunits can also form amyloid fibrils adopting a cross-α structure. Some AMPs in the study exhibited a combination of cross-α and cross-β amyloid fibrils, adding complexity to structure prediction. The results showed that the AlphaFold2 ColabFold models favored α-helical structures in the tested amyloids, successfully predicting the presence of α-helical mated sheets and a hydrophobic core resembling the cross-α configuration. This implies that the AI-based algorithms prefer assemblies of the monomeric state, which was frequently predicted as helical, or capture an α-helical membrane-active form of toxic peptides, which is triggered upon interaction with lipid membranes.     

Switching in amyloid structure within individual fibrils: Implication for strain adaptation, species barrier and strain classification

Amyloid fibrils are highly ordered crystal-like structures. It is generally assumed that individual amyloid fibrils consist of conformationally uniform cross-β-sheet structures that enable the amyloids to replicate their individual conformations via a template-dependent mechanism. Recent studies revealed that amyloids are capable of accommodating a global conformational switch from one amyloid strain to another within individual fibrils. The current review highlights the high adaptation potential of amyloid structures and discusses the implication of these findings for several emerging issues including prion strain adaptation (i.e. gradual change in strain structure). It also proposes that the catalytic activity of an amyloid structure should be separated from its templating effect, and raises the question of strain classification according to their promiscuous or species-specific nature.     

Biology of amyloid: structure, function, and regulation

Amyloids are highly ordered cross-β sheet protein aggregates associated with many diseases including Alzheimer's disease, but also with biological functions such as hormone storage. The cross-β sheet entity comprising an indefinitely repeating intermolecular β sheet motif is unique among protein folds. It grows by recruitment of the corresponding amyloid protein, while its repetitiveness can translate what would be a nonspecific activity as monomer into a potent one through cooperativity. Furthermore, the one-dimensional crystal-like repeat in the amyloid provides a structural framework for polymorphisms. This review summarizes the recent high-resolution structural studies of amyloid fibrils in light of their biological activities. We discuss how the unique properties of amyloids gives rise to many activities and further speculate about currently undocumented biological roles for the amyloid entity. In particular, we propose that amyloids could have existed in a prebiotic world, and may have been the first functional protein fold in living cells.     

Stem-forming regions that are essential for the amyloidogenesis of prion proteins

Prion diseases represent fatal neurodegenerative disorders caused by the aggregation of prion proteins. With regard to the formation of the amyloidogenic cross-β-structure, the initial mechanism in the conversion to a β-structure is critically important. To explore the core regions forming a stem of the amyloid, we designed and prepared a series of peptides comprised of two native sequences linked by a turn-inducing dipeptide moiety and examined their ability to produce amyloids. A sequence alignment of the peptides bearing the ability to form amyloid structures revealed that paired strands consisting of VNITI (residues 180-184) and VTTTT (residues 189-193) are the core regions responsible for initiating the formation of cross-β-structures and for further ordered aggregation. In addition, most of the causative mutations responsible for inherited prion diseases were found to be located in these stem-forming regions on helix H2 and their counterpart on helix H3. Moreover, the volume effect of the nonstem domain, which contains ~200 residues, was deduced to be a determinant of the nature of the association such as oligomerization, because the stem-forming domain is only a small part of a prion protein. Taken together, we conclude that the mechanism underlying the initial stage of amyloidogenesis is the exposure of a newly formed intramolecular β-sheet to a solvent through the partial transition of a native structure from an α-helix to a β-structure. Our results also demonstrate that prion diseases caused by major prion proteins except the prions of some fungi such as yeast are inherent only in mammals, as evidenced by a comparison of the corresponding sequences to the stem-forming regions among different animals.     

Secondary structure switching in Cro protein evolution

We report the solution structure of the Cro protein from bacteriophage P22. Comparisons of its sequence and structure to those of lambda Cro strongly suggest an alpha-to-beta secondary structure switching event during Cro evolution. The folds of P22 Cro and lambda Cro share a three alpha helix fragment comprising the N-terminal half of the domain. However, P22 Cro's C terminus folds as two helices, while lambda Cro's folds as a beta hairpin. The all-alpha fold found for P22 Cro appears to be ancestral, since it also occurs in cI proteins, which are anciently duplicated paralogues of Cro. PSI-BLAST and transitive homology analyses strongly suggest that the sequences of P22 Cro and lambda Cro are globally homologous despite encoding different folds. The alpha+beta fold of lambda Cro therefore likely evolved from its all-alpha ancestor by homologous secondary structure switching, rather than by nonhomologous replacement of both sequence and structure.     

Degradation of fungal prion HET-s(218-289) induces formation of a generic amyloid fold

The prion-forming domain of the fungal prion protein HET-s, HET-s(218-289), is known from solid-state NMR studies to have a β-solenoidal structure; the β-solenoid has the cross-β structure characteristic of all amyloids, but is inherently more complex than the generic stacked β-sheets found in studies of small synthetic peptides. At low pH HET-s(218-289) has also been reported to form an alternative structure, which has not been characterized. We have confirmed by x-ray fiber diffraction that HET-s(218-289) adopts a β-solenoidal structure at neutral pH, and shown that at low pH, it forms either a β-solenoid or a stacked β-sheet structure, depending on the integrity of the protein and the conditions of fibrillization. The low pH stacked-sheet structure is usually formed only by proteolyzed HET-s(218-289), but intact HET-s(218-289) can form stacked sheets when seeded with proteolyzed stacked-sheet HET-s(218-289). The polymorphism of HET-s parallels the structural differences between the infectious brain-derived and the much less infectious recombinant mammalian prion protein PrP. Taken together, these observations suggest that the functional or pathological forms of amyloid proteins are more complex than the simple generic stacked-sheet amyloids commonly formed by short peptides.     

Sequence and structure classification of kinases

Kinases are a ubiquitous group of enzymes that catalyze the phosphoryl transfer reaction from a phosphate donor (usually ATP) to a receptor substrate. Although all kinases catalyze essentially the same phosphoryl transfer reaction, they display remarkable diversity in their substrate specificity, structure, and the pathways in which they participate. In order to learn the relationship between structural fold and functional specificities in kinases, we have done a comprehensive survey of all available kinase sequences (>17,000) and classified them into 30 distinct families based on sequence similarities. Of these families, 19, covering nearly 98% of all sequences, fall into seven general structural folds for which three-dimensional structures are known. These fold groups include some of the most widespread protein folds, such as Rossmann fold, ferredoxin fold, ribonuclease H fold, and TIM beta/alpha-barrel. On the basis of this classification system, we examined the shared substrate binding and catalytic mechanisms as well as variations of these mechanisms in the same fold groups. Cases of convergent evolution of identical kinase activities occurring in different folds are discussed.     

Functional Amyloids: Where Supramolecular Amyloid Assembly Controls Biological Activity or Generates New Functionality

Functional amyloids are a rapidly expanding class of fibrillar protein structures, with a core cross-β scaffold, where novel and advantageous biological function is generated by the assembly of the amyloid. The growing number of amyloid structures determined at high resolution reveal how this supramolecular template both accommodates a wide variety of amino acid sequences and also imposes selectivity on the assembly process. The amyloid fibril can no longer be considered a generic aggregate, even when associated with disease and loss of function. In functional amyloids the polymeric β-sheet rich structure provides multiple different examples of unique control mechanisms and structures that are finely tuned to deliver assembly or disassembly in response to physiological or environmental cues. Here we review the range of mechanisms at play in natural, functional amyloids, where tight control of amyloidogenicity is achieved by environmental triggers of conformational change, proteolytic generation of amyloidogenic fragments, or heteromeric seeding and amyloid fibril stability. In the amyloid fibril form, activity can be regulated by pH, ligand binding and higher order protofilament or fibril architectures that impact the arrangement of associated domains and amyloid stability. The growing understanding of the molecular basis for the control of structure and functionality delivered by natural amyloids in nearly all life forms should inform the development of therapies for amyloid-associated diseases and guide the design of innovative biomaterials.     

Computational modelling of diatom silicic acid transporters predicts a conserved fold with implications for their function and evolution

Diatoms are an important group of algae that can produce intricate silicified cell walls (frustules). The complex process of silicification involves a set of enigmatic integral membrane proteins that are thought to actively transport the soluble precursor of biosilica, dissolved silicic acid. Full-length silicic acid transporters are found widely across the diatoms while homologous shorter proteins have now been identified in a range of other organisms. It has been suggested that modern silicic acid transporters arose from the union of such partial sequences. Here, we present a computational study of the silicic acid transporters and related transporter-like sequences to help understand the structure, function and evolution of this class of membrane protein. The AlphaFold software predicts that all of the protein sequences studied here share a common fold in the membrane domain which is entirely different from the predicted folds of non-homologous silicic acid transporters from plants. Substrate docking reveals how conserved polar residues could interact with silicic acid at a central solvent-accessible binding site, consistent with an alternating access mechanism of transport. The structural conservation between these proteins supports a model where modern silicon transporters evolved from smaller ancestral proteins by gene fusion.     

On the evolution of protein folds: are similar motifs in different protein folds the result of convergence, insertion, or relics of an ancient peptide world?

This paper presents and discusses evidence suggesting how the diversity of domain folds in existence today might have evolved from peptide ancestors. We apply a structure similarity detection method to detect instances where localized regions of different protein folds contain highly similar sequences and structures. Results of performing an all-on-all comparison of known structures are described and compared with other recently published findings. The numerous instances of local sequence and structure similarities within different protein folds, together with evidence from proteins containing sequence and structure repeats, argues in favor of the evolution of modern single polypeptide domains from ancient short peptide ancestors (antecedent domain segments (ADSs)). In this model, ancient protein structures were formed by self-assembling aggregates of short polypeptides. Subsequently, and perhaps concomitantly with the evolution of higher fidelity DNA replication and repair systems, single polypeptide domains arose from the fusion of ADSs genes. Thus modern protein domains may have a polyphyletic origin.     

A search method for homologs of small proteins. Ubiquitin-like proteins in prokaryotic cells?

The question of protein homology versus analogy arises when proteins share a common function or a common structural fold without any statistically significant amino acid sequence similarity. Even though two or more proteins do not have similar sequences but share a common fold and the same or closely related function, they are assumed to be homologs, descendant from a common ancestor. The problem of homolog identification is compounded in the case of proteins of 100 or less amino acids. This is due to a limited number of basic single domain folds and to a likelihood of identifying by chance sequence similarity. The latter arises from two conditions: first, any search of the currently very large protein database is likely to identify short regions of chance match; secondly, a direct sequence comparison among a small set of short proteins sharing a similar fold can detect many similar patterns of hydrophobicity even if proteins do not descend from a common ancestor. In an effort to identify distant homologs of the many ubiquitin proteins, we have developed a combined structure and sequence similarity approach that attempts to overcome the above limitations of homolog identification. This approach results in the identification of 90 probable ubiquitin-related proteins, including examples from the two prokaryotic domains of life, Archaea and Bacteria.     

Robustness of predictions of extremely thermally stable proteins in ancient organisms

A number of studies have addressed the environmental temperatures experienced by ancient life. Computational studies using a nonhomogeneous evolution model have estimated ancestral G + C contents of ribosomal RNAs and the amino acid compositions of ancestral proteins, generating hypotheses regarding the mesophilic last universal common ancestor. In contrast, our previous study computationally reconstructed ancestral amino acid sequences of nucleoside diphosphate kinases using a homogeneous model and then empirically resurrected the ancestral proteins. The thermal stabilities of these ancestral proteins were equivalent to or greater than those of extant homologous thermophilic proteins, supporting the thermophilic universal ancestor theory. In this study, we reinferred ancestral sequences using a dataset from which hyperthermophilic sequences were excluded. We also reinferred ancestral sequences using a nonhomogeneous evolution model. The newly reconstructed ancestral proteins are still thermally stable, further supporting the hypothesis that the ancient organisms contained thermally stable proteins and therefore that they were thermophilic.     

Computational exploration of the network of sequence flow between protein structures

We investigate small sequence adjustments (of one or a few amino acids) that induce large conformational transitions between distinct and stable folds of proteins. Such transitions are intriguing from evolutionary and protein‐design perspectives. They make it possible to search for ancient protein structures or to design protein switches that flip between folds and functions. A network of sequence flow between protein folds is computed for representative structures of the Protein Data Bank. The computed network is dense, on an average each structure is connected to tens of other folds. Proteins that attract sequences from a higher than expected number of neighboring folds are more likely to be enzymes and alpha/beta fold. The large number of connections between folds may reflect the need of enzymes to adjust their structures for alternative substrates. The network of the Cro family is discussed, and we speculate that capacity is an important factor (but not the only one) that determines protein evolution. The experimentally observed flip from all alpha to alpha + beta fold is examined by the network tools. A kinetic model for the transition of sequences between the folds (with only protein stability in mind) is proposed. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

The role of negative selection in protein evolution revealed through the energetics of the native state ensemble

Knowing the determinants of conformational specificity is essential for understanding protein structure, stability, and fold evolution. To address this issue, a novel statistical measure of energetic compatibility between sequence and structure was developed using an experimentally validated model of the energetics of the native state ensemble. This approach successfully matched sequences from a diverse subset of the human proteome to their respective folds. Unexpectedly, significant energetic compatibility between ostensibly unrelated sequences and structures was also observed. Interrogation of these matches revealed a general framework for understanding the origins of conformational specificity within a proteome: specificity is a complex function of both the ability of a sequence to adopt folds other than the native, and ability of a fold to accommodate sequences other than the native. The regional variation in energetic compatibility indicates that the compatibility is dominated by incompatibility of sequence for alternative fold segments, suggesting that evolution of protein sequences has involved substantial negative selection, with certain segments serving as “gatekeepers” that presumably prevent alternative structures. Beyond these global trends, a size dependence exists in the degree to which the energetic compatibility is determined from negative selection, with smaller proteins displaying more negative selection. This partially explains how short sequences can adopt unique folds, despite the higher probability in shorter proteins for small numbers of mutations to increase compatibility with other folds. In providing evolutionary ground rules for the thermodynamic relationship between sequence and fold, this framework imparts valuable insight for rational design of unique folds or fold switches. Proteins 2016; 84:435–447. © 2016 Wiley Periodicals, Inc.     

Thoroughly sampling sequence space: large-scale protein design of structural ensembles

Modeling the inherent flexibility of the protein backbone as part of computational protein design is necessary to capture the behavior of real proteins and is a prerequisite for the accurate exploration of protein sequence space. We present the results of a broad exploration of sequence space, with backbone flexibility, through a novel approach: large-scale protein design to structural ensembles. A distributed computing architecture has allowed us to generate hundreds of thousands of diverse sequences for a set of 253 naturally occurring proteins, allowing exciting insights into the nature of protein sequence space. Designing to a structural ensemble produces a much greater diversity of sequences than previous studies have reported, and homology searches using profiles derived from the designed sequences against the Protein Data Bank show that the relevance and quality of the sequences is not diminished. The designed sequences have greater overall diversity than corresponding natural sequence alignments, and no direct correlations are seen between the diversity of natural sequence alignments and the diversity of the corresponding designed sequences. For structures in the same fold, the sequence entropies of the designed sequences cluster together tightly. This tight clustering of sequence entropies within a fold and the separation of sequence entropy distributions for different folds suggest that the diversity of designed sequences is primarily determined by a structure's overall fold, and that the designability principle postulated from studies of simple models holds in real proteins. This has important implications for experimental protein design and engineering, as well as providing insight into protein evolution.     

Multiplicity of carbohydrate-binding sites in <Emphasis Type="Italic">β</Emphasis>-prism fold lectins: occurrence and possible evolutionary implications

The β-prism II fold lectins of known structure, all from monocots, invariably have three carbohydrate-binding sites in each subunit/domain. Until recently, β-prism I fold lectins of known structure were all from dicots and they exhibited one carbohydrate-binding site per subunit/domain. However, the recently determined structure of the β-prism fold I lectin from banana, a monocot, has two very similar carbohydrate-binding sites. This prompted a detailed analysis of all the sequences appropriate for two-lectin folds and which carry one or more relevant carbohydrate-binding motifs. The very recent observation of a β-prism I fold lectin, griffithsin, with three binding sites in each domain further confirmed the need for such an analysis. The analysis demonstrates substantial diversity in the number of binding sites unrelated to the taxonomical position of the plant source. However, the number of binding sites and the symmetry within the sequence exhibit reasonable correlation. The distribution of the two families of β-prism fold lectins among plants and the number of binding sites in them, appear to suggest that both of them arose through successive gene duplication, fusion and divergent evolution of the same primitive carbohydrate-binding motif involving a Greek key. Analysis with sequences in individual Greek keys as independent units lends further support to this conclusion. It would seem that the preponderance of three carbohydrate-binding sites per domain in monocot lectins, particularly those with the β-prism II fold, is related to the role of plant lectins in defence.     

More than the sum of their parts: on the evolution of proteins from peptides

Despite their seemingly endless diversity, proteins adopt a limited number of structural forms. It has been estimated that 80% of proteins will be found to adopt one of only about 400 folds, most of which are already known. These folds are largely formed by a limited 'vocabulary' of recurring supersecondary structure elements, often by repetition of the same element and, increasingly, elements similar in both structure and sequence are discovered. This suggests that modern proteins evolved by fusion and recombination from a more ancient peptide world and that many of the core folds observed today may contain homologous building blocks. The peptides forming these building blocks would not in themselves have had the ability to fold, but would have emerged as cofactors supporting RNA-based replication and catalysis (the 'RNA world'). Their association into larger structures and eventual fusion into polypeptide chains would have allowed them to become independent of their RNA scaffold, leading to the evolution of a novel type of macromolecule: the folded protein.     

Cell surface proteins in archaeal and bacterial genomes comprising "LVIVD", "RIVW" and "LGxL" tandem sequence repeats are predicted to fold as beta-propeller

Proteins that share even low sequence homologies are known to adopt similar folds. The beta-propeller structural motif is one such example. Identifying sequences that adopt a beta-propeller fold is useful to annotate protein structure and function. Often, tandem sequence repeats provide the necessary signal for identifying beta-propellers in proteins. In our recent analysis to identify cell surface proteins in archaeal and bacterial genomes, we identified some proteins that contain novel tandem repeats "LVIVD", "RIVW" and "LGxL". In this work, based on protein fold predictions and three-dimensional comparative modeling methods, we predicted that these repeat types fold as beta-propeller. Further, the evolutionary trace analysis of all proteins constituting amino acid sequence repeats in beta-propellers suggest that the novel repeats have diverged from a common ancestor.     

Smooth Functional Transition along a Mutational Pathway with an Abrupt Protein Fold Switch

Recent protein design experiments have demonstrated that proteins can migrate between folds through the accumulation of substitution mutations without visiting disordered or nonfunctional points in sequence space. To explore the biophysical mechanism underlying such transitions we use a three-letter continuous protein model with seven atoms per amino acid to provide realistic sequence-structure and sequence-function mappings through explicit simulation of the folding and interaction of model sequences. We start from two 16-amino-acid sequences folding into an α-helix and a β-hairpin, respectively, each of which has a preferred binding partner with 35 amino acids. We identify a mutational pathway between the two folds, which features a sharp fold switch. By contrast, we find that the transition in function is smooth. Moreover, the switch in preferred binding partner does not coincide with the fold switch. Discovery of new folds in evolution might therefore be facilitated by following fitness slopes in sequence space underpinned by binding-induced conformational switching.