The in vivo and in vitro aggregation properties of globular proteins correlate with their conformational stability: the SH3 case

Factors influencing amyloid fibril formation

Protein aggregation is a complex process with several stages that lead to the formation of complex structures and shapes with a broad variability in stability and toxicity. The aggregation process is affected by various factors and environmental conditions that disrupt the protein's original state, including internal factors like mutations, expression levels, and polypeptide chain truncation, as well as external factors, such as dense molecular surroundings, post-translation modifications, and interactions with other proteins, nucleic acids, small molecules, metal ions, chaperones, and lipid membranes. During the aggregation process, the biological activity of an aggregating protein may be reduced or eliminated, whereas the resulting aggregates may have the potential to be immunogenic, or they may have other undesirable properties. Finding the cause(s) of protein aggregation and controlling it to an acceptable level is among the most crucial topics of research in academia and biopharmaceutical companies. This chapter aims to review intrinsic pathways of protein aggregation and potential extrinsic variables that influence this process.

Assessment of Therapeutic Antibody Developability by Combinations of In Vitro and In Silico Methods

Although antibodies have become the fastest-growing class of therapeutics on the market, it is still challenging to develop them for therapeutic applications, which often require these molecules to withstand stresses that are not present in vivo. We define developability as the likelihood of an antibody candidate with suitable functionality to be developed into a manufacturable, stable, safe, and effective drug that can be formulated to high concentrations while retaining a long shelf life. The implementation of reliable developability assessments from the early stages of antibody discovery enables flagging and deselection of potentially problematic candidates, while focussing available resources on the development of the most promising ones. Currently, however, thorough developability assessment requires multiple in vitro assays, which makes it labor intensive and time consuming to implement at early stages. Furthermore, accurate in vitro analysis at the early stage is compromised by the high number of potential candidates that are often prepared at low quantities and purity. Recent improvements in the performance of computational predictors of developability potential are beginning to change this scenario. Many computational methods only require the knowledge of the amino acid sequences and can be used to identify possible developability issues or to rank available candidates according to a range of biophysical properties. Here, we describe how the implementation of in silico tools into antibody discovery pipelines is increasingly offering time- and cost-effective alternatives to in vitro experimental screening, thus streamlining the drug development process. We discuss in particular the biophysical and biochemical properties that underpin developability potential and their trade-offs, review various in vitro assays to measure such properties or parameters that are predictive of developability, and give an overview of the growing number of in silico tools available to predict properties important for antibody development, including the CamSol method developed in our laboratory.

Dynamic protein structures in normal function and pathologic misfolding in systemic amyloidosis

Dynamic and disordered regions in native proteins are often critical for their function, particularly in ligand binding and signaling. In certain proteins, however, such regions can contribute to misfolding and pathologic deposition as amyloid fibrils in vivo. For example, dynamic and disordered regions can promote amyloid formation by destabilizing the native structure, by directly triggering the aggregation, by promoting protein condensation, or by acting as sites of early proteolytic cleavage that favor a release of aggregation-prone fragments or facilitate fibril maturation. At the same time, enhanced dynamics in the native protein state accelerates proteolytic degradation that counteracts amyloid accumulation in vivo. Therefore, the functional need for dynamic protein regions must be balanced against their inherently labile nature. How exactly this balance is achieved and how is it shifted upon amyloidogenic mutations or post-translational modifications? To illustrate possible scenarios, here we review the beneficial and pathologic roles of dynamic and disordered regions in the native states of three families of human plasma proteins that form amyloid precursors in systemic amyloidoses: immunoglobulin light chain, apolipoproteins, and serum amyloid A. Analysis of structure, stability and local dynamics of these diverse proteins and their amyloidogenic variants exemplifies how disordered/dynamic regions can provide a functional advantage as well as an Achilles heel in pathologic amyloid formation.

Cryptic amyloidogenic regions in intrinsically disordered proteins: Function and disease association

The amyloid conformation is considered a fundamental state of proteins and the propensity to populate it a generic property of polypeptides. Multiple proteome-wide analyses addressed the presence of amyloidogenic regions in proteins, nurturing our understanding of their nature and biological implications. However, these analyses focused on highly aggregation-prone and hydrophobic stretches that are only marginally found in intrinsically disordered regions (IDRs). Here, we explore the prevalence of cryptic amyloidogenic regions (CARs) of polar nature in IDRs. CARs are widespread in IDRs and associated with IDPs function, with particular involvement in protein–protein interactions, but their presence is also connected to a risk of malfunction. By exploring this function/malfunction dichotomy, we speculate that ancestral CARs might have evolved into functional interacting regions playing a significant role in protein evolution at the origins of life.

PON-Sol2: Prediction of Effects of Variants on Protein Solubility

Genetic variations have a multitude of effects on proteins. A substantial number of variations affect protein-solvent interactions, either aggregation or solubility. Aggregation is often related to structural alterations, whereas solubilizable proteins in the solid phase can be made again soluble by dilution. Solubility is a central protein property and when reduced can lead to diseases. We developed a prediction method, PON-Sol2, to identify amino acid substitutions that increase, decrease, or have no effect on the protein solubility. The method is a machine learning tool utilizing gradient boosting algorithm and was trained on a large dataset of variants with different outcomes after the selection of features among a large number of tested properties. The method is fast and has high performance. The normalized correct prediction rate for three states is 0.656, and the normalized GC2 score is 0.312 in 10-fold cross-validation. The corresponding numbers in the blind test were 0.545 and 0.157. The performance was superior in comparison to previous methods. The PON-Sol2 predictor is freely available. It can be used to predict the solubility effects of variants for any organism, even in large-scale projects.

Amyloid-Like Peptide Aggregates

Self-assembly of proteins and peptides into the amyloid fold is a widespread phenomenon in the natural world. The structural hallmark of self-assembly into amyloid fibrillar assemblies is the cross-beta motif, which conveys distinct morphological and mechanical properties. The amyloid fibril formation has contrasting results depending on the organism, in the sense that it can bestow an organism with the advantages of mechanical strength and improved functionality or, on the contrary, could give rise to pathological states. In this chapter we review the existing information on amyloid-like peptide aggregates, which could either be derived from protein sequences, but also could be rationally or de novo designed in order to self-assemble into amyloid fibrils under physiological conditions. Moreover, the development of self-assembled fibrillar biomaterials that are tailored for the desired properties towards applications in biomedical or environmental areas is extensively analyzed. We also review computational studies predicting the amyloid propensity of the natural amino acid sequences and the structure of amyloids, as well as designing novel functional amyloid materials.

HDX and In Silico Docking Reveal that Excipients Stabilize G-CSF via a Combination of Preferential Exclusion and Specific Hotspot Interactions

Assuring the stability of therapeutic proteins is a major challenge in the biopharmaceutical industry, and a better molecular understanding of the mechanisms through which formulations influence their stability is an ongoing priority. While the preferential exclusion effects of excipients are well known, the additional presence and impact of specific protein-excipient interactions have proven to be more elusive to identify and characterize. We have taken a combined approach of in silico molecular docking and hydrogen deuterium exchange-mass spectrometry (HDX-MS) to characterize the interactions between granulocyte colony-stimulating factor (G-CSF), and some common excipients. These interactions were related to their influence on the thermal-melting temperatures (T_m) for the nonreversible unfolding of G-CSF in liquid formulations. The residue-level interaction sites predicted in silico correlated well with those identified experimentally and highlighted the potential impact of specific excipient interactions on the T_m of G-CSF.

Nucleic acids in inclusion bodies obtained from E. coli cells expressing human interferon-gamma

BACKGROUND

Inclusion bodies (IBs) are protein aggregates in recombinant bacterial cells containing mainly the target recombinant protein. Although it has been shown that IBs contain functional proteins along with protein aggregates, their direct application as pharmaceuticals is hindered by their heterogeneity and hazardous contaminants with bacterial origin. Therefore, together with the production of soluble species, IBs remain the main source for manufacture of recombinant proteins with medical application. The quality and composition of the IBs affect the refolding yield and further purification of the recombinant protein. The knowledge whether nucleic acids are genuine components or concomitant impurities of the IBs is a prerequisite for the understanding of the IBs formation and for development of optimized protocols for recombinant protein refolding and purification. IBs isolated from Escherichia coli overexpressing human interferon-gamma (hIFNγ), a protein with therapeutic application, were used as a model.

RESULTS

IBs were isolated from E. coli LE392 cells transformed with a hIFNγ expressing plasmid under standard conditions and further purified by centrifugation on a sucrose cushion, followed by several steps of sonication and washings with non-denaturing concentrations of urea. The efficiency of the purification was estimated by SDS-PAGE gel electrophoresis and parallel microbiological testing for the presence of residual intact bacteria. Phenol/chloroform extraction showed that the highly purified IBs contain both DNA and RNA. The latter were studied by UV spectroscopy and agarose gel electrophoresis combined with enzymatic treatment and hybridization. DNA was observed as a diffuse fraction mainly in the range of 250 to 1000 bp. RNA isolated by TRIzol^® also demonstrated a substantial molecular heterogeneity. Hybridization with ³²P-labelled oligonucleotides showed that the IBs contain rRNA and are enriched of hIFNγ mRNA.

CONCLUSIONS

The results presented in this study indicate that the nucleic acids might be intrinsic components rather than co-precipitated impurities in the IBs. We assume that the nucleic acids are active participants in the aggregation of recombinant proteins and formation of the IBs that originate from the transcription and translation machinery of the microbial cell factory. Further studies are needed to ascertain this notion.

Aggrescan3D (A3D) 2.0: prediction and engineering of protein solubility

Protein aggregation is a hallmark of a growing number of human disorders and constitutes a major bottleneck in the manufacturing of therapeutic proteins. Therefore, there is a strong need of in-silico methods that can anticipate the aggregative properties of protein variants linked to disease and assist the engineering of soluble protein-based drugs. A few years ago, we developed a method for structure-based prediction of aggregation properties that takes into account the dynamic fluctuations of proteins. The method has been made available as the Aggrescan3D (A3D) web server and applied in numerous studies of protein structure-aggregation relationship. Here, we present a major update of the A3D web server to version 2.0. The new features include: extension of dynamic calculations to significantly larger and multimeric proteins, simultaneous prediction of changes in protein solubility and stability upon mutation, rapid screening for functional protein variants with improved solubility, a REST-ful service to incorporate A3D calculations in automatic pipelines, and a new, enhanced web server interface. A3D 2.0 is freely available at: http://biocomp.chem.uw.edu.pl/A3D2/

Computational Assessment of Bacterial Protein Structures Indicates a Selection Against Aggregation

The aggregation of proteins compromises cell fitness, either because it titrates functional proteins into non-productive inclusions or because it results in the formation of toxic assemblies. Accordingly, computational proteome-wide analyses suggest that prevention of aggregation upon misfolding plays a key role in sequence evolution. Most proteins spend their lifetimes in a folded state; therefore, it is conceivable that, in addition to sequences, protein structures would have also evolved to minimize the risk of aggregation in their natural environments. By exploiting the AGGRESCAN3D structure-based approach to predict the aggregation propensity of >600 Escherichia coli proteins, we show that the structural aggregation propensity of globular proteins is connected with their abundance, length, essentiality, subcellular location and quaternary structure. These data suggest that the avoidance of protein aggregation has contributed to shape the structural properties of proteins in bacterial cells.

Selective Stabilization and Destabilization of Protein Domains in Tissue-Type Plasminogen Activator Using Formulation Excipients

Multidomain biotherapeutic proteins present additional behavioral and analytical challenges for the optimization of their kinetic stability by formulation. Tissue-type plasminogen activator (tPA) comprises six protein domains that exhibit a complex and pH-dependent thermal unfolding profile, due to partially independent domain unfolding. Here we have used tPA as a model for evaluating the relationships between various thermal unfolding and aggregation parameters in multidomain proteins. We show that changes in the thermal unfolding profile of tPA were parametrized by the overall thermal midpoint transition temperature, T_m, and the Van't Hoff entropy for unfolding, Δ S_vh, which is a measure of unfolding cooperativity. The kinetics of degradation at 45 °C, leading to aggregation, were measured as rates of monomer and activity loss. These two rates were found to be coincident at all pH. Aggregation accelerated at pH 4 due to the early unfolding of the serine protease N-terminal domain (SP-N), whereas at pH 5-8, the fraction unfolded at 45 °C ( f₄₅) was <1%, resulting in a baseline rate of aggregation from the native ensemble. We used a Design of Experiments (DoE) approach to evaluate how formulation excipients impact and control the thermal unfolding profile for tPA and found that the relative stability of each of the tPA domains was dependent on the formulation. Therefore, the optimization of formulations for complex multidomain proteins such as tPA may need to be multiobjective, with careful selection of the desired attributes that improve stability. As aggregation rates (ln v) correlated well to T_m ( R² = 0.77) and Δ S_vh ( R² = 0.71) but not T_agg ( R² = 0.01), we analyzed how formulation excipients and pH would be able to optimize T_m and Δ S_vh. Formulation excipient behaviors were found to group according to their combined impact on T_m and Δ S_vh. The effects of each excipient were often selectively stabilizing or destabilizing to specific tPA domains and changed the stability of particular domains relative to the others. The types of mechanism by which this could occur might involve specific interactions with the protein surface, or otherwise effects that are mediated via the solvent as a result of the different surface hydrophobicities and polarities of each domain.

Third generation antibody discovery methods: in silico rational design

Owing to their outstanding performances in molecular recognition, antibodies are extensively used in research and applications in molecular biology, biotechnology and medicine. Recent advances in experimental and computational methods are making it possible to complement well-established in vivo (first generation) and in vitro (second generation) methods of antibody discovery with novel in silico (third generation) approaches. Here we describe the principles of computational antibody design and review the state of the art in this field. We then present Modular, a method that implements the rational design of antibodies in a modular manner, and describe the opportunities offered by this approach.

An Evaluation of the Potential of NMR Spectroscopy and Computational Modelling Methods to Inform Biopharmaceutical Formulations

Protein-based therapeutics are considered to be one of the most important classes of pharmaceuticals on the market. The growing need to prolong stability of high protein concentrations in liquid form has proven to be challenging. Therefore, significant effort is being made to design formulations which can enable the storage of these highly concentrated protein therapies for up to 2 years. Currently, the excipient selection approach involves empirical high-throughput screening, but does not reveal details on aggregation mechanisms or the molecular-level effects of the formulations under storage conditions. Computational modelling approaches have the potential to elucidate such mechanisms, and rapidly screen in silico prior to experimental testing. Nuclear Magnetic Resonance (NMR) spectroscopy can also provide complementary insights into excipient–protein interactions. This review will highlight the underpinning principles of molecular modelling and NMR spectroscopy. It will also discuss the advancements in the applications of computational and NMR approaches in investigating excipient–protein interactions.

Global Protein Stabilization Does Not Suffice to Prevent Amyloid Fibril Formation

Mutations or cellular conditions that destabilize the native protein conformation promote the population of partially unfolded conformations, which in many cases assemble into insoluble amyloid fibrils, a process associated with multiple human pathologies. Therefore, stabilization of protein structures is seen as an efficient way to prevent misfolding and subsequent aggregation. This has been suggested to be the underlying reason why proteins living in harsh environments, such as the extracellular space, have evolved disulfide bonds. The effect of protein disulfides on the thermodynamics and kinetics of folding has been extensively studied, but much less is known on its effect on aggregation reactions. Here, we designed a single point mutation that introduces a disulfide bond in the all-α FF domain, a protein that, despite being devoid of preformed β-sheets, forms β-sheet-rich amyloid fibrils. The novel and unique covalent bond in the FF domain dramatically increases its thermodynamic stability and folding speed. Nevertheless, these optimized properties cannot counteract the inherent aggregation propensity of the protein, thus indicating that a high global protein stabilization does not suffice to prevent amyloid formation unless it contributes to hide from exposure the specific regions that nucleate the aggregation reaction.

Aggregating sequences that occur in many proteins constitute weak spots of bacterial proteostasis

Aggregation is a sequence-specific process, nucleated by short aggregation-prone regions (APRs) that can be exploited to induce aggregation of proteins containing the same APR. Here, we find that most APRs are unique within a proteome, but that a small minority of APRs occur in many proteins. When aggregation is nucleated in bacteria by such frequently occurring APRs, it leads to massive and lethal inclusion body formation containing a large number of proteins. Buildup of bacterial resistance against these peptides is slow. In addition, the approach is effective against drug-resistant clinical isolates of Escherichiacoli and Acinetobacterbaumannii, reducing bacterial load in a murine bladder infection model. Our results indicate that redundant APRs are weak points of bacterial protein homeostasis and that targeting these may be an attractive antibacterial strategy.

Aggregation is sequence-specific and nucleated by short aggregating protein segments (APR). Here authors use a multidisciplinary approach to show that in E.coli some frequently occurring APRs lead to protein aggregation and ultimately bacterial cell death, which could serve as antibacterial strategy.

A single cysteine post-translational oxidation suffices to compromise globular proteins kinetic stability and promote amyloid formation

Oxidatively modified forms of proteins accumulate during aging. Oxidized protein conformers might act as intermediates in the formation of amyloids in age-related disorders. However, it is not known whether this amyloidogenic conversion requires an extensive protein oxidative damage or it can be promoted just by a discrete, localized post-translational modification of certain residues. Here, we demonstrate that the irreversible oxidation of a single free Cys suffices to severely perturb the folding energy landscape of a stable globular protein, compromise its kinetic stability, and lead to the formation of amyloids under physiological conditions. Experiments and simulations converge to indicate that this specific oxidation-promoted protein aggregation requires only local unfolding. Indeed, a large scale analysis indicates that many cellular proteins are at risk of undergoing this kind of deleterious transition; explaining how oxidative stress can impact cell proteostasis and subsequently lead to the onset of pathological states.

•

The population of aggregation-prone states by natural proteins does not require their extensive oxidation.

•

A single residue irreversible oxidation suffices to promote the formation of amyloid fibrils.

•

Under oxidative stress, many cellular proteins are at risk of aggregating into toxic species.

Crystallographic studies on protein misfolding: Domain swapping and amyloid formation in the SH3 domain

Oligomerization by 3D domain swapping is found in a variety of proteins of diverse size, fold and function. In the early 1960s this phenomenon was postulated for the oligomers of ribonuclease A, but it was not until the 1990s that X-ray diffraction provided the first experimental evidence of this special manner of oligomerization. Nowadays, structural information has allowed the identification of these swapped oligomers in over one hundred proteins. Although the functional relevance of this phenomenon is not clear, this alternative folding of protomers into intertwined oligomers has been related to amyloid formation. Studies on proteins that develop 3D domain swapping might provide some clues on the early stages of amyloid formation. The SH3 domain is a small modular domain that has been used as a model to study the basis of protein folding. Among SH3 domains, the c-Src-SH3 domain emerges as a helpful model to study 3D domain swapping and amyloid formation.

In vivo amyloid aggregation kinetics tracked by time-lapse confocal microscopy in real-time

Amyloid polymerization underlies an increasing number of human diseases. Despite this process having been studied extensively in vitro, aggregation is a difficult process to track in vivo due to methodological limitations and the slow kinetics of aggregation reactions in cells and tissues. Herein we exploit the amyloid properties of the inclusions bodies (IBs) formed by amyloidogenic proteins in bacteria to address the kinetics of in vivo amyloid aggregation. To this aim we used time-lapse confocal microscopy and a fusion of the amyloid-beta peptide (A β42) with a fluorescent reporter. This strategy allowed us to follow the intracellular kinetics of amyloid-like aggregation in real-time and to discriminate between variants exhibiting different in vivo aggregation propensity. Overall, the approach opens the possibility to assess the impact of point mutations as well as potential anti-aggregation drugs in the process of amyloid formation in living cells.

Aggregation propensity of neuronal receptors: potential implications in neurodegenerative disorders

Misfolding and aggregation of proteins in tissues is linked to the onset of a diverse set of human neurodegenerative disorders, including Alzheimer's and Parkinson's diseases. In these pathologies proteins usually aggregate into highly ordered and β-sheet enriched amyloid fibrils. However, the formation of these toxic structures is not restricted to a reduced set of polypeptides but rather an intrinsic property of proteins. This suggests that the number of proteins involved in conformational disorders might be much larger than previously thought. The propensity of a protein to form amyloid assemblies is imprinted in its sequence and can be read using computational approaches. Here, we exploit four of these algorithms to analyze the presence of aggregation-prone regions in the sequence and structure of the extracellular domains of several neuroreceptors, with the idea of identifying patches that can interact anomalously with other aggregation-prone molecules such as the amyloid-β peptide or promote their self-assembly. The number of amyloidogenic regions in these domains is rather low but they are significantly exposed to solvent and therefore are suitable for interactions. We find a significant overlap between aggregation-prone regions and receptors interfaces and/or ligand-binding sites, which illustrates an unavoidable competition between the formation of functional native interactions and that of dangerous amyloid-like contacts leading to disease.

The lower limits for protein stability and foldability in primary hyperoxaluria type I

Mutational effects on protein stability and foldability are important to understand conformational diseases and protein evolution. In this work, we perform a comprehensive investigation on the energetic basis underlying mutational effects on the stability of human alanine:glyoxylate aminotransferase (AGT). We study twenty two variants whose kinetic stabilities span over eleven orders of magnitude and are classified into two groups: i) ten naturally-occurring variants, including the most common mutations causing primary hyperoxaluria type I (PH1); and ii) twelve consensus variants obtained by sequence-alignment statistics. We show that AGT dimer stability determines denaturation rates, and mutations modulate stability by changes in the effective thermodynamic stability, the aggregation propensity of partially/globally unfolded states and subtle energetic changes in the rate-limiting denaturation step. In combination with our previous expression analyses in eukaryotic cells, we propose the existence of two lower limits for AGT stability, one linked to optimal folding efficiency (close to the major allele stability) and the other setting a minimal efficiency compatible with glyoxylate detoxification in vivo (close to the minor allele stability). These lower limits could explain the high prevalence of misfolding as a disease mechanism in PH1 and support the use of pharmacological ligands aimed to increase AGT stability as therapies for this disease.

Cytosolic selection systems to study protein stability

Here we describe biosensors that provide readouts for protein stability in the cytosolic compartment of prokaryotes. These biosensors consist of tripartite sandwich fusions that link the in vitro stability or aggregation susceptibility of guest proteins to the in vivo resistance of host cells to the antibiotics kanamycin, spectinomycin, and nourseothricin. These selectable markers confer antibiotic resistance in a wide range of hosts and are easily quantifiable. We show that mutations within guest proteins that affect their stability alter the antibiotic resistances of the cells expressing the biosensors in a manner that is related to the in vitro stabilities of the mutant guest proteins. In addition, we find that polyglutamine tracts of increasing length are associated with an increased tendency to form amyloids in vivo and, in our sandwich fusion system, with decreased resistance to aminoglycoside antibiotics. We demonstrate that our approach allows the in vivo analysis of protein stability in the cytosolic compartment without the need for prior structural and functional knowledge.

Influence of pH control in the formation of inclusion bodies during production of recombinant sphingomyelinase-D in Escherichia coli

BACKGROUND

Inclusion bodies (IBs) are aggregated proteins that form clusters when protein is overexpressed in heterologous expression systems. IBs have been considered as non-usable proteins, but recently they are being used as functional materials, catalytic particles, drug delivery agents, immunogenic structures, and as a raw material in recombinant therapeutic protein purification. However, few studies have been made to understand how culture conditions affect the protein aggregation and the physicochemical characteristics that lead them to cluster. The objective of our research was to understand how pH affects the physicochemical properties of IBs formed by the recombinant sphingomyelinase-D of tick expressed in E. coli BL21-Gold (DE3) by evaluating two pH culture strategies.

RESULTS

Uncontrolled pH culture conditions favored recombinant sphingomyelinase-D aggregation and IB formation. The IBs of sphingomyelinase-D produced under controlled pH at 7.5 and after 24 h were smaller (<500 nm) than those produced under uncontrolled pH conditions (>500 nm). Furthermore, the composition, conformation and β-structure formation of the aggregates were different. Under controlled pH conditions in comparison to uncontrolled conditions, the produced IBs presented higher resistance to denaturants and proteinase-K degradation, presented β-structure, but apparently as time passes the IBs become compacted and less sensitive to amyloid dye binding.

CONCLUSIONS

The manipulation of the pH has an impact on IB formation and their physicochemical characteristics. Particularly, uncontrolled pH conditions favored the protein aggregation and sphingomyelinase-D IB formation. The evidence may lead to find methodologies for bioprocesses to obtain biomaterials with particular characteristics, extending the application possibilities of the inclusion bodies.

Fluorescent dye ProteoStat to detect and discriminate intracellular amyloid-like aggregates in Escherichia coli

The formation of amyloid aggregates is linked to the onset of an increasing number of human disorders. Thus, there is an increasing need for methodologies able to provide insights into protein deposition and its modulation. Many approaches exist to study amyloids in vitro, but the techniques available for the study of amyloid aggregation in cells are still limited and non-specific. In this study we developed a methodology for the detection of amyloid-like aggregates inside cells that discriminates these ordered assemblies from other intracellular aggregates. We chose bacteria as model system, since the inclusion bodies formed by amyloid proteins in the cytosol of bacteria resemble toxic amyloids both structurally and functionally. Using confocal microscopy, fluorescence spectroscopy, and flow cytometry, we show that the recently developed red fluorescent dye ProteoStat can detect the presence of intracellular amyloid-like deposits in living bacterial cells with high specificity, even when the target proteins are expressed at low levels. This methodology allows quantitation of the intracellular amyloid content, shows the potential to replace in vitro screenings in the search for therapeutic anti-amyloidogenic compounds, and might be useful for identifying conditions that prevent the aggregation of therapeutic recombinant proteins.

Association between foldability and aggregation propensity in small disulfide-rich proteins

AIMS

Disulfide-rich domains (DRDs) are small proteins whose native structure is stabilized by the presence of covalent disulfide bonds. These domains are versatile and can perform a wide range of functions. Many of these domains readily unfold on disulfide bond reduction, suggesting that in the absence of covalent bonding they might display significant disorder.

RESULTS

Here, we analyzed the degree of disorder in 97 domains representative of the different DRDs families and demonstrate that, in terms of sequence, many of them can be classified as intrinsically disordered proteins (IDPs) or contain predicted disordered regions. The analysis of the aggregation propensity of these domains indicates that, similar to IDPs, their sequences are more soluble and have less aggregating regions than those of other globular domains, suggesting that they might have evolved to avoid aggregation after protein synthesis and before they can attain its compact and covalently linked native structure.

INNOVATION AND CONCLUSION

DRDs, which resemble IDPs in the reduced state and become globular when their disulfide bonds are formed, illustrate the link between protein folding and aggregation propensities and how these two properties cannot be easily dissociated, determining the main traits of the folding routes followed by these small proteins to attain their native oxidized states.

Mapping the structure of amyloid nucleation precursors by protein engineering kinetic analysis

Understanding the early molecular mechanisms governing amyloid aggregation is crucial to learn how to prevent it. Here, we used a site-directed mutagenesis approach to explore the molecular mechanism of nucleation of amyloid structure in the N47A Spc-SH3 domain. The changes in the native state stability produced by a series of mutations on each structural element of the domain were uncorrelated with the changes in the aggregation rates, although the overall aggregation mechanism was not altered. Analysis of the thioflavin T initial rates based on a simple kinetic model allowed us to extract thermodynamic magnitudes of the precursor states of nucleation and map the regions of the protein participating in the structure of the amyloidogenic precursors. This structure differs from that of the folding transition state of the SH3 domains, strongly suggesting that the regions of the conformational landscape leading to amyloid formation are divergent from those leading to the native fold.

Protein aggregation propensity is a crucial determinant of intracellular inclusion formation and quality control degradation

Protein aggregation is linked to many pathological conditions, including several neurodegenerative diseases. The aggregation propensities of proteins are thought to be controlled to a large extent by the physicochemical properties encoded in the primary sequence. We have previously exploited a set of amyloid β peptide (Aβ42) variants exhibiting a continuous gradient of intrinsic aggregation propensities to demonstrate that this rule applies in vivo in bacteria. In the present work we have characterized the behavior of these Aβ42 mutants when expressed in yeast. In contrast to bacteria, the intrinsic aggregation propensity is gated by yeast, in such a way that this property correlates with the formation of intracellular inclusions only above a specific aggregation threshold. Proteins displaying solubility levels above this threshold escape the inclusion formation pathway. In addition, the most aggregation-prone variants are selectively cleared by the yeast quality control degradation machinery. Thus, both inclusion formation and proteolysis target the same aggregation-prone variants and cooperate to minimize the presence of these potentially dangerous species in the cytosol. The demonstration that sorting to these pathways in eukaryotes is strongly influenced by protein primary sequence should facilitate the development of rational approaches to predict and hopefully prevent in vivo protein deposition.

Microcanonical thermostatistics of coarse-grained proteins with amyloidogenic propensity

The formation of fibrillar aggregates seems to be a common characteristic of polypeptide chains, although the observation of these aggregates may depend on appropriate experimental conditions. Partially folded intermediates seem to have an important role in the generation of protein aggregates, and a mechanism for this fibril formation considers that these intermediates also correspond to metastable states with respect to the fibrillar ones. Here, using a coarse-grained (CG) off-lattice model, we carry out a comparative analysis of the thermodynamic aspects characterizing the folding transition with respect to the propensity for aggregation of four different systems: two isoforms of the amyloid β-protein, the Src SH3 domain, and the human prion proteins (hPrP). Microcanonical analysis of the data obtained from replica exchange method is conducted to evaluate the free-energy barrier and latent heat in these models. The simulations of the amyloid β isoforms and Src SH3 domain indicated that the folding process described by this CG model is related to a negative specific heat, a phenomenon that can only be verified in the microcanonical ensemble in first-order phase transitions. The CG simulation of the hPrP heteropolymer yielded a continuous folding transition. The absence of a free-energy barrier and latent heat favors the presence of partially unfolded conformations, and in this context, this thermodynamic aspect could explain the reason why the hPrP heteropolymer is more aggregation-prone than the other heteropolymers considered in this study. We introduced the hydrophobic radius of gyration as an order parameter and found that it can be used to obtain reliable information about the hydrophobic packing and the transition temperatures in the folding process.

Thioflavin-S staining coupled to flow cytometry. A screening tool to detect in vivo protein aggregation

Amyloid deposits are associated with an increasing number of human disorders, including Alzheimer's and Parkinson's diseases. Recent studies provide compelling evidence for the existence of amyloid-like conformations in the insoluble bacterial inclusion bodies (IBs) produced during the recombinant expression of amyloidogenic proteins. This makes prokaryotic cells a physiologically relevant system to study the mechanisms of in vivo amyloid deposition. We show here that the application of flow cytometry to detect Thioflavin-S (Th-S) fluorescence provides a fast, robust, quantitative, non-invasive method to screen for the presence of in vivo intracellular amyloid-like aggregates in bacteria, with potential application in the analysis of the impact of genetic mutations or chemical compounds on the aggregation of disease-associated polypeptides.

Thermodynamic stability, unfolding kinetics, and aggregation of the N-terminal actin-binding domains of utrophin and dystrophin

Muscular dystrophy (MD) is the most common genetic lethal disorder in children. Mutations in dystrophin trigger the most common form of MD, Duchenne, and its allelic variant Becker MD. Utrophin is the closest homologue and has been shown to compensate for the loss of dystrophin in human disease animal models. However, the structural and functional similarities and differences between utrophin and dystrophin are less understood. Both proteins interact with actin through their N-terminal actin-binding domain (N-ABD). In this study, we examined the thermodynamic stability and aggregation of utrophin N-ABD and compared with that of dystrophin. Our results show that utrophin N-ABD has spectroscopic properties similar to dystrophin N-ABD. However, utrophin N-ABD has decreased denaturant and thermal stability, unfolds faster, and is correspondingly more susceptible to proteolysis, which might account for its decreased in vivo half-life compared to dystrophin. In addition, utrophin N-ABD aggregates to a lesser extent compared with dystrophin N-ABD, contrary to the general behavior of proteins in which decreased stability enhances protein aggregation. Despite these differences in stability and aggregation, both proteins exhibit deleterious effects of mutations. When utrophin N-ABD mutations analogous in position to the dystrophin disease-causing mutations were generated, they behaved similarly to dystrophin mutants in terms of decreased stability and the formation of cross-β aggregates, indicating a possible role for utrophin mutations in disease mechanisms.

AGGRESCAN: method, application, and perspectives for drug design

Protein aggregation underlies the development of an increasing number of conformational human diseases of growing incidence, such as Alzheimer's and Parkinson's diseases. Furthermore, the accumulation of recombinant proteins as intracellular aggregates represents a critical obstacle for the biotechnological production of polypeptides. Also, ordered protein aggregates constitute novel and versatile nanobiomaterials. Consequently, there is an increasing interest in the development of methods able to forecast the aggregation properties of polypeptides in order to modulate their intrinsic solubility. In this context, we have developed AGGRESCAN, a simple and fast algorithm that predicts aggregation-prone segments in protein sequences, compares the aggregation properties of different proteins or protein sets and analyses the effect of mutations on protein aggregation propensities.

Sequence-based prediction of protein solubility

In order to investigate the relationship between the thermodynamics and kinetics of protein aggregation, we compared the solubility of proteins with their aggregation rates. We found a significant correlation between these two quantities by considering a database of protein solubility values measured using an in vitro reconstituted translation system containing about 70% of Escherichia coli proteins. The existence of such correlation suggests that the thermodynamic stability of the native states of proteins relative to the aggregate states is closely linked with the kinetic barriers that separate them. In order to create the possibility of conducting computational studies at the proteome level to investigate further this concept, we developed a method of predicting the solubility of proteins based on their physicochemical properties.

Deciphering the role of the thermodynamic and kinetic stabilities of SH3 domains on their aggregation inside bacteria

The formation of insoluble deposits by globular proteins underlies the onset of many human diseases. Recent studies suggest a relationship between the thermodynamic stability of proteins and their in vivo aggregation. However, it has been argued that, in the cell, the occurrence of irreversible aggregation might shift the system from equilibrium, in such a way that it could be the rate of unfolding and associated kinetic stability instead of the conformational stability that controls protein deposition. This is an important but difficult to decipher question, because kinetic and thermodynamic stabilities appear usually correlated. Here we address this issue by comparing the in vitro folding kinetics and stability features of a set of non-natural SH3 domains with their aggregation properties when expressed in bacteria. In addition, we compare the in vitro stability of the isolated domains with their effective stability in conditions that mimic the cytosolic environment. Overall, the data argue in favor of a thermodynamic rather than a kinetic control of the intracellular aggregation propensities of small globular proteins in which folding and unfolding velocities largely exceed aggregation rates. These results have implications regarding the evolution of proteins.

Influence of the stability of a fused protein and its distance to the amyloidogenic segment on fibril formation

Conversion of native proteins into amyloid fibrils is irreversible and therefore it is difficult to study the interdependence of conformational stability and fibrillation by thermodynamic analyses. Here we approached this problem by fusing amyloidogenic poly-alanine segments derived from the N-terminal domain of the nuclear poly (A) binding protein PABPN1 with a well studied, reversibly unfolding protein, CspB from Bacillus subtilis. Earlier studies had indicated that CspB could maintain its folded structure in fibrils, when it was separated from the amyloidogenic segment by a long linker. When CspB is directly fused with the amyloidogenic segment, it unfolds because its N-terminal chain region becomes integrated into the fibrillar core, as shown by protease mapping experiments. Spacers of either 3 or 16 residues between CspB and the amyloidogenic segment were not sufficient to prevent this loss of CspB structure. Since the low thermodynamic stability of CspB (ΔG(D) = 12.4 kJ/mol) might be responsible for unfolding and integration of CspB into fibrils, fusions with a CspB mutant with enhanced thermodynamic stability (ΔG(D) = 26.9 kJ/mol) were studied. This strongly stabilized CspB remained folded and prevented fibril formation in all fusions. Our data show that the conformational stability of a linked, independently structured protein domain can control fibril formation.

Protein folding and aggregation in bacteria

Proteins might experience many conformational changes and interactions during their lifetimes, from their synthesis at ribosomes to their controlled degradation. Because, in most cases, only folded proteins are functional, protein folding in bacteria is tightly controlled genetically, transcriptionally, and at the protein sequence level. In addition, important cellular machinery assists the folding of polypeptides to avoid misfolding and ensure the attainment of functional structures. When these redundant protective strategies are overcome, misfolded polypeptides are recruited into insoluble inclusion bodies. The protein embedded in these intracellular deposits might display different conformations including functional and beta-sheet-rich structures. The latter assemblies are similar to the amyloid fibrils characteristic of several human neurodegenerative diseases. Interestingly, bacteria exploit the same structural principles for functional properties such as adhesion or cytotoxicity. Overall, this review illustrates how prokaryotic organisms might provide the bedrock on which to understand the complexity of protein folding and aggregation in the cell.

Optimizing protein stability in vivo

Identifying mutations that stabilize proteins is challenging because most substitutions are destabilizing. In addition to being of immense practical utility, the ability to evolve protein stability in vivo may indicate how evolution has formed today's protein sequences. Here we describe a genetic selection that directly links the in vivo stability of proteins to antibiotic resistance. It allows the identification of stabilizing mutations within proteins. The large majority of mutants selected for improved antibiotic resistance are stabilized both thermodynamically and kinetically, indicating that similar principles govern stability in vivo and in vitro. The approach requires no prior structural or functional knowledge and allows selection for stability without a need to maintain function. Mutations that enhance thermodynamic stability of the protein Im7 map overwhelmingly to surface residues involved in binding to colicin E7, showing how the evolutionary pressures that drive Im7-E7 complex formation have compromised the stability of the isolated Im7 protein.

A single mutation in an SH3 domain increases amyloid aggregation by accelerating nucleation, but not by destabilizing thermodynamically the native state

We investigated the relationship between thermodynamic stability and amyloid aggregation propensity for a set of single mutants of the alpha-spectrin SH3 domain (Spc-SH3). Whilst mutations destabilizing the domain at position 56 did not enhance fibrillation, the N47A mutation increased the rate of amyloid fibril formation by 10-fold. Even under conditions of identical thermodynamic stability, the aggregation rate was much higher for the N47A mutant than for the WT domain. We conclude that the N47A mutation does not change the apparent mechanism of fibrillation or the morphology of the amyloid fibrils, and that its amyloidogenic property is due to its effect upon the rate of the conformational events leading to nucleation and not to its overall destabilizing effect.