Structural Reorganization of α-Synuclein at Low pH Observed by NMR and REMD Simulations

α-Synuclein is an intrinsically disordered protein that appears in aggregated forms in the brains of patients with Parkinson's disease. The conversion from monomer to aggregate is complex, and aggregation rates are sensitive to changes in amino acid sequence and environmental conditions. It has previously been observed that α-synuclein aggregates faster at low pH than at neutral pH. Here, we combine NMR spectroscopy and molecular simulations to characterize α-synuclein conformational ensembles at both neutral and low pH in order to understand how the altered charge distribution at low pH changes the structural properties of these ensembles and leads to an increase in aggregation rate. The N-terminus, which has a small positive charge at neutral pH due to a balance of positively and negatively charged amino acid residues, is very positively charged at low pH. Conversely, the acidic C-terminus is highly negatively charged at neutral pH and becomes essentially neutral and hydrophobic at low pH. Our NMR experiments and replica exchange molecular dynamics simulations indicate that there is a significant structural reorganization within the low-pH ensemble relative to that at neutral pH in terms of long-range contacts, hydrodynamic radius, and the amount of heterogeneity within the conformational ensembles. At neutral pH, there is a very heterogeneous ensemble with transient contacts between the N-terminus and the non-amyloid β component (NAC); however, at low pH, there is a more homogeneous ensemble that exhibits strong contacts between the NAC and the C-terminus. At both pH values, transient contacts between the N- and C-termini are observed, the NAC region shows similar exposure to solvent, and the entire protein shows similar propensities to secondary structure. Based on the comparison of the neutral- and low-pH conformational ensembles, we propose that exposure of the NAC region to solvent and the secondary-structure propensity are not factors that account for differences in propensity to aggregate in this context. Instead, the comparison of the neutral- and low-pH ensembles suggests that the change in long-range interactions between the low- and neutral-pH ensembles, the compaction of the C-terminal region at low pH, and the uneven distribution of charges across the sequence are key to faster aggregation.     

Residual structure, backbone dynamics, and interactions within the synuclein family

The human synuclein protein family includes alpha-synuclein, which has been linked to both familial and sporadic Parkinson's disease, and the highly homologous beta and gamma-synuclein. Mutations in alpha-synuclein cause autosomal dominant early onset Parkinson's, and the protein is found deposited in a fibrillar form in hereditary and idiopathic forms of the disease. No genetic link between beta and gamma-synuclein, and any neurodegenerative disease has been established, and it is generally considered that these proteins are not highly pathogenic. In addition, beta and gamma-synuclein are reported to aggregate less readily than alpha-synuclein in vitro. Indeed, beta-synuclein has been reported to protect against alpha-synuclein aggregation in vitro, as well as alpha-synuclein-mediated toxicity in vivo. Earlier, we compared the structural properties of the highly helical states adopted by all three synucleins in association with detergent micelles in an attempt to delineate the basis for functional differences between the three proteins. Here, we report a comparison of the structural and dynamic properties of the free states of all three proteins in order to shed light on differences that may help to explain their different propensities to aggregate, which in turn may underlie their differing contributions to the etiology of Parkinson's disease. We find that gamma-synuclein closely resembles alpha-synuclein in its free-state residual secondary structure, consistent with the more similar propensities of the two proteins to aggregate in vitro. beta-Synuclein, however, differs significantly from alpha-synuclein, exhibiting a lower predisposition towards helical structure in the second half of its lipid-binding domain, and a higher preference for extended structures in its C-terminal tail. Both beta and gamma-synuclein show less extensive transient long-range structure than that observed in alpha-synuclein. These results raise questions regarding the role of secondary structure propensities and transient long-range contacts in directing synuclein aggregation reactions.     

Solution conformation, backbone dynamics and lipid interactions of the intrinsically unstructured malaria surface protein MSP2

Merozoite surface protein 2 (MSP2), one of the most abundant proteins on the surface of the merozoite stage of Plasmodium falciparum, is a potential component of a malaria vaccine, having shown some efficacy in a clinical trial in Papua New Guinea. MSP2 is a GPI-anchored protein consisting of conserved N- and C-terminal domains and a variable central region. Previous studies have shown that it is an intrinsically unstructured protein with a high propensity for fibril formation, in which the conserved N-terminal domain has a key role. Secondary structure predictions suggest that MSP2 contains long stretches of random coil with very little α-helix or β-strand. Circular dichroism spectroscopy confirms this prediction under physiological conditions (pH 7.4) and in more acidic solutions (pH 6.2 and 3.4). Pulsed field gradient NMR diffusion measurements showed that MSP2 under physiological conditions has a large effective hydrodynamic radius consistent with an intrinsic pre-molten globule state, as defined by Uversky. This was supported by sedimentation velocity studies in the analytical ultracentrifuge. NMR resonance assignments have been obtained for FC27 MSP2, allowing the residual secondary structure and backbone dynamics to be defined. There is some motional restriction in the conserved C-terminal region in the vicinity of an intramolecular disulfide bond. Two other regions show motional restrictions, both of which display helical structure propensities. One of these helical regions is within the conserved N-terminal domain, which adopts essentially the same conformation in full-length MSP2 as in corresponding peptide fragments. We see no evidence of long-range interactions in the full-length protein. MSP2 associates with lipid micelles, but predominantly through the N-terminal region rather than the C terminus, which is GPI-anchored to the membrane in the parasite.     

Helix periodicity, topology, and dynamics of membrane-associated alpha-synuclein

The protein alpha-Synuclein (aS) is a synaptic vesicle-associated regulator of synaptic strength and dopamine homeostasis with a pathological role in Parkinson's disease. The normal function of aS depends on a membrane-associated conformation that is adopted upon binding to negatively charged lipid surfaces. Previously we found that the membrane-binding domain of aS is helical and suggested that it may exhibit an unusual structural periodicity. Here we present a study of the periodicity, topology, and dynamics of detergent micelle-bound aS using paramagnetic spin labels embedded in the micelle or attached to the protein. We show that the helical region of aS completes three full turns every 11 residues, demonstrating the proposed 11/3 periodicity. We also find that the membrane-binding domain is partially buried in the micelle surface and bends toward the hydrophobic interior, but does not traverse the micelle. Deeper submersion of certain regions within the micelle, including the unique lysine-free sixth 11-residue repeat, is observed and may be functionally important. There are no long-range tertiary contacts within this domain, indicating a highly extended configuration. The backbone dynamics of the micelle-bound region are relatively uniform with a slight decrease in flexibility observed toward the C-terminal end. These results clarify the topological features of aS bound to membrane-mimicking detergent micelles, with implications for aS function and pathology.     

Effect of Pseudorepeat Rearrangement on α-Synuclein Misfolding, Vesicle Binding, and Micelle Binding

The pathological and physiological hallmarks of the protein α-synuclein (aS) are its misfolding into cytotoxic aggregates and its binding to synaptic vesicles, respectively. Both events are mediated by seven 11-residue amphiphilic pseudorepeats and, most generally, involve a transition from intrinsically unstructured conformations to structured conformations. Based on aS interactions with aggregation-inhibiting small molecules, an aS variant termed shuffled α-synuclein (SaS), wherein the first six pseudorepeats had been rearranged, was introduced. Here, the effects of this rearrangement on misfolding, vesicle binding, and micelle binding are examined in reference to aS and β-synuclein to study the sequence characteristics underlying these processes. Fibrillization correlates with the distinct clustering of residues with high β-sheet propensities, while vesicle affinities depend on the mode of pseudorepeat interchange and loss. In the presence of micelles, the pseudorepeat region of SaS adopts an essentially continuous helix, whereas aS and β-synuclein encounter a distinct helix break, indicating that a more homogeneous distribution of surfactant affinities in SaS prevented the formation of an extensive helix break in the micelle-bound state. By demonstrating the importance of the distribution of β-sheet propensities and by revealing inhomogeneous aS surfactant affinities, the present study provides novel insights into two central themes of synuclein biology.     

E46K Parkinson’s-Linked Mutation Enhances C-Terminal-to-N-Terminal Contacts in α-Synuclein

Parkinson's disease (PD) is associated with the deposition of fibrillar aggregates of the protein α-synuclein (αS) in neurons. Intramolecular contacts between the acidic C-terminal tail of αS and its N-terminal region have been proposed to regulate αS aggregation, and two originally described PD mutations, A30P and A53T, reportedly reduce such contacts. We find that the most recently discovered PD-linked αS mutation E46K, which also accelerates the aggregation of the protein, does not interfere with C-terminal-to-N-terminal contacts and instead enhances such contacts. Furthermore, we do not observe a substantial reduction in such contacts in the two previously characterized mutants. Our results suggest that C-terminal-to-N-terminal contacts in αS are not strongly protective against aggregation, and that the dominant mechanism by which PD-linked mutations facilitate αS aggregation may be altering the physicochemical properties of the protein such as net charge (E46K) and secondary structure propensity (A30P and A53T).     

The role of C-terminal part of ghrelin in pharmacokinetic profile and biological activity in rats

Ghrelin is an endogenous ligand for growth hormone secretagogue receptor 1a (GHS-R1a), and consists of 28 amino acid residues with octanoyl modification at Ser³. The previous studies have revealed that N-terminal part of ghrelin including modified Ser³ is the active core for the activation of GHS-R1a. On the other hand, the role of C-terminal (8-28) region in ghrelin has not been clarified yet. In the present study, we prepared human ghrelin, C-terminal truncated ghrelin derivatives and anamorelin, a small molecular GHS compound which supposedly mimics the N-terminal active core, and examined GHS-R1a agonist activity in vitro, pharmacokinetic (PK) profile and growth hormone (GH) releasing activity in rats. All compounds demonstrated potent GHS-R1a agonist activities in vitro. Although the lack of C-terminal two amino acids did not modify PK profile and GH releasing activity, the deletion of C-terminal 8 and 20 amino acids affected them, and ghrelin(1-7)-Lys-NH₂ exhibited very short plasma half-life and low GH releasing activity in vivo. In rat plasma, ghrelin(1-7)-Lys-NH₂ was degraded more rapidly than ghrelin, suggesting that C-terminal part of ghrelin protected octanoylation of Ser³ from plasma esterases. Subdiaphragmatic vagotomy significantly attenuated GH response to ghrelin but not to anamorelin. These results suggest that the C-terminal part of ghrelin has an important role in the biological activity in vivo. We also found that ghrelin stimulated GH release mainly via a vagal nerve pathway but anamorelin augmented GH release possibly by directly acting on brain in rats.     

Homodimerization propensity of the intrinsically disordered N-terminal domain of Ultraspiracle from Aedes aegypti

The mosquito Aedes aegypti is the principal vector of dengue, one of the most devastating arthropod-borne viral infections in humans. The isoform specific A/B region, called the N-terminal domain (NTD), is hypervariable in sequence and length and is poorly conserved within the Ultraspiracle (Usp) family. The Usp protein together with ecdysteroid receptor (EcR) forms a heterodimeric complex. Up until now, there has been little data on the molecular properties of the isolated Usp-NTD. Here, we describe the biochemical and biophysical properties of the recombinant NTD of the Usp isoform B (aaUsp-NTD) from A. aegypti. These results, in combination with in silico bioinformatics approaches, indicate that aaUsp-NTD exhibits properties of an intrinsically disordered protein (IDP). We also present the first experimental evidence describing the dimerization propensity of the isolated NTD of Usp. These characteristics also appear for other members of the Usp family in different species, for example, in the Usp-NTD from Drosophila melanogaster and Bombyx mori. However, aaUsp-NTD exhibits the strongest homodimerization potential. We postulate that the unique dimerization of the NTD might be important for Usp function by providing an additional platform for interactions, in addition to the nuclear receptor superfamily dimerization via DNA binding domains and ligand binding domains that has already been extensively documented. Furthermore, the unique NTD–NTD interaction that was observed might contribute new insight into the dimerization propensities of nuclear receptors.     

Sensitivity of secondary structure propensities to sequence differences between alpha- and gamma-synuclein: implications for fibrillation

The synucleins are a family of intrinsically disordered proteins involved in various human diseases. alpha-Synuclein has been extensively characterized due to its role in Parkinson's disease where it forms intracellular aggregates, while gamma-synuclein is overexpressed in a majority of late-stage breast cancers. Despite fairly strong sequence similarity between the amyloid-forming regions of alpha- and gamma-synuclein, gamma-synuclein has only a weak propensity to form amyloid fibrils. We hypothesize that the different fibrillation tendencies of alpha- and gamma-synuclein may be related to differences in structural propensities. Here we have measured chemical shifts for gamma-synuclein and compared them to previously published shifts for alpha-synuclein. In order to facilitate direct comparison, we have implemented a simple new technique for re-referencing chemical shifts that we have found to be highly effective for both disordered and folded proteins. In addition, we have developed a new method that combines different chemical shifts into a single residue-specific secondary structure propensity (SSP) score. We observe significant differences between alpha- and gamma-synuclein secondary structure propensities. Most interestingly, gamma-synuclein has an increased alpha-helical propensity in the amyloid-forming region that is critical for alpha-synuclein fibrillation, suggesting that increased structural stability in this region may protect against gamma-synuclein aggregation. This comparison of residue-specific secondary structure propensities between intrinsically disordered homologs highlights the sensitivity of transient structure to sequence changes, which we suggest may have been exploited as an evolutionary mechanism for fast modulation of protein structure and, hence, function.     

Solution structure and backbone dynamics of human Raf-1 kinase inhibitor protein

Human Raf-1 kinase inhibitor protein (hRKIP) is a small multi-functional protein of 187 residues. It contains a conserved pocket, which binds a wide range of ligands from various small molecules to distinct proteins. To provide a structural basis for the ligand diversity of RKIP, we herein determined the solution structure of hRKIP, and analyzed its structural dynamics. In solution, hRKIP mainly comprises two antiparallel β sheets, two α helices and two 3₁₀ helices. NMR dynamic analysis reveals that the overall structure of hRKIP is rigid, but its C-terminal helix which is close to the ligand-binding site is mobile. In addition, residues around the ligand-binding pocket exhibit significant conformational exchange on the μs–ms timescale. Conformational flexibility may allow the ligand-binding pocket and the C-terminal helix to adopt various conformations to interact with different substrates. This work may shed light on the underlying molecular mechanisms of how hRKIP recognizes and binds diverse substrate ligands.     

Molecular basis of dynamic relocalization of Dictyostelium myosin IB

Class I myosins have a single heavy chain comprising an N-terminal motor domain with actin-activated ATPase activity and a C-terminal globular tail with a basic region that binds to acidic phospholipids. These myosins contribute to the formation of actin-rich protrusions such as pseudopodia, but regulation of the dynamic localization to these structures is not understood. Previously, we found that Acanthamoeba myosin IC binds to acidic phospholipids in vitro through a short sequence of basic and hydrophobic amino acids, BH site, based on the charge density of the phospholipids. The tail of Dictyostelium myosin IB (DMIB) also contains a BH site. We now report that the BH site is essential for DMIB binding to the plasma membrane and describe the molecular basis of the dynamic relocalization of DMIB in live cells. Endogenous DMIB is localized uniformly on the plasma membrane of resting cells, at active protrusions and cell-cell contacts of randomly moving cells, and at the front of motile polarized cells. The BH site is required for association of DMIB with the plasma membrane at all stages where it colocalizes with phosphoinositide bisphosphate/phosphoinositide trisphosphate (PIP(2)/PIP(3)). The charge-based specificity of the BH site allows for in vivo specificity of DMIB for PIP(2)/PIP(3) similar to the PH domain-based specificity of other class I myosins. However, DMIB-head is required for relocalization of DMIB to the front of migrating cells. Motor activity is not essential, but the actin binding site in the head is important. Thus, dynamic relocalization of DMIB is determined principally by the local PIP(2)/PIP(3) concentration in the plasma membrane and cytoplasmic F-actin.     

Insights into nucleic acid conformational dynamics from massively parallel stochastic simulations

The helical hairpin is one of the most ubiquitous and elementary secondary structural motifs in nucleic acids, capable of serving functional roles and participating in long-range tertiary contacts. Yet the self-assembly of these structures has not been well-characterized at the atomic level. With this in mind, the dynamics of nucleic acid hairpin formation and disruption have been studied using a novel computational tool: large-scale, parallel, atomistic molecular dynamics simulation employing an inhomogeneous distributed computer consisting of more than 40,000 processors. Using multiple methodologies, over 500 micro s of atomistic simulation time has been collected for a large ensemble of hairpins (sequence 5'-GGGC[GCAA]GCCU-3'), allowing characterization of rare events not previously observable in simulation. From uncoupled ensemble dynamics simulations in unperturbed folding conditions, we report on 1), competing pathways between the folded and unfolded regions of the conformational space; 2), observed nonnative stacking and basepairing traps; and 3), a helix unwinding-rewinding mode that is differentiated from the unfolding and folding dynamics. A heterogeneous transition state ensemble is characterized structurally through calculations of conformer-specific folding probabilities and a multiplexed replica exchange stochastic dynamics algorithm is used to derive an approximate folding landscape. A comparison between the observed folding mechanism and that of a peptide beta-hairpin analog suggests that although native topology defines the character of the folding landscape, the statistical weighting of potential folding pathways is determined by the chemical nature of the polymer.     

Single proline residues can dictate the oxidative folding pathways of cysteine-rich peptides

The cysteine-rich N and C-terminal domains of minicollagen-1 from Hydra nematocysts fold with excesses of oxidized/reduced glutathione (10:1) into globular structures with distinct cystine frameworks despite their identical cysteine sequence pattern. An additional main difference is the cis conformation of a conserved proline residue in the N-terminal and the trans conformation of this residue in the C-terminal domain. Comparative analysis of the oxidative folding revealed for the C-terminal domain a fast and highly cooperative formation of a single disulfide isomer. Conversely, oxidation of the N-terminal domain proceeds mainly via an intermediate that results from the fast quasi-stochastic disulfide formation according to the proximity rule. The rate of conversion of the bead-like isomer into the globular end-product is largely dominated by the trans-to-cis isomerization of the critical proline residue as well assessed by its replacement with (4R)- and (4S)-fluoroproline known to exhibit distinct propensities for the trans and cis conformation, respectively. Independently, whether the trans or cis conformation is favored by these substitutions, both analogues retain sufficient sequence-encoded information to fold almost quantitatively into the identical cystine framework and thus spatial structure of the parent peptide with the critical proline residue as cis isomer, but at rates significantly lower for the (4R) than for the (4S)-fluoroproline analogue. Correspondingly, other sequence-encoded structural elements have to act as a driving force for these unidirectional folding pathways despite the rather simple sequence composition consisting only of aliphatic residues, some proline and only one aromatic residue (tyrosine) in the core parts of the C and N-terminal domains. The two cysteine-rich domains of minicollagen-1 may well represent ideal targets for ab initio structure calculations in order to learn more about the elementary information encoded in such primordial molecules.     

Structural and functional characterization of the Geobacillus copper nitrite reductase: Involvement of the unique N-terminal region in the interprotein electron transfer with its redox partner

The crystal structures of copper-containing nitrite reductase (CuNiR) from the thermophilic Gram-positive bacterium Geobacillus kaustophilus HTA426 and the amino (N)-terminal 68 residue-deleted mutant were determined at resolutions of 1.3 Å and 1.8 Å, respectively. Both structures show a striking resemblance with the overall structure of the well-known CuNiRs composed of two Greek key β-barrel domains; however, a remarkable structural difference was found in the N-terminal region. The unique region has one β-strand and one α-helix extended to the northern surface of the type-1 copper site. The superposition of the Geobacillus CuNiR model on the electron-transfer complex structure of CuNiR with the redox partner cytochrome c₅₅₁ in other denitrifier system led us to infer that this region contributes to the transient binding with the partner protein during the interprotein electron transfer reaction in the Geobacillus system. Furthermore, electron-transfer kinetics experiments using N-terminal residue-deleted mutant and the redox partner, Geobacillus cytochrome c₅₅₁, were carried out. These structural and kinetics studies demonstrate that the region is directly involved in the specific partner recognition.     

Terminal regions of mouse nebulin: sequence analysis and complementary localization with N-RAP

The regions of mouse nebulin extending from the ends of the super repeats to the C-terminus and N-terminus were cloned and sequenced. Comparison of the mouse sequence with the previously published human sequence shows that the terminal regions of nebulin are highly conserved. The four phosphorylation motifs and SH3 domain found at the C-terminus of mouse nebulin are identical to those found in human nebulin, with the exception of four conservative substitutions. The modules linking this C-terminal region to the super repeats have deletions relative to both fetal and adult human nebulins that correspond to integral numbers of modules, making the mouse C-terminal simple repeat region among the shortest observed to date. The N-terminal region and the C-terminal modules were expressed in Escherichia coli and used for antibody production. Immunofluorescent labeling of these regions of nebulin in isolated myofibrils demonstrates that they are located near the center of the sarcomere and near the Z-line, respectively. Immunogold labeling with antibodies raised against the N-terminal nebulin sequence localizes this region in the A-band near the tips of the thin filaments. Nebulin localization is complementary to that of N-RAP, another muscle-specific protein containing nebulin-like super repeats; nebulin is exclusively found in the sarcomeres, while N-RAP is confined to the terminal bundles of actin filaments at the myotendinous junction. Cell Motil. Cytoskeleton 3:211-222, 2000 Published 2000 Wiley-Liss, Inc.     

Competition between Intramolecular and Intermolecular Interactions in an Amyloid-Forming Protein

Despite much progress in understanding the folding and the aggregation processes of proteins, the rules defining their interplay have yet to be fully defined. This problem is of particular importance since many diseases are initiated by protein unfolding and hence the propensity to aggregate competes with intramolecular collapse and other folding events. Here, we describe the roles of intramolecular and intermolecular interactions in defining the length of the lag time and the apparent rate of elongation of the 100-residue protein human β₂-microglobulin at pH 2.5, commencing from an acid-denatured state that lacks persistent structure but contains significant non-random hydrophobic interactions. Using a combination of site-directed mutagenesis, quantitative kinetic analysis and computational methods, we show that only a single region of about 10 residues in length, determines the rate of fibril formation, despite the fact that other regions exhibit a significant intrinsic propensity for aggregation. We rationalise these results by analysing the effect of incorporating the conformational properties of acid-unfolded β₂-microglobulin and its variants at pH 2.5 as measured by NMR spectroscopy into the Zyggregator aggregation prediction algorithm. These results demonstrate that residual structure in the precursor state modulates the intrinsic propensity of the polypeptide chain to aggregate and that the algorithm developed here allows the key regions for aggregation to be more clearly identified and the rates of their self-association to be predicted. Given the common propensity of unfolded chains to form non-random intramolecular interactions as monomers and to self-assemble subsequently into amyloid fibrils, the approach developed should find widespread utility for the prediction of regions important in amyloid formation and their rates of self-assembly.     

Functional importance of the conserved N-terminal domain of the mitochondrial replicative DNA helicase

The mitochondrial replicative DNA helicase is an essential cellular protein that shows high similarity with the bifunctional primase-helicase of bacteriophage T7, the gene 4 protein (T7 gp4). The N-terminal primase domain of T7 gp4 comprises seven conserved sequence motifs, I, II, III, IV, V, VI, and an RNA polymerase basic domain. The putative primase domain of metazoan mitochondrial DNA helicases has diverged from T7 gp4 and in particular, the primase domain of vertebrates lacks motif I, which comprises a zinc binding domain. Interestingly, motif I is conserved in insect mtDNA helicases. Here, we evaluate the effects of overexpression in Drosophila cell culture of variants carrying mutations in conserved amino acids in the N-terminal region, including the zinc binding domain. Overexpression of alanine substitution mutants of conserved amino acids in motifs I, IV, V and VI and the RNA polymerase basic domain results in increased mtDNA copy number as is observed with overexpression of the wild type enzyme. In contrast, overexpression of three N-terminal mutants W282L, R301Q and P302L that are analogous to human autosomal dominant progressive external ophthalmoplegia mutations results in mitochondrial DNA depletion, and in the case of R301Q, a dominant negative cellular phenotype. Thus whereas our data suggest lack of a DNA primase activity in Drosophila mitochondrial DNA helicase, they show that specific N-terminal amino acid residues that map close to the central linker region likely play a physiological role in the C-terminal helicase function of the protein.     

Low-Level Expression of a Folding-Incompetent Protein in Escherichia coli: Search for the Molecular Determinants of Protein Aggregation In Vivo

Aggregation of peptides and proteins into insoluble amyloid fibrils or related intracellular inclusions is the hallmark of many degenerative diseases, including Alzheimer's disease, Parkinson's disease, and various forms of amyloidosis. In spite of the considerable progress carried out in vitro in elucidating the molecular determinants of the conversion of purified and isolated proteins into amyloid fibrils, very little is known on factors governing this process in the complex environment of living organisms. Taking advantage of increasing evidence that bacterial inclusion bodies consist of amyloid-like aggregates, we have expressed in Escherichia coli both wild type and 21 single-point mutants of the N-terminal domain of the E. coli protein HypF. All variants were expressed as folding-incompetent units in a controlled manner, at low and comparable levels. Their solubilities were measured by quantifying the protein amount contained in the soluble and insoluble fractions by Western blot analysis. A significant negative correlation was found between the solubility of the variants in E. coli and their intrinsic propensity to form amyloid fibrils, predicted using an algorithm previously validated experimentally in vitro on a number of unfolded peptides and proteins, and considering hydrophobicity, β-sheet propensity, and charge as major sequence determinants of the aggregation process. These findings show that the physicochemical parameters previously recognized to govern amyloid formation by fully or partially unfolded proteins are largely applicable in vivo and pave the way for the molecular exploration of a process as complex as protein aggregation in living organisms.