Trapping the on-pathway folding intermediate of Im7 at equilibrium

The four-helical protein Im7 folds via a rapidly formed on-pathway intermediate (k(UI)=3000 s(-1) at pH 7.0, 10 degrees C) that contains three (helices I, II and IV) of the four native alpha-helices. The relatively slow (k(IN)=300 s(-1)) conversion of this intermediate into the native structure is driven by the folding and docking of the six residue helix III onto the developing hydrophobic core. Here, we describe the structural properties of four Im7* variants designed to trap the protein in the intermediate state by disrupting the stabilising interactions formed between helix III and the rest of the protein structure. In two of these variants (I54A and L53AI54A), hydrophobic residues within helix III have been mutated to alanine, whilst in the other two mutants the sequence encompassing the native helix III was replaced by a glycine linker, three (H3G3) or six (H3G6) residues in length. All four variants were shown to be monomeric, as judged by analytical ultracentrifugation, and highly helical as measured by far-UV CD. In addition, all the variants denature co-operatively and have a stability (DeltaG(UF)) and buried hydrophobic surface area (M(UF)) similar to those of the on-pathway kinetic intermediate. Structural characterisation of these variants using 1-anilino-8-napthalene sulphonic acid (ANS) binding, near-UV CD and 1D (1)H NMR demonstrate further that the trapped intermediate ensemble is highly structured with little exposed hydrophobic surface area. Interestingly, however, the structural properties of the variants I54A and L53AI54A differ in detail from those of H3G3 and H3G6. In particular, the single tryptophan residue, located near the end of helix IV, and distant from helix III, is in a distinct environment in the two sets of mutants as judged by fluorescence, near-UV CD and the sensitivity of tryptophan fluorescence to iodide quenching. Overall, the results confirm previous kinetic analysis that demonstrated the hierarchical folding of Im7 via an on-pathway intermediate, and show that this species is a highly helical ensemble with a well-formed hydrophobic core. By contrast with the native state, however, the intermediate ensemble is flexible enough to change in response to mutation, its structural properties being tailored by residues in the sequence encompassing the native helix III.     

The mechanism of proton translocation in respiratory complex I from molecular dynamics

Abstract

Respiratory complex I, the biggest enzyme of respiratory chain, plays a key role in energy production by the mitochondrial respiratory chain and has been implicated in many human neurodegenerative diseases. Recently, the crystal structure of respiratory complex I is reported. We perform 50?ns molecular dynamics simulations on the membrane domain of respiratory complex I under two hypothetical states (oxidized state and reduced state). We find that the density of water molecules in the trans-membrane domain under reduced state is bigger than that under oxidized state. The connecting elements (helix HL and β-hairpins-helix element) fluctuate stronger under reduced state than that under oxidized state, causing more internal water molecules and facilitating the proton conduction. The conformational changes of helix HL and the crucial charged residue Glu in TM5 play key roles in the mechanism of proton translocation. Our results illustrate the dynamic behavior and the potential mechanism of respiratory complex I, which provides the structural basis for drug design of respiratory complex I.     

Homology modeling of the human microsomal glucose 6-phosphate transporter explains the mutations that cause the glycogen storage disease type Ib

Glycogen storage disease type Ib is caused by mutations in the glucose 6-phosphate transporter (G6PT) in the endoplasmic reticulum membrane in liver and kidney. Twenty-eight missense and two deletion mutations that cause the disease were previously shown to reduce or abolish the transporter's activity. However, the mechanisms by which these mutations impair transport remain unknown. On the basis of the recently determined crystal structure of its Escherichia coli homologue, the glycerol 3-phosphate transporter, we built a three-dimensional structural model of human G6PT by homology modeling. G6PT is proposed to consist of 12 transmembrane alpha-helices that are divided into N- and C-terminal domains, with the substrate-translocation pore located between the two domains and the substrate-binding site formed by R28 and K240 at the domain interface. The disease-causing mutations were found to occur at four types of positions: (I) in the substrate-translocation pore, (II) at the N-/C-terminal domain interface, (III) in the interior of the N- and C-terminal domains, and (IV) on the protein surface. Whereas class I mutations affect substrate binding directly, class II mutations, mostly involving changes in side chain size, charge, or both, hinder the conformational change required for substrate translocation. On the other hand, class III and class IV mutations, often introducing a charged residue into a helix bundle or at the protein-lipid interface, probably destabilize the protein. These results also suggest that G6PT operates by a similar antiport mechanism as its E. coli homologue, namely, the substrate binds at the N- and C-terminal domain interface and is then transported across the membrane via a rocker-switch type of movement of the two domains.     

Equilibrium phi-analysis of a molten globule: the 1-149 apoflavodoxin fragment

The apoflavodoxin fragment comprising residues 1-149 that can be obtained by chemical cleavage of the C-terminal alpha-helix of the full-length protein is known to populate a molten globule conformation that displays a cooperative behaviour and experiences two-state urea and thermal denaturation. Here, we have used a recombinant form of this fragment to investigate molten globule energetics and to derive structural information by equilibrium Phi-analysis. We have characterized 15 mutant fragments designed to probe the persistence of native interactions in the molten globule and compared their conformational stability to that of the equivalent full-length apoflavodoxin mutants. According to our data, most of the mutations analysed modify the stability of the molten globule fragment following the trend observed when the same mutations are implemented in the full-length protein. However, the changes in stability observed in the molten globule are much smaller and the Phi-values calculated are (with a single exception) below 0.4. This is consistent with an overall and significant debilitation of the native structure. Nevertheless, the fact that the molten globule fragment can be stabilised using as a guide the native structure of the full-length protein (by increasing helix propensity, optimising charge interactions and filling small cavities) suggests that the overall structure of the molten globule is still quite close to native, in spite of the lowered stability observed.     

Benchmarking of TASSER in the ab initio limit

A significant number of protein sequences in a given proteome have no obvious evolutionarily related protein in the database of solved protein structures, the PDB. Under these conditions, ab initio or template-free modeling methods are the sole means of predicting protein structure. To assess its expected performance on proteomes, the TASSER structure prediction algorithm is benchmarked in the ab initio limit on a representative set of 1129 nonhomologous sequences ranging from 40 to 200 residues that cover the PDB at 30% sequence identity and which adopt alpha, alpha + beta, and beta secondary structures. For sequences in the 40-100 (100-200) residue range, as assessed by their root mean square deviation from native, RMSD, the best of the top five ranked models of TASSER has a global fold that is significantly close to the native structure for 25% (16%) of the sequences, and with a correct identification of the structure of the protein core for 59% (36%). In the absence of a native structure, the structural similarity among the top five ranked models is a moderately reliable predictor of folding accuracy. If we classify the sequences according to their secondary structure content, then 64% (36%) of alpha, 43% (24%) of alpha + beta, and 20% (12%) of beta sequences in the 40-100 (100-200) residue range have a significant TM-score (TM-score > or = 0.4). TASSER performs best on helical proteins because there are less secondary structural elements to arrange in a helical protein than in a beta protein of equal length, since the average length of a helix is longer than that of a strand. In addition, helical proteins have shorter loops and dangling tails. If we exclude these flexible fragments, then TASSER has similar accuracy for sequences containing the same number of secondary structural elements, irrespective of whether they are helices and/or strands. Thus, it is the effective configurational entropy of the protein that dictates the average likelihood of correctly arranging the secondary structure elements.     

Strand 6B deformation and residues exposure towards N-terminal end of helix B during proteinase inhibition by Serpins

Serine Protease inhibitors (Serpins) like antithrombin, antitrypsin, neuroserpin, antichymotrypsin, protein C-inhibitor and plasminogen activator inhibitor is involved in important biological functions like blood coagulation, fibrinolysis, inflammation, cell migration and complement activation. Serpins native state is metastable, which undergoes transformation to a more stable state during the process of protease inhibition. Serpins are prone to conformation defects, however little is known about the factors and mechanisms which promote its conformational change and misfolding. Helix B region in serpins is with several point mutations which result in pathological conditions due to polymerization. Helix B analysis for residue burial and cavity was undertaken to understand its role in serpin structure function. A structural overlap and an accessible surface area analysis showed the deformation of strand 6B and exposure of helix B at N-terminal end in cleaved conformation but not in the native and latent conformation of various inhibitory serpins. A cleaved polymer like conformation of antitrypsin also showed deformation of s6B and helix B exposure. Cavity analysis showed that helix B residues were part of the largest cavity in most of the serpins in the native state which increase in size during the transformation to cleaved and latent states. These data for the first time show the importance of strand 6B deformation and exposure of helix B in smooth insertion of the reactive center loop during serpin inhibition and indicate that helix B exposure due to variants may increase its polymer propensity. ABBREVIATIONS: serpin -serine protease inhibitors RCL -reactive center loop ASA -accessible surface area.     

Small libraries of protein fragments model native protein structures accurately

Prediction of protein structure depends on the accuracy and complexity of the models used. Here, we represent the polypeptide chain by a sequence of rigid fragments that are concatenated without any degrees of freedom. Fragments chosen from a library of representative fragments are fit to the native structure using a greedy build-up method. This gives a one-dimensional representation of native protein three-dimensional structure whose quality depends on the nature of the library. We use a novel clustering method to construct libraries that differ in the fragment length (four to seven residues) and number of representative fragments they contain (25-300). Each library is characterized by the quality of fit (accuracy) and the number of allowed states per residue (complexity). We find that the accuracy depends on the complexity and varies from 2.9A for a 2.7-state model on the basis of fragments of length 7-0.76A for a 15-state model on the basis of fragments of length 5. Our goal is to find representations that are both accurate and economical (low complexity). The models defined here are substantially better in this regard: with ten states per residue we approximate native protein structure to 1A compared to over 20 states per residue needed previously.For the same complexity, we find that longer fragments provide better fits. Unfortunately, libraries of longer fragments must be much larger (for ten states per residue, a seven-residue library is 100 times larger than a five-residue library). As the number of known protein native structures increases, it will be possible to construct larger libraries to better exploit this correlation between neighboring residues. Our fragment libraries, which offer a wide range of optimal fragments suited to different accuracies of fit, may prove to be useful for generating better decoy sets for ab initio protein folding and for generating accurate loop conformations in homology modeling.     

Structural instability of the prion protein upon M205S/R mutations revealed by molecular dynamics simulations

The point mutations M205S and M205R have been demonstrated to severely disturb the folding and maturation process of the cellular prion protein (PrP(C)). These disturbances have been interpreted as consequences of mutation-induced structural changes in PrP, which are suggested to involve helix 1 and its attachment to helix 3, because the mutated residue M205 of helix 3 is located at the interface of these two helices. Furthermore, current models of the prion protein scrapie (PrP(Sc)), which is the pathogenic isoform of PrP(C) in prion diseases, imply that helix 1 disappears during refolding of PrP(C) into PrP(Sc). Based on molecular-dynamics simulations of wild-type and mutant PrP(C) in aqueous solution, we show here that the native PrP(C) structure becomes strongly distorted within a few nanoseconds, once the point mutations M205S and M205R have been applied. In the case of M205R, this distortion is characterized by a motion of helix 1 away from the hydrophobic core into the aqueous environment and a subsequent structural decay. Together with experimental evidence on model peptides, this decay suggests that the hydrophobic attachment of helix 1 to helix 3 at M205 is required for its correct folding into its stable native structure.     

Protein hydrogen exchange mechanism: local fluctuations

Experiments were done to study the dynamic structural motions that determine protein hydrogen exchange (HX) behavior. The replacement of a solvent-exposed lysine residue with glycine (Lys8Gly) in a helix of recombinant cytochrome c does not perturb the native structure, but it entropically potentiates main-chain flexibility and thus can promote local distortional motions and large-scale unfolding. The mutation accelerates amide hydrogen exchange of the mutated residue by about 50-fold, neighboring residues in the same helix by less, and residues elsewhere in the protein not at all, except for Leu98, which registers the change in global stability. The pattern of HX changes shows that the coupled structural distortions that dominate exchange can be several residues in extent, but they expose to exchange only one amide NH at a time. This "local fluctuation" mode of hydrogen exchange may be generally recognized by disparate near-neighbor rates and a low dependence on destabilants (denaturant, temperature, pressure). In contrast, concerted unfolding reactions expose multiple neighboring amide NHs with very similar computed protection factors, and they show marked destabilant sensitivity. In both modes, ionic hydrogen exchange catalysts attack from the bulk solvent without diffusing through the protein matrix.     

Structural and dynamic characterization of an unfolded state of poplar apo-plastocyanin formed under nondenaturing conditions

Plastocyanin is a predominantly beta-sheet protein containing a type I copper center. The conformational ensemble of a denatured state of apo-plastocyanin formed in solution under conditions of low salt and neutral pH has been investigated by multidimensional heteronuclear NMR spectroscopy. Chemical shift assignments were obtained by using three-dimensional triple-resonance NMR experiments to trace through-bond heteronuclear connectivities along the backbone and side chains. The (3)J(HN,Halpha) coupling constants, (15)N-edited proton-proton nuclear Overhauser effects (NOEs), and (15)N relaxation parameters were also measured for the purpose of structural and dynamic characterization. Most of the residues corresponding to beta-strands in the folded protein exhibit small upfield shifts of the (13)C(alpha) and (13)CO resonances relative to random coil values, suggesting a slight preference for backbone dihedral angles in the beta region of (phi,psi) space. This is further supported by the presence of strong sequential d(alphaN)(i, i + 1) NOEs throughout the sequence. The few d(NN)(i, i + 1) proton NOEs that are observed are mostly in regions that form loops in the native plastocyanin structure. No medium or long-range NOEs were observed. A short sequence, between residues 59 and 63, was found to populate a nonnative helical conformation in the unfolded state, as indicated by the shift of the (13)C(alpha), (13)CO, and (1)H(alpha) resonances relative to random coil values and by the decreased values of the (3)J(HN,Halpha) coupling constants. The (15)N relaxation parameters indicate restriction of motions on a nanosecond timescale in this region. Intriguingly, this helical conformation is present in a sequence that is close to but not in the same location as the single short helix in the native folded protein. The results are consistent with earlier NMR studies of peptide fragments of plastocyanin and confirm that the regions of the sequence that form beta-strands in the native protein spontaneously populate the beta-region of (phi,psi) space under folding conditions, even in the absence of stabilizing tertiary interactions. We conclude that the state of apo-plastocyanin present under nondenaturing conditions is a noncompact unfolded state with some evidence of nativelike and nonnative local structuring that may be initiation sites for folding of the protein.     

Glutamate-41 of Vibrio harveyi acyl carrier protein is essential for fatty acid synthase but not acyl-ACP synthetase activity

Bacterial acyl carrier protein (ACP) is a small, acidic, and highly conserved protein that supplies acyl groups for biosynthesis of a variety of lipid products. Recent modelling studies predict that residues primarily in helix II of Escherichia coli ACP (Glu-41, Ala-45) are involved in its interaction with the condensing enzyme FabH of fatty acid synthase. Using recombinant Vibrio harveyi ACP as a template for site-directed mutagenesis, we have shown that an acidic residue at position 41 is essential for V. harveyi fatty acid synthase (but not acyl-ACP synthetase) activity. In contrast, various replacements of Ala-45 were tolerated by both enzymes. None of the mutations introduced dramatic structural changes based on circular dichroism and native gel electrophoresis. These results confirm that Glu-41 of ACP is a critical residue for fatty acid synthase, but not for all enzymes that utilize ACP as a substrate.     

Modulation of the structural integrity of helix F in apomyoglobin by single amino acid replacements

The conformational features of native and mutant forms of sperm-whale apomyoglobin (apoMb) at neutral pH were probed by limited proteolysis experiments utilizing up to eight proteases of different substrate specificities. It was shown that all proteases selectively cleave apoMb at the level of chain segment 82-94 (HEAELKPLAQSHA), encompassing helix F in the X-ray structure of the holo form of the native protein; for example, thermolysin cleaves the Pro 88-Leu 89 peptide bond. These results indicate that helix F is highly flexible or largely disrupted in apoMb. Because helix F contains the helix-breaking Pro 88 residue, we propose that helix F is kept in place in the native holo protein by a variety of helix-heme stabilizing interactions. To modulate the stability of helix F, the Pro88Ala and Pro88Gly mutants were prepared by site-directed mutagenesis, and their conformational properties investigated by both far-UV circular dichroism spectroscopy and limited proteolysis. The helix content of the Pro88Ala mutant was somewhat enhanced with respect to that of both native and Pro88Gly mutant, as expected from the fact that Ala is the strongest helix inducer among the 20 amino acid residues. The rate of limited proteolysis of the three apoMb variants by thermolysin and proteinase K was in the order native > Pro88Gly > Pro88Ala, in agreement with the scale of helix propensity of Ala, Gly, and Pro. The possible role of the flexible/unfolded chain segment 82-94 for the function and fate of apoMb at the cellular level is discussed.     

Molecular organization of 19 S calf thyroglobulin

Controlled freezing and thawing of 19 S calf thyroglobulin resulted in a specific and reproducible breakdown of the protein. Beside the elementary chain (300,000 Da), new, discrete bands are revealed by gel electrophoresis in sodium dodecyl sulfate under reducing conditions. These new species consist of a major peptide of 100,000 Da and several faster-migrating bands. Most of these polypeptides were purified by preparative gel electrophoresis and individually digested in formic acid with CNBr. A comparative gel electrophoresis under denaturating and reducing conditions of (i) the fragments obtained from the native protein, (ii) the electrophoretically purified elementary chain, (iii) the 100,000-Da peptide, and (iv) a smaller fragment (of about 50,000) was performed. It revealed a very close homology among the peptide maps of the intact 19 S, the elementary chain, and the 100,000-Da peptide. Furthermore, it was shown that the digestion products of the smaller fragment, present in the peptide map of the native protein, were absent in both the elementary chain and the 100,000-Da species. These results support the idea that calf thyroglobulin, even though it has an apparently complex molecular organization, contains structural motifs which are repeated in the elementary chain.     

NMR characterization of the conformational fluctuations of the human lymphocyte function‐associated antigen‐1 I‐domain

Lymphocyte function‐associated antigen‐1 (LFA‐1) is an integrin protein that transmits information across the plasma membrane through the so‐called inside‐out and outside‐in signaling mechanisms. To investigate these mechanisms, we carried out an NMR analysis of the dynamics of the LFA‐1 I‐domain, which has enabled us to characterize the motions of this domain on a broad range of timescales. We studied first the internal motions on the nanosecond timescale by spin relaxation measurements and model‐free analysis. We then extended this analysis to the millisecond timescale motions by measuring ¹⁵N‐¹H residual dipolar couplings of the backbone amide groups. We analyzed these results in the context of the three major conformational states of the I‐domain using their corresponding X‐ray crystallographic structures. Our results highlight the importance of the low‐frequency motions of the LFA‐1 I‐domain in the inactive apo‐state. We found in particular that α‐helix 7 is in a position in the apo‐closed state that cannot be fully described by any of the existing X‐ray structures, as it appears to be in dynamic exchange between different conformations. This type of motion seems to represent an inherent property of the LFA‐1 I‐domain and might be relevant for controlling the access to the allosteric binding pocket, as well as for the downward displacement of α‐helix 7 that is required for the activation of LFA‐1.     

On the role of some conserved and nonconserved amino acid residues in the transitional state and intermediate of apomyoglobin folding

The contributions of some amino acid residues in the A, B, G, and H helices to the formation of the folding nucleus and folding intermediate of apomyoglobin were estimated. The effects of point substitutions of Ala for hydrophobic amino acid residues on the structural stability of the native (N) protein and its folding intermediate (I), as well as on the folding/unfolding rates for four mutant apomyoglobin forms, were studied. The equilibrium and kinetic studies of the folding/unfolding rates of these mutant proteins in a wide range of urea concentrations demonstrated that their native state was considerably destabilized as compared with the wild-type protein, whereas the stability of the intermediate state changed moderately. It was shown that the amino acid residues in the A, G, and H helices contributed insignificantly to the stabilization of the apomyoglobin folding nucleus in the rate-limiting I ? N transition, taking place after the formation of the intermediate, whereas the residue of the B helix was of great importance in the formation of the folding nucleus in this transition.     

Effects of alanine substitutions in alpha-helices of sperm whale myoglobin on protein stability.

The peptide backbones in folded native proteins contain distinctive secondary structures, alpha-helices, beta-sheets, and turns, with significant frequency. One question that arises in folding is how the stability of this secondary structure relates to that of the protein as a whole. To address this question, we substituted the alpha-helix-stabilizing alanine side chain at 16 selected sites in the sequence of sperm whale myoglobin, 12 at helical sites on the surface of the protein, and 4 at obviously internal sites. Substitution of alanine for bulky side chains at internal sites destabilizes the protein, as expected if packing interactions are disrupted. Alanine substitutions do not uniformly stabilize the protein, either in capping positions near the ends of helices or at mid-helical sites near the surface of myoglobin. When corrected for the extent of exposure of each side chain replaced by alanine at a mid-helix position, alanine replacement still has no clear effect in stabilizing the native structure. Thus linkage between the stabilization of secondary structure and tertiary structure in myoglobin cannot be demonstrated, probably because of the relatively small free energy differences between side chains in stabilizing isolated helix. By contrast, about 80% of the variance in free energy observed can be accounted for by the loss in buried surface area of the native residue substituted by alanine. The differential free energy of helix stabilization does not account for any additional variation.     

A conformational equilibrium in a protein fragment caused by two consecutive capping boxes: 1H-, 13C-NMR, and mutational analysis.

The conformational properties of an 18 residues peptide spanning the entire sequence, L1KTPA5QFDAD10ELRAA15MKG, of the first helix (A-helix) of domain 2 of annexin I, were thoroughly investigated. This fragment exhibits several singular features, and in particular, two successive potential capping boxes, T3xxQ6 and D8xxE11. The former corresponds to the native hydrogen bond network stabilizing the alpha helix N-terminus in the protein; the latter is a non-native capping box able to break the helix at residue D8, and is observed in the domain 2 partially folded state. Using 2D-NMR techniques, we showed that two main populations of conformers coexist in aqueous solution. The first corresponds to a single helix extending from T3 to K17. The second corresponds to a broken helix at residue Ds. Four mutants, T3A, F7A, D8A, and E11A, were designed to further analyze the role of key amino acids in the equilibrium between the two ensembles of conformers. The sensitivity of NMR parameters to account for the variations in the populations of conformers was evaluated for each peptide. Our data show the delta13Calpha chemical shift to be the most relevant parameter. We used it to estimate the population ratio in the various peptides between the two main ensembles of conformers, the full helix and the broken helix. For the WT, E11A, and F7A peptides, these ratios are respectively 35/65, 60/40, 60/40. Our results were compared to the data obtained from helix/coil transition algorithms.     

Stability and folding kinetics of a ubiquitin mutant with a strong propensity for nonnative beta-hairpin conformation in the unfolded state

A F45W mutant of yeast ubiquitin has been used as a model system to examine the effects of nonnative local interactions on protein folding and stability. Mutating the native TLTGK G-bulged type I turn in the N-terminal beta-hairpin to NPDG stabilizes a nonnative beta-strand alignment in the isolated peptide fragment. However, NMR structural analysis of the native and mutant proteins shows that the NPDG mutant is forced to adopt the native beta-strand alignment and an unfavorable type I NPDG turn. The mutant is significantly less stable (approximately 9 kJ mol(-1)) and folds 30 times slower than the native sequence, demonstrating that local interactions can modulate protein stability and that attainment of a nativelike beta-hairpin conformation in the transition state ensemble is frustrated by the turn mutations. Surprising, alcoholic cosolvents [5-10% (v/v) TFE] are shown to accelerate the folding rate of the NPDG mutant. We conclude, backed-up by NMR data on the peptide fragments, that even though nonnative states in the denatured ensemble are highly populated and their stability further enhanced in the presence of cosolvents, the simultaneous increase in the proportion of nativelike secondary structure (hairpin or helix), in rapid equilibrium with nonnative states, is sufficient to accelerate the folding process. It is evident that modulating local interactions and increasing nonnative secondary structure propensities can change protein stability and folding kinetics. However, nonlocal contacts formed in the global cooperative folding event appear to determine structural specificity.     

A highly amyloidogenic region of hen lysozyme

Amyloid fibrils obtained after incubating hen egg-white lysozyme (HEWL) at pH 2.0 and 65 degrees C for extended periods of time have been found to consist predominantly of fragments of the protein corresponding to residues 49-100, 49-101, 53-100 and 53-101, derived largely from the partial acid hydrolysis of Asp-X peptide bonds. These internal fragments of HEWL encompass part of the beta-domain and all the residues forming the C-helix in the native protein, and contain two internal disulfide bridges Cys64-Cys80 and Cys76-Cys94. The complementary protein fragments, including helices A, B and D of the native protein, are not significantly incorporated into the network of fibrils, but remain largely soluble, in agreement with their predicted lower propensities to aggregate. Further analysis of the properties of different regions of HEWL to form amyloid fibrils was carried out by studying fragments produced by limited proteolysis of the protein by pepsin. Here, we show that only fragment 57-107, but not fragment 1-38/108-129, is able to generate well-defined amyloid fibrils under the conditions used. This finding is of particular importance, as the beta-domain and C-helix of the highly homologous human lysozyme have been shown to unfold locally in the amyloidogenic variant D67H, which is associated with the familial cases of systemic amyloidosis linked to lysozyme deposition. The identification of the highly amyloidogenic character of this region of the polypeptide chain provides strong support for the involvement of partially unfolded species in the initiation of the aggregation events that lead to amyloid deposition in clinical disease.     

Sequence and structural determinants of Cu, Zn superoxide dismutase aggregation

Diverse point mutations in the enzyme Cu, Zn superoxide dismutase (SOD1) are linked to its aggregation in the familial form of the disease amyotrophic lateral sclerosis. The disease-associated mutations are known to destabilize the protein, but the structural basis of the aggregation of the destabilized protein and the structure of aggregates are not well understood. Here, we investigate in silico the sequence and structural determinants of SOD1 aggregation: (1) We identify sequence fragments in SOD1 that have a high aggregation propensity, using only the sequence of SOD1, and (2) we perform molecular dynamics simulations of the SOD1 dimer folding and misfolding. In both cases, we identify identical regions of the protein as having high propensity to form intermolecular interactions. These regions correspond to the N- and C-termini, and two crossover loops and two beta-strands in the Greek-key native fold of SOD1. Our results suggest that the high aggregation propensity of mutant SOD1 may result from a synergy of two factors: the presence of highly amyloidogenic sequence fragments ("hot spots"), and the presence of these fragments in regions of the protein that are structurally most likely to form intermolecular contacts under destabilizing conditions. Therefore, we postulate that the balance between the self-association of aggregation-prone sequences and the specific structural context of these sequences in the native state determines the aggregation propensity of proteins.