Slow formation of aggregation‐resistant β‐sheet folding intermediates

Protein folding has been studied extensively for decades, yet our ability to predict how proteins reach their native state from a mechanistic perspective is still rudimentary at best, limiting our understanding of folding‐related processes in vivo and our ability to manipulate proteins in vitro. Here, we investigate the in vitro refolding mechanism of a large β‐helix protein, pertactin, which has an extended, elongated shape. At 55 kDa, this single domain, all‐β‐sheet protein allows detailed analysis of the formation of β‐sheet structure in larger proteins. Using a combination of fluorescence and far‐UV circular dichroism spectroscopy, we show that the pertactin β‐helix refolds remarkably slowly, with multiexponential kinetics. Surprisingly, despite the slow refolding rates, large size, and β‐sheet‐rich topology, pertactin refolding is reversible and not complicated by off‐pathway aggregation. The slow pertactin refolding rate is not limited by proline isomerization, and 30% of secondary structure formation occurs within the rate‐limiting step. Furthermore, site‐specific labeling experiments indicate that the β‐helix refolds in a multistep but concerted process involving the entire protein, rather than via initial formation of the stable core substructure observed in equilibrium titrations. Hence pertactin provides a valuable system for studying the refolding properties of larger, β‐sheet‐rich proteins, and raises intriguing questions regarding the prevention of aggregation during the prolonged population of partially folded, β‐sheet‐rich refolding intermediates. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Molecular dynamics simulations of an engineered T4 lysozyme exclude helix to sheet transition,and provide insights into long distance,intra‐protein switchable motion

An engineered variant of T4 lysozyme serves as a model for studying induced remote conformational changes in a full protein context. The design involves a duplicated surface helix, flanked by two loops, that switches between two different conformations spanning about 20 Å. Molecular dynamics simulations of the engineered protein, up to 1 μs, rule out α‐helix to β‐sheet transitions within the duplicated helix as suggested by others. These simulations highlight how the use of different force fields can lead to radical differences in the structure of the protein. In addition, Markov state modeling and transition path theory were employed to map a 6.6 μs simulation for possible early intermediate states and to provide insights into the onset of the switching motion. The putative intermediates involve the folding of one helical turn in the C‐terminal loop through energy driven, sequential rearrangement of nearby salt bridges around the key residue Arg63. These results provide a first step towards understanding the energetics and dynamics of a rather complicated intra‐protein motion.     

Determination of the secondary structural elements of chicken liver fatty acid binding protein by two‐dimensional homonuclear NMR

A conformational study in solution of the fatty acid binding protein from chicken liver is presented. The nearly complete sequence‐specific ¹H resonance assignment was achieved from homonuclear two‐dimensional nmr experiments using a sample of native protein. The principal elements of secondary structure were identified: 10 antiparallel β‐strands and one helical segment followed by a turn comprising 5 residues. These elements correspond closely with those of the crystal structure of the related protein, and two new secondary structural features obtained from the nmr data are the β‐sheet conformation between the first and the last β‐strand in the protein sequence, as well as a helical loop at the N‐terminus of the polypeptide chain. © 1999 John Wiley & Sons, Inc. Biopoly 50: 1–11, 1999     

Helix‐sheet packing in proteins

The three‐dimensional structure of a protein is organized around the packing of its secondary structure elements. Although much is known about the packing geometry observed between α‐helices and between β‐sheets, there has been little progress on characterizing helix–sheet interactions. We present an analysis of the conformation of αβ₂ motifs in proteins, corresponding to all occurrences of helices in contact with two strands that are hydrogen bonded. The geometry of the αβ₂ motif is characterized by the azimuthal angle θ between the helix axis and an average vector representing the two strands, the elevation angle ψ between the helix axis and the plane containing the two strands, and the distance D between the helix and the strands. We observe that the helix tends to align to the two strands, with a preference for an antiparallel orientation if the two strands are parallel; this preference is diminished for other topologies of the β‐sheet. Side‐chain packing at the interface between the helix and the strands is mostly hydrophobic, with a preference for aliphatic amino acids in the strand and aromatic amino acids in the helix. From the knowledge of the geometry and amino acid propensities of αβ₂ motifs in proteins, we have derived different statistical potentials that are shown to be efficient in picking native‐like conformations among a set of non‐native conformations in well‐known decoy datasets. The information on the geometry of αβ₂ motifs as well as the related statistical potentials have applications in the field of protein structure prediction. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Non‐native α‐helix formation is not necessary for folding of lipocalin: Comparison of burst‐phase folding between tear lipocalin and β‐lactoglobulin

Tear lipocalin and β‐lactoglobulin are members of the lipocalin superfamily. They have similar tertiary structures but unusually low overall sequence similarity. Non‐native helical structures are formed during the early stage of β‐lactoglobulin folding. To address whether the non‐native helix formation is found in the folding of other lipocalin superfamily proteins, the folding kinetics of a tear lipocalin variant were investigated by stopped‐flow methods measuring the time‐dependent changes in circular dichroism (CD) spectrum and small‐angle X‐ray scattering (SAXS). CD spectrum showed that extensive secondary structures are not formed during a burst‐phase (within a measurement dead time). The SAXS data showed that the radius of gyration becomes much smaller than in the unfolded state during the burst‐phase, indicating that the molecule is collapsed during an early stage of folding. Therefore, non‐native helix formation is not general for folding of all lipocalin family members. The non‐native helix content in the burst‐phase folding appears to depend on helical propensities of the amino acid sequence. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

A thermodynamic approach to the conformational preferences of the 180–195 segment derived from the human prion protein α2‐helix

On consideration that intrinsic structural weakness could affect the segment spanning the α2‐helical residues 173–195 of the PrP, we have investigated the conformational stabilities of some synthetic Ala‐scanned analogs of the peptide derived from the 180–195 C‐terminal sequence, using a novel approach whose theoretical basis originates from protein thermodynamics. Even though a quantitative comparison among peptides could not be assessed to rank them according to the effect caused by single amino acid substitution, as a general trend, all peptides invariably showed an appreciable preference for an α‐type organization, consistently with the fact that the wild‐type sequence is organized as an α‐helix in the native protein. Moreover, the substitution of whatever single amino acid in the wild‐type sequence reduced the gap between the α‐ and the β‐propensity, invariably enhancing the latter, but in any case this gap was larger than that evaluated for the full‐length α2‐helix‐derived peptide. It appears that the low β‐conformation propensity of the 180–195 region depends on the simultaneous presence of all of the Ala‐scanned residues, indirectly confirming that the N‐terminal 173–179 segment could play a major role in determining the chameleon conformational behavior of the entire 173–195 region in the PrP. Copyright © 2008 European Peptide Society and John Wiley & Sons, Ltd.     

Conformational analysis of the full‐length M2 protein of the influenza A virus using solid‐state NMR

The influenza A M2 protein forms a proton channel for virus infection and mediates virus assembly and budding. While extensive structural information is known about the transmembrane helix and an adjacent amphipathic helix, the conformation of the N‐terminal ectodomain and the C‐terminal cytoplasmic tail remains largely unknown. Using two‐dimensional (2D) magic‐angle‐spinning solid‐state NMR, we have investigated the secondary structure and dynamics of full‐length M2 (M2FL) and found them to depend on the membrane composition. In 2D ¹³C DARR correlation spectra, 1,2‐dimyristoyl‐sn‐glycero‐3‐phosphocholine (DMPC)‐bound M2FL exhibits several peaks at β‐sheet chemical shifts, which result from water‐exposed extramembrane residues. In contrast, M2FL bound to cholesterol‐containing membranes gives predominantly α‐helical chemical shifts. Two‐dimensional J‐INADEQUATE spectra and variable‐temperature ¹³C spectra indicate that DMPC‐bound M2FL is highly dynamic while the cholesterol‐containing membranes significantly immobilize the protein at physiological temperature. Chemical‐shift prediction for various secondary‐structure models suggests that the β‐strand is located at the N‐terminus of the DMPC‐bound protein, while the cytoplasmic domain is unstructured. This prediction is confirmed by the 2D DARR spectrum of the ectodomain‐truncated M2(21–97), which no longer exhibits β‐sheet chemical shifts in the DMPC‐bound state. We propose that the M2 conformational change results from the influence of cholesterol, and the increased helicity of M2FL in cholesterol‐rich membranes may be relevant for M2 interaction with the matrix protein M1 during virus assembly and budding. The successful determination of the β‐strand location suggests that chemical‐shift prediction is a promising approach for obtaining structural information of disordered proteins before resonance assignment.     

Temperature‐induced partially unfolded state of hUBF HMG Box‐5: Conformational and dynamic investigations of the Box‐5 thermal intermediate ensemble

Human upstream binding factor (hUBF) HMG Box‐5 is a highly conserved protein domain, containing 84 amino acids and belonging to the family of the nonspecific DNA‐binding HMG boxes. Its native structure adopts a twisted L shape, which consists of three α‐helices and two hydrophobic cores: the major wing and the minor wing. In this article, we report a reversible three‐state thermal unfolding equilibrium of hUBF HMG Box‐5, which is investigated by differential scanning calorimetry (DSC), circular dichroism spectroscopy, fluorescence spectroscopy, and NMR spectroscopy. DSC data show that Box‐5 unfolds reversibly in two separate stages. Spectroscopic analyses suggest that different structural elements exhibit noncooperative transitions during the unfolding process and that the major form of the Box‐5 thermal intermediate ensemble at 55°C shows partially unfolded characteristics. Compared with previous thermal stability studies of other boxes, it appears that Box‐5 possesses a more stable major wing and two well separated subdomains. NMR chemical shift index and sequential ¹H^N_i‐¹H^N_i+1 NOE analyses indicate that helices 1 and 2 are native‐like in the thermal intermediate ensemble, while helix 3 is partially unfolded. Detailed NMR relaxation dynamics are compared between the native state and the intermediate ensemble. Our results implicate a fluid helix‐turn‐helix folding model of Box‐5, where helices 1 and 2 potentially form the helix 1‐turn‐helix 2 motif in the intermediate, while helix 3 is consolidated only as two hydrophobic cores form to stabilize the native structure. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

BuildBeta—A system for automatically constructing beta sheets

We describe a method that can thoroughly sample a protein conformational space given the protein primary sequence of amino acids and secondary structure predictions. Specifically, we target proteins with β‐sheets because they are particularly challenging for ab initio protein structure prediction because of the complexity of sampling long‐range strand pairings. Using some basic packing principles, inverse kinematics (IK), and β‐pairing scores, this method creates all possible β‐sheet arrangements including those that have the correct packing of β‐strands. It uses the IK algorithms of ProteinShop to move α‐helices and β‐strands as rigid bodies by rotating the dihedral angles in the coil regions. Our results show that our approach produces structures that are within 4–6 Å RMSD of the native one regardless of the protein size and β‐sheet topology although this number may increase if the protein has long loops or complex α‐helical regions. Proteins 2010. © Published 2009 Wiley‐Liss, Inc.     

Folding of Top7 in unbiased all‐atom Monte Carlo simulations

For computational studies of protein folding, proteins with both helical and β‐sheet secondary structure elements are very challenging, as they expose subtle biases of the physical models. Here, we present reproducible folding of a 92 residue α/β protein (residues 3–94 of Top7, PDB ID: 1QYS) in computer simulations starting from random initial conformations using a transferable physical model which has been previously shown to describe the folding and thermodynamic properties of about 20 other smaller proteins of different folds. Top7 is a de novo designed protein with two α‐helices and a five stranded β‐sheet. Experimentally, it is known to be unusually stable for its size, and its folding transition distinctly deviates from the two‐state behavior commonly seen in natural single domain proteins. In our all‐atom implicit solvent parallel tempering Monte Carlo simulations, Top7 shows a rapid transition to a group of states with high native‐like secondary structure, and a much slower subsequent transition to the native state with a root mean square deviation of about 3.5 Å from the experimentally determined structure. Consistent with experiments, we find Top7 to be thermally extremely stable, although the simulations also find a large number of very stable non‐native states with high native‐like secondary structure. Proteins 2013; 81:1446–1456. © 2013 Wiley Periodicals, Inc.     

Conformational dynamics is more important than helical propensity for the folding of the all α‐helical protein Im7

Im7 folds via an on‐pathway intermediate that contains three of the four native α‐helices. The missing helix, helix III, is the shortest and its failure to be formed until late in the pathway is related to frustration in the structure. Im7H3M3, a 94‐residue variant of the 87‐residue Im7 in which helix III is the longest of the four native helices, also folds via an intermediate. To investigate the structural basis for this we calculated the frustration in the structure of Im7H3M3 and used NMR to investigate its dynamics. We found that the native state of Im7H3M3 is highly frustrated and in equilibrium with an intermediate state that lacks helix III, similar to Im7. Model‐free analysis identified residues with chemical exchange contributions to their relaxation that aligned with the residues predicted to have highly frustrated interactions, also like Im7. Finally, we determined properties of urea‐denatured Im7H3M3 and identified four clusters of interacting residues that corresponded to the α‐helices of the native protein. In Im7 the cluster sizes were related to the lengths of the α‐helices with cluster III being the smallest but in Im7H3M3 cluster III was also the smallest, despite this region forming the longest helix in the native state. These results suggest that the conformational properties of the urea‐denatured states promote formation of a three‐helix intermediate in which the residues that form helix III remain non‐helical. Thus it appears that features of the native structure are formed early in folding linked to collapse of the unfolded state.     

Geometrical principles of homomeric β‐barrels and β‐helices: Application to modeling amyloid protofilaments

Examples of homomeric β‐helices and β‐barrels have recently emerged. Here we generalize the theory for the shear number in β‐barrels to encompass β‐helices and homomeric structures. We introduce the concept of the “β‐strip,” the set of parallel or antiparallel neighboring strands, from which the whole helix can be generated giving it n‐fold rotational symmetry. In this context, the shear number is interpreted as the sum around the helix of the fixed register shift between neighboring identical β‐strips. Using this approach, we have derived relationships between helical width, pitch, angle between strand direction and helical axis, mass per length, register shift, and number of strands. The validity and unifying power of the method is demonstrated with known structures including α‐hemolysin, T4 phage spike, cylindrin, and the HET‐s(218‐289) prion. From reported dimensions measured by X‐ray fiber diffraction on amyloid fibrils, the relationships can be used to predict the register shift and the number of strands within amyloid protofilaments. This was used to construct models of transthyretin and Alzheimer β(40) amyloid protofilaments that comprise a single strip of in‐register β‐strands folded into a “β‐strip helix.” Results suggest both stabilization of an individual β‐strip helix and growth by addition of further β‐strip helices can involve the same pair of sequence segments associating with β‐sheet hydrogen bonding at the same register shift. This process would be aided by a repeat sequence. Hence, understanding how the register shift (as the distance between repeat sequences) relates to helical dimensions will be useful for nanotube design.     

Two‐photon excited lasing of Coumarin 307 for lysozyme amyloid fibrils detection

Amyloid fibrils are a well‐recognized hallmark of neurodegeneration. A common approach to detect amyloid fibrils is staining with organic molecules and monitoring optical properties using fluorescence spectroscopy. However, the structural diversity of amyloids necessitates new sensitive methods and probes that can be reliably used to characterize them. Here, Coumarin 307 is applied for lysozyme fibrils detection by observation of laser action in the process of two‐photon excited stimulated emission. It is shown that the lasing threshold and spectrum significantly depend on the adopted structure (α‐helix or β‐sheet) of the lysozyme protein, whereas fluorescence spectrum is insensitive to the protein structure. The applications of coherent stimulated emission light that can be emitted deep inside a scattering medium can be particularly promising for imaging and therapeutic purposes in the neurodegeneration field. Two‐photon excitation with the near‐infrared light, which allows the deepest penetration of tissues, is an important advantage of the method.     

An alternative structural isoform in amyloid‐like aggregates formed from thermally denatured human γD‐crystallin

The eye lens protein γD‐crystallin contributes to cataract formation in the lens. In vitro experiments show that γD‐crystallin has a high propensity to form amyloid fibers when denatured, and that denaturation by acid or UV‐B photodamage results in its C‐terminal domain forming the β‐sheet core of amyloid fibers. Here, we show that thermal denaturation results in sheet‐like aggregates that contain cross‐linked oligomers of the protein, according to transmission electron microscopy and SDS‐PAGE. We use two‐dimensional infrared spectroscopy to show that these aggregates have an amyloid‐like secondary structure with extended β‐sheets, and use isotope dilution experiments to show that each protein contributes approximately one β‐strand to each β‐sheet in the aggregates. Using segmental ¹³C labeling, we show that the organization of the protein's two domains in thermally induced aggregates results in a previously unobserved structure in which both the N‐terminal and C‐terminal domains contribute to β‐sheets. We propose a model for the structural organization of the aggregates and attribute the recruitment of the N‐terminal domain into the fiber structure to intermolecular cross linking.     

Mechanism of formation of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin‐binding protein G from Streptococcus. IV. Implication for the mechanism of folding of the parent protein

A 34‐residue α/β peptide [IG(28–61)], derived from the C‐terminal part of the B3 domain of the immunoglobulin binding protein G from Streptoccocus, was studied using CD and NMR spectroscopy at various temperatures and by differential scanning calorimetry. It was found that the C‐terminal part (a 16‐residue‐long fragment) of this peptide, which corresponds to the sequence of the β‐hairpin in the native structure, forms structure similar to the β‐hairpin only at T = 313 K, and the structure is stabilized by non‐native long‐range hydrophobic interactions (Val47–Val59). On the other hand, the N‐terminal part of IG(28–61), which corresponds to the middle α‐helix in the native structure, is unstructured at low temperature (283 K) and forms an α‐helix‐like structure at 305 K, and only one helical turn is observed at 313 K. At all temperatures at which NMR experiments were performed (283, 305, and 313 K), we do not observe any long‐range connectivities which would have supported packing between the C‐terminal (β‐hairpin) and the N‐terminal (α‐helix) parts of the sequence. Such interactions are absent, in contrast to the folding pathway of the B domain of protein G, proposed recently by Kmiecik and Kolinski (Biophys J 2008, 94, 726–736), based on Monte‐Carlo dynamics studies. Alternative folding mechanisms are proposed and discussed. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 469–480, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Cofactor‐induced reversible folding of Flavodoxin‐4 from Lactobacillus acidophilus

Flavodoxins in combination with the flavin mononucleotide (FMN) cofactor play important roles for electron transport in prokaryotes. Here, novel insights into the FMN‐binding mechanism to flavodoxins‐4 were obtained from the NMR structures of the apo‐protein from Lactobacillus acidophilus (YP_193882.1) and comparison of its complex with FMN. Extensive reversible conformational changes were observed upon FMN binding and release. The NMR structure of the FMN complex is in agreement with the crystal structure (PDB ID: 3EDO ) and exhibits the characteristic flavodoxin fold, with a central five‐stranded parallel β–sheet and five α‐helices forming an α/β‐sandwich architecture. The structure differs from other flavoproteins in that helix α2 is oriented perpendicular to the β‐sheet and covers the FMN‐binding site. This helix reversibly unfolds upon removal of the FMN ligand, which represents a unique structural rearrangement among flavodoxins.     

Kinetically trapped metastable intermediate of a disulfide‐deficient mutant of the starch‐binding domain of glucoamylase

Refolding of a thermally unfolded disulfide‐deficient mutant of the starch‐binding domain of glucoamylase was investigated using differential scanning calorimetry, isothermal titration calorimetry, CD, and ¹H NMR. When the protein solution was rapidly cooled from a higher temperature, a kinetic intermediate was formed during refolding. The intermediate was unexpectedly stable compared with typical folding intermediates that have short half‐lives. It was shown that this intermediate contained substantial secondary structure and tertiary packing and had the same binding ability with β‐cyclodextrin as the native state, suggesting that the intermediate is highly‐ordered and native‐like on the whole. These characteristics differ from those of partially folded intermediates such as molten globule states. Far‐UV CD spectra showed that the secondary structure was once disrupted during the transition from the intermediate to the native state. These results suggest that the intermediate could be an off‐pathway type, possibly a misfolded state, that has to undergo unfolding on its way to the native state.     

NMR‐based characterization of a refolding intermediate of β2‐microglobulin labeled using a wheat germ cell‐free system

In patients with dialysis‐related amyloidosis, β2‐microglobulin (β2‐m) is a major structural component of amyloid fibrils. It has been suggested that the partial unfolding of β2‐m is a prerequisite to the formation of amyloid fibrils, and that the folding intermediate trapped by the non‐native trans‐Pro32 isomer leads to the formation of amyloid fibrils. Although clarifying the structure of this refolding intermediate by high resolution NMR spectroscopy is important, this has been made difficult by the limited lifetime of the intermediate. Here, we studied the structure of the refolding intermediate using a combination of amino acid selective labeling with wheat germ cell‐free protein synthesis and NMR techniques. The HSQC spectra of β2‐ms labeled selectively at either phenylalanine, leucine, or valine enabled us to monitor the structures of the refolding intermediate. The results suggested that the refolding intermediate has an overall fold and cores similar to the native structure, but contains disordered structures around Pro32. The fluctuation of the β‐sheet regions especially the last half of the βB strand and the first half of the βE strand, both suggested to be important for amyloidogenicity, may transform β2‐m into an amyloidogenic structure.     

PackHelix: a tool for helix-sheet packing during protein structure prediction

The three‐dimensional structure of a protein is organized around the packing of its secondary structure elements. Predicting the topology and constructing the geometry of structural motifs involving α‐helices and/or β‐strands are therefore key steps for accurate prediction of protein structure. While many efforts have focused on how to pack helices and on how to sample exhaustively the topologies and geometries of multiple strands forming a β‐sheet in a protein, there has been little progress on generating native‐like packings of helices on sheets. We describe a method that can generate the packing of multiple helices on a given β‐sheet for αβα sandwich type protein folds. This method mines the results of a statistical analysis of the conformations of αβ₂ motifs in protein structures to provide input values for the geometric attributes of the packing of a helix on a sheet. It then proceeds with a geometric builder that generates multiple arrangements of the helices on the sheet of interest by sampling through these values and performing consistency checks that guarantee proper loop geometry between the helices and the strands, minimal number of collisions between the helices, and proper formation of a hydrophobic core. The method is implemented as a module of ProteinShop. Our results show that it produces structures that are within 4–6 Å RMSD of the native one, regardless of the number of helices that need to be packed, though this number may increase if the protein has several helices between two consecutive strands in the sequence that pack on the sheet formed by these two strands. Proteins 2011; Published 2011 Wiley‐Liss, Inc.     

Quaternary conformational stability: The effect of reversible self‐association on the fibrillation of two insulin analogs

Under conditions relevant to the manufacturing of insulin (e.g., pH 3, room temperature), biosynthetic human insulin (BHI), and Lispro insulin (Lispro) require a nucleation step to initiate aggregation. However, upon seeding with preformed aggregates, both insulins rapidly aggregate into nonnative fibrils. Far ultraviolet circular dichroism (far‐UV CD) and second derivative Fourier transform infrared (2D‐FTIR) spectroscopic analyses show that the fibrillation process involves a change in protein secondary structure from α‐helical in native insulin to predominantly β‐sheet in the nonnative fibrils. After seeding, Lispro aggregates faster than BHI, likely because of a reduced propensity to reversibly self‐associate. Composition gradient multi‐angle light scattering (CG‐MALS) analyses show that Lispro is more monomeric than BHI, whereas their conformational stabilities measured by denaturant‐induced unfolding are statistically indistinguishable. For both BHI and Lispro, as the protein concentration increases, the apparent first‐order rate constant for soluble protein loss decreases. To explain these phenomena, we propose an aggregation model that assumes fibril growth through monomer addition with competitive inhibition by insulin dimers. Biotechnol. Bioeng. 2011;108: 2359–2370. © 2011 Wiley Periodicals, Inc.