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

Exploration of structural features of monomeric helical peptides designed with a genetic algorithm

A genetic algorithm (GA)-based strategy to dissect the determinants of peptide folding into alpha-helix was developed. The structural information of helical peptides was obtained with respect to patterns of sequence variability. In many previously reported studies the intrinsic alpha-helical propensities of amino acids although sequence-dependent are apparently independent of the amino acid position. In this research, monomeric helical peptides selected from possible sequences produced by a GA-chemical synthesis were analyzed to identify possible influential structural features. These hexadeca-peptides were obtained after four successive generations. A total of 128 synthetic peptides were evaluated via circular dichroism (CD) measurements in aqueous solution, while the mean ellipticity at 222 nm confirmed the monomeric state of the peptides. The results presented here show that our GA-based strategy may be useful in the design of proteins with increased alpha-helix content.

Mimicking Photosynthesis in a Computationally Designed Synthetic Metalloprotein

While advances in protein design have made possible the construction of protein architectures with nativelike properties and predictable structures and function, there are as of yet no examples of functional, protein-based, solar energy conversion systems. This communication describes the design and characterization of an artificial reaction center (RC) protein that closely resembles the function of the natural photosynthetic RC. The synthetic protein, designed by the protein design program CORE, participates in multiple reduction/oxidation cycles with exogenous acceptors/donors following photoexcitation. The designed metalloprotein, aRC, consists of a tetrahelical bundle functionalized with two bis-histidine bound metal cofactors: a Ru(bpy)2 moiety and a heme group. Two distinct bis-histidine binding sites were engineered for each of these metal centers. Photoexcitation of aRC results in rapid ET from the RuII complex to the heme group (kET >/= 5 x 1010 s-1) yielding a long-lived (70 ns) charge-separated state (CSS), RuIII/FeII. This long-lived CSS participates in subsequent ET reactions with exogenous donors and acceptors in multiple photocycles, thus mimicking the basic function of native photosynthetic RCs. This study illustrates the successful design and construction of a protein-based functional charge separation device using a combination of automated computational protein design and knowledge of the engineering principles of biological electron tunneling extracted from natural electron-transfer systems. To our knowledge, this represents the first example of a functional protein-based artificial reaction center.

Self-assembling of the porphyrin-linked acyclic penta- and heptapeptides in aqueous trifluoroethanol

The conjugates of porphyrin with links to the acyclic penta- and heptapeptides were synthesized to mimic natural multiple porphyrin systems. The linear penta- and heptapeptide with hydrophilic/hydrophobic alternative sequences took a random structure in aqueous trifluoroethanol (TFE). However, these polypeptides took a beta-sheet structure in the same solvent when the N-terminal Cys linked to the porphyrin, suggesting that the conjugates self-assembled via the intermolecular hydrophobic interaction between the porphyrins. The circular dichroism (CD) spectra, UV/vis spectra, size exclusion chromatography (SEC), and (1)H NMR spectroscopy supported the self-assembling. In the self-assembled structure of the pentapeptide linking porphyrin at the p-phenyl position (9), the porphyrins were involved in two porphyrin-porphyrin interactions, i.e., the side-by-side interaction between the neighboring polypeptide chains and the face-to-face interaction between the first and the third peptide chains. The CD spectra of 9 showed two sets of Cotton effects probably arising from these two interactions. The UV/vis spectra also supported the above interpretation, showing multiple absorptions in the longwave and shortwave shifted regions. The SEC analyses showed the assembled structure of the conjugates. The (1)H NMR signals of the porphyrin rings of 9 were hardly observed in D(2)O-CD(3)OD because of the shortened spin-spin relaxation time T(2)().

Monosaccharide templates for de novo designed 4-alpha-helix bundle proteins: template effects in carboproteins

De novo design and total chemical synthesis of proteins provide powerful approaches to critically test our understanding of protein folding, structure, and stability. The 4-alpha-helix bundle is a frequently studied structure in which four amphiphilic alpha-helical peptide strands form a hydrophobic core. Assembly of protein models on a template has been suggested as a way to reduce the entropy of folding. We have previously developed the concept of carbohydrates as templates in the de novo design of protein models termed 'carboproteins'. Here we present the chemical synthesis of three 8.1 kDa 4-alpha-helix bundles by oxime ligation of tetra-aminooxyacetyl functionalized D-galacto-, D-gluco-, and D-altropyranoside templates with an amphiphilic C-terminal hexadecapeptide aldehyde sequence. CD spectroscopy indicated that the choice of template has an effect on the overall structure of the carboprotein, as the altro-based carboprotein was found to be more alpha-helical than the corresponding galacto- and gluco-carboproteins. However, an influence on stability could not be detected in the present experiments, as the three carboproteins gave similar free energy of foldings (deltaG(F)H2O) and melting points in chemical and thermal denaturation experiments.

De novo design of foldable proteins with smooth folding funnel: automated negative design and experimental verification

De novo sequence design of foldable proteins provides a way of investigating principles of protein architecture. We performed fully automated sequence design for a target structure having a three-helix bundle topology and synthesized the designed sequences. Our design principle is different from the conventional approach, in that instead of optimizing interactions within the target structure, we design the global shape of the protein folding funnel. This includes automated implementation of negative design by explicitly requiring higher free energy of the denatured state. The designed sequences do not have significant similarity to those of any natural proteins. The NMR and CD spectroscopic data indicated that one designed sequence has a well-defined three-dimensional structure as well as alpha-helical content consistent with the target.

Noncovalent binding of a reaction intermediate by a designed helix-loop-helix motif-implications for catalyst design

In our search for a catalyst for the transamination reaction of aspartic acid to form oxaloacetate, twenty-five forty-two-residue sequences were designed to fold into helix-loop-helix dimers and form binding sites for the key intermediate along the reaction pathway, the aldimine. This intermediate is formed from aspartic acid and the cofactor pyridoxal phosphate. The design of the binding sites followed a strategy in which exclusively noncovalent forces were used for binding the aldimine. Histidine residues were incorporated to catalyse the rate-limiting 1,3 proton transfer reaction that converts the aldimine into the ketimine, an intermediate that is subsequently hydrolysed to form oxaloacetate and pyridoxamine phosphate. The two most efficient catalysts, T-4 and T-16, selected from the pool of sequences by a simple screening procedure, were shown by CD and NMR spectroscopies to bind the aldimine intermediate with dissociation constants in the millimolar range. The mean residue ellipticity of T-4 in aqueous solution at pH 7.4 and a concentration of 0.75 mM was -18500 deg x cm(2) dmol(-1). Upon addition of 6 mm l-aspartic acid and 1.5 mM pyridoxal phosphate to form the aldimine, the mean residue ellipticity changed to -19900 deg x cm(2) dmol(-1). The corresponding mean residue ellipticities of T-16 were -21200 deg x cm(2) dmol(-1) and -24000 deg x cm(2) dmol(-1). These results show that the helical content increased in the presence of the aldimine, and that the folded polypeptides bound the aldimine. The (1)H NMR relaxation time of the imine CH proton of the aldimine was affected by the presence of T-4 as was the (31)P NMR resonance linewidth. The catalytic efficiencies of T-4 and T-16 were compared to that of imidazole and found to be more than three orders of magnitude larger. The designed binding sites were thus shown to be capable of binding the aldimine in close proximity to His residues, by noncovalent forces, into conformations that proved to be catalytically active. The results show for the first time the design of well-defined catalytic sites that bind a reaction intermediate with enzyme-like affinities under equilibrium conditions and represent an important advance in de novo catalyst design.

Lanthanide-binding helix-turn-helix peptides: solution structure of a designed metallonuclease

A designed lanthanide-binding chimeric peptide based on the strikingly similar geometries of the EF-hand and helix-turn-helix (HTH) motifs was investigated by NMR and CD spectroscopy and found to retain the same overall solution structure of the parental motifs. CD spectroscopy showed that the 33-mer peptide P3W folds on binding lanthanides, with an increase in alpha-helicity from 20% in the absence of metal to 38% and 35% in the presence of excess Eu(III) and La(III) ions, respectively. The conditional binding affinities of P3W for La(III) (5.9 +/- 0.3 microM) and for Eu(III) (6.2 +/- 0.3 microM) (pH 7.8, 5 mM Tris) were determined by tryptophan fluorescence titration. The La(III) complex of peptide P3, which differs from P3W by only one Trp-to-His substitution, has much less signal dispersion in the proton NMR spectra than LaP3W, indicating that the Trp residue is a critical hydrophobic anchor for maintaining a well-folded helix-turn-helix structure. A chemical-shift index analysis indicates the metallopeptide has a helix-loop-helix secondary structure. A structure calculated by using nuclear Overhauser effect and other NMR constraints reveals that P3W not only has a tightly folded metal-binding loop but also retains the alpha-alpha corner supersecondary structure of the parental motifs. Although the solution structure is undefined at both the N and C termini, the NMR structure confirms the successful incorporation of a metal-binding loop into a HTH sequence.

Rapid cooperative two-state folding of a miniature alpha-beta protein and design of a thermostable variant

There is a great deal of interest in developing small stably folded miniature proteins. A limited number of these molecules have been described, however they typically have not been characterized in depth. In particular, almost no detailed studies of the thermodynamics and folding kinetics of these proteins have been reported. Here we describe detailed studies of the thermodynamics and kinetics of folding of a 39 residue mixed alpha-beta protein (NTL9(1-39)) derived from the N-terminal domain of the ribosomal protein L9. The protein folds cooperatively and rapidly in a two-state fashion to a native state typical of those found for normal globular proteins. At pH 5.4 in 20mM sodium acetate, 100mM NaCl the temperature of maximum stability is 6 degrees C, the t(m) is 65.3 degrees C, deltaH degrees (t(m)) is between 24.6 kcalmol(-1) and 26.3 kcalmol(-1), and deltaC(p) degrees is 0.38 kcalmol(-1)deg(-1). The thermodynamic parameters are in the range expected on the basis of per residue values determined from databases of globular proteins. H/2H exchange measurements reveal a set of amides that exchange via global unfolding, exactly as expected for a normal cooperatively folded globular protein. Kinetic measurements show that folding is two-state folding. The folding rate is 640 s(-1) and the value of deltaG degrees calculated from the folding and unfolding rates is in excellent agreement with the equilibrium value. A designed thermostable variant, generated by mutating K12 to M, was characterized and found to have a t(m) of 82 degrees C. Equilibrium and kinetic measurements demonstrate that its folding is cooperative and two-state.

Analysis of the binding energies of testosterone, 5alpha-dihydrotestosterone, androstenedione and dehydroepiandrosterone sulfate with an antitestosterone antibody

Molecular dynamics simulations and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) free energy calculations were used to study the binding of testosterone (TES), 5alpha-dihydrotestosterone (5ADHT), androstenedione (AND), and dehydroepiandrosterone sulfate (DHEAS) to the monoclonal antitestosterone antibody 3-C(4)F(5). The relative binding free energy of TES and AND was also calculated with free energy perturbation (FEP) simulations. The antibody 3-C(4)F(5) has a relatively high affinity (3 x 10(8) M(-1)) and on overall good binding profile for testosterone but its cross-reactivity with DHEAS has been the main reason for the failure to use this antibody in clinical immunoassays. The relative binding free energies obtained with the MM-PBSA method were 1.5 kcal/mol for 5ADHT, 3.8 kcal/mol for AND, and 4.3 kcal/mol for DHEAS, as compared to TES. When a water molecule of the ligand binding site, observed in the antibody-TES crystal structure, was explicitly included in MM-PBSA calculations, the relative binding energies were 3.4, 4.9, and 5.4 kcal/mol for 5ADHT, AND, and DHEAS, respectively. The calculated numbers are in correct order but larger than the corresponding experimental energies of 1.3, 1.5, and 2.6 kcal/mol, respectively. The fact that the MM-PBSA method reproduced the relative binding free energies of DHEAS, a steroid having a negatively charged sulfate group, and the neutrally charged TES, 5ADHT, and AND in satisfactory agreement with experiment shows the robustness of the method in predicting relative binding affinities. The 800-ps FEP simulations predicted that the antibody 3-C(4)F(5) binds TES 1.3 kcal/mol tighter than AND. Computational mutagenesis of selected amino acid residues of the ligand binding site revealed that the lower affinities of AND and DHEAS as compared to TES are due to a combined effect of several residues, each contributing a small fraction to the tighter binding of TES. An exception to this is Tyr99H, whose mutation to Ala lowered the binding of DHEAS 0.7 kcal/mol more than the binding of TES. This is probably due to the hydrogen bonding interaction formed between the OH group of Tyr99H and the sulfate group of DHEAS. Computational mutagensis data also showed that the affinity of the steroids to the antitestosterone antibody 3-C(4)F(5) would be enhanced if Trp47H were repositioned so that it would make more extensive contacts with the bound ligands. In addition, the binding of steroids to antitestosterone, antiprogesterone, and antiestradiol antibodies is discussed.

Hydrogen exchange, core modules, and new designed proteins

A strategy for design of new proteins that mimic folding properties of native proteins is based on peptides modeled on the slow exchange cores of natural proteins. We have synthesized peptides, called core modules, that correspond to the elements of secondary structure that carry the very slowest exchanging amides in a protein. The expectation is that, if soluble in water, core modules will form conformational ensembles that favor native-like structure. Core modules modeled on natural bovine pancreatic trypsin inhibitor have been shown by NMR studies to meet this expectation. The next step toward production of a native state mimic is to further shift the conformational bias of a core module toward more ordered structure by promoting module-module interactions that are mutually stabilizing. For this, two core modules were incorporated into a single molecule by means of a long cross-link. From a panel of several two-module peptides, one very promising lead emerged; it is called BetaCore. BetaCore is monomeric in water and forms a new fold composed of a four-stranded, antiparallel beta-sheet. The single, dominant conformation of BetaCore is characterized by various NMR experiments. Here we compare the individual core module to the two-module BetaCore and discuss the progressive stabilization of intramodule structure and the formation of new intermodule interactions.

Computational de novo design, and characterization of an A(2)B(2) diiron protein

Diiron proteins are found throughout nature and have a diverse range of functions; proteins in this class include methane monooxygenase, ribonucleotide reductase, Delta(9)-acyl carrier protein desaturase, rubrerythrin, hemerythrin, and the ferritins. Although each of these proteins has a very different overall fold, in every case the diiron active site is situated within a four-helix bundle. Additionally, nearly all of these proteins have a conserved Glu-Xxx-Xxx-His motif on two of the four helices with the Glu and His residues ligating the iron atoms. Intriguingly, subtle differences in the active site can result in a wide variety of functions. To probe the structural basis for this diversity, we designed an A(2)B(2) heterotetrameric four-helix bundle with an active site similar to those found in the naturally occurring diiron proteins. A novel computational approach was developed for the design, which considers the energy of not only the desired fold but also alternatively folded structures. Circular dichroism spectroscopy, analytical ultracentrifugation, and thermal unfolding studies indicate that the A and B peptides specifically associate to form an A(2)B(2) heterotetramer. Further, the protein binds Zn(II) and Co(II) in the expected manner and shows ferroxidase activity under single turnover conditions.

Role of stabilization centers in 4 helix bundle proteins

Stabilization centers (SCs) were shown to play an important role in preventing decay of three-dimensional protein structures. These residue clusters, stabilized by cooperative long-range interactions, were proposed to serve as anchoring points for arranging secondary structure elements. In all-alpha proteins, SC elements appear less frequently than in all-beta, alpha/beta, and alpha+beta proteins suggesting that tertiary structure formation of all-alpha proteins is governed by different principles than in other protein classes. Here we analyzed the relation between the formation of stabilization centers and the inter-axial angles (Omega) of alpha-helices in 4 helix bundle proteins. In the distance range, where dipoles have dominant effect on the helix pair arrangement, those helix pairs, where residues from both helices participate in SC elements, appear as parallel more frequently than those helices where no SC elements are present. For SC containing helix pairs, the energetic difference between the parallel and anti-parallel states decreases considerably from 1.1 kcal/mol to 0.4 kcal/mol. Although the observed effect is weak for more distant helices, a competition between the SC element formation and the optimal dipole-dipole interaction of alpha-helices is proposed as a mechanism for tertiary structure formation in 4 helix bundle proteins. The SC-forming potential of different arrangements as well as the pitfalls of the SC definition are also discussed.

BetaCore, a designed water soluble four-stranded antiparallel beta-sheet protein

BetaCore is a designed approximately 50-residue protein in which two BPTI-derived core modules, CM I and CM II, are connected by a 22-atom cross-link. At low temperature and pH 3, homo- and heteronuclear NMR data report a dominant folded ('f') conformation with well-dispersed chemical shifts, i, i+1 periodicity, numerous long-range NOEs, and slowed amide hydrogen isotope exchange patterns that is a four-stranded antiparallel beta-sheet with nonsymmetrical and specific association of CM I and CM II. BetaCore 'f' conformations undergo reversible, global, moderately cooperative, non-two-state thermal transitions to an equilibrium ensemble of unfolded 'u' conformations. There is a significant energy barrier between 'f' and 'u' conformations. This is the first designed four-stranded antiparallel beta-sheet that folds in water.

A hierarchic approach to the design of hexameric helical barrels

The design of large macromolecular assemblies is an endeavor with implications for protein engineering as well as nanotechnology. A hierarchic approach was used to design an antiparallel hexameric, tubular assembly of helices. In previous studies, a domain-swapped, dimeric three-helix bundle was designed from first principles. In the crystal lattice, three dimers associate around a 3-fold rotational axis to form a hexameric assembly. Although this hexameric assembly was not observed in solution, it was possible to stabilize its formation by changing three polar residues per monomer to hydrophobic (two Phe and one Trp) residues. Molecular models based on the crystallographic coordinates of DSD (PDB accession code 1G6U) show that these side-chains pack in the central cavity (the "supercore") of the hexameric bundle. Analytical ultracentrifugation, fluorescence spectroscopy, CD spectroscopy, and guanidine-HCl denaturation were used to determine the assembly of the hexamer. To probe the requirements for stabilizing the hexamer, we systematically varied the polarity and steric bulk of one of the Phe residues in the supercore of the hexamer. Depending on the nature of this side-chain, it is possible to modulate the stability of the hexamer in a predictable manner. This family of hexameric proteins may provide a useful framework for the construction of proteins that change their oligomeric states in response to binding of small molecules.

Toward the synthesis of artificial proteins: the discovery of an amphiphilic helical peptoid assembly

While nature exploits folded biopolymers to achieve molecular recognition and catalysis, comparable abiological heteropolymer systems have been difficult to create. We synthesized and identified abiological peptoid heteroploymers capable of binding a dye. Using combinatorial synthesis, we constructed a library of 3400 amphiphilic 15-mer peptoids on an ultra-high-capacity beaded support. Individual macrobeads, each containing a single peptoid sequence, were arrayed into plates, cleaved, and screened in aqueous solution to locate dye binding heteropolymer assemblies. Resynthesis and characterization demonstrated the formation of defined helical assemblies as judged by size-exclusion chromatography, circular dichroism, and analytical ultracentrifugation. Inspired by nature's process of sequence variation and natural selection, we identified rare abiological sequence-specific heteropolymers that begin to mimic the structure and functional properties of their biological counterparts.

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.

The metal site as a template for the metalloprotein structure formation

Achieving a thorough explanation of the behavior of metal sites in the formation of native metalloprotein structures is an exciting challenge in the biochemistry of metallobiomacromolecules. This study presents a personal insight into the subject. It is proposed that a metal center and its exogenous ligand compose a template. A template may impose a clear stereochemical preference on the loose peptide chains, and organize them into natural stereospecificity via the metal-ligand interaction, a long-range and strong interaction. Therefore, the stable peptide conformation induced by the template effect surrounding a template polyhedron could be called a template-mediated structural motif (TMSM).

Computational protein design

The field of computational protein design is reaching its adolescence. Protein design algorithms have been applied to design or engineer proteins that fold, fold faster, catalyze, catalyze faster, signal, and adopt preferred conformational states. Further developments of scoring functions, sampling strategies, and optimization methods will expand the range of applicability of computational protein design to larger and more varied systems, with greater incidence of success. Developments in this field are beginning to have significant impact on biotechnology and chemical biology.

Toward beta-peptide tertiary structure: self-association of an amphiphilic 14-helix in aqueous solution

A major frontier in foldamer research is creation of unnatural oligomers that adopt discrete tertiary structures; at present, only biopolymers are known to fold into such compact conformations. We report an initial step toward helix-bundle tertiary structure in the beta-peptide realm by showing that a 10-residue beta-peptide designed to adopt an amphiphilic helical conformation forms small soluble aggregates in water. Sedimentation equilibrium data indicate that the aggregated state falls in the tetramer-hexamer size range. [structure: see text]

Energetic analysis of binding of progesterone and 5 beta-androstane-3,17-dione to anti-progesterone antibody DB3 using molecular dynamics and free energy calculations

Molecular dynamics simulations and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) free energy calculations were used to study the energetics of the binding of progesterone (PRG) and 5 beta-androstane-3,17-dione (5AD) to anti-PRG antibody DB3. Although the two steroids bind to DB3 in different orientations, their binding affinities are of the same magnitude, 1 nM for PRG and 8 nM for 5AD. The calculated relative binding free energy of the steroids, 8.8 kJ/mol, is in fair agreement with the experimental energy, 5.4 kJ/mol. In addition, computational alanine scanning was applied to study the role of selected amino acid residues of the ligand-binding site on the steroid cross-reactivity. The electrostatic and van der Waals components of the total binding free energies were found to favour more the binding of PRG, whereas solvation energies were more favourable for the binding of 5AD. The differences in the free energy components are due to the binding of the A rings of the steroids to different binding pockets: PRG is bound to a pocket in which electrostatic antibody-steroid interactions are dominating, whereas 5AD is bound to a pocket in which van der Waals and hydrophobic interactions dominate.

Searching sequence space for protein catalysts

Genetic selection was used to explore the probability of finding enzymes in protein sequence space. Large degenerate libraries were prepared by replacing all secondary structure units in a dimeric, helical bundle chorismate mutase with simple binary-patterned modules based on a limited set of four polar and four nonpolar residues. Two-stage in vivo selection yielded catalytically active variants possessing biophysical and kinetic properties typical of the natural enzyme even though approximately 80% of the protein originates from the simplified modules and >90% of the protein consists of only eight different amino acids. This study provides a quantitative assessment of the number of sequences compatible with a given fold and implicates previously unidentified residues needed to form a functional active site. Given the extremely low incidence of enzymes in completely unbiased libraries, strategies that combine chemical information with genetic selection, like the one used here, may be generally useful in designing novel protein scaffolds with tailored activities.

Emerging principles of de novo catalyst design

Considerable progress has been made in the understanding of how to exploit hydrophobic and charge-charge interactions in forming binding sites for peptides and small molecules in folded polypeptide catalysts. This knowledge has enabled the introduction of feedback and control functions into catalytic cycles and the construction of folded polypeptide catalysts that follow saturation kinetics. Major advances have also been made in the design of metalloproteins and metallopeptides, especially with regards to understanding redox potential control.

Total synthesis of cytochrome b562 by native chemical ligation using a removable auxiliary

We have completed the total chemical synthesis of cytochrome b562 and an axial ligand analogue, [SeMet(7)]cyt b562, by thioester-mediated chemical ligation of unprotected peptide segments. A novel auxiliary-mediated native chemical ligation that enables peptide ligation to be applied to protein sequences lacking cysteine was used. A cleavable thiol-containing auxiliary group, 1-phenyl-2-mercaptoethyl, was added to the alpha-amino group of one peptide segment to facilitate amide bond-forming ligation. The amine-linked 1-phenyl-2-mercaptoethyl auxiliary was stable to anhydrous hydrogen fluoride used to cleave and deprotect peptides after solid-phase peptide synthesis. Following native chemical ligation with a thioester-containing segment, the auxiliary group was cleanly removed from the newly formed amide bond by treatment with anhydrous hydrogen fluoride, yielding a full-length unmodified polypeptide product. The resulting polypeptide was reconstituted with heme and folded to form the functional protein molecule. Synthetic wild-type cyt b562 exhibited spectroscopic and electrochemical properties identical to the recombinant protein, whereas the engineered [SeMet(7)]cyt b562 analogue protein was spectroscopically and functionally distinct, with a reduction potential shifted by approximately 45 mV. The use of the 1-phenyl-2-mercaptoethyl removable auxiliary reported here will greatly expand the applicability of total protein synthesis by native chemical ligation of unprotected peptide segments.