Redesign of a four-helix bundle protein by phage display coupled with proteolysis and structural characterization by NMR and X-ray crystallography

To test whether it is practical to use phage display coupled with proteolysis for protein design, we used this approach to convert a partially unfolded four-helix bundle protein, apocytochrome b(562), to a stably folded four-helix bundle protein. Four residues expected to form a hydrophobic core were mutated. One residue was changed to Trp to provide a fluorescence probe for studying the protein's physical properties and to partially fill the void left by the heme. The other three positions were randomly mutated. In addition, another residue in the region to be redesigned was substituted with Arg to provide a specific cutting site for protease Arg-c. This library of mutants was displayed on the surface of phage and challenged with protease Arg-c to select stably folded proteins. The consensus sequence that emerged from the selection included hydrophobic residues at only one of the three positions and non-hydrophobic residues at the other two. Nevertheless, the selected proteins were thermodynamically very stable. The structure of a selected protein was characterized using multi-dimensional NMR. All four helices were formed in the structure. Further, site-directed mutagenesis was used to change one of the two non-hydrophobic residues to a hydrophobic residue, which increased the stability of the protein, indicating that the selection result was not based solely on the protein's global stability and that local structural characteristics may also govern the selection. This conclusion is supported by the crystal structure of another mutant that has two hydrophobic residues substituted for the two non-hydrophobic residues. These results suggest that the hydrophobic interactions in the core are not sufficient to dictate the selection and that the location of the cutting site of the protease also influences the selection of structures.     

Core-directed protein design. II. Rescue of a multiply mutated and destabilized variant of ubiquitin.

We have applied the method described in the preceding paper [Finucane, M. D., et al. (1999) Biochemistry 38, 11604-11612], namely, stability-based selection using phage display, to explore the sequence requirements for packing in the hydrophobic core of ubiquitin. In contrast to the parent protein, which was a structurally compromised mutant, the selected variants could be overexpressed and purified in yields for structural studies. In particular, CD and NMR measurements showed that the selectants folded correctly to stable native-like structures. These points demonstrate the utility of our core-directed method for stabilizing and redesigning proteins. In addition and in contrast to foregoing studies on other proteins, which suggest that hydrophobic cores permit substitutions provided that hydrophobicity and core volumes are generally conserved, we find that the core of ubiquitin is surprisingly intolerant of amino acid substitutions; variants that survived our selection showed a clear consensus for the wild-type sequence. It is probable that our results differed from those from other groups for two reasons. First, ubiquitin may be unusual in that it has strict sequence requirements for its structure and stability. We discuss this result in light of sequence conservation in the eukaryotic ubiquitins and proteins of the ubiquitin structural superfamily. Second, our mutants were selected solely on the basis of stability, in contrast to the other studies that rely on function-based selection. The latter may lead to proteins that are more plastic and tolerant of substitutions.     

The role of aromatic residues in the hydrophobic core of the villin headpiece subdomain

Small autonomously folding proteins are of interest as model systems to study protein folding, as the same molecule can be used for both experimental and computational approaches. The question remains as to how well these minimized peptide model systems represent larger native proteins. For example, is the core of a minimized protein tolerant to mutation like larger proteins are? Also, do minimized proteins use special strategies for specifying and stabilizing their folded structure? Here we examine these questions in the 35‐residue autonomously folding villin headpiece subdomain (VHP subdomain). Specifically, we focus on a cluster of three conserved phenylalanine (F) residues F47, F51, and F58, that form most of the hydrophobic core. These three residues are oriented such that they may provide stabilizing aromatic–aromatic interactions that could be critical for specifying the fold. Circular dichroism and 1D‐NMR spectroscopy show that point mutations that individually replace any of these three residues with leucine were destabilized, but retained the native VHP subdomain fold. In pair‐wise replacements, the double mutant that retains F58 can adopt the native fold, while the two double mutants that lack F58 cannot. The folding of the double mutant that retains F58 demonstrates that aromatic–aromatic interactions within the aromatic cluster are not essential for specifying the VHP subdomain fold. The ability of the VHP subdomain to tolerate mutations within its hydrophobic core indicates that the information specifying the three dimensional structure is distributed throughout the sequence, as observed in larger proteins. Thus, the VHP subdomain is a legitimate model for larger, native proteins.     

De novo design of the hydrophobic core of ubiquitin.

We have previously reported the development and evaluation of a computational program to assist in the design of hydrophobic cores of proteins. In an effort to investigate the role of core packing in protein structure, we have used this program, referred to as Repacking of Cores (ROC), to design several variants of the protein ubiquitin. Nine ubiquitin variants containing from three to eight hydrophobic core mutations were constructed, purified, and characterized in terms of their stability and their ability to adopt a uniquely folded native-like conformation. In general, designed ubiquitin variants are more stable than control variants in which the hydrophobic core was chosen randomly. However, in contrast to previous results with 434 cro, all designs are destabilized relative to the wild-type (WT) protein. This raises the possibility that beta-sheet structures have more stringent packing requirements than alpha-helical proteins. A more striking observation is that all variants, including random controls, adopt fairly well-defined conformations, regardless of their stability. This result supports conclusions from the cro studies that non-core residues contribute significantly to the conformational uniqueness of these proteins while core packing largely affects protein stability and has less impact on the nature or uniqueness of the fold. Concurrent with the above work, we used stability data on the nine ubiquitin variants to evaluate and improve the predictive ability of our core packing algorithm. Additional versions of the program were generated that differ in potential function parameters and sampling of side chain conformers. Reasonable correlations between experimental and predicted stabilities suggest the program will be useful in future studies to design variants with stabilities closer to that of the native protein. Taken together, the present study provides further clarification of the role of specific packing interactions in protein structure and stability, and demonstrates the benefit of using systematic computational methods to predict core packing arrangements for the design of proteins.     

Hydrophobic clusters in protein structures

Globular proteins fold such that the hydrophobic groups are packed inside forming hydrophobic clusters, and the hydrophilic groups are present on the surface. In this article we analyze clusters of hydrophobic groups of atoms in 781 protein structures selected from the PDB. Our analysis showed that every structure consists of two types of clusters: at least one large cluster that forms the hydrophobic core and probably dictates the protein fold; and numerous smaller clusters, which might be involved in the stabilization of the fold. We also analyzed the preference of the hydrophobic groups in each of the amino acids toward forming hydrophobic clusters. We find that hydrophobic groups from the hydrophilic amino acids also contribute toward cluster formation. Proteins 2008. © 2008 Wiley‐Liss, Inc.     

Frequencies of amino acid strings in globular protein sequences indicate suppression of blocks of consecutive hydrophobic residues

Patterns of hydrophobic and hydrophilic residues play a major role in protein folding and function. Long, predominantly hydrophobic strings of 20-22 amino acids each are associated with transmembrane helices and have been used to identify such sequences. Much less attention has been paid to hydrophobic sequences within globular proteins. In prior work on computer simulations of the competition between on-pathway folding and off-pathway aggregate formation, we found that long sequences of consecutive hydrophobic residues promoted aggregation within the model, even controlling for overall hydrophobic content. We report here on an analysis of the frequencies of different lengths of contiguous blocks of hydrophobic residues in a database of amino acid sequences of proteins of known structure. Sequences of three or more consecutive hydrophobic residues are found to be significantly less common in actual globular proteins than would be predicted if residues were selected independently. The result may reflect selection against long blocks of hydrophobic residues within globular proteins relative to what would be expected if residue hydrophobicities were independent of those of nearby residues in the sequence.     

Using physical features of protein core packing to distinguish real proteins from decoys

The ability to consistently distinguish real protein structures from computationally generated model decoys is not yet a solved problem. One route to distinguish real protein structures from decoys is to delineate the important physical features that specify a real protein. For example, it has long been appreciated that the hydrophobic cores of proteins contribute significantly to their stability. We used two sources to obtain datasets of decoys to compare with real protein structures: submissions to the biennial Critical Assessment of protein Structure Prediction competition, in which researchers attempt to predict the structure of a protein only knowing its amino acid sequence, and also decoys generated by 3DRobot, which have user‐specified global root‐mean‐squared deviations from experimentally determined structures. Our analysis revealed that both sets of decoys possess cores that do not recapitulate the key features that define real protein cores. In particular, the model structures appear more densely packed (because of energetically unfavorable atomic overlaps), contain too few residues in the core, and have improper distributions of hydrophobic residues throughout the structure. Based on these observations, we developed a feed‐forward neural network, which incorporates key physical features of protein cores, to predict how well a computational model recapitulates the real protein structure without knowledge of the structure of the target sequence. By identifying the important features of protein structure, our method is able to rank decoy structures with similar accuracy to that obtained by state‐of‐the‐art methods that incorporate many additional features. The small number of physical features makes our model interpretable, emphasizing the importance of protein packing and hydrophobicity in protein structure prediction.     

Folding specificity induced by loop stiffness

To test the importance of loop stiffness in restricting the heterogeneity of transition state ensemble, we relaxed the distal loop of 10 unstable redesigned hydrophobic core mutants of alpha-spectrin SH3 domain. This was achieved by replacing Asp48 by Gly at the tip of the distal hairpin. Although the change was local in nature, the effect on stabilization was not uniform across the core mutants tested. There is an inverse rough correlation between the stabilization and the increase in buried hydrophobic volume, with respect to the wild type. Interestingly enough, proteins that although unstable are properly folded become molten globule-like after relaxation of the distal loop. These results highlight the importance of stiffness in restricting the conformational heterogeneity of a protein during the folding reaction. An interplay between unspecific hydrophobic interactions and constraint induced by polar interactions, or in this case local stiffness, is essential to achieve a well-ordered folded structure.     

N-Terminus-Mediated Solution Structure of Dimerization Domain of PRC1

Simple SummaryThe solution structure of the N-terminal domain of Protein Regulator of Cytokinesis 1 (PRC1) was determined, and compared with the previously published crystal structure, significant differences were found. Extensive analyses were carried out to find the true reason for the differences between the solution and crystal structures, we discovered that this might be related to the conformation of residue M1, which is buried in the protein core of the solution structure, while situated outside of the hydrophobic core in the crystal structure. In this study, we have carried out a series of examinations using various methods and confirmed that the N terminal conformation is the key point in better describing the structure of PRC1 dimerization domain under solution conditions.AbstractMicrotubule-associated proteins (MAPs) are essential for the accurate division of a cell into two daughter cells. These proteins target specific microtubules to be incorporated into the spindle midzone, which comprises a special array of microtubules that initiate cytokinesis during anaphase. A representative member of the MAPs is Protein Regulator of Cytokinesis 1 (PRC1), which self-multimerizes to cross-link microtubules, the malfunction of which might result in cancerous cells. The importance of PRC1 multimerization makes it a popular target for structural studies. The available crystal structure of PRC1 has low resolution (>3 Å) and accuracy, limiting a better understanding of the structure-related functions of PRC1. Therefore, we used NMR spectroscopy to better determine the structure of the dimerization domain of PRC1. The NMR structure shows that the PRC1 N terminus is crucial to the overall structure integrity, but the crystal structure bespeaks otherwise. We systematically addressed the role of the N terminus by generating a series of mutants in which N-terminal residues methionine (Met1) and arginine (Arg2) were either deleted, extended or substituted with other rationally selected amino acids. Each mutant was subsequently analyzed by NMR spectroscopy or fluorescence thermal shift assays for its structural or thermal stability; we found that N-terminal perturbations indeed affected the overall protein structure and that the solution structure better reflects the conformation of PRC1 under solution conditions. These results reveal that the structure of PRC1 is governed by its N terminus through hydrophobic interactions with other core residues, such hitherto unidentified N-terminal conformations might shed light on the structure–function relationships of PRC1 or other proteins. Therefore, our study is of major importance in terms of identifying a novel structural feature and can further the progress of protein folding and protein engineering.     

The Hydrophobic Core Sequence Modulates the Neurotoxic and Secondary Structure Properties of the Prion Peptide 106-126

The neurodegeneration seen in spongiform encephalopathies is believed to be mediated by protease-resistant forms of the prion protein (PrP). A peptide encompassing residues 106-126 of human PrP has been shown to be neurotoxic in vitro. The neurotoxicity of PrP106-126 appears to be dependent upon its adoption of an aggregated fibril structure. To examine the role of the hydrophobic core, AGAAAAGA, on PrP106-126 toxicity, we performed structure-activity analyses by substituting two or more hydrophobic residues for the hydrophilic serine residue to decrease its hydrophobicity. A peptide with a deleted alanine was also synthesized. We found all the peptides except the deletion mutant were no longer toxic on mouse cerebellar neuronal cultures. Circular dichroism analysis showed that the nontoxic PrP peptides had a marked decrease in beta-sheet structure. In addition, the mutants had alterations in aggregability as measured by turbidity, Congo red binding, and fibril staining using electron microscopy. These data show that the hydrophobic core sequence is important for PrP106-126 toxicity probably by influencing its assembly into a neurotoxic structure. The hydrophobic sequence may similarly affect aggregation and toxicity observed in prion diseases.     

Simulating the minimum core for hydrophobic collapse in globular proteins.

To investigate the nature of hydrophobic collapse considered to be the driving force in protein folding, we have simulated aqueous solutions of two model hydrophobic solutes, methane and isobutylene. Using a novel methodology for determining contacts, we can precisely follow hydrophobic aggregation as it proceeds through three stages: dispersed, transition, and collapsed. Theoretical modeling of the cluster formation observed by simulation indicates that this aggregation is cooperative and that the simulations favor the formation of a single cluster midway through the transition stage. This defines a minimum solute hydrophobic core volume. We compare this with protein hydrophobic core volumes determined from solved crystal structures. Our analysis shows that the solute core volume roughly estimates the minimum core size required for independent hydrophobic stabilization of a protein and defines a limiting concentration of nonpolar residues that can cause hydrophobic collapse. These results suggest that the physical forces driving aggregation of hydrophobic molecules in water is indeed responsible for protein folding.     

The relationship between conservation, thermodynamic stability, and function in the SH3 domain hydrophobic core

To investigate the relationships between sequence conservation, protein stability, and protein function, we have measured the thermodynamic stability, folding kinetics, and in vitro peptide-binding activity of a large number of single-site substitutions in the hydrophobic core of the Fyn SH3 domain. Comparison of these data to that derived from an analysis of a large alignment of SH3 domain sequences revealed a very good correlation between the distinct pattern of conservation observed at each core position and the thermodynamic stability of mutants. Conservation was also found to correlate well with the unfolding rates of mutants, but not to the folding rates, suggesting that evolution selects more strongly for optimal native state packing interactions than for maximal folding rates. Structural analysis suggests that residue-residue core packing interactions are very similar in all SH3 domains, which provides an explanation for the correlation between conservation and mutant stability effects studied in a single SH3 domain. We also demonstrate a correlation between stability and the in vivo activity of mutants, and between conservation and activity. However, the relationship between conservation and activity was very strong only for the three most conserved hydrophobic core positions. The weaker correlation between activity and conservation seen at the other seven core positions indicates that maintenance of protein stability is the dominant selective pressure at these positions. In general, the pattern of conservation at hydrophobic core positions appears to arise from conserved packing constraints, and can be effectively utilized to predict the destabilizing effects of amino acid substitutions.     

Designing proteins from the inside out

Globular proteins are characterized by the specific and tight packing of hydrophobic side-chains in the so-called "hydrophobic core." Formation of the core is key in folding, stabilization, and conformational specificity. The critical role of hydrophobic cores in maintaining the highly ordered structures present in natural proteins justifies the tremendous efforts devoted to their redesign. Both experimental and computational combinatorial-based approaches have been reported in the last years as powerful protein design tools. These manage to explore large regions of the sequence/conformational space, allowing the search for alternative protein core arrangements displaying native-like properties. The overall results obtained from core design projects have contributed significantly to our present knowledge of protein folding and function. In addition, core design has worked as a benchmark for the development of ambitious protein design projects that nowadays are allowing the de novo design of novel protein structures and functions.     

Mutational analysis of the channel and loop sequences of human ferritin H-chain.

Human ferritin H-chain mutants were obtained by engineering the recombinant protein expressed by Escherichia coli. The mutagenesis were directed to the C-terminal sequence forming the hydrophobic channel, to the hydrophilic channel and to the loop sequence. The mutants were analysed for extent of expression, for stability, for capacity to incorporate iron and for kinetics of iron uptake and iron oxidation. Of the 22 mutants analysed only two with deletions of single residues in the loop sequence and one with deletion of the last 28 amino acid residues did not assemble into ferritin-like proteins. The other mutants assembled correctly and showed similar chemical/physical properties to the wild-type; they included duplication of an 18-amino acid-residue stretch, deletion of the last 22 and the last seven residues and various mutations of single amino acid residues. Two mutants with extensive alteration in the C-terminal sequence had a diminished thermostability associated with incapability to incorporate iron though they still catalysed iron oxidation. The mutants with alterations of the sequence around the hydrophilic channel showed diminished iron uptake and oxidation kinetics, together with a slightly larger apparent molecular size. The results indicate (i) that two of the sequences are important for ferritin assembly/stability, (ii) that the presence of the hydrophobic channel is essential for formation of the iron core and (iii) that the sites of iron interaction and the path of iron penetration into ferritin remain unidentified.     

A tripartite structure of the signals that determine protein insertion into the endoplasmic reticulum membrane

Multilineage colony stimulating factor is a secretory protein with a cleavable signal sequence that is unusually long and hydrophobic. Using molecular cloning techniques we exchanged sequences NH2- or COOH-terminally flanking the hydrophobic signal sequence. Such modified fusion proteins still inserted into the membrane but their signal sequence was not cleaved. Instead the proteins were now anchored in the membrane by the formerly cleaved signal sequence (signal-anchor sequence). They exposed the NH2 terminus on the exoplasmic and the COOH terminus on the cytoplasmic side of the membrane. We conclude from our results that hydrophilic sequences flanking the hydrophobic core of a signal sequence can determine cleavage by signal peptidase and insertion into the membrane. It appears that negatively charged amino acid residues close to the NH2 terminal side of the hydrophobic segment are compatible with translocation of this segment across the membrane. A tripartite structure is proposed for signal-anchor sequences: a hydrophobic core region that mediates targeting to and insertion into the ER membrane and flanking hydrophilic segments that determine the orientation of the protein in the membrane.     

Hydrophobic core formation in protein complex of cathepsin

The “fuzzy oil drop” model assumes that the idealized hydrophobic core in a protein body can be described by a 3D Gauss function. The structure of the 1ICF protein (cathepsin), which participates in the proteolysis process and has cysteine-type peptidase activity, has been analyzed on the basis of the “fuzzy oil drop” model. The authors have determined the contribution of individual exon fragments to the creation of a common hydrophobic core and assessed the involvement of each chain in this process, depending on the number of complexed chains. Quantitative assessment of exons, chains, dimers, and the whole complex suggest that each of these units plays a different role in shaping the protein’s hydrophobic core.     

Effect of subdomain interactions on methyl group dynamics in the hydrophobic core of villin headpiece protein

Thermostable villin headpiece protein (HP67) consists of the N‐terminal subdomain (residues 10–41) and the autonomously folding C‐terminal subdomain (residues 42–76) which pack against each other to form a structure with a unified hydrophobic core. The X‐ray structures of the isolated C‐terminal subdomain (HP36) and its counterpart in HP67 are very similar for the hydrophobic core residues. However, fine rearrangements of the free energy landscape are expected to occur because of the interactions between the two subdomains. We detect and characterize these changes by comparing the µs‐ms time scale dynamics of the methyl‐bearing side chains in isolated HP36 and in HP67. Specifically, we probe three hydrophobic side chains at the interface of the two subdomains (L42, V50, and L75) as well as at two residues far from the interface (L61 and L69). Solid‐state deuteron NMR techniques are combined with computational modeling for the detailed characterization of motional modes in terms of their kinetic and thermodynamic parameters. The effect of interdomain interactions on side chain dynamics is seen for all residues but L75. Thus, changes in dynamics because of subdomain interactions are not confined to the site of perturbation. One of the main results is a two‐ to threefold increase in the value of the activation energies for the rotameric mode of motions in HP67 compared with HP36. Detailed analysis of configurational entropies and heat capacities complement the kinetic view of the degree of the disorder in the folded state.     

Selection and analyses of variants of a designed protein suggest importance of hydrophobicity of partially buried sidechains for protein stability at high temperatures

Computationally designed proteins of high stability provide specimen in addition to natural proteins for the study of sequence‐structure stability relationships at the very high end of protein stability spectrum. The melting temperature of E_1r26, a protein we previously designed using the A Backbone‐based Amino aCid Usage Survey (ABACUS) sequence design program, is above 110 °C, more than 50 °C higher than that of the natural thioredoxin protein whose backbone (PDB ID 1R26 ) has been used as the design target. Using an experimental selection approach, we obtained variants of E_1r26 that remain folded but are of reduced stability, including one whose unfolding temperature and denaturing guanidine concentration are similar to those of 1r26. The mutant unfolds with a certain degree of cooperativity. Its structure solved by X‐ray crystallography agrees with that of 1r26 by a root mean square deviation of 1.3 Å, adding supports to the accuracy of the ABACUS method. Analyses of intermediate mutants indicate that the substitution of two partially buried hydrophobic residues (isoleucine and leucine) by polar residues (threonine and serine, respectively) are responsible for the dramatic change in the unfolding temperature. It is suggested that the effects of mutations located in rigid secondary structure regions, but not those in loops, may be well predicted through ABACUS mutation energy analysis. The results also suggest that hydrophobic effects involving intermediately buried sidechains can be critically important for protein stability at high temperatures.     

De novo design of the hydrophobic cores of proteins.

We have developed and experimentally tested a novel computational approach for the de novo design of hydrophobic cores. A pair of computer programs has been written, the first of which creates a "custom" rotamer library for potential hydrophobic residues, based on the backbone structure of the protein of interest. The second program uses a genetic algorithm to globally optimize for a low energy core sequence and structure, using the custom rotamer library as input. Success of the programs in predicting the sequences of native proteins indicates that they should be effective tools for protein design. Using these programs, we have designed and engineered several variants of the phage 434 cro protein, containing five, seven, or eight sequence changes in the hydrophobic core. As controls, we have produced a variant consisting of a randomly generated core with six sequence changes but equal volume relative to the native core and a variant with a "minimalist" core containing predominantly leucine residues. Two of the designs, including one with eight core sequence changes, have thermal stabilities comparable to the native protein, whereas the third design and the minimalist protein are significantly destabilized. The randomly designed control is completely unfolded under equivalent conditions. These results suggest that rational de novo design of hydrophobic cores is feasible, and stress the importance of specific packing interactions for the stability of proteins. A surprising aspect of the results is that all of the variants display highly cooperative thermal denaturation curves and reasonably dispersed NMR spectra. This suggests that the non-core residues of a protein play a significant role in determining the uniqueness of the folded structure.     

Probing phenylalanine environments in oligomeric structures with pentafluorophenylalanine and cyclohexylalanine

Stabilization of protein structures and protein-protein interactions are critical in the engineering of industrially useful enzymes and in the design of pharmaceutically valuable ligands. Hydrophobic interactions involving phenylalanine residues play crucial roles in protein stability and protein-protein/peptide interactions. To establish an effective method to explore the hydrophobic environments of phenylalanine residues, we present a strategy that uses pentafluorophenylalanine (F5Phe) and cyclohexylalanine (Cha). In this study, substitution of F5Phe or Cha for three Phe residues at positions 328, 338, and 341 in the tetramerization domain of the tumor suppressor protein p53 was performed. These residues are located at the interfaces of p53-p53 interactions and are important in the stabilization of the tetrameric structure. The stability of the p53 tetrameric structure did not change significantly when F5Phe-containing peptides at positions Phe328 or Phe338 were used. In contrast, the substitution of Cha for Phe341 in the hydrophobic core enhanced the stability of the tetrameric structure with a T(m) value of 100 degrees C. Phe328 and Phe338 interact with each other through pi-interactions, whereas Phe341 is buried in the surrounding alkyl side-chains of the hydrophobic core of the p53 tetramerization domain. Furthermore, high pressure-assisted denaturation analysis indicated improvement in the occupancy of the hydrophobic core. Considerable stabilization of the p53 tetramer was achieved by filling the identified cavity in the hydrophobic core of the p53 tetramer. The results indicate the status of the Phe residues, indicating that the "pair substitution" of Cha and F5Phe is highly suitable for probing the environments of Phe residues.