Design of a monomeric 23-residue polypeptide with defined tertiary structure

A Stable Miniature Protein with Oxaloacetate Decarboxylase Activity

An 18-residue miniature enzyme, Apoxaldie-1, has been designed, based on the known structure of the neurotoxic peptide apamin. Three lysine residues were introduced on the solvent-exposed face of the apamin alpha-helix to serve as an active site for decarboxylation of oxaloacetate. The oxidised form of Apoxaldie-1, in which two disulfide bonds stabilise the alpha-helix, formed spontaneously. CD spectroscopy measurements revealed that, in its oxidised form, Apoxaldie-1 adopted a stably folded structure, which was lost upon reduction of the disulfide bonds. Despite its small size and the absence of a designed binding pocket, Apoxaldie-1 displayed saturation kinetics in its oxidised form and catalysed the decarboxylation of oxaloacetate at a rate that was almost four orders of magnitude faster than that observed with n-butylamine. This rivals the performance of the best synthetic oxaloacetate decarboxylases reported to date. Unlike those, however, Apoxaldie-1 displayed significant stability. It maintained its secondary structure at temperatures in excess of 75 degrees C, in the presence of high concentrations of guanidinium chloride and at pH values as low as 2.2. Apamin-based catalysts have potential for the generation of miniature peptides that display activity under nonphysiological conditions.

The crystal structure of human endonuclease VIII-like 1 (NEIL1) reveals a zincless finger motif required for glycosylase activity

In prokaryotes, two DNA glycosylases recognize and excise oxidized pyrimidines: endonuclease III (Nth) and endonuclease VIII (Nei). The oxidized purine 8-oxoguanine, on the other hand, is recognized by Fpg (also known as MutM), a glycosylase that belongs to the same family as Nei. The recent availability of the human genome sequence allowed the identification of three human homologs of Escherichia coli Nei. We report here the crystal structure of a human Nei-like (NEIL) enzyme, NEIL1. The structure of NEIL1 exhibits the same overall fold as E. coli Nei, albeit with an unexpected twist. Sequence alignments had predicted that NEIL1 would lack a zinc finger, and it was therefore expected to use a different DNA-binding motif instead. Our structure revealed that, to the contrary, NEIL1 contains a structural motif composed of two antiparallel beta-strands that mimics the antiparallel beta-hairpin zinc finger found in other Fpg/Nei family members but lacks the loops that harbor the zinc-binding residues and, therefore, does not coordinate zinc. This "zincless finger" appears to be required for NEIL1 activity, because mutating a very highly conserved arginine within this motif greatly reduces the glycosylase activity of the enzyme.

Designing proteins from the inside out

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

Simulations of the role of water in the protein-folding mechanism

There are many unresolved questions regarding the role of water in protein folding. Does water merely induce hydrophobic forces, or does the discrete nature of water play a structural role in folding? Are the nonadditive aspects of water important in determining the folding mechanism? To help to address these questions, we have performed simulations of the folding of a model protein (BBA5) in explicit solvent. Starting 10,000 independent trajectories from a fully unfolded conformation, we have observed numerous folding events, making this work a comprehensive study of the kinetics of protein folding starting from the unfolded state and reaching the folded state and with an explicit solvation model and experimentally validated rates. Indeed, both the raw TIP3P folding rate (4.5 +/- 2.5 micros) and the diffusion-constant corrected rate (7.5 +/- 4.2 micros) are in strong agreement with the experimentally observed rate of 7.5 +/- 3.5 micros. To address the role of water in folding, the mechanism is compared with that predicted from implicit solvation simulations. An examination of solvent density near hydrophobic groups during folding suggests that in the case of BBA5, there are water-induced effects not captured by implicit solvation models, including signs of a "concurrent mechanism" of core collapse and desolvation.

De novo designed cyclic-peptide heme complexes

The structural characterization of de novo designed metalloproteins together with determination of chemical reactivity can provide a detailed understanding of the relationship between protein structure and functional properties. Toward this goal, we have prepared a series of cyclic peptides that bind to water-soluble metalloporphyrins (FeIII and CoIII). Neutral and positively charged histidine-containing peptides bind with a high affinity, whereas anionic peptides bind only weakly to the negatively charged metalloporphyrin. Additionally, it was found that the peptide becomes helical only in the presence of the metalloporphyrin. CD experiments confirm that the metalloporphyrin binds specific cyclic peptides with high affinity and with isodichroic behavior. Thermal unfolding experiments show that the complex has "native-like" properties. Finally, NMR spectroscopy produced well dispersed spectra and experimental restraints that provide a high-resolution solution structure of the complexed peptide.

Mechanism-based protein design: attempted "nucleation-condensation" approach to a possible minimal helix-bundle protein

In an intended mechanism-based de novo approach, a 22-mer peptide was so designed as to make it both a stereochemically nucleatable and hydrophobically condensable minimal globular protein. Framework-like nucleation of a triple-helix bundle was targeted by employing as folding nucleators composite beta-turns that could both nucleate helices and place them in close juxtaposition for possible interhelical interaction. To promote the targeted triple-helix bundle to condense as a globular protein, an amphipathic sequence pattern was adopted for possible hydrophobic interhelical interaction. A predominantly helicogenic 22-mer amphipathic peptide was thus designed, punctuating it with composite type II'-III and type II-Asx type beta-turns as the helix nucleators cum chain reversal elements. The peptide made by solid-phase synthesis was shown by NMR and CD to be a nascent and distorted triple-helix bundle in a trifluoroethanol (TFE)-water mixture, but more or less a random coil in water. A fold nucleation effect is evident in the TFE-water mixture, but apparently the hydrophobic effect cannot sustain the peptide conformational order in water. A lack of synergy between folding nucleation and hydrophobic condensation of the peptide is possible. Indeed, a mismatch between the sequential H,P pattern of the peptide and its nascent-type globular fold in a TFE-water mixture is evident based on a simulated annealing study guided by NMR.

Ferrocenyl-Bearing Cyclopseudopeptides as Redox-Switchable Cation Receptors

A family of ferrocenyl-bearing cyclopseudopeptides (1-10) designed for redox-switchable receptors of cations was synthesized. Circular dichroism (CD) and cyclic valtammetry (CV) studies of cation binding properties in both the reduced (K1) and oxidized (K2) forms revealed that the binding preference is mainly governed by the charges and radius of the guest cation as well as by the suitability of the host to accommodate the guest. Particularly worth mentioning is the fact that some synthesized cyclopseudopeptides showed high binding affinity and selectivity toward alkaline-earth ions. For example, the K1 of compound 2 binding with Ca2+ is 4.37 x 10(6) mol(-1) x L and its Ca2+/K+ selectivity is 3.1 x 10(5):1, both values are much greater than those of an excellent natural ionophore, valinomycin (1 x 10(5) mol(-1) x L and 0.33:1, respectively). The linear relationship between the shifts of half-wave potentials (deltaE(1/2)) and the radius/charge [r/(+)] ratios suggests that the sensitivity of electrochemical responses to cation complexation be dominated by repulsion factors between the redox center and the incoming cation guest.

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.

Ab initio folding simulation of the Trp-cage mini-protein approaches NMR resolution

Here, we report a 100 ns molecular dynamics simulation of the folding process of a recently designed autonomous-folding mini-protein designated as tc5b with a new AMBER force field parameter set developed based on condensed-phase quantum mechanical calculations and a Generalized Born continuum solvent model. Starting from its fully extended conformation, our simulation has produced a final structure resembling that of NMR native structure to within 1A main-chain root mean square deviation. Remarkably, the simulated structure stayed in the native state for most part of the simulation after it reached the state. Of greater significance is that our simulation has not only reached the correct main-chain conformation, but also a very high degree of accuracy in side-chain packing conformation. This feat has traditionally been a challenge for ab initio simulation studies. In addition to characterization of the trajectory, comparison of our results to experimental data is also presented. Analysis of the trajectory suggests that the rate-limiting step of folding of this mini-protein is the packing of the Trp side-chain.

Reinvestigation of the proposed folding and self-association of the Neuropeptide Head Activator

The Neuropeptide Head Activator (HA), pGlu-Pro-Pro-Gly-Gly-Ser-Lys-Val-Ile-Leu-Phe (pGlu is pyroglutamic acid), is involved in head-specific growth and differentiation processes in the freshwater coelenterate Hydra attenuata. Peptides of identical sequence have also been isolated from higher-organism tissues such as human and bovine hypothalamus. Early studies by molecular sieve chromatography suggested that HA dimerizes with high affinity (K(d) approximately 1 nM). This dimerization was proposed to occur via antiparallel beta-sheet formation between the Lys(7)-Phe(11) segments in each HA molecule. We conducted biophysical studies on synthetic HA in order to gain insight into its structure and aggregation tendencies. We found by analytical ultracentrifugation that HA is monomeric at low millimolar concentrations. Studies by (1)H-NMR revealed that HA did not adopt any significant secondary structure in solution. We found no NOEs that would support the proposed dimer structure. We probed the propensity of the Lys(7)-Phe(11) fragment to form antiparallel beta-sheet by designing peptides in which two such fragments are joined by a two-residue linker. These peptides were intended to form stable beta-hairpin structures with cross-strand interactions that mimic those of the proposed HA dimer interface. We found that the HA-derived fragments may be induced to form intramolecular beta-sheet, albeit only weakly, when linked by the highly beta-hairpin-promoting D-Pro-Gly turn, but not when linked by the more flexible Gly-Gly unit. These findings suggest that the postulated mode of HA dimerization and the proposed propensity of the molecule to form discrete aggregates with high affinity are incorrect.

D-amino acid residues in peptides and proteins

We have investigated the D-amino acid residues present in Protein Data Bank (PDB) entries, categorizing them into "real" D-residues and artifacts. In polypeptide chains of more than 20 residues, only a single instance of a "real" D-residue, other than those deliberately designed or engineered, was found. This example was the result of a slow chemical epimerization process. Another 12 designed D-residues were found in these longer polypeptide chains. Smaller peptides of 20 or fewer residues contained 479 "real" D-residues, the majority in various gramicidin, actinomycin, or cyclosporin structures. We found 148 PDB entries with "real" D-residues and a further 186, in which all apparent D-residues are artifacts. Investigating the (phi, psi) preferences of the "real" D-residues, we found that the region around (-60 degrees, -45 degrees ) was almost completely unoccupied, even though it is not formally disallowed. We link the low propensity to occupy this region with the alpha-helix destabilizing properties of D-residues.

Absolute comparison of simulated and experimental protein-folding dynamics

Protein folding is difficult to simulate with classical molecular dynamics. Secondary structure motifs such as alpha-helices and beta-hairpins can form in 0.1-10 micros (ref. 1), whereas small proteins have been shown to fold completely in tens of microseconds. The longest folding simulation to date is a single 1- micro s simulation of the villin headpiece; however, such single runs may miss many features of the folding process as it is a heterogeneous reaction involving an ensemble of transition states. Here, we have used a distributed computing implementation to produce tens of thousands of 5-20-ns trajectories (700 micros) to simulate mutants of the designed mini-protein BBA5. The fast relaxation dynamics these predict were compared with the results of laser temperature-jump experiments. Our computational predictions are in excellent agreement with the experimentally determined mean folding times and equilibrium constants. The rapid folding of BBA5 is due to the swift formation of secondary structure. The convergence of experimentally and computationally accessible timescales will allow the comparison of absolute quantities characterizing in vitro and in silico (computed) protein folding.

Contribution of individual zinc ligands to metal binding and peptide folding of zinc finger peptides

Little is known about the contribution of individual zinc-ligating amino acid residues for coupling between zinc binding and protein folding in zinc finger domains. To understand such roles of each zinc ligand, four zinc finger mutant peptides corresponding to the second zinc finger domain of Sp1 were synthesized. In the mutant peptides, glycine was substituted for one of four zinc ligands. Their metal binding and folding properties were spectroscopically characterized and compared to those of the native zinc finger peptide. In particular, the electronic charge-transfer and d-d bands of the Co(II)-substituted peptide complexes were used to examine the metal coordination number and geometry. Fluorescence emission studies revealed that the mutant peptides are capable of binding zinc despite removing one ligand. Circular dichroism results clearly showed the induction of an alpha-helix by zinc binding. In addition, the structures of certain mutant zinc finger peptides were simulated by molecular dynamics calculation. The information indicates that His23 and the hydrophobic core formed between the alpha-helix and the beta-sheet play an essential role in alpha-helix induction. This report demonstrates that each ligand does not contribute equally to alpha-helix formation and coordination geometry in the zinc finger peptide.

The neuropeptide Y monomer in solution is not folded in the pancreatic-polypeptide fold

Fluorescence-labelled analogs of NPY, a 36-amino acid peptide amide, were synthesized by solid-phase peptide synthesis and used for fluorescence-resonance energy transfer studies to investigate the conformation. Energy-transfer efficiency measurements in different media at the concentration of 10 microM are in agreement with a model of the NPY structure proposed by NMR studies (performed at millimolar concentration) in which the C-terminal part of the molecule adopts an alpha-helical conformation while the N-terminal part is flexible. According to this model, the alpha-helix is stabilized by intermolecular hydrophobic interactions because of the formation of dimers. The decrease of the peptide concentration causes a shift of the dimerization equilibrium toward the monomeric form. Energy-transfer efficiency measurements performed at lower concentrations do not support the hypothesis of the folding back of the N-terminal tail onto the C-terminal alpha-helix to yield the so-called "PP-fold" conformation. This structure is observed in the crystal structure of avian pancreatic polypeptide, a member of the NPY peptide hormone family, and it has been considered to be the bioactive one. Our results complete the structural characterization of NPY in solution at concentration ranges in which NMR experiments are not feasible. Furthermore, these results open the way to study the conformation of the receptor-bound ligand.

Emergence of highly designable protein-backbone conformations in an off-lattice model

Despite the variety of protein sizes, shapes, and backbone configurations found in nature, the design of novel protein folds remains an open problem. Within simple lattice models it has been shown that all structures are not equally suitable for design. Rather, certain structures are distinguished by unusually high designability: the number of amino acid sequences for which they represent the unique lowest energy state; sequences associated with such structures possess both robustness to mutation and thermodynamic stability. Here we report that highly designable backbone conformations also emerge in a realistic off-lattice model. The highly designable conformations of a chain of 23 amino acids are identified and found to be remarkably insensitive to model parameters. Although some of these conformations correspond closely to known natural protein folds, such as the zinc finger and the helix-turn-helix motifs, others do not resemble known folds and may be candidates for novel fold design.

Improving coiled-coil stability by optimizing ionic interactions

Alpha-helical coiled coils are a common protein oligomerization motif stabilized mainly by hydrophobic interactions occurring along the coiled-coil interface. We have recently designed and solved the structure of a two-heptad repeat coiled-coil peptide that is stabilized further by a complex network of inter- and intrahelical salt-bridges in addition to the hydrophobic interactions. Here, we extend and improve the de novo design of this two heptad-repeat peptide by four newly designed peptides characterized by different types of ionic interactions. The contribution of these different types of ionic interactions to coiled-coil stability are analyzed by CD spectroscopy and analytical ultracentrifugation. We show that all peptides are highly alpha-helical and two of them are 100% dimeric under physiological conditions. Furthermore, we have solved the X-ray structure of the most stable of these peptides and the rational design principles are verified by comparing this structure to the structure of the parent peptide. We show that by combining the most favorable inter- and intrahelical salt-bridge arrangements it is possible to design coiled-coil oligomerization domains with improved stability properties.

A new zinc binding fold underlines the versatility of zinc binding modules in protein evolution

Many different zinc binding modules have been identified. Their abundance and variety suggests that the formation of zinc binding folds might be relatively common. We have determined the structure of CH1(1), a 27-residue peptide derived from the first cysteine/histidine-rich region (CH1) of CREB binding protein (CBP). This peptide forms a highly ordered zinc-dependent fold that is distinct from known folds. The structure differs from a subsequently determined structure of a larger region from the CH3 region of CBP, and the CH1(1) fold probably represents a nonphysiologically active form. Despite this, the fold is thermostable and tolerant to both multiple alanine mutations and changes in the zinc-ligand spacing. Our data support the idea that zinc binding domains may arise frequently. Additionally, such structures may prove useful as scaffolds for protein design, given their stability and robustness.

Infinite pleated beta -sheet formed by the beta-hairpin Boc-beta-Phe-beta-Phe-D-Pro-Gly-beta-Phe-beta-Phe-OMe

A beta-hairpin conformation and extended beta-pleated sheet assembly have been characterized by single crystal x-ray diffraction for the synthetic peptide t-butoxycarbonyl--beta-Phe-beta-Phe-D-Pro-Gly-beta-Phe-beta-Phe-methyl ester [beta-Phe: (S)-beta(3) homophenylalanine]. The centrally located D-Pro-Gly segment nucleates a chain reversal in a type II' beta-turn conformation. Two intramolecular cross-strand hydrogen bonds stabilize the peptide fold. Intermolecular NH...O[double bond]C hydrogen bonds (two on each side of the hairpin) connect the hairpins into an infinitely extended beta-sheet. The beta-residues cause all C[double bond]O groups to point in the same direction, resulting in a "polar" sheet by the unidirectional alignment of NH...O[double bond]C hydrogen bonds. In contrast, beta-sheets formed by alpha-residues have alternating directions for the hydrogen bonds, thus resulting in an "apolar" sheet. The crystallographic parameters for C(53)H(66)N(6)O(9) x CH(3)OH are: space group P2(1), a = 9.854(2) A, b = 10.643(2) A, c = 25.296(4) A, beta = 100.39(2) degrees, Z = 2, agreement factor R(1) = 0.065 for 3,706 data observed >4 sigma(F) and a resolution of 0.90 A.

Rationally designed mutations convert de novo amyloid-like fibrils into monomeric beta-sheet proteins

Amyloid fibrils are associated with a variety of neurodegenerative maladies including Alzheimer's disease and the prion diseases. The structures of amyloid fibrils are composed of beta-strands oriented orthogonal to the fibril axis ("cross beta" structure). We previously reported the design and characterization of a combinatorial library of de novo beta-sheet proteins that self-assemble into fibrillar structures resembling amyloid. The libraries were designed by using a "binary code" strategy, in which the locations of polar and nonpolar residues are specified explicitly, but the identities of these residues are not specified and are varied combinatorially. The initial libraries were designed to encode proteins containing amphiphilic beta-strands separated by reverse turns. Each beta-strand was designed to be seven residues long, with polar (open circle) and nonpolar (shaded circle) amino acids arranged with an alternating periodicity ([see text]). The initial design specified the identical polar/nonpolar pattern for all of the beta-strands; no strand was explicitly designated to form the edges of the resulting beta-sheets. With all beta-strands preferring to occupy interior (as opposed to edge) locations, intermolecular oligomerization was favored, and the proteins assembled into amyloid-like fibrils. To assess whether explicit design of edge-favoring strands might tip the balance in favor of monomeric beta-sheet proteins, we have now redesigned the first and/or last beta-strands of several sequences from the original library. In the redesigned beta-strands, the binary pattern is changed from [see text] (K denotes lysine). The presence of a lysine on the nonpolar face of a beta-strand should disfavor fibrillar structures because such structures would bury an uncompensated charge. The nonpolar right arrow lysine mutations, therefore, would be expected to favor monomeric structures in which the [see text] sequences form edge strands with the charged lysine side chain accessible to solvent. To test this hypothesis, we constructed several second generation sequences in which the central nonpolar residue of either the N-terminal beta-strand or the C-terminal beta-strand (or both) is changed to lysine. Characterization of the redesigned proteins shows that they form monomeric beta-sheet proteins.

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.

Oligomeric beta(beta)(alpha) miniprotein motifs: pivotal role of single hinge residue in determining the oligomeric state

The role of a single glycine hinge residue in the structure of BBAT1, a beta(beta)(alpha) peptide that forms a discrete homotrimeric structure in solution, was evaluated with 11 new peptide sequences which differ only in the identity of the residue at the hinge position. The integrity of the structure and oligomeric state of the peptides was evaluated by using a combination of analytical ultracentrifugation and circular dichroism spectroscopy. Initially, it was discovered that the glycine hinge adopts backbone dihedral angles favored in D-amino acids and that incorporation of D-alanine at the hinge position stabilizes the trimer species. Subsequently, the effect of the side chains of different D-amino acids at the hinge position was evaluated. While incorporation of polar amino acids led to a destabilization of the oligomeric form of the peptide, only peptides including D-Ser or D-Asp at the hinge position were able to achieve a discrete trimer species. Incorporation of hydrophobic amino acids D-Leu and D-Phe led to oligomerization beyond a trimer to a tetrameric form. The dramatic differences among the thermodynamic stabilities and oligomeric states of these peptides illustrates the pivotal role of the hinge residue in the oligomerization of the beta(beta)(alpha) peptides.

Length-dependent stability and strand length limits in antiparallel beta -sheet secondary structure

Designed peptides that fold autonomously to specific conformations in aqueous solution are useful for elucidating protein secondary structural preferences. For example, autonomously folding model systems have been essential for establishing the relationship between alpha-helix length and alpha-helix stability, which would be impossible to probe with alpha-helices embedded in folded proteins. Here, we use designed peptides to examine the effect of strand length on antiparallel beta-sheet stability. alpha-Helices become more stable as they grow longer. Our data show that a two-stranded beta-sheet ("beta-hairpin") becomes more stable when the strands are lengthened from five to seven residues, but that further strand lengthening to nine residues does not lead to further beta-hairpin stabilization for several extension sequences examined. (In one case, all-threonine extension, there may be an additional stabilization on strand lengthening from seven to nine residues.) These results suggest that there may be an intrinsic limit to strand length for most sequences in antiparallel beta-sheet secondary structure.

Rapid boundary element solvation electrostatics calculations in folding simulations: successful folding of a 23-residue peptide

Solvation effects play a profound role in the energetics of protein folding. While a continuum dielectric model of solvation may provide a sufficiently accurate estimate of the solvation effects, until now this model was too computationally expensive and unstable for folding simulations. Here we proposed a fast yet accurate and robust implementation of the boundary element solution of the Poisson equation, the REBEL algorithm. Using our earlier double-energy scheme, we included for the first time the mathematically rigorous continuous REBEL solvation term in our Biased Probability Monte Carlo (BPMC) simulations of the peptide folding. The free energy of a 23-residue beta beta alpha-peptide was then globally optimized with and without the solvation electrostatics contribution. An ensemble of beta beta alpha conformations was found at and near the global minimum of the energy function with the REBEL electrostatic solvation term. Much poorer correspondence to the native solution structure was found in the "control" simulations with a traditional method to account for solvation via a distance-dependent dielectric constant. Each simulation took less than 40 h (21 h without electrostatic solvation calculation) on a single Alpha 677 MHz CPU and involved more than 40,000 solvation energy evaluations. This work demonstrates for the first time that such a simulation can be performed in a realistic time frame. The proposed procedure may eliminate the energy evaluation accuracy bottleneck in folding simulations.

Interstrand side chain--side chain interactions in a designed beta-hairpin: significance of both lateral and diagonal pairings

The contributions of interstrand side chain-side chain contacts to beta-sheet stability have been examined with an autonomously folding beta-hairpin model system. RYVEV(D)PGOKILQ-NH2 ((D)P = D-proline, O = ornithine) has previously been shown to adopt a beta-hairpin conformation in aqueous solution, with a two-residue loop at D-Pro-Gly. In the present study, side chains that display interstrand NOEs (Tyr-2, Lys-9, and Leu-11) are mutated to alanine or serine, and the conformational impact of the mutations is assessed. In the beta-hairpin conformation Tyr-2 and Leu-11 are directly across from one another (non-hydrogen bonded pair). This "lateral" juxtaposition of two hydrophobic side chains appears to contribute to beta-hairpin conformational stability, which is consistent with results from other beta-sheet model studies and with statistical analyses of interstrand residue contacts in protein crystal structures. Interaction between the side chains of Tyr-2 and Lys-9 also stabilizes the beta-hairpin conformation. Tyr-2/Lys-9 is a "diagonal" interstrand juxtaposition because these residues are not directly across from one another in terms of the hydrogen bonding registry between the strands. This diagonal interaction arises from the right-handed twist that is commonly observed among beta-sheets. Evidence of diagonal side chain-side chain contacts has been observed in other autonomously folding beta-sheet model systems, but we are not aware of other efforts to determine whether a diagonal interaction contributes to beta-sheet stability.

Effects of turn residues in directing the formation of the beta-sheet and in the stability of the beta-sheet

The designed peptide (denoted 20-mer, sequence VFITS(D)PGKTYTEV(D)PGOKILQ) has been shown to form a three-strand antiparallel beta-sheet. It is generally believed that the (D)Pro-Gly segment has the propensity to adopt a type II' beta-turn, thereby promoting the formation of this beta-sheet. Here, we replaced (D)Pro-Gly with Asp-Gly, which should favor a type I' turn, to examine the influence of different type of turns on the stability of the beta-sheet. Contrary to our expectation, the mutant peptide, denoted P6D, forms a five-residue type I turn plus a beta-bulge between the first two strands due to a one amino-acid frameshift in the hydrogen bonding network and side-chain inversion of the first beta-strand. In contrast, the same kind of substitution at (D)Pro-14 in the double mutant, denoted P6DP14D, does not yield the same effect. These observations suggest that the SDGK sequence disfavors the type I' conformation while the VDGO sequence favors a type I' turn, and that the frameshift in the first strand provides a way for the peptide to accommodate a disfavored turn sequence by protruding a bulge in the formation of the beta-hairpin. Thus, different types of turns can affect the stability of a beta-structure.

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.

Crystal structure of a synthetic cyclodecapeptide for template-assembled synthetic protein design

The structural prototype of a new generation of regioselectively addressable functionalized templates (RAFTs) for use in protein de novo design has been synthesized and crystallized. The structure of the aromatically substituted cyclodecapeptide was determined by X-ray diffraction; it consists of an antiparallel beta sheet spanned by heterochirally induced type IIprime prime or minute beta turns, similar to that observed in gramicidin S. The three-dimensional structure of the artificial template was also examined by an NMR spectroscopic analysis in solution and shown to be compatible with a beta-sheet plane suitable for accommodating secondary functional peptide fragments for the synthesis of template-assembled synthetic proteins (TASPs).

Molecular carpentry: piecing together helices and hairpins in designed peptides

The design of a peptide that contains two distinct elements of secondary structure, helix and beta-hairpin, is described. Two designed 17-residue peptides: Boc-Val-Ala-Leu-Aib-Val-Ala-Leu-Gly-Gly-Leu-Phe-Val-D-Pro-Gly-Leu-Phe-Val-OMe (I) and Boc-Leu-Aib-Val-Ala-Leu-Aib-Val-Gly-Gly-Leu-Val-Val-D-Pro-Gly-Leu-Val-Val-OMe (II) have been conformationally characterized by NMR spectroscopy. Peptides I and II contain a seven-residue helical module at the N terminus and a eight-residue beta-hairpin module at the C terminus, which are connected by a conformationally flexible Gly-Gly segment. The choice of the secondary-structure modules is based upon prior crystallographic and spectroscopic analysis of the individual modules. Analysis of 500 MHz 1H NMR data, recorded as solutions in methanol, suggests that the observed pattern of chemical shifts, 3JHN CalphaH values, temperature coefficients of the NH chemical shifts, and backbone inter-residue nuclear Overhauser effects favor helical structures for residues 1-7 and beta-hairpin structures for residues 10-17. The spectroscopic data are compatible with termination of the helical segment by formation of a Schellman motif; this restricts Gly(8) to a left-handed alpha-helical conformation. Gly(9) is the only residue with multiple conformational possibilities in phi,psi space. Possible orientations of the two secondary-structure modules are considered. This study validates the use of stereochemically rigid peptide modules as prefabricated elements in the construction of synthetic protein mimics.

Tryptophan zippers: stable, monomeric beta -hairpins

A structural motif, the tryptophan zipper (trpzip), greatly stabilizes the beta-hairpin conformation in short peptides. Peptides (12 or 16 aa in length) with four different turn sequences are monomeric and fold cooperatively in water, as has been observed previously for some hairpin peptides. However, the folding free energies of the trpzips exceed substantially those of all previously reported beta-hairpins and even those of some larger designed proteins. NMR structures of three of the trpzip peptides reveal exceptionally well-defined beta-hairpin conformations stabilized by cross-strand pairs of indole rings. The trpzips are the smallest peptides to adopt an unique tertiary fold without requiring metal binding, unusual amino acids, or disulfide crosslinks.

Oligomerization of uniquely folded mini-protein motifs: development of a homotrimeric betabetaalpha peptide

The discovery of a discretely folded homotrimeric betabetaalpha motif (BBAT1) was recently reported (J. Am. Chem. Soc. 2001, 123, 1002-1003). Herein the design, synthesis, and analysis of a small library of peptides which led to the isolation of BBAT1 is described. betabetaalpha peptides based on the monomeric sequence of BBA5 (Folding Des. 1998, 120, 95-103) were synthesized to include the anthranilic acid/nitrotyrosine fluorescence quenching pair to rapidly detect interpeptide association. In the first generation of peptides synthesized, truncations in the loop region connecting the beta-hairpin to the alpha-helix revealed that a two-residue deletion in the loop promoted an interpeptide association as detected by fluorescence quenching. An additional library of 22 loop-truncated betabetaalpha peptides was subsequently synthesized to include a variety of sequence mutations in an effort to enhance the observed peptide-peptide binding. From the fluorescence quenching screen, peptide B2 was found to possess the strongest fluorescence-quenching response, indicative of a strong peptide-peptide association. Due the poor solubility of peptide B2, the S-methylated cysteine at position 9 in the loop was substituted with a glycine to generate peptide BBAT1 which possessed greatly improved water solubility and formed discrete trimers. The successful design of this oligomeric betabetaalpha structure will likely aid the design of more complex alpha-beta superstructures and further our understanding of the factors controlling protein-protein interactions at alpha-beta protein interfaces.

Rescuing a destabilized protein fold through backbone cyclization

We describe the physicochemical characterization of various circular and linear forms of the approximately 60 residue N-terminal Src homology 3 (SH3) domain from the murine c-Crk adapter protein. Structural, dynamic, thermodynamic, kinetic and biochemical studies reveal that backbone circularization does not prevent the adoption of the natural folded structure in any of the circular proteins. Both the folding and unfolding rate of the protein increased slightly upon circularization. Circularization did not lead to a significant thermodynamic stabilization of the full-length protein, suggesting that destabilizing enthalpic effects (e.g. strain) negate the expected favorable entropic contribution to overall stability. In contrast, we find circularization results in a dramatic stabilization of a truncated version of the SH3 domain lacking a key glutamate residue. The ability to rescue the destabilized mutant indicates that circularization may be a useful tool in protein engineering programs geared towards generating minimized proteins.

Peptide hybrids containing alpha - and beta-amino acids: structure of a decapeptide beta-hairpin with two facing beta-phenylalanine residues

A beta-hairpin conformation has been characterized in crystals of the decapeptide t-butoxycarbonyl-Leu-Val-beta Phe-Val-(D)Pro-Gly-Leu-beta Phe-Val-Val-methyl ester [beta Phe; (S)-beta(3) homophenylalanine] by x-ray diffraction. The polypeptide chain reversal is nucleated by the centrally positioned (D)Pro-Gly segment, which adopts a type-I' beta-turn conformation. Four intramolecular cross-strand hydrogen bonds stabilize the peptide fold. The beta Phe(3) and beta Phe(8) residues occupy facing positions on the hairpin, with the side chains projecting on opposite faces of the beta-sheet. At the site of insertion of beta-residues, the polarity of the peptide units along each strand reverses, as compared with the alpha-peptide segments. In this analog, a small segment of a polar sheet is observed, where adjacent CO and NH groups line up in opposite directions in each strand. In the crystal, an extended beta-sheet is formed by hydrogen bonding between strands of antiparallel pairs of beta-hairpins. The crystallographic parameters for C(65)H(102)N(10)O(13) x 3H(2)O are: space group P2(1)2(1)2(1); a = 19.059(8) A, b = 19.470(2) A, c = 21.077(2) A; Z = 4; agreement factor R(1) = 9.12% for 3,984 data observed >4 sigma(F) and a resolution of 0.90 A.

Designed beta-hairpin peptides with defined tight turn stereochemistry

The conformational analysis of two synthetic octapeptides, Boc-Leu-Val-Val-D-Pro-L-Ala-Leu-Val-Val-OMe (1) and Boc-Leu-Val-Val-D-Pro-D-Ala-Leu-Val-Val-OMe (2) has been carried out in order to investigate the effect of beta-turn stereochemistry on designed beta-hairpin structures. Five hundred megahertz (1)H NMR studies establish that both peptides 1 and 2 adopt predominantly beta-hairpin conformations in methanol solution. Specific nuclear Overhauser effects provide evidence for a type II' beta-turn conformation for the D-Pro-L-Ala segment in 1, while the NMR data suggest that the type I' D-Pro-D-Ala beta-turn conformation predominates in peptide 2. Evidence for a minor conformation in peptide 2, in slow exchange on the NMR time scale, is also presented. Interstrand registry is demonstrated in both peptides 1 and 2. The crystal structure of 1 reveals two independent molecules in the crystallographic asymmetric unit, both of which adopt beta-hairpin conformations nucleated by D-Pro-L-Ala type II' beta-turns and are stabilized by three cross-strand hydrogen bonds. CD spectra for peptides 1 and 2 show marked differences, presumably as a consequence of the superposition of spectral bands arising from both beta-turn and beta-strand conformations.

Large anhydrous polyalanine ions: evidence for extended helices and onset of a more compact state

Ion mobility measurements and molecular modeling calculations have been used to examine the conformations of large multiply charged polyalanine peptides. Two series of [Ala(n)+3H](3+) conformations which do not interconvert during the 10 to 30 ms experimental timescales are observed: a family of elongated structures for n = 18 to 39 and a series of more compact conformations for n = 24 to 41. The more compact state becomes the dominant conformer type for n > 32. Molecular modeling studies and comparisons of calculated collision cross sections with experiment indicate that the elongated ions have extended helical conformations. We suggest that the more compact state corresponds to a new conformer type: a folded hinged helix-coil state in which helical and coil regions have similar physical dimensions. The competition between extended and compact states is rationalized by considering differences in charge stabilization and entropy.

The role of the hinge region of the luteinizing hormone receptor in hormone interaction and signal generation

Luteinizing hormone receptor, a G protein-coupled receptor, consists of two halves, the N-terminal extracellular hormone binding domain (exodomain) and the C-terminal membrane-associated, signal-generating domain (endodomain). The exodomain has seven to nine Leu-rich repeats, which are generally thought to form a 1/3 donut-like structure and interact with human choriogonadotropin (hCG). The resulting hCG-exodomain complex adjusts the structure and its association with the endodomain, which results in signal generation in the endodomain. It is unclear whether the rigid 1/3 donut structure could provide the agility and versatility of this dynamic action. In addition, there is no clue as to where the endodomain contact point (the signal modulator) in the exodomain is. To address these issues, the exodomain was examined by Ala scan and multiple substitutions, while receptor peptides were used for photoaffinity labeling and affinity cross-linking. Our results show that the C-flanking sequence (hinge region), Thr(250)-Gln(268), of the Leu-rich repeats (LRRs) specifically interacts with hCG, preferentially hCGalpha. This interaction is inhibited by exoloop 2 of the endodomain but not by exoloops 1 and 3, suggesting an intimate relationship between Thr(250)-Gln(268), exoloop 2, and hCG. Taken together, our observations in this article suggest a new paradigm that the LRRs contact the front of hCG, while both flanking regions of the LRRs interact with the sides of hCG. This would trap hCG in the 1/3 donut structure of the LRRs and enhance the binding affinity. In addition, mutations of conserved Ser(255) in the sequence can constitutively activate the receptor. This provides a clue for the signal modulator in the exodomain. In contrast, a phenyl or phenolic group is necessary at conserved Tyr(253) for targeting the receptor to the surface.