Characterizing a partially folded intermediate of the villin headpiece domain under non-denaturing conditions: contribution of His41 to the pH-dependent stability of the N-terminal subdomain

The contribution of interactions involving the imidazole ring of His41 to the pH-dependent stability of the villin headpiece (HP67) N-terminal subdomain has been investigated by nuclear magnetic resonance (NMR) spin relaxation. NMR-derived backbone N-H order parameters (S2) for wild-type (WT) HP67 and H41Y HP67 indicate that reduced conformational flexibility of the N-terminal subdomain in WT HP67 is due to intramolecular interactions with the His41 imidazole ring. These interactions, together with desolvation effects, contribute to significantly depress the pKa of the buried imidazole ring in the native state. 15N R1rho relaxation dispersion data indicate that WT HP67 populates a partially folded intermediate state that is 10.9 kJ mol(-1) higher in free energy than the native state under non-denaturing conditions at neutral pH. The partially folded intermediate is characterized as having an unfolded N-terminal subdomain while the C-terminal subdomain retains a native-like fold. Although the majority of the residues in the N-terminal subdomain sample a random-coil distribution of conformations, deviations of backbone amide 1H and 15N chemical shifts from canonical random-coil values for residues within 5A of the His41 imidazole ring indicate that a significant degree of residual structure is maintained in the partially folded ensemble. The pH-dependence of exchange broadening is consistent with a linear three-state exchange model whereby unfolding of the N-terminal subdomain is coupled to titration of His41 in the partially folded intermediate with a pKa,I=5.69+/-0.07. Although maintenance of residual interactions with the imidazole ring in the unfolded N-terminal subdomain appears to reduce pKa,I compared to model histidine compounds, protonation of His41 disrupts these interactions and reduces the difference in free energy between the native state and partially folded intermediate under acidic conditions. In addition, chemical shift changes for residues Lys70-Phe76 in the C-terminal subdomain suggest that the HP67 actin binding site is disrupted upon unfolding of the N-terminal subdomain, providing a potential mechanism for regulating the villin-dependent bundling of actin filaments.     

Free-energy landscape of the villin headpiece in an all-atom force field

We investigate the landscape of the internal free-energy of the 36 amino acid villin headpiece with a modified basin hopping method in the all-atom force field PFF01, which was previously used to predictively fold several helical proteins with atomic resolution. We identify near native conformations of the protein as the global optimum of the force field. More than half of the twenty best simulations started from random initial conditions converge to the folding funnel of the native conformation, but several competing low-energy metastable conformations were observed. From 76,000 independently generated conformations we derived a decoy tree which illustrates the topological structure of the entire low-energy part of the free-energy landscape and characterizes the ensemble of metastable conformations. These emerge as similar in secondary content, but differ in tertiary arrangement.     

Empirical solvation models in the context of conformational energy searches: application to bovine pancreatic trypsin inhibitor.

Continuum solvation models that estimate free energies of solvation as a function of solvent accessible surface area are computationally simple enough to be useful for predicting protein conformation. The behavior of three such solvation models has been examined by applying them to the minimization of the conformational energy of bovine pancreatic trypsin inhibitor. The models differ only with regard to how the constants of proportionality between free energy and surface area were derived. Each model was derived by fitting to experimentally measured equilibrium solution properties. For two models, the solution property was free energy of hydration. For the third, the property was NMR coupling constants. The purpose of this study is to determine the effect of applying these solvation models to the nonequilibrium conformations of a protein arising in the course of global searches for conformational energy minima. Two approaches were used: (1) local energy minimization of an ensemble of conformations similar to the equilibrium conformation and (2) global search trajectories using Monte Carlo plus minimization starting from a single conformation similar to the equilibrium conformation. For the two models derived from free energy measurements, it was found that both the global searches and local minimizations yielded conformations more similar to the X-ray crystallographic structures than did searches or local minimizations carried out in the absence of a solvation component of the conformational energy. The model derived from NMR coupling constants behaved similarly to the other models in the context of a global search trajectory. For one of the models derived from measured free energies of hydration, it was found that minimization of an ensemble of near-equilibrium conformations yielded a new ensemble in which the conformation most similar to the X-ray determined structure PTI4 had the lowest total free energy. Despite the simplicity of the continuum solvation models, the final conformation generated in the trajectories for each of the models exhibited some of the characteristics that have been reported for conformations obtained from molecular dynamics simulations in the presence of a bath of explicit water molecules. They have smaller root mean square (rms) deviations from the experimentally determined conformation, fewer incorrect hydrogen bonds, and slightly larger radii of gyration than do conformations derived from search trajectories carried out in the absence of solvent.     

Native-like mean structure in the unfolded ensemble of small proteins

The nature of the unfolded state plays a great role in our understanding of proteins. However, accurately studying the unfolded state with computer simulation is difficult, due to its complexity and the great deal of sampling required. Using a supercluster of over 10,000 processors we have performed close to 800 micros of molecular dynamics simulation in atomistic detail of the folded and unfolded states of three polypeptides from a range of structural classes: the all-alpha villin headpiece molecule, the beta hairpin tryptophan zipper, and a designed alpha-beta zinc finger mimic. A comparison between the folded and the unfolded ensembles reveals that, even though virtually none of the individual members of the unfolded ensemble exhibits native-like features, the mean unfolded structure (averaged over the entire unfolded ensemble) has a native-like geometry. This suggests several novel implications for protein folding and structure prediction as well as new interpretations for experiments which find structure in ensemble-averaged measurements.     

Conformational characteristics of unstructured peptides: alpha-synuclein

We have performed replica-exchange molecular dynamics simulations on 41 residue peptides containing NAC region of alpha-synuclein in various force fields and solvent conditions. Alpha-synuclein is known to be the major cause of Parkinson's disease by amyloid-like aggregation, and one of the natively unfolded proteins. To investigate conformational characteristics of intrinsically unstructured peptides, we carried out structural analysis by introducing 'representative structure' for ensemble of structures occurring during the overall trajectory. Representative structures may be defined by using either coordinate averaging or distance averaging. When applied to the natively folded proteins such as villin headpiece, structural analysis based on representative structure was found to yield consistent results with those obtained from conventional analysis. Individual conformations obtained from the simulations of NAC peptide for various conditions show flexible structures close to random coil. Secondary structure contents and free energy surfaces showed dependency on solvent conditions, which may be interpreted as another manifestation of structural diversity. It is found that representative structures can provide useful information about structural characteristics of intrinsically unstructured proteins.     

Heterogeneity even at the speed limit of folding: large-scale molecular dynamics study of a fast-folding variant of the villin headpiece

We have performed molecular dynamics simulations on a set of nine unfolded conformations of the fastest-folding protein yet discovered, a variant of the villin headpiece subdomain (HP-35 NleNle). The simulations were generated using a new distributed computing method, yielding hundreds of trajectories each on a time scale comparable to the experimental folding time, despite the large (10,000 atom) size of the simulation system. This strategy eliminates the need to assume a two-state kinetic model or to build a Markov state model. The relaxation to the folded state at 300 K from the unfolded configurations (generated by simulation at 373 K) was monitored by a method intended to reflect the experimental observable (quenching of tryptophan by histidine). We also monitored the relaxation to the native state by directly comparing structural snapshots with the native state. The rate of relaxation to the native state and the number of resolvable kinetic time scales both depend upon starting structure. Moreover, starting structures with folding rates most similar to experiment show some native-like structure in the N-terminal helix (helix 1) and the phenylalanine residues constituting the hydrophobic core, suggesting that these elements may exist in the experimentally relevant unfolded state. Our large-scale simulation data reveal kinetic complexity not resolved in the experimental data. Based on these findings, we propose additional experiments to further probe the kinetics of villin folding.     

Coupling Between Folding and Ionization Equilibria: Effects of pH on the Conformational Preferences of Polypeptides

A new approach to the conformational study of polypeptides is presented. It considers explicitly the coupling between the conformation of the molecule and the ionization equilibria at a given pH value. Calculations of the solvation free energy and free energy of ionization of a 17-residue polypeptide are carried out using a fast multigrid boundary element method (MBE). The MBE method uses an adaptive tessellation of the molecular surface by boundary elements with non-regular size to solve the Poisson equation rapidly, and with a high degree of accuracy. The MBE method is integrated into the ECEPP (Empirical Conformational Energy Program for Peptides) algorithm to compute the coupling between the ionization state and the conformation of the molecule.This approach has been applied to study the conformational preference of a short polypeptide for which the available NMR and CD experimental data indicate that conformations containing a right-handed α-helical segment are energetically more favorable at low values of pH. The results of calculations using the present method agree quite well with experiments, in contrast to previous applications with standard techniques (using pre-assigned charges at each pH) that were not able to reproduce the experimental findings. Also, it is shown how the coupling to the conformation leads to different degrees of ionization of a given type of residue, for example glutamic acid, at different positions in the amino acid sequence, at any given pH. The results of this study provide a sound basis to discuss the origin of the stability of polypeptide conformations, and its dependence on the environmental conditions.     

Multistate folding of the villin headpiece domain

The villin headpiece (HP67) is a 67 residue, monomeric protein derived from the C-terminal domain of villin. Wild-type HP67 (WT HP67) is the smallest fragment of villin that retains strong in vitro actin-binding activity. WT HP67 is made up of two subdomains, which form a tightly packed interface. The C-terminal subdomain of WT HP67, denoted HP35, is rich in helical structure, folds in isolation, and has been widely used as a model system for folding studies. In contrast, very little is known about the folding of the intact villin headpiece domain. Here, NMR, CD and H/2H amide exchange measurements are used to follow the pH, thermal and urea-induced unfolding of WT HP67 and a mutant (HP67 H41Y) in which a buried conserved histidine in the N-terminal subdomain, His41, has been mutated to Tyr. Although most small proteins display two-state equilibrium unfolding, the results presented here demonstrate that unfolding of the villin headpiece is a multistate process. The presence of a folded N-terminal subdomain is shown to stabilize the C-terminal subdomain, increasing the midpoints of the thermal and urea-induced unfolding transitions and increasing protection factors for H/2H exchange. Histidine 41 has been shown to act as a pH-dependent switch in wild-type HP67: the N-terminal subdomain is unfolded when His41 is protonated, while the C-terminal subdomain remains folded irrespective of the protonation state of His41. Mutation of His41 to Tyr eliminates the segmental pH-dependent unfolding of the headpiece. The mutation stabilizes both domains, but folding is still multistate, indicating that His41 is not solely responsible for the unusual equilibrium unfolding behavior of villin headpiece domain.     

Rapid generation of a representative ensemble of N-glycan conformations

Glycosylated proteins are ubiquitous components of extracellular matrices and cellular surfaces where their oligosaccharide moieties are implicated in a wide range of cell-cell and cell-matrix recognition events. Glycans constitute highly flexible molecules. Only a small number of glycan X-ray structures is available for which sufficient electron density for an entire oligosaccharide chain has been observed. An unambiguous structure determination based on NMR-derived geometric constraints alone is often not possible. Time consuming computational approaches such as Monte Carlo calculations and molecular dynamics simulations have been widely used to explore the conformational space accessible to complex carbohydrates. The generation of a comprehensive data base for N-glycan fragments based on long time molecular dynamics simulations is presented. The fragments are chosen in such a way that the effects of branched N-glycan structures are taken into account. The prediction database constitutes the basis of a procedure to generate a complete set of all possible conformations for a given N-glycan. The constructed conformations are ranked according to their energy content. The resulting conformations are in reasonable agreement with experimental data. A web interface has been established (http://www.dkfz.de/spec/glydict/), which enables to input any N-glycan of interest and to receive an ensemble of generated conformations within a few minutes.     

Efficient high level expression of peptides and proteins as fusion proteins with the N-terminal domain of L9: application to the villin headpiece helical subdomain

The efficient expression of small to midsize polypeptides and small marginally stable proteins can be difficult. A new protein fusion system is developed to allow the expression of peptides and small proteins. The polypeptide of interest is linked via a Factor Xa cleavage sequence to the C-terminus of the N-terminal domain of the ribosomal protein L9 (NTL9). NTL9 is a small (56 residue) basic protein. The C-terminus of the protein is part of an alpha-helix which extends away from the globular structure thus additional domains can be fused without altering the fold of NTL9. NTL9 expresses at high levels, is extremely soluble, and remains fully folded over a wide temperature and pH range. The protein has a high net positive charge, facilitating purification of fusion proteins by ion exchange chromatography. NTL9 fusions can also be easily purified by reverse phase HPLC. As a test case we demonstrate the high level expression of a small, 36 residue, three helix bundle, the villin headpiece subdomain. This protein is widely used as a model system for folding studies and the development of a simple expression system should facilitate experimental studies of the subdomain. The yield of purified fusion protein is 70 mg/L of culture and the yield of purified villin headpiece subdomain is 24 mg/L of culture. We also demonstrate the use of the fusion system to express a smaller marginally folded peptide fragment of the villin headpiece domain.     

Interplay of charge distribution and conformation in peptides: comparison of theory and experiment

We assessed the correlation between charge distribution and conformation of flexible peptides by comparing the theoretically calculated potentiometric-titration curves of two model peptides, Ac-Lys5-NHMe (a model of poly-L-lysine) and Ac-Lys-Ala11-Lys-Gly2-Tyr-NH2 (P1) in water and methanol, with the experimental curves. The calculation procedure consisted of three steps: (i) global conformational search of the peptide under study using the electrostatically driven Monte Carlo (EDMC) method with the empirical conformational energy program for peptides (ECEPP)/3 force field plus the surface-hydration (SRFOPT) or the generalized Born surface area (GBSA) solvation model as well as a molecular dynamics method with the assisted model building and energy refinement (AMBER)99/GBSA force field; (ii) reevaluation of the energy in the pH range considered by using the modified Poisson-Boltzmann approach and taking into account all possible protonation microstates of each conformation, and (iii) calculation of the average degree of protonation of the peptide at a given pH value by Boltzmann averaging over conformations. For Ac-Lys5-NHMe, the computed titration curve agrees qualitatively with the experimental curve of poly-L-lysine in 95% methanol. The experimental titration curves of peptide P1 in water and methanol indicate a remarkable downshift of the first pK(a) value compared to the values for reference compounds (n-butylamine and phenol, respectively), suggesting the presence of a hydrogen bond between the tyrosine hydroxyl oxygen and the H(epsilon) proton of a protonated lysine side chain. The theoretical titration curves agree well with the experimental curves, if conformations with such hydrogen bonds constitute a significant part of the ensemble; otherwise, the theory predicts too small a downward pH shift.     

A phage display-based method for determination of relative affinities of mutants. Application of the actin-binding motifs in thymosin beta 4 and the villin headpiece

We propose phage display combined with enzyme-linked immunosorbent assay as a tool for the systematic analysis of protein-protein interactions by investigating the binding behavior of variants to a partner protein. Via enzyme-linked immunosorbent assay we determine both the amount of fusion protein presented at the phage surface and the amount of complex formed, the ratio of which is proportional to the affinity. Hence this method enables us to calculate the relative affinities of a large number of mutants. As model systems, we investigated actin-binding motifs conserved in a number of proteins binding monomeric or filamentous actin. The hexapeptide motifs LKKTET, present in thymosin beta4, and LKKEKG, present in the villin headpiece, were mutated, and the variants were analyzed. Study of the positional tolerance allows postulating that the motifs, although similar in primary structures adopt different conformations when bound to actin. In addition, our data show that the second and the fourth amino acid of the thymosin beta4 motif and the first three residues of the villin headpiece motif are most important for actin binding. The latter result challenges the charged crown hypothesis for the villin headpiece filamentous actin interaction.     

Use of MM-PB/SA in estimating the free energies of proteins: application to native, intermediates, and unfolded villin headpiece

We investigated the stability of three different ensembles of the 36-mer villin headpiece subdomain, the native, a compact folding intermediate, and the random coil. Structures were taken from a 1-micros molecular dynamics folding simulation and a 100-ns control simulation on the native structure. Our approach for each conformation is to first determine the solute internal energy from the molecular mechanics potential and then to add the change resulting from solvation (DeltaG(solv)). Explicit water was used to run the simulation, and a continuum model was used to estimate DeltaG(solv) with the finite difference Poisson-Boltzmann model accounting for the polarization part and a linearly surface area-dependent term for the non-polar part. We leave out the solute vibrational entropy from these values but demonstrate that there is no statistical difference among the native, folding intermediate, and random coil ensembles. We find the native ensemble to be approximately 26 kcal/mol more stable than the folding intermediate and approximately 39 kcal/mol more stable than the random coil ensemble. With an experimental estimate for the free energy of denaturation equal to 3 kcal/mol, we approximate the non-native degeneracy to lie between 10(16) and 10.(25) We also present a possible scheme for the mechanism of folding, first-order exponential decay of a putative transition state, with an estimate for the t(1/2) of folding of approximately 1 micros.     

Parallel tempering simulations of HP-36

We report results from all-atom Monte Carlo simulations of the 36-residue villin headpiece subdomain HP-36. Protein-solvent interactions are approximated by an implicit solvent model. The parallel tempering is used to overcome the problem of slow convergence in low-temperature protein simulations. Our results show that this technique allows one to sample native-like structures of small proteins and points out the need for improved energy functions.     

Preferred conformations of a linear RGD tripeptide.

Conformational study of RGD tripeptides in the nonhydrated and hydrated states was carried out using an empirical potential function ECEPP/3 and the hydration shell model in order to investigate preferred conformations and factors responsible for their stability. RGD tripeptides in the nonhydrated and hydrated states can be interpreted as existing as an ensemble of feasible conformations rather than as a single dominant conformation from the analysis of distributions of backbone conformations, hydrogen bonds and beta-turns. The different distributions of conformations for the neutral and zwitterionic RGD tripeptides in both states may indicate that the conformation of the RGD tripeptide is liable to depend on solvent polarity and pH values. beta-Turn populations for the neutral tripeptide in both states are reasonably consistent with NMR measurements on linear RGD-containing peptides. The degradation of RGD tripeptide seems to be attributed mainly to the hydrogen bonds between the Asp side-chain and the backbone of Asp residue or C-terminal NHMe group, rather than to the flexible backbones of Gly and Asp residues.     

Simulation of folding of a small alpha-helical protein in atomistic detail using worldwide-distributed computing

By employing thousands of PCs and new worldwide-distributed computing techniques, we have simulated in atomistic detail the folding of a fast-folding 36-residue alpha-helical protein from the villin headpiece. The total simulated time exceeds 300 micros, orders of magnitude more than previous simulations of a molecule of this size. Starting from an extended state, we obtained an ensemble of folded structures, which is on average 1.7A and 1.9A away from the native state in C(alpha) distance-based root-mean-square deviation (dRMS) and C(beta) dRMS sense, respectively. The folding mechanism of villin is most consistent with the hydrophobic collapse view of folding: the molecule collapses non-specifically very quickly ( approximately 20ns), which greatly reduces the size of the conformational space that needs to be explored in search of the native state. The conformational search in the collapsed state appears to be rate-limited by the formation of the aromatic core: in a significant fraction of our simulations, the C-terminal phenylalanine residue packs improperly with the rest of the hydrophobic core. We suggest that the breaking of this interaction may be the rate-determining step in the course of folding. On the basis of our simulations we estimate the folding rate of villin to be approximately 5micros. By analyzing the average features of the folded ensemble obtained by simulation, we see that the mean folded structure is more similar to the native fold than any individual folded structure. This finding highlights the need for simulating ensembles of molecules and averaging the results in an experiment-like fashion if meaningful comparison between simulation and experiment is to be attempted. Moreover, our results demonstrate that (1) the computational methodology exists to simulate the multi-microsecond regime using distributed computing and (2) that potential sets used to describe interatomic interactions may be sufficiently accurate to reach the folded state, at least for small proteins. We conclude with a comparison between our results and current protein-folding theory.     

Quantitative Determination of Site-Specific Conformational Distributions in an Unfolded Protein by Solid-State Nuclear Magnetic Resonance

Solid-state nuclear magnetic resonance (NMR) techniques are used to investigate the structure of the 35-residue villin headpiece subdomain (HP35) in folded, partially denatured, and fully denatured states. Experiments are carried out in frozen glycerol/water solutions, with chemical denaturation by guanidine hydrochloride (GdnHCl). Without GdnHCl, two-dimensional solid-state ¹³C NMR spectra of samples prepared with uniform ¹³C labeling of selected residues show relatively sharp cross-peaks at chemical shifts that are consistent with the known three-helix bundle structure of folded HP35. At high GdnHCl concentrations, most cross-peaks broaden and shift, qualitatively indicating disruption of the folded structure and development of static conformational disorder in the frozen denatured state. Conformational distributions at one residue in each helical segment are probed quantitatively with three solid-state NMR techniques that provide independent constraints on backbone ? and ψ torsion angles in samples with sequential pairs of carbonyl ¹³C labels. Without GdnHCl, the combined data are well fit by α-helical conformations. At [GdnHCl] = 4.5 M, corresponding to the approximate denaturation midpoint, the combined data are well fit by a combination of α-helical and partially extended conformations at each site, but with a site-dependent population ratio. At [GdnHCl] = 7.0 M, corresponding to the fully denatured state, the combined data are well fit by a combination of partially extended and polyproline II conformations, again with a site-dependent population ratio. Two entirely different models for conformational distributions lead to nearly the same best-fit distributions, demonstrating the robustness of these conclusions. This work represents the first quantitative investigation of site-specific conformational distributions in partially folded and unfolded states of a protein by solid-state NMR.     

Continuum solvent studies of the stability of RNA hairpin loops and helices

We apply continuum solvent models to investigate the relative stability of various conformational forms for two RNA sequences, GGAC(UUCG)GUCC and GGUG(UGAA)CACC. In the first part, we compare alternate hairpin conformations to explore the reliability of these models to discriminate between different local conformations. A second part looks at the hairpin-duplex conversion for the UUCG sequence, identifying major contributors to the thermodynamics of a much large scale transition. Structures were taken as snapshots from multi-nanosecond molecular dynamics simulations computed in a consistent fashion using explicit solvent and with long-range electrostatics accounted for using the Particle-Mesh Ewald procedure. The electrostatic contribution to solvation energies were computed using both a finite-difference Poisson-Boltzmann (PB) model and a pairwise Generalized Born model; non-electrostatic contributions were estimated with a surface-area dependent term. To these solvation free energies were added the mean solute internal energies (determined from a molecular mechanics potential) and estimates of the solute entropy (from a harmonic analysis). Consistent with experiment and with earlier solvated molecular dynamics simulations, the UUCG hairpin was found to prefer conformers close to a recent NMR structure determination in preference to those from an earlier NMR study. Similarly, results for the UGAA hairpin favored an NMR-derived structure over that to be expected for a generic GNRA hairpin loop. Experimental free energies are not known for the hairpin/duplex conversion, but must be close to zero since hairpins are seen in solution and duplexes in crystals; out calculations find a value near zero and illustrate the expected interplay of solvation, salt effects and entropy in affecting this equilibrium.     

Determination of conformational equilibrium of peptides in solution by NMR spectroscopy and theoretical conformational analysis: application to the calibration of mean-field solvation models.

Peptides occur in solution as ensembles of conformations rather than in a fixed conformation. The existing energy functions are usually inadequate to predict the conformational equilibrium in solution, because of failure to account properly for solvation, if the solvent is not considered explicitly (which is usually prohibitively expensive). NMR data are therefore widely incorporated into theoretical conformational analysis. Because of conformational flexibility, restrained molecular dynamics (with restraints derived from NMR data), which is usually applied to determine protein conformation is of limited use in the case of peptides. Instead, (a) the restraints are averaged within predefined time windows during molecular dynamics (MD) simulations (time averaging), (b) multiple-copy MD simulations are carried out and the restraints are averaged over the copies (ensemble averaging), or (c) a representative ensemble of sterically feasible conformations is generated and the weights of the conformations are then fitted so that the computed average observables match the experimental data (weight fitting). All these approaches are briefly discussed in this article. If an adequate force field is used, conformations with large statistical weights obtained from the weight-fitting procedure should also have low energies, which can be implemented in force field calibration. Such a procedure is particularly attractive regarding the parameterization of the solvation energy in nonaqueous solvents, e.g., dimethyl sulfoxide, for which thermodynamic solvation data are scarce. A method for calibration of solvation parameters in dimethyl sulfoxide, which is based on this principle was recently proposed by C. Baysal and H. Meirovitch (Journal of the American Chemical Society, 1998, Vol. 120, pp. 800--812), in which the energy gap between the conformations compatible with NMR data and the alternative conformations is maximized. In this work we propose an alternative method based on the principle that the best-fitting statistical weights of conformations should match the Boltzmann weights computed with the force field applied. Preliminary results obtained using three test peptides of varying conformational mobility: H-Ser(1)-Pro(2)-Lys(3)-Leu(4)-OH, Ac-Tyr(1)-D-Phe(2)-Ser(3)-Pro(4)-Lys(5)-Leu(6)-NH(2), and cyclo(Tyr(1)-D-Phe(2)-Ser(3)-Pro(4)-Lys(5)-Leu(6)) are presented.