Frequencies of amino acid strings in globular protein sequences indicate suppression of blocks of consecutive hydrophobic residues

Order of amino acid recruitment into the genetic code resolved by last universal common ancestor's protein domains

The current "consensus" order in which amino acids were added to the genetic code is based on potentially biased criteria, such as the absence of sulfur-containing amino acids from the Urey-Miller experiment which lacked sulfur. More broadly, abiotic abundance might not reflect biotic abundance in the organisms in which the genetic code evolved. Here, we instead identify which protein domains date to the last universal common ancestor (LUCA) and then infer the order of recruitment from deviations of their ancestrally reconstructed amino acid frequencies from the still-ancient post-LUCA controls. We find that smaller amino acids were added to the code earlier, with no additional predictive power in the previous consensus order. Metal-binding (cysteine and histidine) and sulfur-containing (cysteine and methionine) amino acids were added to the genetic code much earlier than previously thought. Methionine and histidine were added to the code earlier than expected from their molecular weights and glutamine later. Early methionine availability is compatible with inferred early use of S-adenosylmethionine and early histidine with its purine-like structure and the demand for metal binding. Even more ancient protein sequences-those that had already diversified into multiple distinct copies prior to LUCA-have significantly higher frequencies of aromatic amino acids (tryptophan, tyrosine, phenylalanine, and histidine) and lower frequencies of valine and glutamic acid than single-copy LUCA sequences. If at least some of these sequences predate the current code, then their distinct enrichment patterns provide hints about earlier, alternative genetic codes.

Discovery of an on-pathway protein folding intermediate illuminates the kinetic competition between folding and misfolding

Our current understanding of protein folding is based predominantly on studies of small (<150 aa) proteins that refold reversibly from a chemically denatured state. As protein length increases, the competition between off-pathway misfolding and on-pathway folding likewise increases, creating a more complex energy landscape. Little is known about how intermediates populated during the folding of larger proteins affect navigation of this more complex landscape. Previously, we reported extremely slow folding rates for the 539 aa β-helical passenger domain of pertactin (P.69T), including conditions that favor the formation of a kinetically trapped, off-pathway partially folded state (PFS). The existence of an on-pathway intermediate for P.69T folding was speculated but its characterization remained elusive. In this work, we exploited the extremely slow kinetics of PFS unfolding to develop a double-jump "denaturant challenge" assay. With this assay, we identified a transient unfolding intermediate, PFS*, that adopts a similar structure to PFS, including C-terminal folded structure and a disordered N-terminus, yet unfolds much more quickly than PFS. Additional experiments revealed that PFS* also functions as an on-pathway intermediate for P.69T folding. Collectively, these results support a two-step, C-to-N-terminal model for P.69T folding: folding initiates in the C-terminus with the rate-limiting formation of the transient on-pathway PFS* intermediate, which sits at the junction of the kinetic competition between folding and misfolding. Notably, processive folding from C-to-N-terminus also occurs during C-to-N-terminal translocation of P.69T across the bacterial outer membrane. These results illuminate the crucial role of kinetics when navigating a complex energy landscape for protein folding.

Distribution and solvent exposure of Hsp70 chaperone binding sites across the Escherichia coli proteome

Many proteins must interact with molecular chaperones to achieve their native state in the cell. Yet, how chaperone binding-site characteristics affect the folding process is poorly understood. The ubiquitous Hsp70 chaperone system prevents client-protein aggregation by holding unfolded conformations and by unfolding misfolded states. Hsp70 binding sites of client proteins comprise a nonpolar core surrounded by positively charged residues. However, a detailed analysis of Hsp70 binding sites on a proteome-wide scale is still lacking. Further, it is not known whether proteins undergo some degree of folding while chaperone bound. Here, we begin to address the above questions by identifying Hsp70 binding sites in 2258 Escherichia coli (E. coli) proteins. We find that most proteins bear at least one Hsp70 binding site and that the number of Hsp70 binding sites is directly proportional to protein size. Aggregation propensity upon release from the ribosome correlates with number of Hsp70 binding sites only in the case of large proteins. Interestingly, Hsp70 binding sites are more solvent-exposed than other nonpolar sites, in protein native states. Our findings show that the majority of E. coli proteins are systematically enabled to interact with Hsp70 even if this interaction only takes place during a fraction of the protein lifetime. In addition, our data suggest that some conformational sampling may take place within Hsp70-bound states, due to the solvent exposure of some chaperone binding sites in native proteins. In all, we propose that Hsp70-chaperone-binding traits have evolved to favor Hsp70-assisted protein folding devoid of aggregation.

The Role of Proteolysis in Amyloidosis

Amyloidoses are a group of diseases associated with deposits of amyloid fibrils in different tissues. So far, 36 different types of amyloidosis are known, each due to the misfolding and accumulation of a specific protein. Amyloid deposits can be found in several organs, including the heart, brain, kidneys, and spleen, and can affect single or multiple organs. Generally, amyloid-forming proteins become prone to aggregate due to genetic mutations, acquired environmental factors, excessive concentration, or post-translational modifications. Interestingly, amyloid aggregates are often composed of proteolytic fragments, derived from the degradation of precursor proteins by yet unidentified proteases, which display higher amyloidogenic tendency compared to precursor proteins, thus representing an important mechanism in the onset of amyloid-based diseases. In the present review, we summarize the current knowledge on the proteolytic susceptibility of three of the main human amyloidogenic proteins, i.e., transthyretin, β-amyloid precursor protein, and α-synuclein, in the onset of amyloidosis. We also highlight the role that proteolytic enzymes can play in the crosstalk between intestinal inflammation and amyloid-based diseases.

Proteome-wide landscape of solubility limits in a bacterial cell

Proteins are prone to aggregate when expressed above their solubility limits. Aggregation may occur rapidly, potentially as early as proteins emerge from the ribosome, or slowly, following synthesis. However, in vivo data on aggregation rates are scarce. Here, we classified the Escherichia coli proteome into rapidly and slowly aggregating proteins using an in vivo image-based screen coupled with machine learning. We find that the majority (70%) of cytosolic proteins that become insoluble upon overexpression have relatively low rates of aggregation and are unlikely to aggregate co-translationally. Remarkably, such proteins exhibit higher folding rates compared to rapidly aggregating proteins, potentially implying that they aggregate after reaching their folded states. Furthermore, we find that a substantial fraction (~ 35%) of the proteome remain soluble at concentrations much higher than those found naturally, indicating a large margin of safety to tolerate gene expression changes. We show that high disorder content and low surface stickiness are major determinants of high solubility and are favored in abundant bacterial proteins. Overall, our study provides a global view of aggregation rates and hence solubility limits of proteins in a bacterial cell.

Large-scale design and refinement of stable proteins using sequence-only models

Engineered proteins generally must possess a stable structure in order to achieve their designed function. Stable designs, however, are astronomically rare within the space of all possible amino acid sequences. As a consequence, many designs must be tested computationally and experimentally in order to find stable ones, which is expensive in terms of time and resources. Here we use a high-throughput, low-fidelity assay to experimentally evaluate the stability of approximately 200,000 novel proteins. These include a wide range of sequence perturbations, providing a baseline for future work in the field. We build a neural network model that predicts protein stability given only sequences of amino acids, and compare its performance to the assayed values. We also report another network model that is able to generate the amino acid sequences of novel stable proteins given requested secondary sequences. Finally, we show that the predictive model-despite weaknesses including a noisy data set-can be used to substantially increase the stability of both expert-designed and model-generated proteins.

Obtaining Tertiary Protein Structures by the ab Initio Interpretation of Small Angle X-ray Scattering Data

Small angle X-ray scattering (SAXS) is an important tool for investigating the structure of proteins in solution. We present a novel ab initio method representing polypeptide chains as discrete curves used to derive a meaningful three-dimensional model from only the primary sequence and SAXS data. High resolution structures were used to generate probability density functions for each common secondary structural element found in proteins, which are used to place realistic restraints on the model curve’s geometry. This is coupled with a novel explicit hydration shell model in order to derive physically meaningful three-dimensional models by optimizing against experimental SAXS data. The efficacy of this model is verified on an established benchmark protein set, and then it is used to predict the lysozyme structure using only its primary sequence and SAXS data. The method is used to generate a biologically plausible model of the coiled-coil component of the human synaptonemal complex central element protein.

A topology-based investigation of protein interaction sites using Hydrophobic Cluster Analysis

Hydrophobic clusters, as defined by Hydrophobic Cluster Analysis (HCA), are conditioned binary patterns, made of hydrophobic and non-hydrophobic positions, whose limits fit well those of regular secondary structures. They were proved to be useful for predicting secondary structures in proteins from the only information of a single amino acid sequence and have permitted to assess, in a comprehensive way, the leading role of binary patterns in secondary structure preference towards a particular state. Here, we considered the available experimental 3D structures of protein globular domains to enlarge our previously reported hydrophobic cluster database (HCDB), almost doubling the number of hydrophobic cluster species (each species being defined by a unique binary pattern) that represent the most frequent structural bricks encountered within protein globular domains. We then used this updated HCDB to show that the hydrophobic amino acids of discordant clusters, i.e. those less abundant clusters for which the observed secondary structure is in disagreement with the binary pattern preference of the species to which they belong, are more exposed to solvent and are more involved in protein interfaces than the hydrophobic amino acids of concordant clusters. As amino acid composition differs between concordant/discordant clusters, considering binary patterns may be used to gain novel insights into key features of protein globular domain cores and surfaces. It can also provide useful information on possible conformational plasticity, including disorder to order transitions.

A Shift in Aggregation Avoidance Strategy Marks a Long-Term Direction to Protein Evolution

To detect a direction to evolution, without the pitfalls of reconstructing ancestral states, we need to compare "more evolved" to "less evolved" entities. But because all extant species have the same common ancestor, none are chronologically more evolved than any other. However, different gene families were born at different times, allowing us to compare young protein-coding genes to those that are older and hence have been evolving for longer. To be retained during evolution, a protein must not only have a function, but must also avoid toxic dysfunction such as protein aggregation. There is conflict between the two requirements: hydrophobic amino acids form the cores of protein folds, but also promote aggregation. Young genes avoid strongly hydrophobic amino acids, which is presumably the simplest solution to the aggregation problem. Here we show that young genes' few hydrophobic residues are clustered near one another along the primary sequence, presumably to assist folding. The higher aggregation risk created by the higher hydrophobicity of older genes is counteracted by more subtle effects in the ordering of the amino acids, including a reduction in the clustering of hydrophobic residues until they eventually become more interspersed than if distributed randomly. This interspersion has previously been reported to be a general property of proteins, but here we find that it is restricted to old genes. Quantitatively, the index of dispersion delineates a gradual trend, i.e., a decrease in the clustering of hydrophobic amino acids over billions of years.

Extending human IgG half-life using structure-guided design

Engineering of antibodies for improved pharmacokinetics through enhanced binding to the neonatal Fc receptor (FcRn) has been demonstrated in transgenic mice, non-human primates and humans. Traditionally, such approaches have largely relied on random mutagenesis and display formats, which fail to address related critical attributes of the antibody, such as effector functions or biophysical stability. We have developed a structure- and network-based framework to interrogate the engagement of IgG with multiple Fc receptors (FcRn, C1q, TRIM21, FcγRI, FcγRIIa/b, FcγRIIIa) simultaneously. Using this framework, we identified features that govern Fc-FcRn interactions and identified multiple distinct pathways for enhancing FcRn binding in a pH-specific manner. Network analysis provided a novel lens to study the allosteric impact of half-life-enhancing Fc mutations on FcγR engagement, which occurs distal to the FcRn binding site. Applying these principles, we engineered a panel of unique Fc variants that enhance FcRn binding while maintaining robust biophysical properties and wild type-like binding to activating receptors. An antibody harboring representative Fc designs demonstrates a half-life improvement of > 9 fold in transgenic mice and > 3.5 fold in cynomolgus monkeys, and maintains robust effector functions such as antibody-dependent cell-mediated cytotoxicity and complement-dependent cytotoxicity.

Distribution of dipeptides in different protein structural classes: an effort to find new similarities

Finding any regularity in the sequences of proteins and determining their correlation with structural features are of great interest for an understanding of molecular biology. We statistically analyzed the relative frequencies of all 400 possible dipeptides in a data set containing randomly selected proteins of different defined structural classes including all-alpha, all-beta, alpha + beta and alpha/beta families. We found that the distribution of dipeptides is not the same for different structural classes, and some of them are significantly far from a random distribution. A tendency of a given amino acid to localize in the first or second position of a dipeptide depending on the structural class of protein was also found. Interestingly, some amino acids may be substituted for each other in the first or second positions of specific dipeptides in each structural class. This finding apparently contrasts with the routine expectation from the viewpoint of amino acid properties, as classically understood.

Peptide binding landscapes: Specificity and homophilicity across sequence space in a lattice model

Peptide aggregation frequently involves sequences with strong homophilic binding character, i.e., sequences that self-assemble with like species in a crowded cellular environment, in the face of a multitude of other peptides or proteins as potential heterophilic binding partners. What kinds of sequences display a strong tendency towards homophilic binding and self-assembly, and what are the origins of this behavior? Here, we consider how sequence specificity in oligomerization processes plays out in a simple two-dimensional (2D) lattice statistical-thermodynamic peptide model that permits exhaustive examination of the entire sequence and configurational landscapes. We find that sequences with strong self-specificities have either alternating hydrophobic and hydrophilic residues or short patches of hydrophobic residues, both which minimize intramolecular hydrophobic interactions in part due to the constraints of the 2D lattice. We also find that these specificities are highly sensitive to entropic and free energetic features of the unbound conformational state, such that direct binding interaction energies alone do not capture the complete behavior. These results suggest that the ability of particular peptide sequences to self-assemble and aggregate in a many-protein environment reflects a precise balance of direct binding interactions and behavior in the unbound (monomeric) state.

Antibody Aggregation: Insights from Sequence and Structure

Monoclonal antibodies (mAbs) are the fastest-growing biological therapeutics with important applications ranging from cancers, autoimmunity diseases and metabolic disorders to emerging infectious diseases. Aggregation of mAbs continues to be a major problem in their developability. Antibody aggregation could be triggered by partial unfolding of its domains, leading to monomer-monomer association followed by nucleation and growth. Although the aggregation propensities of antibodies and antibody-based proteins can be affected by the external experimental conditions, they are strongly dependent on the intrinsic antibody properties as determined by their sequences and structures. In this review, we describe how the unfolding and aggregation susceptibilities of IgG could be related to their cognate sequences and structures. The impact of antibody domain structures on thermostability and aggregation propensities, and effective strategies to reduce aggregation are discussed. Finally, the aggregation of antibody-drug conjugates (ADCs) as related to their sequence/structure, linker payload, conjugation chemistry and drug-antibody ratio (DAR) is reviewed.

The Lepidopteran endoribonuclease-U domain protein P102 displays dramatically reduced enzymatic activity and forms functional amyloids

Hemocytes of Heliothis virescens (F.) (Lepidoptera, Noctuidae) larvae produce a protein, P102, with a putative endoribonuclease-U domain. In previous works we have shown that P102 is involved in Lepidopteran immune response by forming amyloid fibrils, which catalyze and localize melanin deposition around non-self intruders during encapsulation, preventing harmful systemic spreading. Here we demonstrate that P102 belongs to a new class of proteins that, at least in Lepidoptera, has a diminished endoribonuclease-U activity probably due to the lack of two out of five catalytically essential residues. We show that the P102 homolog from Trichoplusia ni (Lepidoptera, Noctuidae) displays catalytic site residues identical to P102, a residual endoribonuclease-U activity and the ability to form functional amyloids. On the basis of these results as well as sequence and structural analyses, we hypothesize that all the Lepidoptera endoribonuclease-U orthologs with catalytic site residues identical to P102 form a subfamily with similar function.

Structural insights into Bacillus thuringiensis Cry, Cyt and parasporin toxins

Since the first X-ray structure of Cry3Aa was revealed in 1991, numerous structures of B. thuringiensis toxins have been determined and published. In recent years, functional studies on the mode of action and resistance mechanism have been proposed, which notably promoted the developments of biological insecticides and insect-resistant transgenic crops. With the exploration of known pore-forming toxins (PFTs) structures, similarities between PFTs and B. thuringiensis toxins have provided great insights into receptor binding interactions and conformational changes from water-soluble to membrane pore-forming state of B. thuringiensis toxins. This review mainly focuses on the latest discoveries of the toxin working mechanism, with the emphasis on structural related progress. Based on the structural features, B. thuringiensis Cry, Cyt and parasporin toxins could be divided into three categories: three-domain type α-PFTs, Cyt toxin type β-PFTs and aerolysin type β-PFTs. Structures from each group are elucidated and discussed in relation to the latest data, respectively.

Edge strand engineering prevents native-like aggregation in Sulfolobus solfataricus acylphosphatase

β-proteins are constantly threatened by the risk of aggregation because β-sheets are inherently structured for edge-to-edge interactions. To avoid native-like aggregation, evolution has resulted in a set of strategies that prevent intermolecular β-interactions. Acylphosphatase from Sulfolobus solfataricus (Sso AcP) represents a suitable model for the study of such a process. Under conditions promoting aggregation, Sso AcP acquires a native-like conformational state whereby an unstructured N-terminal segment interacts with the edge β-strand B4 of an adjacent Sso AcP molecule. Because B4 is poorly protected against aggregation, this interaction triggers the aggregation cascade without the need for unfolding. Recently, three single Sso AcP mutants (V84D, Y86E and V84P) were designed to engineer additional protection against aggregation in B4 and were observed to successfully impair native-like aggregation in all three variants at the expense of a lower stability. To understand the structural basis of the reduced aggregation propensity and lower stability, the crystal structures of the Sso AcP variants were determined in the present study. Structural analysis reveals that the V84D and Y86E mutations exert protection by the insertion of an edge negative charge. A conformationally less regular B4 underlies protection against aggregation in the V84P mutant. The thermodynamic basis of instability is discussed. Moreover, kinetic experiments indicate that aggregation of the three mutants is not native-like and is independent of the interaction between B4 and the unstructured N-terminal segment. The reported data rationalize previous evidence regarding Sso AcP native-like aggregation and provide a basis for the design of aggregation-free proteins.

DATABASE

The atomic coordinates and related experimental data for the Sso AcP mutants V84P, V84D, ΔN11 Y86E have been deposited in the Protein Data Bank under accession numbers 4OJ3, 4OJG and 4OJH, respectively.

STRUCTURED DIGITAL ABSTRACT

• Sso AcP and Sso AcP bind by fluorescence technology (View interaction).

Chimera-induced folding: implications for amyloidosis

The discoveries that non-native proteins have a role in amyloidosis and that multiple protein misfolding diseases can occur concurrently suggest that cross-seeding of amyloidogenic proteins may be central to misfolding. To study this process, a synthetic chimeric amyloidogenic protein (YEHK21-YE8) composed of two components, one that readily folds to form fibrils (YEHK21) and one that does not (YE8), was designed. Secondary structural conformational changes during YEHK21-YE8 aggregation demonstrate that, under the appropriate conditions, YEHK21 is able to induce fibril formation of YE8. The unambiguous demonstration of the induction of folding and fibrillation within a single molecule illuminates the factors controlling this process and hence suggests the importance of those factors in amyloidogenic diseases.

Learning To Fold Proteins Using Energy Landscape Theory

This review is a tutorial for scientists interested in the problem of protein structure prediction, particularly those interested in using coarse-grained molecular dynamics models that are optimized using lessons learned from the energy landscape theory of protein folding. We also present a review of the results of the AMH/AMC/AMW/AWSEM family of coarse-grained molecular dynamics protein folding models to illustrate the points covered in the first part of the article. Accurate coarse-grained structure prediction models can be used to investigate a wide range of conceptual and mechanistic issues outside of protein structure prediction; specifically, the paper concludes by reviewing how AWSEM has in recent years been able to elucidate questions related to the unusual kinetic behavior of artificially designed proteins, multidomain protein misfolding, and the initial stages of protein aggregation.

Frustration in the energy landscapes of multidomain protein misfolding

Frustration from strong interdomain interactions can make misfolding a more severe problem in multidomain proteins than in single-domain proteins. On the basis of bioinformatic surveys, it has been suggested that lowering the sequence identity between neighboring domains is one of nature's solutions to the multidomain misfolding problem. We investigate folding of multidomain proteins using the associative-memory, water-mediated, structure and energy model (AWSEM), a predictive coarse-grained protein force field. We find that reducing sequence identity not only decreases the formation of domain-swapped contacts but also decreases the formation of strong self-recognition contacts between β-strands with high hydrophobic content. The ensembles of misfolded structures that result from forming these amyloid-like interactions are energetically disfavored compared with the native state, but entropically favored. Therefore, these ensembles are more stable than the native ensemble under denaturing conditions, such as high temperature. Domain-swapped contacts compete with self-recognition contacts in forming various trapped states, and point mutations can shift the balance between the two types of interaction. We predict that multidomain proteins that lack these specific strong interdomain interactions should fold reliably.

Structural basis of fluorescence quenching in caspase activatable-GFP

Apoptosis is critical for organismal homeostasis and a wide variety of diseases. Caspases are the ultimate executors of the apoptotic programmed cell death pathway. As caspases play such a central role in apoptosis, there is significant demand for technologies to monitor caspase function. We recently developed a caspase activatable-GFP (CA-GFP) reporter. CA-GFP is unique due to its "dark" state, where chromophore maturation of the GFP is inhibited by the presence of a C-terminal peptide. Here we show that chromophore maturation is prevented because CA-GFP does not fold into the robust β-barrel of GFP until the peptide has been cleaved by active caspase. Both CA-GFP and GFP₁₋₁₀ , a split form of GFP lacking the 11th strand, have similar secondary structure, different from mature GFP. A similar susceptibility to proteolytic digestion indicates that this shared structure is not the robust, fully formed GFP β-barrel. We have developed a model that suggests that as CA-GFP is translated in vivo it follows the same folding path as wild-type GFP; however, the presence of the appended peptide does not allow CA-GFP to form the barrel of the fully matured GFP. CA-GFP is therefore held in a "pro-folding" intermediate state until the peptide is released, allowing it to continue folding into the mature barrel geometry. This new understanding of the structural basis of the dark state of the CA-GFP reporter will enable manipulation of this mechanism in the development of reporter systems for any number of cellular processes involving proteases and potentially other enzymes.

Cellular strategies for regulating functional and nonfunctional protein aggregation

Growing evidence suggests that aggregation-prone proteins are both harmful and functional for a cell. How do cellular systems balance the detrimental and beneficial effect of protein aggregation? We reveal that aggregation-prone proteins are subject to differential transcriptional, translational, and degradation control compared to nonaggregation-prone proteins, which leads to their decreased synthesis, low abundance, and high turnover. Genetic modulators that enhance the aggregation phenotype are enriched in genes that influence expression homeostasis. Moreover, genes encoding aggregation-prone proteins are more likely to be harmful when overexpressed. The trends are evolutionarily conserved and suggest a strategy whereby cellular mechanisms specifically modulate the availability of aggregation-prone proteins to (1) keep concentrations below the critical ones required for aggregation and (2) shift the equilibrium between the monomeric and oligomeric/aggregate form, as explained by Le Chatelier’s principle. This strategy may prevent formation of undesirable aggregates and keep functional assemblies/aggregates under control.

Amyloidogenic properties of a D/N mutated 12 amino acid fragment of the C-terminal domain of the Cholesteryl-Ester Transfer Protein (CETP)

The cholesteryl-ester transfer protein (CETP) facilitates the transfer of cholesterol esters and triglycerides between lipoproteins in plasma where the critical site for its function is situated in the C-terminal domain. Our group has previously shown that this domain presents conformational changes in a non-lipid environment when the mutation D(470)N is introduced. Using a series of peptides derived from this C-terminal domain, the present study shows that these changes favor the induction of a secondary β-structure as characterized by spectroscopic analysis and fluorescence techniques. From this type of secondary structure, the formation of peptide aggregates and fibrillar structures with amyloid characteristics induced cytotoxicity in microglial cells in culture. These supramolecular structures promote cell cytotoxicity through the formation of reactive oxygen species (ROS) and change the balance of a series of proteins that control the process of endocytosis, similar to that observed when β-amyloid fibrils are employed. Therefore, a fine balance between the highly dynamic secondary structure of the C-terminal domain of CETP, the net charge, and the physicochemical characteristics of the surrounding microenvironment define the type of secondary structure acquired. Changes in this balance might favor misfolding in this region, which would alter the lipid transfer capacity conducted by CETP, favoring its propensity to substitute its physiological function.

Interactions in the native state of monellin, which play a protective role against aggregation

A series of recent studies have provided initial evidence about the role of specific intra-molecular interactions in maintaining proteins in their soluble state and in protecting them from aggregation. Here we show that the amino acid sequence of the protein monellin contains two aggregation-prone regions that are prevented from initiating aggregation by multiple non-covalent interactions that favor their burial within the folded state of the protein. By investigating the behavior of single-chain monellin and a series of five of its mutational variants using a variety of biochemical, biophysical and computational techniques, we found that weakening of the non-covalent interaction that stabilizes the native state of the protein leads to an enhanced aggregation propensity. The lag time for fibrillation was found to correlate with the apparent midpoint of thermal denaturation for the series of mutational variants, thus showing that a reduced thermal stability is associated with an increased aggregation tendency. We rationalize these findings by showing that the increase in the aggregation propensity upon mutation can be predicted in a quantitative manner through the increase in the exposure to solvent of the amyloidogenic regions of the sequence caused by the destabilization of the native state. Our findings, which are further discussed in terms of the structure of monellin and the perturbation by the amino acid substitutions of the contact surface between the two subdomains that compose the folded state of monellin, provide a detailed description of the specific intra-molecular interactions that prevent aggregation by stabilizing the native state of a protein.

Role of large hydrophobic residues in proteins

Large Hydrophobic Residues (LHR) such as phenylalanine, isoleucine, leucine, methionine and valine play an important role in protein structure and activity. We describe the role of LHR in complete set of protein sequences in 15 different species. That is the distribution of LHR in different proteins of different species is reported. It is observed that the proteins prefer to have 27% of large hydrophobic residues in total and all along the sequence. It is also observed that proteins accumulate more LHR in its active sites. A window analysis on these protein sequences shows that the 27% of LHR is more frequent at window length of 45 amino acids. The influenza virus and P. falciparum show a random distribution of LHR in its proteins compared to other model organisms.

The effect of loops on the structural organization of alpha-helical membrane proteins

Loops connecting the transmembrane (TM) alpha-helices in membrane proteins are expected to affect the structural organization of the thereby connected helices and the helical bundles as a whole. This effect, which has been largely ignored previously, is studied here by analyzing the x-ray structures of 41 alpha-helical membrane proteins. First we define the loop flexibility ratio, R, and find that 53% of the loops are stretched, where a stretched loop constrains the distance between the two connected helices. The significance of this constraining effect is supported by experiments carried out with bacteriorhodopsin and rhodopsin, in which cutting or eliminating their (predominately stretched) loops has led to a decrease in protein stability, and for rhodopsin, in most cases, also to the destruction of the structure. We show that for nonstretched loops in the extramembranous regions, the fraction of hydrophobic residues is comparable to that for soluble proteins; furthermore (as is also the case for soluble proteins), the hydrophobic residues in these regions are preferentially buried. This is expected to lead to the compact structural organization of the loops, which is transferred to the TM helices, causing them to assemble. We argue that a soluble protein complexed with a membrane protein similarly promotes compactness; other properties of such complexes are also studied. We calculate complementary attractive interactions between helices, including hydrogen bonds and van der Waals interactions of sequential motifs, such as GXXXG. The relative and combined effects of all these factors on the association of the TM helices are discussed and protein structures with only a few of these factors are analyzed. Our study emphasizes the need for classifying membrane proteins into groups according to structural organization. This classification should be considered when procedures for structural analysis or prediction are developed and applied. Detailed analysis of each structure is provided at http://flan.blm.cs.cmu.edu/memloop/

Folding by numbers: primary sequence statistics and their use in studying protein folding

The exponential growth over the past several decades in the quantity of both primary sequence data available and the number of protein structures determined has provided a wealth of information describing the relationship between protein primary sequence and tertiary structure. This growing repository of data has served as a prime source for statistical analysis, where underlying relationships between patterns of amino acids and protein structure can be uncovered. Here, we survey the main statistical approaches that have been used for identifying patterns within protein sequences, and discuss sequence pattern research as it relates to both secondary and tertiary protein structure. Limitations to statistical analyses are discussed, and a context for their role within the field of protein folding is given. We conclude by describing a novel statistical study of residue patterning in β-strands, which finds that hydrophobic (i,i+2) pairing in β-strands occurs more often than expected at locations near strand termini. Interpretations involving β-sheet nucleation and growth are discussed.

Amyloid peptides and proteins in review

Amyloids are filamentous protein deposits ranging in size from nanometres to microns and composed of aggregated peptide beta-sheets formed from parallel or anti-parallel alignments of peptide beta-strands. Amyloid-forming proteins have attracted a great deal of recent attention because of their association with over 30 diseases, notably neurodegenerative conditions like Alzheimer's, Huntington's, Parkinson's, Creutzfeldt-Jacob and prion disorders, but also systemic diseases such as amyotrophic lateral sclerosis (Lou Gehrig's disease) and type II diabetes. These diseases are all thought to involve important conformational changes in proteins, sometimes termed misfolding, that usually produce beta-sheet structures with a strong tendency to aggregate into water-insoluble fibrous polymers. Reasons for such conformational changes in vivo are still unclear. Intermediate aggregated state(s), rather than precipitated insoluble polymeric aggregates, have recently been implicated in cellular toxicity and may be the source of aberrant pathology in amyloid diseases. Numerous in vitro studies of short and medium length peptides that form amyloids have provided some clues to amyloid formation, with an alpha-helix to beta-sheet folding transition sometimes implicated as an intermediary step leading to amyloid formation. More recently, quite a few non-pathological amyloidogenic proteins have also been identified and physiological properties have been ascribed, challenging previous implications that amyloids were always disease causing. This article summarises a great deal of current knowledge on the occurrence, structure, folding pathways, chemistry and biology associated with amyloidogenic peptides and proteins and highlights some key factors that have been found to influence amyloidogenesis.

Prevention of amyloid-like aggregation as a driving force of protein evolution

Uncontrolled protein aggregation is a constant challenge in all compartments of living organisms. The failure of a peptide or protein to remain soluble often results in pathology. So far, more than 40 human diseases have been associated with the formation of extracellular fibrillar aggregates - known as amyloid fibrils - or structurally related intracellular deposits. It is well known that molecular chaperones and elaborate quality control mechanisms exist in the cell to counteract aggregation. However, an increasing number of reports during the past few years indicate that proteins have also evolved structural and sequence-based strategies to prevent aggregation. This review describes these strategies and the selection pressures that exist on protein sequences to combat their uncontrolled aggregation. We will describe the different types of mechanism evolved by proteins that adopt different conformational states including normally folded proteins, intrinsically disordered polypeptide chains, elastomeric systems and multimodular proteins.

Generic cell dysfunction in neurodegenerative disorders: role of surfaces in early protein misfolding, aggregation, and aggregate cytotoxicity

Recent knowledge supports the idea that early protein aggregates share basic structural features and are responsible for cytotoxicity underlying neurodegeneration; in most cases, early aggregate cytotoxicity apparently proceeds through similar molecular mechanisms and results in similar biochemical modifications. Data suggest that aggregate cytotoxicity may be considered a generic property of the oligomers preceding fibril appearance. Oligomers can interact with cell membranes, impairing their structural organization and destroying their selective ion permeability, eventually culminating with cell death. This process can be influenced by the physicochemical features and aggregation state of amyloids as well as by the physical and biochemical features of cell surfaces. The roles of synthetic and biological surfaces in affecting protein folding and misfolding, in speeding up aggregate nucleation, and as targets of aggregate toxicity is gaining consideration. Recent research has highlighted the involvement of surfaces as protein-misfolding chaperones and aggregation catalysts and their effects in these phenomena.

The distribution of residues in a polypeptide sequence is a determinant of aggregation optimized by evolution

It has been shown that the propensity of a protein to form amyloid-like fibrils can be predicted with high accuracy from the knowledge of its amino acid sequence. It has also been suggested, however, that some regions of the sequences are more important than others in determining the aggregation process. Here, we have addressed this issue by constructing a set of "sequence scrambled" variants of the first 29 residues of horse heart apomyoglobin (apoMb(1-29)), in which the sequence was modified while maintaining the same amino acid composition. The clustering of the most amyloidogenic residues in one region of the sequence was found to cause a marked increase of the elongation rate (k(agg)) and a remarkable shortening of the lag phase (t(lag)) of the fibril growth, as determined by far-UV circular dichroism and thioflavin T fluorescence. We also show that taking explicitly into consideration the presence of aggregation-promoting regions in the predictive methods results in a quantitative agreement between the theoretical and observed k(agg) and t(lag) values of the apoMb(1-29) variants. These results, together with a comparison between homologous segments from the family of globins, indicate the existence of a negative selection against the clustering of highly amyloidogenic residues in one or few regions of polypeptide sequences.

Protein aggregation in silico

Protein aggregation is a challenge to the successful manufacture of protein therapeutics; it can impose severe limitations on purification yields and compromise formulation stability. Advances in computer power, and the wealth of computational studies pertaining to protein folding, have facilitated the development of molecular simulation as a tool to investigate protein misfolding and aggregation. Here, we highlight the successes of protein aggregation studies carried out in silico, with a particular emphasis on studies related to biotechnology. To conclude, we discuss future prospects for the field, and identify several biotechnology-related problems that would benefit from molecular simulation.

Insights into protein structure and function from disorder-complexity space

Intrinsically disordered proteins have a wide variety of important functional roles. However, the relationship between sequence and function in these proteins is significantly different than that for well-folded proteins. In a previous work, we showed that the propensity to be disordered can be recognized based on sequence composition alone. Here that analysis is furthered by examining the relationship of disorder propensity to sequence complexity, where the metrics for these two properties depend only on composition. The distributions of 40 amino acid peptides from both ordered and disordered proteins are graphed in this disorder-complexity space. An analysis of Swiss-Prot shows that most peptides have high complexity and relatively low disorder. However, there are also an appreciable number of low complexity-high disorder peptides in the database. In contrast, there are no low complexity-low disorder peptides. A similar analysis for peptides in the PDB reveals a much narrower distribution, with few peptides of low complexity and high disorder. In this case, the bounds of the disorder-complexity distribution are well defined and might be used to evaluate the likelihood that a peptide can be crystallized with current methods. The disorder-complexity distributions of individual proteins and sets of proteins grouped by function are also examined. Among individual proteins, there is an enormous variety of distributions that in some cases can be rationalized with regard to function. Groups of functionally related proteins are found to have distributions that are similar within each group but show notable differences between groups. Finally, a pattern matching algorithm is used to search for proteins with particular disorder-complexity distributions. The results suggest that this approach might be used to identify relationships between otherwise dissimilar proteins.

Protein misfolding, functional amyloid, and human disease

Peptides or proteins convert under some conditions from their soluble forms into highly ordered fibrillar aggregates. Such transitions can give rise to pathological conditions ranging from neurodegenerative disorders to systemic amyloidoses. In this review, we identify the diseases known to be associated with formation of fibrillar aggregates and the specific peptides and proteins involved in each case. We describe, in addition, that living organisms can take advantage of the inherent ability of proteins to form such structures to generate novel and diverse biological functions. We review recent advances toward the elucidation of the structures of amyloid fibrils and the mechanisms of their formation at a molecular level. Finally, we discuss the relative importance of the common main-chain and side-chain interactions in determining the propensities of proteins to aggregate and describe some of the evidence that the oligomeric fibril precursors are the primary origins of pathological behavior.

Synthetic non-peptide mimetics of alpha-helices

Proteins in nature fold into native conformations in which combinations of peripherally projected aliphatic, aromatic and ionic functionalities direct a wide range of properties. Alpha-helices, one of the most common protein secondary structures, serve as important recognition regions on protein surfaces for numerous protein-protein, protein-DNA and protein-RNA interactions. These interactions are characterized by conserved structural features within the alpha-helical domain. Rational design of structural mimetics of these domains with synthetic small molecules has proven an effective means to modulate such protein functions. In this tutorial review we discuss strategies that utilize synthetic small-molecule antagonists to selectively target essential protein-protein interactions involved in certain diseases. We also evaluate some of the protein-protein interactions that have been or are potential targets for alpha-helix mimetics.

Diminishing of aggregation for bovine liver catalase through acidic residues modification

The tendency of proteins to aggregate is an important problem in biotechnology and the pharmaceutical industry. Because proteins in the aggregated state generally do not have the same biological activity as proteins in the native state. In order to prevent aggregation, it is essential to know the effective parameters in anti-aggregation mechanism. Using a chemical protein modification approach, UV-vis and fluorescence spectroscopies and circular dichroism spectropolarimetry, this study investigates the parameters involved in anti-aggregation mechanism of bovine liver catalase. Our findings clearly indicate that the modified bovine liver catalase provides better protection than the native enzyme against thermal aggregation. It seems that a decrease in hydrophobicity resulting in chemical modification plays an important role in preventing aggregation.

Point mutations in protein globular domains: contributions from function, stability and misfolding

Several contrasting hypotheses have been formulated about the influence of functional and conformational properties, like stability and avoidance of misfolding, on the evolution of protein globular domains. Selection at functional sites has been suggested to be detrimental to stability or coupled to it. Avoidance of misfolding may be achieved by discarding misfolding-prone sequences or by maintaining a stable native state and thus destabilizing partially or fully unfolded states from which misfolding can take place. We have performed a hierarchical analysis of a large database of point mutations to dissect the relative contributions of function, stability and misfolding in the evolution of natural sequences. We show that at catalytic sites, selection for function overrules selection for stability but find no evidence for an anticorrelation between function and stability. Selection for stability plays a secondary role at binding sites, but is not fully coupled to selection for function. Remarkably, we did not find a selective pressure against misfolding-prone sequences in globular proteins at the level of individual positions. We suggest that such a selection would compromise native-state stability due to a correlation between the stabilities of native and misfolded states. Stabilization of the native state is the most frequent way in which natural proteins avoid misfolding.

Entropic criteria for protein folding derived from recurrences: six residues patch as the basic protein word

Some research has suggested that patches of six constitute an important amino acid window length in proteins for conveying information. We present database evidence that supports this conjecture, as well as additional recurrence-based data that characterization and quantification of these words affect the folding/aggregation features of proteins. Other indirect evidence is presented and discussed.

High polar content of long buried blocks of sequence in protein domains suggests selection against amyloidogenic non-polar sequences

Native protein structures achieve stability in part by burying hydrophobic side-chains. About 75% of all amino acid residues buried in protein interiors are non-polar. Buried residues are not uniformly distributed in protein sequences, but sometimes cluster as contiguous polypeptide stretches that run through the interior of protein domain structures. Such regions have an intrinsically high local sequence density of non-polar residues, creating a potential problem: local non-polar sequences also promote protein misfolding and aggregation into non-native structures such as the amyloid fibrils in Alzheimer's disease. Here we show that long buried blocks of sequence in protein domains of known structure have, on average, a lower content of non-polar amino acids (about 70%) than do isolated buried residues (about 80%). This trend is observed both in small and in large protein domains and is independent of secondary structure. Long, completely non-polar buried stretches containing many large side-chains are particularly avoided. Aspartate residues that are incorporated in long buried stretches were found to make fewer polar interactions than those in short stretches, hinting that they may be destabilizing to the native state. We suggest that evolutionary pressure is acting on non-native properties, causing buried polar residues to be placed at positions where they would break up aggregation-prone non-polar sequences, perhaps even at some cost to native state stability.

Hydrophobic collapse in (in silico) protein folding

A model of hydrophobic collapse, which is treated as the driving force for protein folding, is presented. This model is the superposition of three models commonly used in protein structure prediction: (1) 'oil-drop' model introduced by Kauzmann, (2) a lattice model introduced to decrease the number of degrees of freedom for structural changes and (3) a model of the formation of hydrophobic core as a key feature in driving the folding of proteins. These three models together helped to develop the idea of a fuzzy-oil-drop as a model for an external force field of hydrophobic character mimicking the hydrophobicity-differentiated environment for hydrophobic collapse. All amino acids in the polypeptide interact pair-wise during the folding process (energy minimization procedure) and interact with the external hydrophobic force field defined by a three-dimensional Gaussian function. The value of the Gaussian function usually interpreted as a probability distribution is treated as a normalized hydrophobicity distribution, with its maximum in the center of the ellipsoid and decreasing proportionally with the distance versus the center. The fuzzy-oil-drop is elastic and changes its shape and size during the simulated folding procedure.

Frequencies of hydrophobic and hydrophilic runs and alternations in proteins of known structure

Patterns of alternation of hydrophobic and polar residues are a profound aspect of amino acid sequences, but a feature not easily interpreted for soluble proteins. Here we report statistics of hydrophobicity patterns in proteins of known structure in a current protein database as compared with results from earlier, more limited structure sets. Previous studies indicated that long hydrophobic runs, common in membrane proteins, are underrepresented in soluble proteins. Long runs of hydrophobic residues remain significantly underrepresented in soluble proteins, with none longer than 16 residues observed. These long runs most commonly occur as buried alpha helices, with extended hydrophobic strands less common. Avoiding aggregation of partially folded intermediates during intracellular folding remains a viable explanation for the rarity of long hydrophobic runs in soluble proteins. Comparison between database editions reveals robustness of statistics on aqueous proteins despite an approximately twofold increase in nonredundant sequences. The expanded database does now allow us to explain several deviations of hydrophobicity statistics from models of random sequence in terms of requirements of specific secondary structure elements. Comparison to prior membrane-bound protein sequences, however, shows significant qualitative changes, with the average hydrophobicity and frequency of long runs of hydrophobic residues noticeably increasing between the database editions. These results suggest that the aqueous proteins of solved structure may represent an essentially complete sample of the universe of aqueous sequences, while the membrane proteins of known structure are not yet representative of the universe of membrane-associated proteins, even by relatively simple measures of hydrophobic patterns.

Untangling influences of hydrophobicity on protein sequences and structures

We perform a statistical analysis of solvent accessibility and hydrophobicity profiles of a representative set of proteins. The joint probability distribution is well fitted to a multivariable Gaussian, which takes a relatively simple form when expressed in terms of the Fourier transforms of the profiles. This allows us to quantify the asymmetric manner by which these profiles influence each other. For example, the alpha-helix periodicity in sequence hydrophobicity is dictated by the solvent accessibility of structures, and not vice versa, possibly indicating the faster evolution of sequences compared to structures. The decorrelated hydrophobicity and solvent accessibility profiles show distinct behaviors at long periods, where sequence hydrophobicity fluctuates less, while solvent accessibility fluctuates more than average. The correlations between the two profiles can be interpreted as the Boltzmann weight of the solvation energy at room temperature, consistent with earlier observations.

Influence of denatured and intermediate states of folding on protein aggregation

We simulate the aggregation thermodynamics and kinetics of proteins L and G, each of which self-assembles to the same alpha/beta [corrected] topology through distinct folding mechanisms. We find that the aggregation kinetics of both proteins at an experimentally relevant concentration exhibit both fast and slow aggregation pathways, although a greater proportion of protein G aggregation events are slow relative to those of found for protein L. These kinetic differences are correlated with the amount and distribution of intrachain contacts formed in the denatured state ensemble (DSE), or an intermediate state ensemble (ISE) if it exists, as well as the folding timescales of the two proteins. Protein G aggregates more slowly than protein L due to its rapidly formed folding intermediate, which exhibits native intrachain contacts spread across the protein, suggesting that certain early folding intermediates may be selected for by evolution due to their protective role against unwanted aggregation. Protein L shows only localized native structure in the DSE with timescales of folding that are commensurate with the aggregation timescale, leaving it vulnerable to domain swapping or nonnative interactions with other chains that increase the aggregation rate. Folding experiments that characterize the structural signatures of the DSE, ISE, or the transition state ensemble (TSE) under nonaggregating conditions should be able to predict regions where interchain contacts will be made in the aggregate, and to predict slower aggregation rates for proteins with contacts that are dispersed across the fold. Since proteins L and G can both form amyloid fibrils, this work also provides mechanistic and structural insight into the formation of prefibrillar species.