De novo design of a monomeric three-stranded antiparallel beta-sheet

Here we describe the NMR conformational study of a 20-residue linear peptide designed to fold into a monomeric three-stranded antiparallel beta-sheet in aqueous solution. Experimental and statistical data on amino acid beta-turn and beta-sheet propensities, cross-strand side-chain interactions, solubility criteria, and our previous experience with beta-hairpins were considered for a rational selection of the peptide sequence. Sedimentation equilibrium measurements and NMR dilution experiments provide evidence that the peptide is monomeric. Analysis of 1H and 13C-NMR parameters of the peptide, in particular NOEs and chemical shifts, and comparison with data obtained for two 12-residue peptides encompassing the N- and C-segments of the designed sequence indicates that the 20-residue peptide folds into the expected conformation. Assuming a two-state model, the exchange kinetics between the beta-sheet and the unfolded peptide molecules is in a suitable range to estimate the folding rate on the basis of the NMR linewidths of several resonances. The time constant for the coil-beta-sheet transition is of the order of several microseconds in the designed peptide. Future designs based on this peptide system are expected to contribute greatly to our knowledge of the many factors involved in beta-sheet formation and stability.     

Folding and stability of the three-stranded beta-sheet peptide Betanova: insights from molecular dynamics simulations

The dynamics of the three-stranded beta-sheet peptide Betanova has been studied at four different temperatures (280, 300, 350, and 450 K by molecular dynamics simulation techniques, in explicit water. Two 20-ns simulations at 280 K indicate that the peptide remains very flexible under "folding" conditions sampling a range of conformations that together satisfy the nuclear magnetic resonance (NMR)-derived experimental constraints. Two simulations at 300 K (above the experimental folding temperature) of 20 ns each show partial formation of "native"-like structure, which also satisfies most of the NOE constraints at 280 K. At higher temperature, the presence of compact states, in which a series of hydrophobic contacts remain present, are observed. This is consistent with experimental observations regarding the role of hydrophobic contacts in determining the peptide's stability and in initiating the formation of turns and loops. A set of different structures is shown to satisfy NMR-derived distance restraints and a possible mechanism for the folding of the peptide into the NMR-determined structure is proposed.     

Folding of beta/alpha-unit scrambled forms of S. cerevisiae triosephosphate isomerase: Evidence for autonomy of substructure formation and plasticity of hydrophobic and hydrogen bonding interactions in core of (beta/alpha)8-barrel

We show that residues at the interfaces of protein-protein complexes have higher side-chain energy than other surface residues. Eight different sets of protein complexes were analyzed. For each protein pair, the complex structure was used to identify the interface residues in the unbound monomer structures. Side-chain energy was calculated for each surface residue in the unbound monomer using our previously developed scoring function.1 The mean energy was calculated for the interface residues and the other surface residues. In 15 of the 16 monomers, the mean energy of the interface residues was higher than that of other surface residues. By decomposing the scoring function, we found that the energy term of the buried surface area of non-hydrogen-bonded hydrophilic atoms is the most important factor contributing to the high energy of the interface regions. In spite of lacking hydrophilic residues, the interface regions were found to be rich in buried non-hydrogen-bonded hydrophilic atoms. Although the calculation results could be affected by the inaccuracy of the scoring function, patch analysis of side-chain energy on the surface of an isolated protein may be helpful in identifying the possible protein-protein interface. A patch was defined as 20 residues surrounding the central residue on the protein surface, and patch energy was calculated as the mean value of the side-chain energy of all residues in the patch. In 12 of the studied monomers, the patch with the highest energy overlaps with the observed interface. The results are more remarkable when only three residues with the highest energy in a patch are averaged to derive the patch energy. All three highest-energy residues of the top energy patch belong to interfacial residues in four of the eight small protomers. We also found that the residue with the highest energy score on the surface of a small protomer is very possibly the key interaction residue.     

Removal of the N-terminal hexapeptide from human beta2-microglobulin facilitates protein aggregation and fibril formation

The solution structure and stability of N-terminally truncated beta2-microglobulin (deltaN6beta2-m), the major modification in ex vivo fibrils, have been investigated by a variety of biophysical techniques. The results show that deltaN6beta2-m has a free energy of stabilization that is reduced by 2.5 kcal/mol compared to the intact protein. Hydrogen exchange of a mixture of the truncated and full-length proteins at microM concentrations at pH 6.5 monitored by electrospray mass spectrometry reveals that deltaN6beta2-m is significantly less protected than its wild-type counterpart. Analysis of deltaN6beta2-m by NMR shows that this loss of protection occurs in beta strands I, III, and part of II. At mM concentration gel filtration analysis shows that deltaN6beta2-m forms a series of oligomers, including trimers and tetramers, and NMR analysis indicates that strand V is involved in intermolecular interactions that stabilize this association. The truncated species of beta2-microglobulin was found to have a higher tendency to self-associate than the intact molecule, and unlike wild-type protein, is able to form amyloid fibrils at physiological pH. Limited proteolysis experiments and analysis by mass spectrometry support the conformational modifications identified by NMR and suggest that deltaN6beta2-m could be a key intermediate of a proteolytic pathway of beta2-microglobulin. Overall, the data suggest that removal of the six residues from the N-terminus of beta2-microglobulin has a major effect on the stability of the overall fold. Part of the tertiary structure is preserved substantially by the disulfide bridge between Cys25 and Cys80, but the pairing between beta-strands far removed from this constrain is greatly perturbed.     

Folding thermodynamics of three beta-sheet peptides: a model study

We study the folding thermodynamics of a beta-hairpin and two three-stranded beta-sheet peptides using a simplified sequence-based all-atom model, in which folding is driven mainly by backbone hydrogen bonding and effective hydrophobic attraction. The native populations obtained for these three sequences are in good agreement with experimental data. We also show that the apparent native population depends on which observable is studied; the hydrophobicity energy and the number of native hydrogen bonds give different results. The magnitude of this dependence matches well with the results obtained in two different experiments on the beta-hairpin.     

Synthesis and secondary structure of oligo(alpha-isobutyl beta,L-aspartate)s

The oligo(beta-peptide)s, hexa(alpha-isobutyl beta,L-aspartate) (Hex-AIBLA) and octa(alpha-isobutyl beta,L-aspartate) (Oct-AIBLA), were synthesized in solution by using standard coupling methods. Powder x-ray diffraction showed that the octamer crystallized in the two helical crystal forms known to exist in the homologous poly(beta-peptide), whereas the hexamer seemed to adopt an extended conformation. Both CD and 1H-NMR spectra of Oct-AIBLA in MeOH revealed the presence of a regularly folded conformation in this solvent, presumably the 3(14) helix. The helix-to-coil transition of Oct-AIBLA was observed to take place upon heating in both MeOH and CHCl3, in the second case associated with a not-well-defined aggregation-disaggregation process. The spectroscopic evidence obtained on the presence of folded structures in Hex-AIBLA were much weaker than for the octamer.     

Terminal sidechain packing of a designed β‐hairpin influences conformation and stability

While end capping in α‐helices is well understood, the concept of capping a β‐hairpin is a relatively recent development; to date, favorable Coulombic interactions are the only example of sidechains at the termini influencing the overall stability of a β‐hairpin. While cross‐strand hydrophobic residues generally provide hairpin stabilization, particular when flanking the turn region, those remote from this location appear to provide little stabilization. While probing for an optimal residue at a hydrogen bond position near the terminus of a designed β‐hairpin a conservative, hydrophobic, V → I mutation was observed to not only result in a significant change in fold population but also effected major changes in the structuring shifts at numerous sites in the peptide. Mutational studies reveal that there is an interaction between the sidechain at this H‐bonded site and the sidechain at the C‐terminal non‐H‐bonded site of the hairpin. This interaction, which appears to be hydrophobic in character, requires a highly twisted hairpin structure. Modifications at the C‐terminal site, for example an E → A mutation (ΔΔG_U = 6 kJ/mol), have profound affects on fold structure and stability. The data suggests that this may be a case of hairpin end capping by the formation of a hydrophobic cluster. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 557–564, 2009. 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     

Dimerization of beta B2-crystallin: the role of the linker peptide and the N- and C-terminal extensions.

beta B2- and gamma B-crystallins of vertebrate eye lens are 2-domain proteins in which each domain consists of 2 Greek key motifs connected by a linker peptide. Although the folding topologies of beta B2- and gamma B-domains are very similar, gamma B-crystallin is always monomeric, whereas beta B2-crystallin associates to homodimers. It has been suggested that the linker or the protruding N- and C-terminal arms of beta B2-crystallin (not present in gamma B) are a necessary requirement for this association. In order to investigate the role of these segments for dimerization, we constructed two beta B2 mutants. In the first mutant, the linker peptide was replaced with the one from gamma B (beta B2 gamma L). In the second mutant, the N- and C-terminal arms of 15- and 12-residues length were deleted (beta B2 delta NC). The beta B2 gamma L mutant is monomeric, whereas the beta B2 delta NC mutant forms dimers and tetramers that cannot be interconverted without denaturation. The spectral properties of the beta B2 mutants, as well as their stabilities against denaturants, resemble those of wild-type beta B2-crystallin, thus indicating that the overall peptide fold of the subunits is not changed significantly. We conclude that the peptide linker in beta B2-crystallin is necessary for dimerization, whereas the N- and C-terminal arms appear to be involved in preventing the formation of higher homo-oligomers.     

The ensemble folding kinetics of the FBP28 WW domain revealed by an all-atom Monte Carlo simulation in a knowledge-based potential

In this work, we apply a detailed all‐atom model with a transferable knowledge‐based potential to study the folding kinetics of Formin‐Binding protein, FBP28, which is a canonical three‐stranded β‐sheet WW domain. Replica exchange Monte Carlo simulations starting from random coils find native‐like (Cα RMSD of 2.68 Å) lowest energy structure. We also study the folding kinetics of FBP28 WW domain by performing a large number of ab initio Monte Carlo folding simulations. Using these trajectories, we examine the order of formation of two β‐hairpins, the folding mechanism of each individual β‐hairpin, and transition state ensemble (TSE) of FBP28 WW domain and compare our results with experimental data and previous computational studies. To obtain detailed structural information on the folding dynamics viewed as an ensemble process, we perform a clustering analysis procedure based on graph theory. Further, a rigorous P_fold analysis is used to obtain representative samples of the TSEs showing good quantitative agreement between experimental and simulated Φ values. Our analysis shows that the turn structure between first and second β strands is a partially stable structural motif that gets formed before entering the TSE in FBP28 WW domain and there exist two major pathways for the folding of FBP28 WW domain, which differ in the order and mechanism of hairpin formation. Proteins 2011. © 2011 Wiley‐Liss, Inc.     

Beta-hairpin formation in aqueous solution and in the presence of trifluoroethanol: a (1)H and (13)C nuclear magnetic resonance conformational study of designed peptides

In order to check our current knowledge on the principles involved in beta-hairpin formation, we have modified the sequence of a 3:5 beta-hairpin forming peptide with two different purposes, first to increase the stability of the formed 3:5 beta-hairpin, and second to convert the 3:5 beta-hairpin into a 2:2 beta-hairpin. The conformational behavior of the designed peptides was investigated in aqueous solution and in 30% trifluoroethanol (TFE) by analysis of the following nuclear magnetic resonance (NMR) parameters: nuclear Overhauser effect (NOE) data, and C(alpha)H, (13)C(alpha), and (13)C(beta) conformational shifts. From the differences in the ability to adopt beta-hairpin structures in these peptides, we have arrived to the following conclusions: (i) beta-Hairpin population increases with the statistical propensity of residues to occupy each turn position. (ii) The loop length, and in turn, the beta-hairpin type, can be modified as a function of the type of turn favored by the loop sequence. These two conclusions reinforce previous results about the importance of beta-turn sequence in beta-hairpin folding. (iii) Side-chain packing on each face of the beta-sheet may play a major role in beta-hairpin stability; hence simplified analysis in terms of isolated pair interactions and intrinsic beta-sheet propensities is insufficient. (iv) Contributions to beta-hairpin stability of turn and strand sequences are not completely independent. (v) The burial of hydrophobic surface upon beta-hairpin formation that, in turn, depends on side-chain packing also contributes to beta-hairpin stability. (vi) As previously observed, TFE stabilizes beta-hairpin structures, but the extent of the contribution of different factors to beta-hairpin formation is sometimes different in aqueous solution and in 30% TFE.     

A detailed unfolding pathway of a (beta/alpha)8-barrel protein as studied by molecular dynamics simulations

The (beta/alpha)(8)-barrel is the most common protein fold. Similar structural properties for folding intermediates of (beta/alpha)(8)-barrel proteins involved in tryptophan biosynthesis have been reported in a number of experimental studies; these intermediates have the last two beta-strands and three alpha-helices partially unfolded, with other regions of the polypeptide chain native-like in conformation. To investigate the detailed folding/unfolding pathways of these (beta/alpha)(8)-barrel proteins, temperature-induced unfolding simulations of N-(5'-phosphoribosyl)anthranilate isomerase from Escherichia coli were carried out using a special-purpose parallel computer system. Unfolding simulations at five different temperatures showed a sequential unfolding pathway comprised of several events. Early events in unfolding involved disruption of the last two strands and three helices, producing an intermediate ensemble similar to those detected in experimental studies. Then, denaturation of the first two betaalpha units and separation of the sixth strand from the fifth took place independently. The remaining central betaalphabetaalphabeta module persisted the longest during all simulations, suggesting an important role for this module as the incipient folding scaffold. Our simulations also predicted the presence of a nucleation site, onto which several hydrophobic residues condensed forming the foundation for the central betaalphabetaalphabeta module.     

Determination of alpha-helix and beta-sheet stability in the solid state: a solid-state NMR investigation of poly(L-alanine)

The relative stability of alpha-helix and beta-sheet secondary structure in the solid state was investigated using poly(L-alanine) (PLA) as a model system. Protein folding and stability has been well studied in solution, but little is known about solid-state environments, such as the core of a folded protein, where peptide packing interactions are the dominant factor in determining structural stability. (13)C cross-polarization with magic angle spinning (CPMAS) NMR spectroscopy was used to determine the backbone conformation of solid powder samples of 15-kDa and 21.4-kDa PLA before and after various sample treatments. Reprecipitation from helix-inducing solvents traps the alpha-helical conformation of PLA, although the method of reprecipitation also affects the conformational distribution. Grinding converts the secondary structure of PLA to a final steady-state mixture of 55% beta-sheet and 45% alpha-helix at room temperature regardless of the initial secondary structure. Grinding PLA at liquid nitrogen temperatures leads to a similar steady-state mixture with 60% beta-sheet and 40% alpha-helix, indicating that mechanical shear force is sufficient to induce secondary structure interconversion. Cooling the sample in liquid nitrogen or subjecting it to high pressure has no effect on secondary structure. Heating the sample without grinding results in equilibration of secondary structure to 50% alpha-helix/50% beta-sheet at 100 degrees C when starting from a mostly alpha-helical state. No change was observed upon heating a beta-sheet sample, perhaps due to kinetic effects and the different heating rate used in the experiments. These results are consistent with beta-sheet approximately 260 J/mol more stable than alpha-helix in solid-state PLA.     

The importance of surface loops for stabilizing an eightfold beta alpha barrel protein.

An important step in understanding how a protein folds is to determine those regions of the sequence that are critical to both its stability and its folding pathway. We chose phosphoribosyl anthranilate isomerase from Escherichia coli, which is a monomeric representative of the (beta alpha)8 barrel family of proteins, to construct a variant that carries an internal tandem duplication of the fifth beta alpha module. This (beta alpha)9 variant was enzymically active and therefore must have a wild-type (beta alpha)8 core. It had a choice a priori to fold to three different folding frames, which are distinguished by carrying the duplicated segment as an insert into one out of three different loops. Steady-state kinetic constants, the fluorescence properties of a crucial tryptophan residue, and limited proteolysis showed that the stable (beta alpha)9 variant carries the insertion between beta-strand 5 and alpha-helix 5. This preference can be explained by the important role of loops between alpha helices and beta strands in stabilizing the structure of the enzyme.     

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.     

Structures of single‐layer β‐sheet proteins evolved from β‐hairpin repeats

Free‐standing single‐layer β‐sheets are extremely rare in naturally occurring proteins, even though β‐sheet motifs are ubiquitous. Here we report the crystal structures of three homologous, single‐layer, anti‐parallel β‐sheet proteins, comprised of three or four twisted β‐hairpin repeats. The structures reveal that, in addition to the hydrogen bond network characteristic of β‐sheets, additional hydrophobic interactions mediated by small clusters of residues adjacent to the turns likely play a significant role in the structural stability and compensate for the lack of a compact hydrophobic core. These structures enabled identification of a family of secreted proteins that are broadly distributed in bacteria from the human gut microbiome and are putatively involved in the metabolism of complex carbohydrates. A conserved surface patch, rich in solvent‐exposed tyrosine residues, was identified on the concave surface of the β‐sheet. These new modular single‐layer β‐sheet proteins may serve as a new model system for studying folding and design of β‐rich proteins.     

The role of proline‐containing peptide triads in β‐sheet formation: A kinetic study

The design of biomimetic materials through molecular self‐assembly is a growing area of modern nanotechnology. With problems of protein folding, self‐assembly, and sequence–structure relationships as essential in nanotechnology as in biology, the effect of the nucleation of β‐hairpin formation by proline on the folding process has been investigated in model studies. Previously such studies were limited to investigations of the influence of proline on the formation of turns in short peptide sequences. The effect of proline‐based triads on the folding of an 11‐kDa amyloidogenic peptide GH₆[(GA)₃GY(GA)₃GE]₈GAH₆ ( YE8 ) was investigated by selective substitution of the proline‐substituted triads at the γ‐turn sites. The folding and fibrillation of the singly proline‐substituted polypeptides, e.g., GH6? [(GA)3GY(GA)3GE]7(GA)3GY(GA)3PD? GAH6 ( 8PD ), and doubly proline‐substituted polypeptides, e.g., GH6? [(GA)3GY(GA)3GE]3(GA)3GY(GA)3PD[(GA)3GY(GA)3GE]3(GA)3GY(GA)3PD? GAH6 ( 4,8PD ), were directly monitored by circular dichroism and deep UV resonance Raman and fluorescence spectroscopies. These findings were used to identify the essential folding domains, i.e., the minimum number of β‐strands necessary for stable folding. These experimental findings may be especially useful in the design and construction of peptidic materials for a wide range of applications as well as in understanding the mechanisms of folding critical to fibril formation. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 339–350, 2015.     

Folding of intracellular retinol and retinoic acid binding proteins

The folding mechanisms of cellular retinol binding protein II (CRBP II), cellular retinoic acid binding protein I (CRABP I), and cellular retinoic acid binding protein II (CRABP II) were examined. These beta-sheet proteins have very similar structures and higher sequence homologies than most proteins in this diverse family. They have similar stabilities and show completely reversible folding at equilibrium with urea as a denaturant. The unfolding kinetics of these proteins were monitored during folding and unfolding by circular dichroism (CD) and fluorescence. During unfolding, CRABP II showed no intermediates, CRABP I had an intermediate with nativelike secondary structure, and CRBP II had an intermediate that lacked secondary structure. The refolding kinetics of these proteins were more similar. Each protein showed a burst-phase change in intensity by both CD and fluorescence, followed by a single observed phase by both CD and fluorescence and one or two additional refolding phases by fluorescence. The fluorescence spectral properties of the intermediate states were similar and suggested a gradual increase in the amount of native tertiary structure present for each step in a sequential path. However, the rates of folding differed by as much as 3 orders of magnitude and were slower than those expected from the contact order and topology of these proteins. As such, proteins with the same final structure may not follow the same route to the native state.     

Early formation of a beta hairpin during folding of staphylococcal nuclease H124L as detected by pulsed hydrogen exchange

Pulsed hydrogen exchange methods were used to follow the formation of structure during the refolding of acid-denatured staphylococcal nuclease containing a stabilizing Leu substitution at position 124 (H124L SNase). The protection of more than 60 backbone amide protons in uniformly (15)N-labeled H124L SNase was monitored as a function of refolding time by heteronuclear two-dimensional NMR spectroscopy. As found in previous studies of staphylococcal nuclease, partial protection was observed for a subset of amide protons even at the earliest folding time point (10 msec). Protection indicative of marginally stable hydrogen-bonded structure in an early folding intermediate was observed at over 30 amide positions located primarily in the beta-barrel and to a lesser degree in the alpha-helical domain of H124L SNase. To further characterize the folding intermediate, protection factors for individual amide sites were measured by varying the pH of the labeling pulse at a fixed refolding time of 16 msec. Protection factors >5.0 were observed only for amide positions in a beta-hairpin formed by strands 2 and 3 of the beta-barrel domain and a single site near the C-terminus. The results indicate that formation of stable hydrogen-bonded structure in a core region of the beta-sheet is among the earliest structural events in the folding of SNase and may serve as a nucleation site for further structure formation.     

Determination of protein fold class from Raman or Raman optical activity spectra using random forests

Knowledge of the fold class of a protein is valuable because fold class gives an indication of protein function and evolution. Fold class can be accurately determined from a crystal structure or NMR structure, though these methods are expensive, time-consuming, and inapplicable to all proteins. In contrast, vibrational spectra [infra-red, Raman, or Raman optical activity (ROA)] are rapidly obtained for proteins under wide range of biological molecules under diverse experimental and physiological conditions. Here, we show that the fold class of a protein can be determined from Raman or ROA spectra by converting a spectrum into data of 10 cm⁻¹ bin widths and applying the random forest machine learning algorithm. Spectral data from 605 and 1785 cm⁻¹ were analyzed, as well as the amide I, II, and III regions in isolation and in combination. ROA amide II and III data gave the best performance, with 33 of 44 proteins assigned to one of the correct four top-level structural classification of proteins (SCOP) fold class (all α, all β, α and β, and disordered). The method also shows which spectral regions are most valuable in assigning fold class.