Solution structure of alpha t alpha, a helical hairpin peptide of de novo design

A refined pH-dependent coarse-grained model for peptide structure prediction in aqueous solution

Introduction: Peptides carry out diverse biological functions and the knowledge of the conformational ensemble of polypeptides in various experimental conditions is important for biological applications. All fast dedicated softwares perform well in aqueous solution at neutral pH. Methods: In this study, we go one step beyond by combining the Debye-Hückel formalism for charged-charged amino acid interactions and a coarse-grained potential of the amino acids to treat pH and salt variations. Results: Using the PEP-FOLD framework, we show that our approach performs as well as the machine-leaning AlphaFold2 and TrRosetta methods for 15 well-structured sequences, but shows significant improvement in structure prediction of six poly-charged amino acids and two sequences that have no homologous in the Protein Data Bank, expanding the range of possibilities for the understanding of peptide biological roles and the design of candidate therapeutic peptides.

Effect of Mutations on the Global and Site-Specific Stability and Folding of an Elementary Protein Structural Motif

Understanding the folding mechanism of proteins requires detailed knowledge of the roles of individual amino acid residues in stabilization of specific elements and local segments of the native structure. Recently, we have utilized the combination of circular dichroism (CD) and site-specific ¹³C isotopically edited infrared spectroscopy (IR) coupled with the Ising-like model for protein folding to map the thermal unfolding at the residue level of a de novo designed helix-turn-helix motif αtα. Here we use the same methodology to study how the sequence of local thermal unfolding is affected by selected mutations introduced into the most and least stable parts of the motif. Seven different mutants of αtα are screened to find substitutions with the most pronounced effects on the overall stability. Subsequently, thermal unfolding of two mutated αtα sequences is studied with site-specific resolution, using four distinct ¹³C isotopologues of each. The data are analyzed with the Ising-like model, which builds on a previous parametrization for the original αtα sequence and tests different ways of incorporating the amino acid substitution. We show that for both more and less stable mutants only the adjustment of all interaction parameters of the model can yield a satisfactory fit to the experimental data. The stabilizing and destabilizing mutations result, respectively, in a similar increase and decrease of the stability of all probed local segments, irrespective of their position with respect to the mutation site. Consequently, the relative order of their unfolding remains essentially unchanged. These results underline the importance of the interconnectivity of the stabilizing interaction network and cooperativity of the protein structure, which is evident even in a small motif with apparently noncooperative, heterogeneous unfolding. Overall, our findings are consistent with the native structure being the dominant factor in determining the folding mechanism, regardless of the details of its overall or local thermodynamic stabilization.

The Intracellular Loop of the Na+/Ca2+ Exchanger Contains an "Awareness Ribbon"-Shaped Two-Helix Bundle Domain

The Na⁺/Ca²⁺ exchanger (NCX) is a ubiquitous single-chain membrane protein that plays a major role in regulating the intracellular Ca²⁺ homeostasis by the counter transport of Na⁺ and Ca²⁺ across the cell membrane. Other than its prokaryotic counterpart, which contains only the transmembrane domain and is self-sufficient as an active ion transporter, the eukaryotic NCX protein possesses in addition a large intracellular loop that senses intracellular calcium signals and controls the activation of ion transport across the membrane. This provides a necessary layer of regulation for the more complex function of eukaryotic cells. The Ca²⁺ sensor in the intracellular loop is known as the Ca²⁺-binding domain (CBD12). However, how the signaling of the allosteric intracellular Ca²⁺ binding propagates and results in transmembrane ion transportation still lacks a detailed explanation. Further structural and dynamics characterization of the intracellular loop flanking both sides of CBD12 is therefore imperative. Here, we report the identification and characterization of another structured domain that is N-terminal to CBD12 in the intracellular loop using solution nuclear magnetic resonance (NMR) spectroscopy. The atomistic structure of this domain reveals that two tandem long α-helices, connected by a short linker, form a stable crossover two-helix bundle (THB), resembling an "awareness ribbon". Considering the highly conserved amino acid sequence of the THB domain, the detailed structural and dynamics properties of the THB domain will be common among NCXs from different species and will contribute toward the understanding of the regulatory mechanism of eukaryotic Na⁺/Ca²⁺ exchangers.

Folding Simulations of an α-Helical Hairpin Motif αtα with Residue-Specific Force Fields

α-Helical hairpin (two-helix bundle) is a structure motif composed of two interacting helices connected by a turn or a short loop. It is an important model for protein folding studies, filling the gap between isolated α-helix and larger all-α domains. Here, we present, for the first time, successful folding simulations of an α-helical hairpin. Our RSFF1 and RSFF2 force fields give very similar predicted structures of this αtα peptide, which is in good agreement with its NMR structure. Our simulations also give site-specific stability of α-helix formation in good agreement with amide hydrogen exchange experiments. Combining the folding free energy landscapes and analyses of structures sampled in five different ranges of the fraction of native contacts (Q), a folding mechanism of αtα is proposed. The most stable sites of Q9-E15 in helix-1 and E24-A30 in helix-2 close to the loop region act as the folding initiation sites. The formation of interhelix side-chain contacts also initiates near the loop region, but some residues in the central parts of the two helices also form contacts quite early. The two termini fold at a final stage, and the loop region remains flexible during the whole folding process. This mechanism is similar to the "zipping out" pathway of β-hairpin folding.

Sequence, structure, and cooperativity in folding of elementary protein structural motifs

Residue-level unfolding of two helix-turn-helix proteins--one naturally occurring and one de novo designed--is reconstructed from multiple sets of site-specific (13)C isotopically edited infrared (IR) and circular dichroism (CD) data using Ising-like statistical-mechanical models. Several model variants are parameterized to test the importance of sequence-specific interactions (approximated by Miyazawa-Jernigan statistical potentials), local structural flexibility (derived from the ensemble of NMR structures), interhelical hydrogen bonds, and native contacts separated by intervening disordered regions (through the Wako-Saitô-Muñoz-Eaton scheme, which disallows such configurations). The models are optimized by directly simulating experimental observables: CD ellipticity at 222 nm for model proteins and their fragments and (13)C-amide I' bands for multiple isotopologues of each protein. We find that data can be quantitatively reproduced by the model that allows two interacting segments flanking a disordered loop (double sequence approximation) and incorporates flexibility in the native contact maps, but neither sequence-specific interactions nor hydrogen bonds are required. The near-identical free energy profiles as a function of the global order parameter are consistent with expected similar folding kinetics for nearly identical structures. However, the predicted folding mechanism for the two motifs is different, reflecting the order of local stability. We introduce free energy profiles for "experimental" reaction coordinates--namely, the degree of local folding as sensed by site-specific (13)C-edited IR, which highlight folding heterogeneity and contrast its overall, average description with the detailed, local picture.

Immunisation with foamy virus Bet fusion proteins as novel strategy for HIV-1 epitope delivery

The induction of 2F5- and 4E10-like antibodies broadly neutralising HIV-1 and targeting the membrane external proximal region (MPER) of the transmembrane envelope protein gp41 would be a major advancement for the development of a preventive HIV-1 vaccine, but successful attempts remain rare. Recent studies demonstrated that broadly reactive antibodies develop relatively late during infection and after intensive affinity maturation. Therefore, a prolonged antigen delivery might be beneficial to induce them. Replicating foamy viruses which are characterised by apathogenic but persistent infection could represent suitable carrier viruses for this purpose. In order to develop such a system, we modified the accessory foamy virus Bet protein to contain the MPER of gp41, or the MPER linked to the stabilising fusion peptide proximal region of gp41 and analysed here the antigenic and immunogenic properties of such hybrid proteins. The antigens, expressed and purified to homogeneity, were recognised by the monoclonal antibodies 2F5 and 4E10 with nanomolar affinities and induced high levels of antibodies specific to gp41 after immunisation of rats. The antisera also bound to virus particles attached to infected cells, and peptide-based epitope mapping showed that they recognised the 2F5 epitope. Although no HIV-1 neutralising activity was observed, the presented data demonstrate that using the foamy virus Bet for HIV-1 epitope delivery is successfully applicable. Together with the attractive potential for sustained antigen expression after transfer to replicating virus, these results should therefore provide a first basis for the development of chimeric foamy viruses as novel HIV-1 vaccine vectors.

Fundamental processes of protein folding: measuring the energetic balance between helix formation and hydrophobic interactions

Theories of protein folding often consider contributions from three fundamental elements: loops, hydrophobic interactions, and secondary structures. The pathway of protein folding, the rate of folding, and the final folded structure should be predictable if the energetic contributions to folding of these fundamental factors were properly understood. alphatalpha is a helix-turn-helix peptide that was developed by de novo design to provide a model system for the study of these important elements of protein folding. Hydrogen exchange experiments were performed on selectively 15N-labeled alphatalpha and used to calculate the stability of hydrogen bonds within the peptide. The resulting pattern of hydrogen bond stability was analyzed using a version of Lifson-Roig model that was extended to include a statistical parameter for tertiary interactions. This parameter, x, represents the additional statistical weight conferred upon a helical state by a tertiary contact. The hydrogen exchange data is most closely fit by the XHC model with an x parameter of 9.25. Thus the statistical weight of a hydrophobic tertiary contact is approximately 5.8x the statistical weight for helix formation by alanine. The value for the x parameter derived from this study should provide a basis for the understanding of the relationship between hydrophobic cluster formation and secondary structure formation during the early stages of protein folding.

A multibody, whole-residue potential for protein structures, with testing by Monte Carlo simulated annealing

A new multibody, whole-residue potential for protein tertiary structure is described. The potential is based on the local environment surrounding each main-chain alpha carbon (CA), defined as the set of all residues whose CA coordinates lie within a spherical volume of set radius in 3-dimensional (3D) space surrounding that position. It is shown that the relative positions of the CAs in these local environments belong to a set of preferred templates. The templates are derived by cluster analysis of the presently available database of over 3000 protein chains (750,000 residues) having not more than 30% sequence similarity. For each template is derived also a set of residue propensities for each topological position in the template. Using lookup tables of these derived templates, it is then possible to calculate an energy for any conformation of a given protein sequence. The application of the potential to ab initio protein tertiary structure prediction is evaluated by performing Monte Carlo simulated annealing on test protein sequences.

ARTIST: An activated method in internal coordinate space for sampling protein energy landscapes

We present the first applications of an activated method in internal coordinate space for sampling all-atom protein conformations, the activation-relaxation technique for internal coordinate space trajectories (ARTIST). This method differs from all previous internal coordinate-based studies aimed at folding or refining protein structures in that conformational changes result from identifying and crossing well-defined saddle points connecting energy minima. Our simulations of four model proteins containing between 4 and 47 amino acids indicate that this method is efficient for exploring conformational space in both sparsely and densely packed environments, and offers new perspectives for applications ranging from computer-aided drug design to supramolecular assembly.

Tuning the conformation properties of a peptide by glycosylation and phosphorylation

We have deployed the alpha-helical hairpin peptide (alpha-helix/turn/alpha-helix) and used it as a model system to explore how glycosylation and phosphorylation might affect the conformational properties of the peptide. The native conformations of the modified peptides in buffer solution have been compared with that of the wild-type peptide by nuclear magnetic resonance spectroscopy. Circular dichroism spectroscopy was used to probe the effects of an O-linked beta-GlcNAc and a phosphate group on the overall folding stability of the peptide. Finally, the rate of fibrillogenesis was used to infer the effects of these chemical modifications on the alpha-to-beta transition as well as the rate of nucleation of amyloidogenesis.

De novo backbone and sequence design of an idealized alpha/beta-barrel protein: evidence of stable tertiary structure

We have designed, synthesized, and characterized a 216 amino acid residue sequence encoding a putative idealized alpha/beta-barrel protein. The design was elaborated in two steps. First, the idealized backbone was defined with geometric parameters representing our target fold: a central eight parallel-stranded beta-sheet surrounded by eight parallel alpha-helices, connected together with short structural turns on both sides of the barrel. An automated sequence selection algorithm, based on the dead-end elimination theorem, was used to find the optimal amino acid sequence fitting the target structure. A synthetic gene coding for the designed sequence was constructed and the recombinant artificial protein was expressed in bacteria, purified and characterized. Far-UV CD spectra with prominent bands at 222nm and 208nm revealed the presence of alpha-helix secondary structures (50%) in fairly good agreement with the model. A pronounced absorption band in the near-UV CD region, arising from immobilized aromatic side-chains, showed that the artificial protein is folded in solution. Chemical unfolding monitored by tryptophan fluorescence revealed a conformational stability (DeltaG(H2O)) of 35kJ/mol. Thermal unfolding monitored by near-UV CD revealed a cooperative transition with an apparent T(m) of 65 degrees C. Moreover, the artificial protein did not exhibit any affinity for the hydrophobic fluorescent probe 1-anilinonaphthalene-8-sulfonic acid (ANS), providing additional evidence that the artificial barrel is not in the molten globule state, contrary to previously designed artificial alpha/beta-barrels. Finally, 1H NMR spectra of the folded and unfolded proteins provided evidence for specific interactions in the folded protein. Taken together, the results indicate that the de novo designed alpha/beta-barrel protein adopts a stable three-dimensional structure in solution. These encouraging results show that de novo design of an idealized protein structure of more than 200 amino acid residues is now possible, from construction of a particular backbone conformation to determination of an amino acid sequence with an automated sequence selection algorithm.

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.

Sampling activated mechanisms in proteins with the activation-relaxation technique

The activated dynamics of proteins occur on time scales of milliseconds and longer. Standard all-atom molecular dynamics simulations are limited to much shorter times, of the order of tens of nanoseconds. Therefore, many activated mechanisms that are crucial for long-time dynamics will not be observed in such molecular dynamics simulation; different methods are required. Here, we describe in detail the activation-relaxation technique (ART) that generates directly activated mechanisms. The method is defined in the configurational energy landscape and defines moves in a two step fashion: (a) a configuration is first brought from a local minimum to a nearby first-order saddle point (the activation); and (b) the configuration is relaxed to a new metastable state (the relaxation). The method has already been applied to a wide range of problems in condensed matter, including metallic glasses, amorphous semiconductors and silica glass. We review the algorithm in detail, discuss some previously published results and present simulations of activated mechanisms for a two-helix bundle protein using an all-atom energy function.

Evaluation of a fast implicit solvent model for molecular dynamics simulations

A solvation term based on the solvent accessible surface area (SASA) is combined with the CHARMM polar hydrogen force field for the efficient simulation of peptides and small proteins in aqueous solution. Only two atomic solvation parameters are used: one is negative for favoring the direct solvation of polar groups and the other positive for taking into account the hydrophobic effect on apolar groups. To approximate the water screening effects on the intrasolute electrostatic interactions, a distance-dependent dielectric function is used and ionic side chains are neutralized. The use of an analytical approximation of the SASA renders the model extremely efficient (i.e., only about 50% slower than in vacuo simulations). The limitations and range of applicability of the SASA model are assessed by simulations of proteins and structured peptides. For the latter, the present study and results reported elsewhere show that with the SASA model it is possible to sample a significant amount of folding/unfolding transitions, which permit the study of the thermodynamics and kinetics of folding at an atomic level of detail.

Ab initio protein structure prediction using physicochemical potentials and a simplified off-lattice model

This study describes a computational method for ab inito protein structure prediction. Protein conformation has been modeled by using six optimized backbone torsion angles and fixed side chains approximating rotationally averaged real side chains. The approximations aim to keep complexity of the structure description to a minimum without seriously compromising the accuracy of the structural representation. An evolutionary Monte Carlo algorithm has been developed to search through this restricted conformational space to locate low-energy protein structures. A simple physicochemical force field has been developed to assess the energies of different conformations within this structural description. The corresponding residue interaction energies are based on hydrophobic, hydrophilic, steric, and hydrogen-bonding potentials. The search procedure has been used to locate native energy minima from primary sequence alone. The 3-D structures of polypeptides up to 38 residues with both beta and alpha secondary structural elements have been accurately predicted. The search procedure has been found to be highly efficient and follows an energetically and structurally plausible pathway to locate native populations. The simple force field described in the study has been compared with a more complex all-atom model and been found to be similarly effective in predicting the structures of proposed independent folding units. Proteins 2001;43:186-202.

De novo design of helical bundles as models for understanding protein folding and function

De novo protein design has proven to be a powerful tool for understanding protein folding, structure, and function. In this Account, we highlight aspects of our research on the design of dimeric, four-helix bundles. Dimeric, four-helix bundles are found throughout nature, and the history of their design in our laboratory illustrates our hierarchic approach to protein design. This approach has been successfully applied to create a completely native-like protein. Structural and mutational analysis allowed us to explore the determinants of native protein structure. These determinants were then applied to the design of a dinuclear metal-binding protein that can now serve as a model for this important class of proteins.

De novo design and structural characterization of proteins and metalloproteins

De novo protein design has recently emerged as an attractive approach for studying the structure and function of proteins. This approach critically tests our understanding of the principles of protein folding; only in de novo design must one truly confront the issue of how to specify a protein's fold and function. If we truly understand proteins, it should be possible to design receptors, enzymes, and ion channels from scratch. Further, as this understanding evolves and is further refined, it should be possible to design proteins and biomimetic polymers with properties unprecedented in nature.

Retrostructural analysis of metalloproteins: application to the design of a minimal model for diiron proteins

De novo protein design provides an attractive approach for the construction of models to probe the features required for function of complex metalloproteins. The metal-binding sites of many metalloproteins lie between multiple elements of secondary structure, inviting a retrostructural approach to constructing minimal models of their active sites. The backbone geometries comprising the metal-binding sites of zinc fingers, diiron proteins, and rubredoxins may be described to within approximately 1 A rms deviation by using a simple geometric model with only six adjustable parameters. These geometric models provide excellent starting points for the design of metalloproteins, as illustrated in the construction of Due Ferro 1 (DF1), a minimal model for the Glu-Xxx-Xxx-His class of dinuclear metalloproteins. This protein was synthesized and structurally characterized as the di-Zn(II) complex by x-ray crystallography, by using data that extend to 2.5 A. This four-helix bundle protein is comprised of two noncovalently associated helix-loop-helix motifs. The dinuclear center is formed by two bridging Glu and two chelating Glu side chains, as well as two monodentate His ligands. The primary ligands are mostly buried in the protein interior, and their geometries are stabilized by a network of hydrogen bonds to second-shell ligands. In particular, a Tyr residue forms a hydrogen bond to a chelating Glu ligand, similar to a motif found in the diiron-containing R2 subunit of Escherichia coli ribonucleotide reductase and the ferritins. DF1 also binds cobalt and iron ions and should provide an attractive model for a variety of diiron proteins that use oxygen for processes including iron storage, radical formation, and hydrocarbon oxidation.

Structures of the intradiskal loops and amino terminus of the G-protein receptor, rhodopsin

The intradiskal surface of the transmembrane protein, rhodopsin, consists of the amino terminal domain and three loops connecting six of the seven transmembrane helices. This surface corresponds to the extracellular surface of other G-protein receptors. Peptides that represent each of the extramembraneous domains on this surface (three loops and the amino terminus) were synthesized. These peptides also included residues which, based on a hydrophobic plot, could be expected to be part of the transmembrane helix. The structure of each of these peptides in solution was then determined using two-dimensional 1H nuclear magnetic resonance. All peptide domains showed ordered structures in solution. The structures of each of the peptides from intradiskal loops of rhodopsin exhibited a turn in the central region of the peptide. The ends of the peptides show an unwinding of the transmembrane helices to form this turn. The amino terminal domain peptide exhibited alpha-helical regions with breaks and bends at proline residues. This region forms a compact domain. Together, the structures for the loop and amino terminus domains indicate that the intradiskal surface of rhodopsin is ordered. These data further suggest a structural motif for short loops in transmembrane proteins. The ordered structures of these loops, in the absence of the transmembrane helices, indicate that the primary sequences of these loops are sufficient to code for the turn.

A polar, solvent-exposed residue can be essential for native protein structure

BACKGROUND

A large energy gap between the native state and the non-native folded states is required for folding into a unique three-dimensional structure. The features that define this energy gap are not well understood, but can be addressed using de novo protein design. Previously, alpha(2)D, a dimeric four-helix bundle, was designed and shown to adopt a native-like conformation. The high-resolution solution structure revealed that this protein adopted a bisecting U motif. Glu7, a solvent-exposed residue that adopts many conformations in solution, might be involved in defining the unique three-dimensional structure of alpha(2)D.

RESULTS

A variety of hydrophobic and polar residues were substituted for Glu7 and the dynamic and thermodynamic properties of the resulting proteins were characterized by analytical ultracentrifugation, circular dichroism spectroscopy, and nuclear magnetic resonance spectroscopy. The majority of substitutions at this solvent-exposed position had little affect on the ability to fold into a dimeric four-helix bundle. The ability to adopt a unique conformation, however, was profoundly modulated by the residue at this position despite the similar free energies of folding of each variant.

CONCLUSIONS

Although Glu7 is not involved directly in stabilizing the native state of alpha(2)D, it is involved indirectly in specifying the observed fold by modulating the energy gap between the native state and the non-native folded states. These results provide experimental support for hypothetical models arising from lattice simulations of protein folding, and underscore the importance of polar interfacial residues in defining the native conformations of proteins.

Amyloid beta-protein fibrillogenesis. Structure and biological activity of protofibrillar intermediates

Alzheimer's disease is characterized by extensive cerebral amyloid deposition. Amyloid deposits associated with damaged neuropil and blood vessels contain abundant fibrils formed by the amyloid beta-protein (Abeta). Fibrils, both in vitro and in vivo, are neurotoxic. For this reason, substantial effort has been expended to develop therapeutic approaches to control Abeta production and amyloidogenesis. Achievement of the latter goal is facilitated by a rigorous mechanistic understanding of the fibrillogenesis process. Recently, we discovered a novel intermediate in the pathway of Abeta fibril formation, the amyloid protofibril (Walsh, D. M., Lomakin, A., Benedek, G. B., Condron, M. M., and Teplow, D. B. (1997) J. Biol. Chem. 272, 22364-22372). We report here results of studies of the assembly, structure, and biological activity of these polymers. We find that protofibrils: 1) are in equilibrium with low molecular weight Abeta (monomeric or dimeric); 2) have a secondary structure characteristic of amyloid fibrils; 3) appear as beaded chains in rotary shadowed preparations examined electron microscopically; 4) give rise to mature amyloid-like fibrils; and 5) affect the normal metabolism of cultured neurons. The implications of these results for the development of therapies for Alzheimer's disease and for our understanding of fibril assembly are discussed.

Uniquely folded mini-protein motifs

Mini-proteins containing fewer than 40 amino acids provide simple model systems for studying protein folding and stability as well as serving as scaffolds for the rational design of new functional motifs. This article reviews current progress on the design and characterization of discretely folded mini-protein motifs.

Solution structure and dynamics of a de novo designed three-helix bundle protein

Although de novo protein design is an important endeavor with implications for understanding protein folding, until now, structures have been determined for only a few 25- to 30-residue designed miniproteins. Here, the NMR solution structure of a complex 73-residue three-helix bundle protein, alpha3D, is reported. The structure of alpha3D was not based on any natural protein, and yet it shows thermodynamic and spectroscopic properties typical of native proteins. A variety of features contribute to its unique structure, including electrostatics, the packing of a diverse set of hydrophobic side chains, and a loop that incorporates common capping motifs. Thus, it is now possible to design a complex protein with a well defined and predictable three-dimensional structure.

From coiled coils to small globular proteins: design of a native-like three-helix bundle

A monomolecular native-like three-helix bundle has been designed in an iterative process, beginning with a peptide that noncooperatively assembled into an antiparallel three-helix bundle. Three versions of the protein were designed in which specific interactions were incrementally added. The hydrodynamic and spectroscopic properties of the proteins were examined by size exclusion chromatography, sedimentation equilibrium, fluorescence spectroscopy, and NMR. The thermodynamics of folding were evaluated by monitoring the thermal and guanidine-induced unfolding transitions using far UV circular dichroism spectroscopy. The attainment of a unique, native-like state was achieved through the introduction of: (1) helix capping interactions; (2) electrostatic interactions between partially exposed charged residues; (3) a diverse collection of apolar side chains within the hydrophobic core.