Salt effects on hydrophobic‐core formation in folding of a helical miniprotein studied by molecular dynamics simulations

We have investigated effects of salt ions on folding events of a helical miniprotein chicken villin headpiece subdomain HP36. Low concentrations of ions alter electrostatic interactions between charged groups of a protein and can change the populations of conformers. Here, we compare two data sets of folding simulations of HP36 in explicit water solvent with or without ions. For efficient sampling of the conformational space of HP36, the multicanonical replica‐exchange molecular dynamics method was employed. Our analyses suggest that salt alters salt‐bridging nature of the protein at later stages of folding at room temperature. Especially, more nonnative, nonlocal salt bridges are formed at near‐native conformations in pure water. Our analyses also show that such salt‐bridge formation hinders the fully native hydrophobic‐core packing at the final stages of folding. Proteins 2014; 82:933–943. © 2013 Wiley Periodicals, Inc.     

Aromatic and methyl NOEs highlight hydrophobic clustering in the unfolded state of an SH3 domain

The N-terminal SH3 domain of Drosophila drk (drkN SH3 domain) exists in equilibrium between a folded (F(exch)) state and a relatively compact unfolded (U(exch)) state under nondenaturing conditions. Selectively labeled samples of the domain have been analyzed by NOESY NMR experiments to probe residual hydrophobic clustering in the U(exch) state. The labeling strategy included selective protonation of aromatic rings or delta-methyl groups on Ile and Leu residues in a highly deuterated background. Combined with long mixing times, the methods permitted observation of significant numbers of long-range interactions between hydrophobic side chains, providing evidence for multiple conformers involving non-native hydrophobic clusters around the Trp 36 indole. Comparison of these data with previously reported HN-HN NOEs yields structural insight into the diversity of structures within the U(exch) ensemble in the drkN SH3 domain. Many of the HN-HN NOEs are consistent with models containing compact residual nativelike secondary structure and greater exposure of the Trp 36 indole to solvent, similar to kinetic intermediates formed in the hierarchic condensation model of folding. However, the methyl and aromatic NOE data better fit conformations with non-native burial of the Trp indole surrounded by hydrophobic groups and more loosely formed beta-structure; these structural characteristics are more consistent with those of kinetic intermediates formed during the hydrophobic collapse mechanism of folding. This suite of NOE data provides a more complete picture of the structures that span the U(exch) state ensemble, from conformers with non-native structure but long-range contacts to those that are highly nativelike. Together, the results are also consistent with the folding funnel view involving multiple folding pathways for this molecule.     

Backbone-driven collapse in unfolded protein chains

Collapse of unfolded protein chains is an early event in folding. It affects structural properties of intrinsically disordered proteins, which take a considerable fraction of the human proteome. Collapse is generally believed to be driven by hydrophobic forces imposed by the presence of nonpolar amino acid side chains. Contributions from backbone hydrogen bonds to protein folding and stability, however, are controversial. To date, the experimental dissection of side-chain and backbone contributions has not yet been achieved because both types of interactions are integral parts of protein structure. Here, we realized this goal by applying mutagenesis and chemical modification on a set of disordered peptides and proteins. We measured the protein dimensions and kinetics of intra-chain diffusion of modified polypeptides at the level of individual molecules using fluorescence correlation spectroscopy, thereby avoiding artifacts commonly caused by aggregation of unfolded protein material in bulk. We found no contributions from side chains to collapse but, instead, identified backbone interactions as a source sufficient to form globules of native-like dimensions. The presence of backbone hydrogen bonds decreased polypeptide water solubility dramatically and accelerated the nanosecond kinetics of loop closure, in agreement with recent predictions from computer simulation. The presence of side chains, instead, slowed loop closure and modulated the dimensions of intrinsically disordered domains. It appeared that the transient formation of backbone interactions facilitates the diffusive search for productive conformations at the early stage of folding and within intrinsically disordered proteins.     

Side-chain interactions in the folding pathway of a Fyn SH3 domain mutant studied by relaxation dispersion NMR spectroscopy

A major challenge to the study of protein folding is the fact that intermediate states along the reaction pathway are generally unstable and thus difficult to observe. Recently developed NMR relaxation dispersion experiments present an avenue to accessing such states, providing kinetic, thermodynamic, and structural information for intermediates with small (greater than or equal to approximately 1%) populations at equilibrium. We have employed these techniques to study the three-state folding reaction of the G48M Fyn SH3 domain. Using (13)C-, (1)H-, and (15)N-based methods, we have characterized backbone and side-chain interactions in the folded, unfolded, intermediate, and transition states, thereby mapping the energy landscape of the protein. We find that the intermediate, populated to approximately 1%, contains nativelike structure in a central beta-sheet, and is disordered at the amino and carboxy termini. The intermediate is stabilized by side-chain van der Waals contacts, yet (13)C chemical shifts indicate that methyl-containing residues remain disordered. This state has a partially structured backbone and a collapsed yet mobile hydrophobic core and thus closely resembles a molten globule. Nonpolar side-chain contacts are formed in the unfolded-intermediate transition state; these interactions are disrupted in the intermediate-folded transition state, possibly allowing side chains to rearrange as they adopt the native packing configuration. This work illustrates the power of novel relaxation dispersion experiments in characterizing excited states that are "invisible" in even the most sensitive of NMR experiments.     

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.     

On the unyielding hydrophobic core of villin headpiece

Villin headpiece (HP67) is a small, autonomously-folding domain that has become a model system for understanding the fundamental tenets governing protein folding. In this communication, we explore the role that Leu61 plays in the structure and stability of the construct. Deletion of Leu61 results in a completely unfolded protein that cannot be expressed in Escherichia coli. Omission of only the aliphatic leucine side chain (HP67 L61G) perturbed neither the backbone conformation nor the orientation of local hydrophobic side chains. As a result, a large, solvent-exposed hydrophobic pocket, a negative replica of the leucine side-chain, was created on the surface. The loss of the hydrophobic interface between leucine 61 and the hydrophobic pocket destabilized the construct by ~3.3 kcal/mol. Insertion of a single glycine residue immediately before Leu61 (HP67 L61[GL]) was also highly destabilizing and had the effect of altering the backbone conformation (α-helix to π-helix) in order to precisely preserve the wild-type position and conformation of all hydrophobic residues, including Leu61. In addition to demonstrating that the hydrophobic side-chain of Leu61 is critically important for the stability of villin headpiece, our results are consistent with the notion that the precise interactions present within the hydrophobic core, rather than the hydrogen bonds that define the secondary structure, specify a protein's fold.     

Peptide models of folding initiation sites of bovine beta-lactoglobulin: identification of nativelike hydrophobic interactions involving G and H strands

In an attempt to characterize the early folding events in bovine beta-lactoglobulin (BLG), a set of peptides, covering the flexible N-terminal region and the stable C-terminus beta-core, was synthesized and analyzed by circular dichroism and by nuclear magnetic resonance in water, trifluoroethanol (TFE), and sodium dodecyl sulfate (SDS) below and above the critical micellar concentration. The role of local and long-range hydrophobic interactions in guiding the folding has been investigated. For the peptide fragment covering the more flexible N-terminal region of BLG (beta-strands A, B), where both theoretical predictions and kinetic refolding experiments suggested the formation of non-native alpha-helix, no native long-range contacts were identified, and a helical secondary structure was stabilized only in the presence of 25 mM SDS. At variance, in 50% (v/v) TFE, native, long-range hydrophobic interactions were observed in the peptide covering the core region comprising G and H beta-strands. The side chains involved in these interactions form a nativelike hydrophobic cluster, thus suggesting that the GH region may act as the folding initiation site for BLG. This result is reinforced by the identification, in the urea denaturated BLG, of residual structure located at the level of the GH interface, as evidenced by NMR analysis. These results, in excellent agreement with kinetic, thermodynamic, and cold denaturation folding data, once more underline the utmost importance of the GH region for the stability and folding of BLG. Severe aggregation effects prevented the structural analysis of the peptide covering the EFGH region, indicating that this larger segment does not represent an independent folding domain and that the terminal alpha-helix is necessary for stabilizing the BLG folding core.     

A model study of protein nascent chain and cotranslational folding using hydrophobic-polar residues

To study protein nascent chain folding during biosynthesis, we investigate the folding behavior of models of hydrophobic and polar (HP) chains at growing length using both two-dimensional square lattice model and an optimized three-dimensional 4-state discrete off-lattice model. After enumerating all possible sequences and conformations of HP heteropolymers up to length N = 18 and N = 15 in two and three-dimensional space, respectively, we examine changes in adopted structure, stability, and tolerance to single point mutation as the nascent chain grows. In both models, we find that stable model proteins have fewer folded nascent chains during growth, and often will only fold after reaching full length. For the few occasions where partial chains of stable proteins fold, these partial conformations on average are very similar to the corresponding parts of the final conformations at full length. Conversely, we find that sequences with fewer stable nascent chains and sequences with native-like folded nascent chains are more stable. In addition, these stable sequences in general can have many more point mutations and still fold into the same conformation as the wild type sequence. Our results suggest that stable proteins are less likely to be trapped in metastable conformations during biosynthesis, and are more resistant to point-mutations. Our results also imply that less stable proteins will require the assistance of chaperone and other factors during nascent chain folding. Taken together with other reported studies, it seems that cotranslational folding may not be a general mechanism of in vivo protein folding for small proteins, and in vitro folding studies are still relevant for understanding how proteins fold biologically.     

Minimization of cavity size ensures protein stability and folding: structures of Phe46-replaced bovine pancreatic RNase A

The Phe46 residue, located in the hydrophobic core of RNase A, was replaced with other hydrophobic residues, leucine, valine, or alanine, and their X-ray crystallographic structures were determined up to 1.50-1.80 A resolution in an attempt to examine the relationship between structural changes and conformational stability or folding kinetics. The backbone structure of F46L, F46V, and F46A was indistinguishable from that of the wild-type enzyme, retaining the correct active site structure. However, one water molecule was included in the hydrophobic core of F46A, forming two hydrogen bonds with the backbone peptide chain. The side chain of Met29 in F46V and F46A adopted two different conformations in an equal occupancy. A trapped water molecule and two conformations of Met29 represent changes that minimize the cavity volume. Nevertheless, the replacement of Phe46 with the above residues resulted in a marked decrease in both thermal stability and folding reaction. Thus, Phe46 ensures the thermal stability and the rapid and correct folding of RNase A by the role it plays in forming a highly packed, hydrophobic core.     

Identification of formation of initial native structure in onconase from an unfolded state

In the oxidative folding of onconase, the stabilization of intermediates early in the folding process gives rise to efficient formation of its biologically active form. To identify the residues responsible for the initial formation of structured intermediates, the transition from an ensemble of unstructured three-disulfide species, 3S(U), to a single structured three-disulfide intermediate species, des-[30-75] or 3S(F), at pH 8.0 and 25 °C was examined. This transition was first monitored by far-UV circular dichroism spectroscopy at pH 8.0 and 25 °C, showing that it occurs with the formation of secondary structure, presumably because of native interactions. The time dependence of formation of nativelike structure was then followed by nuclear magnetic resonance spectroscopy after we had arrested the transition at different times by lowering the pH to 3 and then acquiring (1)H-(15)N heteronuclear single-quantum coherence spectra at pH 3 and 16 °C to identify amide hydrogens that become part of nativelike structure. H/D exchange was utilized to reduce the intensity of resonances from backbone amide hydrogens not involved in structure, without allowing exchange of backbone amide hydrogens involved in initial structure. Six hydrogen-bonding residues, namely, Tyr38, Lys49, Ser82, Cys90, Glu91, and Ala94, were identified as being involved in the earliest detectable nativelike structure before complete formation of des-[30-75] and are further stabilized later in the formation of this intermediate through S-S/SH interchange. By observing the stabilization of the structures of these residues by their neighboring residues, we have identified the initial, nativelike structural elements formed in this transition, providing details of the initial events in the oxidative folding of onconase.     

Analysis on folding of misgurin using two-dimensional HP model

Misgurin is an antimicrobial peptide from the loach, while the hydrophobic-polar (HP) model is a way to study the folding conformations and native states in peptide and protein although several amino acids cannot be classified either hydrophobic or polar. Practically, the HP model requires extremely intensive computations, thus it has yet to be used widely. In this study, we use the two-dimensional HP model to analyze all possible folding conformations and native states of misgurin with conversion of natural amino acids according to the normalized amino acid hydrophobicity index as well as the shortest benchmark HP sequence. The results show that the conversion of misgurin into HP sequence with glycine as hydrophobic amino acid at pH 2 has 1212 folding conformations with the same native state of minimal energy -6; the conversion of glycine as polar amino acid at pH 2 has 13,386 folding conformations with three native states of minimal energy -5; the conversion of glycine as hydrophobic amino acid at pH 7 has 2538 folding conformations with three native states of minimal energy -5; and the conversion of glycine as polar amino acid at pH 7 has 12,852 folding conformations with three native states of minimal energy -4. Those native states can be ranked according to the normalized amino acid hydrophobicity index. The detailed discussions suggest two ways to modify misgurin.     

Role of side-chains in the cooperative beta-hairpin folding of the short C-terminal fragment derived from streptococcal protein G

A short C-terminal fragment of immunoglobulin-binding domain of streptococcal protein G is known to form nativelike beta-hairpin at physiological conditions. To understand the cooperative folding of the short peptide, eight Ala-substituted mutants of the fragment were investigated with respect to their structural stabilities by analyzing temperature dependence of NMR signals. On comparison of the obtained thermodynamic parameters, we found that the nonpolar residues Tyr45 and Phe52 and the polar residues Asp46 and Thr49 are crucial for the beta-hairpin folding. The results suggest a strong interaction between the nonpolar side chains that participates in a putative hydrophobic cluster and that the polar side chains form a fairly rigid conformation around the loop (46-51). We also investigated the complex formation of the mutants with N-terminal fragment at the variety of temperature to get their thermal unfolding profiles and found that the mutations on the residues Asp46 and Thr49 largely destabilized the complexes, while substitution of Asp47 slightly stabilized the complex. From these results, we deduced that both the hydrophobic cluster formation and the rigidity of the loop (46-51) cooperatively stabilize the beta-hairpin structure of the fragment. These interactions which form a stable beta-hairpin may be the initial structural scaffold which is important in the early folding events of the whole domain.     

Synthesis and characterization of a beta-hairpin peptide that represents a 'core module' of bovine pancreatic trypsin inhibitor (BPTI)

A new strategy for the design and construction of peptide fragments that can achieve defined, nativelike secondary structure is presented. The strategy is based upon the hypothesis that 'core elements' of a protein, synthesized in a single polypeptide molecule, will favor nativelike structure, and that by incorporating a cross-link, nativelike core structure will dominate the ensemble as the more extended conformations are excluded. 'Core elements' are the elements of packed secondary structure that contain the slowest exchanging backbone amide protons in the native protein. The 'core elements' in bovine pancreatic trypsin inhibitor (BPTI) are the two long strands of antiparallel beta-sheet (residues 18-24 and 29-35) and the small beta-bridge (residues 43-44). To test the design strategy, we synthesized an 'oxidized core module', which contains the antiparallel strands connected by a modified reverse turn (A27 replaced by D), a natural disulfide cross-link at the open end of the hairpin, and N- and C-termini blocking groups. A peptide with identical sequence but lacking the disulfide cross-link at the open end was used as the 'reduced core module' control. The conformational behavior of both peptides was examined using (1)H NMR spectroscopy. Chemical shift dispersion, chemical shift deviation from random coil values, sequential and long-range NOEs, and H/D amide exchange rates were compared for the two peptides. We conclude that the ensemble of oxidized and reduced core module conformations samples both nativelike 4:4 and non-native 3:5 beta-hairpin structure, and that the oxidized module samples nativelike structure for a greater fraction of the time than the reduced module.     

Burst-phase expansion of native protein prior to global unfolding in SDS

Although numerous studies have been directed at understanding early folding events through the characterization of folding intermediates, there are few reports on the very late folding events, i.e. on the events taking place on the native side of the folding barrier and on alternative conformations of the folded state. To shed further light on these issues, we have characterized by protein engineering the structure of an expanded but native-like intermediate that accumulates transiently in the unfolding reaction of the small protein S6 in the presence of SDS. The results show that the SDS micelles attack the native protein in the dead-time of the denaturation experiment, causing an expansion of the hydrophobic core prior to the major unfolding transition. We distinguish two forms of the unfolding intermediate that are correlated with the micellar structure. With spherical micelles, the expansion is seen mainly as a weakening of the interactions which anchor the two alpha-helices to the core of the S6 structure. With cylindrical micelles, prevalent at higher SDS concentrations, the expansion is more global and produces a species which closely resembles the transition-state structure for unfolding in GdmCl. Despite the highly weakened core, the micelle-associated intermediate displays cooperative unfolding, indicating a significant structural plasticity of the species on the native side of the folding barrier in the presence of SDS.     

New structural insights into the refolding of ribonuclease T1 as seen by time-resolved Fourier-transform infrared spectroscopy

To get new structural insights into different phases of the renaturation of ribonuclease T1 (RNase T1), the refolding of the thermally unfolded protein was initiated by rapid temperature jumps and detected by time-resolved Fourier-transform infrared spectroscopy. The characteristic spectral changes monitoring the formation of secondary structure and tertiary contacts were followed on a time scale of 10(-3) to 10(3) seconds permitting the characterization of medium and slow folding reactions. Additionally, structural information on the folding events that occurred within the experimental dead time was indirectly accessed by comparative analysis of kinetic and steady-state refolding data. At slightly destabilizing refolding temperatures of 45 degrees C, which is close to the unfolding transition region, no specific secondary or tertiary structure is formed within 180 ms. After this delay all infrared markers bands diagnostic for individual structural elements indicate a strongly cooperative and relatively fast folding, which is not complicated by the accumulation of intermediates. At strongly native folding temperatures of 20 degrees C, a folding species of RNase T1 is detected within the dead time, which already possesses significant amounts of antiparallel beta-sheets, turn structures, and to some degree tertiary contacts. The early formed secondary structure is supposed to comprise the core region of the five-stranded beta-sheet. Despite these nativelike characteristics the subsequent refolding events are strongly heterogeneous and slow. The refolding under strongly native conditions is completed by an extremely slow formation or rearrangement of a locally restricted beta-sheet region accompanied by the further consolidation of turns and denser backbone packing. It is proposed that these late events comprise the final packing of strand 1 (residues 40-42) of the five-stranded beta-sheet against the rest of this beta-sheet system within an otherwise nativelike environment. This conclusion was supported by the comparison of refolding of RNase T1 and its variant W59Y RNase T1 that enabled the assignment of these very late events to the trans-->cis isomerization reaction of the prolyl peptide bond preceding Pro-39.     

Mechanism of protein folding. IV. Forming and breaking of disulfide bonds in bovine pancreatic tripsin inhibitor.

The folding mechanism of bovine pancreatic tripsin inhibitor (BPTI) is explained theoretically on the basis of the island model, where the driving force of folding is hydrophobic interaction. For this purpose, we take a look at the formation and breaking of disulfide bonds during the folding process of BPTI. The intermediate conformations and the native one are successfully obtained, which satisfy the so-called "lampshade" geometrical criterion for the formation of the disulfide bonds. The folding pathway is consistent with the renaturation experiment by Creighton. In addition, an elaborate treatment of side chains of amino acid residues by the software programme CHARMm confirms quantitatively the formation of disulfide bridges.     

Nonrandom structure in the urea-unfolded Escherichia coli outer membrane protein X (OmpX)

On the basis of sequence-specific resonance assignments for the complete polypeptide backbone and most of the amino acid side chains by heteronuclear nuclear magnetic resonance (NMR) spectroscopy, the urea-unfolded form of the outer membrane protein X (OmpX) from Escherichia coli has been structurally characterized. (1)H-(1)H nuclear Overhauser effects (NOEs), dispersion of the chemical shifts, amide proton chemical shift temperature coefficients, amide proton exchange rates, and (15)N[(1)H]-NOEs show that OmpX in 8 M urea at pH 6.5 is globally unfolded, but adopts local nonrandom conformations in the polypeptide segments of residues 73-82 and 137-145. For these two regions, numerous medium-range and longer-range NOEs were observed, which were used as the input for structure calculations of these polypeptide segments with the program DYANA. The segment 73-82 forms a quite regular helical structure, with only loosely constrained amino acid side chains. In the segment 137-145, the tryptophan residue 140 forms the core of a small hydrophobic cluster. Both nonrandom structures are present with an abundance of about 25% of the protein molecules. The sequence-specific NMR assignment and the physicochemical characterization of urea-denatured OmpX presented in this paper are currently used as a platform for investigations of the folding mechanism of this integral membrane protein.     

Crystal structure of a pH-stabilized mutant of villin headpiece

Villin-type headpiece domains are compact F-actin-binding motifs that have been used extensively as a model system to investigate protein folding by both experimental and computational methods. Villin headpiece (HP67) harbors a highly helical, thermostable, and autonomously folding subdomain in the C terminus (HP35), and because of this feature, HP67 is usually considered to be composed of a N- and C-terminal subdomain. Unlike the C-terminal subdomain, the N-terminal subdomain consists mainly of loops and turns, and the folding is dependent upon the presence of the C-terminal subdomain. The pH sensitivity of this subdomain is thought to arise from, at least partially, protonation of H41 buried in the hydrophobic core. Substitution of this histidine with tyrosine, another permissive residue at this position for naturally occurring sequences, increases not only the pH stability of HP67 but also the thermal stability and the cooperativity of thermal unfolding over a wide pH range (0.9-7.5). The crystal structures of wild-type HP67 and the H41Y mutant, determined under the same conditions, indicate that the H41Y substitution causes only localized rearrangement around the mutated residue. The F-actin-binding motif remains essentially the same after the mutation, accounting for the negligible effect of the mutation on F-actin affinity. The hydrogen bond formed between the imidazole ring of H41 and the backbone carbonyl of E14 of HP67 is eliminated by the H41Y mutation, which renders the extreme N terminus of H41Y more mobile; the hydrogen bond formed between the imidazole ring of H41 and the backbone nitrogen of D34 is replaced with that between the hydroxyl group of Y41 and the backbone nitrogen of D34 after the H41Y substitution. The increased hydrophobicity of tyrosine compensates for the loss of hydrogen bonds in the extreme N terminus and accounts for the increased stability and cooperativity of the H41Y mutant.     

Exploring sequence/folding space: folding studies on multiple hydrophobic core mutants of ubiquitin

The stability, dynamic, and structural properties of ubiquitin and two multiple hydrophobic core mutants were studied. One of the mutants (U4) has seven substitutions in the hydrophobic core (M1L, I3L, V5I, I13F, L15V, V17M, and V26L). On average, its side chains are larger than the wild-type, and it can thus be thought of as having an overpacked core. The other mutant (U7) has two substitutions (I3V and I13V). On average, it has smaller side chains than the wild-type, and it can therefore be considered to be underpacked. The three proteins are well-folded and show similar backbone dynamics (T(1), T(2), and HNOE values), indicating that the regular secondary structure extends over the same residue ranges. The crystallographic structure of U4 was determined. The final R(factor) and R(free) are 0.198 and 0.248, respectively, at 2.18 A resolution. The structure of U4 is very similar to wild-type ubiquitin. Remarkably, there are almost no changes in the positions of the C(alpha) atoms along the entire backbone, and the hydrogen-bonding network is maintained. The mutations of the hydrophobic core are accommodated by small movements of side chains in the core of mutated and nonmutated residues. Unfolding and refolding kinetic studies revealed that U4 unfolds with the highest rates; however, its refolding rate constants are very similar to those of the wild-type protein. Conversely, U7 seems to be the most destabilized protein; its refolding rate constant is smaller than the other two proteins. This was confirmed by stopped-flow techniques and by H/D exchange methodologies. This work illustrates the possibility of repacking the hydrophobic core of small proteins and has important implications in the de novo design of stable proteins.     

Energy minimization method using automata network for sequence and side-chain conformation prediction from given backbone geometry

Globular proteins have high packing densities as a result of residue side chains in the core achieving a tight, complementary packing. The internal packing is considered the main determinant of native protein structure. From that point of view, we present here a method of energy minimization using an automata network to predict a set of amino acid sequences and their side-chain conformations from a desired backbone geometry for de novo design of proteins. Using discrete side-chain conformations, that is, rotamers, the sequence generation problem from a given backbone geometry becomes one of combinatorial problems. We focused on the residues composing the interior core region and predicted a set of amino acid Sequences and their side-chain conformations only from a given backbone geometry. The kinds of residues were restricted to six hydrophobic amino acids (Ala, Ile, Met, Leu, Phe, and Val) because the core regions are almost always composed of hydrophobic residues. The obtained sequences were well packed as was the native sequence. The method can be used for automated sequence generation in the de novo design of proteins. © 1994 Wiley-Liss, Inc.