Refinement of the thrombin-bound structure of a hirudin peptide by a restrained electrostatically driven Monte Carlo method.

Energy refinement of the structure of a linear peptide, hirudin56-65, bound to thrombin was carried out using a conformational search method in combination with restrained minimization. Five conformations originated from nmr data and distance geometry calculations having a similar global folding pattern but quite different backbone conformations were used as the starting structures. As a result of this approach, a series of low-energy conformations compatible with a set of upper and lower bounds of interproton distances determined from transferred nuclear Overhauser effects were found. A comparison among the lowest energy conformations of each run showed that the combination of energy refinement plus distance constraints led to a very well-defined structure for both the backbone and the side chains of the last 7 residues of the polypeptide. Furthermore, the low-energy conformations generated with this technique contain a segment of 3(10)-helix involving the last 5 residues at the COOH terminal end.     

Conformational analysis of CCK-B agonists using 1H-nmr and restrained molecular dynamics: Comparison of biologically active Boc-Trp-(N-Me)Nle-Asp-Phe-NH2 and inactive Boc-Trp-(N-Me)Phe-Asp-Phe-NH2

The tetrapeptide Boc-Trp-(N-Me)Nle-Asp-Phe-NH₂ is a potent CCK-B agonist. Replacement in this analogue of the norleucine residue by a phenylalanine, to yield Boc-Trp-(N-Me)Phe-Asp-Phe-NH₂, led to a 740-fold decrease in affinity whereas the same decrease in affinity was not observed in their nonmethylated counterparts. In order to ascertain the conformational preferences of these two N-methylated tetrapeptides, a study by two-dimensional (2D) nmr spectroscopy and molecular modeling was undertaken. The solution conformation of the two peptides was examined by ¹H-nmr in a d₆-DMSO/H₂O (80 : 20) mixture. A cis-trans equilibrium, induced by N-methylation, was observed for both analogues, and the proton spectra of the two retamers were fully characterized in each case. ¹H-¹H distance constraints, derived from 2D nuclear Overhauser effect spectroscopy and rotating frame nuclear Overhauser effect spectroscopy experiments, were used as inputs for subsequent restrained molecular dynamics simulations. Comparisons of the nmr and molecular modeling data point toward distinct conformational preferences for these two peptides with an opposite spatial orientation of the Trp residue, and could explain the large difference in their biological activities. Furthermore, the tridimensional structure of Boc-Trp-(N-Me)Nle-Asp-Phe-NH₂ could serve as a model for the design of nonpeptide CCK-B agonists. © 1994 John Wiley & Sons, Inc.     

Analysis of side-chain conformational distributions in neutrophil peptide-5 NMR structures

The side-chain conformations have been analyzed in the antimicrobial peptide, Neutrophil Peptide-5 (NP-5), whose structure was independently generated from nmr-derived distance constraints using a distance geometry algorithm. The side-chain and peptide dihedral angle distributions in the nmr structures were compared with those constructed from a data base of high-resolution protein crystal structures. The side-chain conformational preferences for NP-5 in solution are significantly different from those observed in the crystal structure data base. These results indicate that the side-chain conformations are quite disordered for many of the residues of NP-5. The absence of a correlation between the width of the conformational distribution and surface accessibility suggests that the disorder may be due to limitations in the structural information extracted from the nmr data rather than to molecular motion. However, it is also observed that the degree of conformational disorder is only weakly correlated with the number of nuclear Overhauser enhancements to a given side chain. Possible reasons for this are discussed. Molecular mechanics refinement of these structures did not significantly change the side-chain populations. Anomolously wide distributions are observed for rotations about the peptide bonds and the disulfide bonds in the NP-5 distance geometry structures, which are improved by the refinement. The very high degree of order observed for the central dihedral angle of the disulfide bond in the high-resolution crystal data base suggests that the rotation about this bond in proteins is determined by the local potential.     

Synthesis and conformational studies by 1H- and 13C-NMR spectroscopy of a novel, sterically constrained analogue of thyrotropin-releasing hormone

A novel thyrotropin-releasing hormone (TRH) analogue, [2,4-MePro3]-TRH (2,4-MePro: 2-carboxy-2,4-methanopyrrolidine), has been synthesized using a rapid solid phase peptide synthesis method, and its conformational properties investigated by one- and two-dimensional (2D) nmr spectroscopy and by proton Overhauser measurements. Following a published approach, calibrated interproton Overhauser effects were used together with distance geometry analysis to deduce that the single conformation of the His-2,4-MePro tertiary amide bond is trans in aqueous solution. This conclusion was corroborated by 2D dipolar-correlated (NOESY) spectroscopy. A preferentially extended conformation is indicated by the nmr data, similar to that of TRH. The phi, psi conformational space of 2,4-MePro is, however, significantly different from that of trans proline and the structural consequences of these differences at the C-terminus are discussed. The distribution of histidine side-chain conformations in the TRH analogue was deduced from coupling constants and from the short-range interaction between the imidazole ring and one of the prochiral faces of the 2,4-MePro side chain.     

Analysis of the 1H-NMR spectra of ferrichrome peptides. II. The amide resonances

The high-resolution ¹H-nmr study of the ferrichrome cyclohexapeptides, in d6 -DMSO solutions, has been extended to the amide NH spectral region. A total of ten diamagnetic analogues of ferrichrome that differ in the coordinated metal ion (Al³⁺, Ga³⁺ or Co³⁺), the primary structure, the nature of the bidentate hydroxamate moiety, or the isotope compositions (¹⁴N, ¹⁵N) have been investigated. The ³J_αNH values reflect regiorous conformational isomorphism throughout the complete suite of analogues, quite independent of the residue occupancy of each site. Totally resolved amide multiplets have been obtained in most cases and the four-line (doublet of doublets) appearances of glycyl NH resonances has been observed for the first times; these data enabled stereospecific assignment and accurate analysis of the NH-CαH proton spin systems. The high resolution was made possible by the use of a suitable spectral deconvolution shceme at 360 MHz. The fine structure, extraordinarily well displayed in the ¹⁵N-peptide spectra, provides a series of parameter values whose consistency has been checked by computer simulation. Since the crystallographic structure for two of the ferric peptides is known to 0.002-Å resolution, a ³J vs θ correspondence could be confidently established. A Karplus curve was derived from the combined x-ray and nmr data: It is suggested that seriously nonplanar amides can exhibit ³J_αNH values higher than predicted by the ferrichrome curve.     

NMR three-dimensional solution structure of the serine protease inhibitor cyclotheonamide A

The nmr solution conformation of cyclotheonamide A (CtA) was determined in aqueous media. The data produced 15 distance and 10 torsional constraints which were used to generate conformations using restrained simulated annealing (SA) and distance geometry/simulated annealing (DG/SA) calculations. Two different calculation protocols were performed to ensure proper sampling of conformational space and even though the torsional restraints were input differently, both calculation methods yielded the same conformation of CtA. In the structure calculations, all solutions of the Karplus equation were sampled simultaneously using the restrained SA protocol and large ranges were used for the dihedral restraints in the DG/SA protocol because all solutions to the Karplus equation could not be sampled simultaneously. The solution conformation was also compared to the solid state x-ray conformations of CtA bound to thrombin and trypsin. The conformation of the residues important for active site binding (d -Phe, h-Arg, and Pro) are nearly identical in aqueous solution and solid state with largest differences at the a-Ala and v-Tyr residues. CtA appears to be preordered in structure and does not undergo a significant conformational change upon binding to the enzyme active site. © 1997 John Wiley & Sons, Inc.     

Conformation of CD4-derived cyclic hexapeptides by nmr and molecular dynamics

Two cyclic hexapeptides, cyclo[Ala¹-D -Ala²-Ser³-Phe⁴-Gly⁵-Ser⁶] and cyclo[Ala¹-Gly²-Ser³-Phe⁴-Gly⁵-Ser⁶], derived from the loop portion of the C′C″ ridge of CD4, were characterized by high-resolution nmr spectroscopy and simulated annealing studies. In DMSO-d₆ both of these peptides display a single conformer on the nmr time scale with two intramolecular H-bond (1 ← 4) stabilized β-turns at positions 2–3 and 5–6. The nmr derived distance constraints were used in simulated annealing calculations to generate the solution structures. These structures adopt energetically comparable conformational substates that are not resolvable on the nmr time scale. In aqueous solution, the H-bond stabilized β-turn conformation for cyclo [Ala-D -Ala-Ser-Phe-Gly-Ser] is no longer the predominant structural form. Structures generated using molecular dynamics simulations with no experimental constraints were compared with those from nmr analysis. The correlation between these two sets of structures allows the use of molecular simulations as a predictive tool for the conformational analysis of small peptides. © 1994 John Wiley & Sons, Inc.     

Accessible conformations of melanin-concentrating hormone: a molecular dynamics approach

Molecular dynamics simulations have been used to search for the accessible conformations of the melanin-concentrating hormone (MCH). The studies have been performed on native MCH and two of its peptide fragments, a cyclic MCH(5-14) fragment and a linear MCH(5-14) fragment. An analysis of the molecular dynamics trajectories of the three peptides indicates that two regions of the peptide have characteristic conformational properties that may be important for the biological activity. One is a region around Gly8, which is conformationally mobile, and the other is around Pro13, which shows unusual rigidity. The molecular dynamics simulation results are discussed in terms of backbone structural features like beta turns, side-chain interactions, and orientations of the disulfide bridge. The results of this analysis are used to suggest new analogues that will modify the conformational features of the peptide and further define the conformational requirements for activity. Finally, the results are related to nmr studies of the peptide and reveal agreements between the experimental nuclear Overhauser effect constraints and some of the accessible conformations obtained from the simulation.     

Comparison of helix-stabilizing effects of alpha,alpha-dialkyl glycines with linear and cycloalkyl side chains

The ability of alpha, alpha-di-n-alkyl glycines with linear and cyclic alkyl side chains to stabilize helical conformations has been compared using a model heptapeptide sequence. The conformations of five synthetic heptapeptides (Boc-Val-Ala-Leu-Xxx-Val-Ala-Leu-OMe, Xxx = Ac8c, Ac7c, Aib, Dpg, and Deg, where Ac8c = 1-aminocyclooctane-1-carboxylic acid, Ac7c = 1-aminocycloheptane-1-carboxylic acid, Aib = alpha-aminoisobutyric acid, Dpg = alpha,alpha-di-n-propyl glycine, Deg = alpha,alpha-di-n-ethyl glycine) have been investigated. In crystals, helical conformations have been demonstrated by x-ray crystallography for the peptides, R-Val-Ala-Leu-Dpg-Val-Ala-Leu-OMe, (R = Boc and acetyl). Solution conformations of the five peptides have been studied by 1H-nmr. In the apolar solvent CDCl3, all five peptides favor helical conformations in which the NH groups of residues 3-7 are shielded from the solvent. Successive NiH<-->Ni + 1H nuclear Overhauser effects over the length of the sequence support a major population of continuous helical conformations. Solvent titration experiments in mixtures of CDCl3/DMSO provide evidence for solvent-dependent conformational transitions that are more pronounced for the Deg and Dpg peptides. Solvent-dependent chemical shift variations and temperature coefficients in DMSO suggest that the conformational distributions in the Deg/Dpg peptides are distinctly different from the Aib/Acnc peptides in a strongly solvating medium. Nuclear Overhauser effects provide additional evidence for the population of extended backbone conformations in the Dpg peptide, while a significant residual population of helical conformations is still detectable in the isomeric Ac7c peptide in DMSO.     

Optimization of the binding properties of a synthetic anion receptor using rational and combinatorial strategies

The solution conformations of tetrameric and hexameric cyclopeptides containing alternating L-proline and 6-aminopicolinic acid subunits strongly depend on solvent polarity. Whereas in polar solvents, such as d6-DMSO, both peptides prefer on average symmetric conformations with converging NH groups, in less polar chloroform intramolecular hydrogen bonds to the peptide NH groups stabilize other, and in the case of the hexapeptide, non-symmetrical conformations. Independent of the solvent, both peptides interact with anions via their NH groups but whereas anion binding requires a cleavage of the intramolecular hydrogen bonds accompanied by a conformational reorganization in chloroform, in polar solvents the peptides are already well preorganized for anion complexation. Complex formation between anions and the cyclic hexapeptide was even detected in highly competitive D2O/CD3OD or H2O/CH3CN mixtures, which was attributed to the special sandwich-type structure of the complexes formed. Stabilizing these 2:1 aggregates by covalently linking two cyclopeptide rings together affords ditopic receptors with a high anion affinity in protic solvents. Complex stability depends on the structure of the linker with which the two receptor moieties are connected and even more potent anion receptors were obtained by a dynamic combinatorial optimization of this linking unit.     

Solution conformational study of Scyliorhinin I analogues with conformational constraints by two-dimensional NMR and theoretical conformational analysis.

Two analogues of Scyliorhinin I (Scyl), a tachykinin with N-MeLeu in position 8 and a 1,5-disubstituted tetrazole ring between positions 7 and 8, introduced in order to generate local conformational constraints, were synthesized using the solid-phase method. Conformational studies in water and DMSO-d6 were performed on these peptides using a combination of the two-dimensional NMR technique and theoretical conformational analysis. The algorithm of conformational search consisted of the following three stages: (i) extensive global conformational analysis in order to find all low-energy conformations; (ii) calculation of the NOE effects and vicinal coupling constants for each of the low energy conformations; (iii) determining the statistical weights of these conformations by means of a nonlinear least-squares procedure, in order to obtain the best fit of the averaged simulated spectrum to the experimental one. In both solvents the three-dimensional structure of the analogues studied can be interpreted only in terms of an ensemble of multiple conformations. For [MeLeu8]Scyl, the C-terminal 6-10 fragment adopts more rigid structure than the N-terminal one. In the case of the analogue with the tetrazole ring in DMSO-d6 the three-dimensional structure is characterized by two dominant conformers with similar geometry of their backbones. They superimpose especially well (RMSD = 0.28 A) in the 6-9 fragments. All conformers calculated in both solvents superimpose in their C-terminal fragments much better than those of the first analogue. The results obtained indicate that the introduction of the tetrazole ring into the Scyl molecule rigidifies its structure significantly more than that of MeLeu.     

Conformational studies on bombesin antagonists: CD and NMR characterization of [Thr6, Leu13 psi(CH2NH) Met14] bombesin (6-14).

The conformational flexibility of the [Thr6, Leu13 psi(CH2NH) Met14] bombesin (6-14) nonapeptide has been studied by CD and one- and two-dimensional (1D and 2D) nmr techniques. The CD and nmr parameters in different solvents and in a micellar environment (SDS) are compared with the data collected for the parent bombesin (BN) and [D-Phe12, Leu14]BN. A preliminary investigation on spantide is also reported. In particular, the results obtained from CD measurements indicate that there is a shift from random coil structures, in aqueous solutions, toward folded structures in apolar media (2,2,2-trifluoroethanol) and in a membrane-mimetic environment (40 mM SDS) for all three peptides, namely BN, [D-Phe12, Leu14]BN, and [Thr6, Leu13 psi(CH2NH) Met14]BN (6-14). Spantide, which also possesses some inhibitory activity against BN but very little sequence similarity, even in water, shows an ordered conformation. Nuclear magnetic resonance parameters such as backbone NH-alpha CH coupling constant values, amidic temperature coefficients, and the presence of only sequential nuclear Overhauser effects have not provided, so far, any clear evidence for a preferential ordered structure in the peptides studied, and this may be due to rapid exchange among different conformers in the nmr time scale.     

Conformational analysis of a IgG1 hinge peptide derivative in solution determined by NMR spectroscopy and refined by restrained molecular dynamics simulations.

The hinge region links the antigen binding Fab part to the constant Fc domain in immunoglobulins. For the hinge peptide derivative [AcThr(OtBu)-Cys-Pro-Pro-Cys-Pro-Ala-ProNH2]2 the assignment of the 1H and 13C resonances was achieved by two-dimensional nmr techniques: total correlation spectroscopy (TOCSY), nuclear Overhauser enhancement spectroscopy (NOESY), rotating frame nuclear Overhauser enhancement spectroscopy (ROESY), heteronuclear multiple quantum coherence (HMQC) transfer, and a HSQC (modified Overbodenhausen experiment) with high resolution in F1, which was several times folded in F1 but still phase correctable. Conformational relevant parameters (78 nuclear Overhauser effect distance restraints, 3JHH for prochiral assignments, temperature gradients) were determined by nmr and served as input data for molecular dynamics (MD) structure refinement. A simulated model compound corresponding to the [Cys-Pro-Pro-Cys]2 core elongated by the peptide chains in the Fab and Fc direction served as a starting structure for the final MD run. The conformation calculated in in vacuo does not agree with the C2 symmetry required from nmr data, but the structure obtained by a water simulation fulfills the requirement. Here the core of the hinge peptide derivative adopts a polyproline II double helix as in the x-ray structure of IgG1. Hence, segments responsible for the internal flexibility are located outside the core as confirmed by the flexibility of the solvent exposed C termini.     

Determination of interspin distances between spin labels attached to insulin: comparison of electron paramagnetic resonance data with the X-ray structure

A method was developed to determine the interspin distances of two or more nitroxide spin labels attached to specific sites in proteins. This method was applied to different conformations of spin-labeled insulins. The electron paramagnetic resonance (EPR) line broadening due to dipolar interaction is determined by fitting simulated EPR powder spectra to experimental data, measured at temperatures below 200 K to freeze the protein motion. The experimental spectra are composed of species with different relative nitroxide orientations and interspin distances because of the flexibility of the spin label side chain and the variety of conformational substates of proteins in frozen solution. Values for the average interspin distance and for the distance distribution width can be determined from the characteristics of the dipolar broadened line shape. The resulting interspin distances determined for crystallized insulins in the R6 and T6 structure agree nicely with structural data obtained by x-ray crystallography and by modeling of the spin-labeled samples. The EPR experiments reveal slight differences between crystal and frozen solution structures of the B-chain amino termini in the R6 and T6 states of hexameric insulins. The study of interspin distances between attached spin labels can be applied to obtain structural information on proteins under conditions where other methods like two-dimensional nuclear magnetic resonance spectroscopy or x-ray crystallography are not applicable.     

Combined effect of the DeltaPhe or DeltaAla residue and the p-nitroanilide group on a didehydropeptides conformation

Two series of dehydropeptides of the general formulae Boc-Gly-X-Phe-p-NA, Boc-Gly-Gly-X-Phe-p-NA, Gly-X-Gly-Phe-p-NA.TFA, and Boc-Gly-X-Gly-Phe-p-NA, with X = Delta(Z)Phe and DeltaAla, were studied with NMR in DMSO and CDCl(3)-DMSO, and with CD in MeOH, MeCN, and TFE. The NMR spectra measured in DMSO suggest that peptides with the DeltaPhe residue next to Phe are folded whereas peptides with Gly between DeltaPhe and Phe are less ordered. NMR spectra of DeltaAla-containing peptides indicate that these peptides are flexible and their conformational equilibria are populated by many different conformations. The CD spectra show that conformational properties of the peptides studied are distinctly influenced by a mutual position of the dehydroamino acid residue and the p-NA group. They indicate that all dehydropeptides with the DeltaPhe residue, Boc-Gly-DeltaAla-Phe-p-NA, and Boc-Gly-Gly-DeltaAla-Phe-p-NA adopt ordered conformations in all solvents studied, presumably of the beta-turn type. The last two peptides exhibit surprising chiroptical properties. Their spectra show exciton coupling-like couplets in the region of the p-NA group absorption. This shape of CD spectra suggests a rigid, chiral conformation with a fixed disposition of the p-NA group. The CD spectra indicate that Boc-Gly-DeltaAla-Gly-Phe-p-NA and Gly-DeltaAla-Gly-Phe-p-NA.TFA are unordered, independently of the solvent.     

Conformational changes due to vicinal glycosylation: The branched α-l-Rhap(1–2)[β-d-Galp(1–3)]-β-d-Glc1-OMe trisaccharide compared with its parent disaccharides*

Conformations of the α-l -Rhap(1-2)-β-d -Glc1-OMe and β-d -Galp(1-3)-β-d -Glc1-OMe disaccharides and the branched title trisaccharide were examined in DMSO-d₆ solution by ¹H-nmr. The distance mapping procedure was based on rotating frame nuclear Overhauser effect (NOE) constraints involving C- and O-linked protons, and hydrogen-bond constraints manifested by the splitting of the OH nmr signals for partially deuteriated samples. An “isotopomer-selected NOE” method for the unequivocal identification of mutually hydrogen-bonded hydroxyl groups was suggested. The length of hydrogen bonds thus detected is considered the only one motionally nonaveraged nmr-derived constraint. Molecular mechanics and molecular dynamics methods were used to model the conformational properties of the studied oligosaccharides. Complex conformational search, relying on a regular Φ,Ψ-grid based scanning of the conformational space of the selected glycosidic linkage, combined with simultaneous modeling of different allowed orientations of the pendant groups and the third, neighboring sugar residue, has been carried out. Energy minimizations were performed for each member of the Φ,Ψ grid generated set of conformations. Conformational clustering has been done to group the minimized conformations into families with similar values of glycosidic torsion angles. Several stable syn and anti conformations were found for the 1→2 and 1→3 bonds in the studied disaccharides. Vicinal glycosylation affected strongly the occupancy of conformational states in both branches of the title trisaccharide. The preferred conformational family of the trisaccharide (with average Φ,Ψ values of 38°, 17° for the 1→2 and 48°, 1° for the 1→3 bond, respectively) was shown by nmr to be stabilized by intramolecular hydrogen bonding between the nonbonded Rha and Gal residues. © 1998 John Wiley & Sons, Inc. Biopoly 46: 417–432, 1998     

Conformational analysis of cyclic peptides in solution

The strategy and tactics of conformational analysis of cyclic peptides in solution is demonstrated by the example of cyclo(-D-Pro-Phe-Thr-Phe-Trp-Phe-). Spin-locked experiments like rotating frame nuclear Overhauser enhancement spectroscopy (ROESY), ROTO, and TOCSY are successfully applied to assign all proton signals and to obtain distance information. A crude conformational model was built using the nmr data. This starting model was refined by restrained molecular dynamics (MD) calculations using ROE derived distances and fixed bond angles as determined from homo- and heteronuclear coupling constants. To mimic the solvent and to reduce artifacts in an in vacuo calculation the charges of the solvent-exposed NH protons were gradually reduced according to the temperature gradients. The thus obtained "conformation" (mean of a 40 ps MD trajectory) shows very close similarity to x-ray structures in an orthorhombic and in two monoclinic crystal modifications of the same compound. The main difference is the breaking of an intermolecular hydrogen bond of the threonine hydroxyl group on dissolution of the crystal and forming an intramolecular hydrogen bond in solution.     

Solution structure of a LewisX analogue by off-resonance 1H NMR spectroscopy without use of an internal distance reference

Summary A recent ¹H NMR method has been applied to the determination of the solution structure and internal dynamics of a synthetic mixed C/O trisaccharide related to sialyl Lewis^x. Varying the rf field offset in ROESY-type experiments enabled the measurement of longitudinal and transverse dipolar cross-relaxation rates with high accuracy. Assuming that for each proton pair the motion could be represented by a single exponential autocorrelation function, it was possible to derive geometrical parameters (r) and dynamic parameters _cp. With this assumption, 224 cross-relaxation rates have been transformed into 30 interproton distance constraints and 30 dipolar correlation times. The distance constraints have been used in a simulated-annealing procedure. This trisaccharide exhibits a structure close to the O-glycosidic analogue, but its flexibility seems highly reduced. On the basis of the determined structure and dynamics, it is shown that no conformational exchange occurs, the molecule existing in the form of a unique family in aqueous solution. In order to assess the quality of the resulting structures and to validate this new experimental procedure of distance extraction, we finally compare these solution structures to the ones obtained using three different sets of distances deduced from three choices of internal reference. It appears that this procedure allows the determination of the most precise and accurate solution.Abbreviations COSY correlation spectroscopy - NOE nuclear Overhauser enhancement - NOESY nuclear Overhauser enhancement spectroscopy; rmsd, root-mean-square deviation - ROESY rotating frame Overhauser enhancement spectroscopy - SLe^x sialyl Lewis^x - TOCSY total correlation spectroscopy     

Precise structural determination of weakly binding peptides by utilizing dihedral angle constraints

Structural determination of target-bound conformations of peptides is of primary importance for the optimization of peptide ligands and peptide–mimetic design. In the structural determination of weakly binding ligands, transferred nuclear Overhauser effect (TrNOE) methods have been widely used. However, not many distance constraints can be obtained from small peptide ligands by TrNOE, especially for peptides bound to a target molecule in an extended conformation. Therefore, for precise structural determination of weakly binding peptides, additional structural constraints are required. Here, we present a strategy to systematically introduce dihedral angle constraints obtained from multiple transferred cross-correlated relaxation experiments and demonstrate precise structures of weakly binding peptides. As a result, we could determine the bioactive conformations of phage-derived peptide ligands and define their core binding motifs.     

Necessary conditions for avoiding incorrect polypeptide folds in conformational search by energy minimization

Low energy conformations have been generated for melittin, pancreatic polypeptide, and ribonuclease S-peptide, both in the vicinity of x-ray structures by energy refinement and by an unconstrained search over the entire conformational space. Since the structural polymorphism of these medium-sized peptides in crystal and solution is moderate, comparing the calculated conformations to x-ray and nmr data provides information on local and global behavior of potential functions. Local analysis includes standardization calculations, which show that models with standard geometry can approximate good resolution x-ray data with less than 0.5 Å rms deviation (RMSD). However, the atomic coordinates are shifted up to 2 Å RMSD by local energy minimization, and thus 2 Å is generally the smallest RMSD value one can target in a conformational search using the same energy evaluation models. The unconstrained search was performed by a buildup-type method based on dynamic programming. To accelerate the generation of structures in the conformational search, we used the ECEPP potential, defined in terms of standard polypeptide geometry. A number of low energy conformations were further refined by relaxing the assumption of standard bond lengths and bond angles through the use of the CHARMM potential, and the hydrophobic folding energies of Eisenberg and McLachlan were calculated. Each conformation is described in terms of the RMSD from the native, hydrogen-bonding structure, solvent-acessible surface area, and the ratio of surfaces corresponding to nonpolar and polar residues. The unconstrained search finds conformations that are different from the native, sometimes substantially, and in addition, have lower conformational energies than the native. The origin of deviations is different for each of the three peptides, but in all examples the refined x-ray structures have lower energies than the calculated incorrect folds when (1) the assumption of standard bond lengths and bond angles is relaxed; (2) a small and constant effective dielectric permittivity (ε < 10) is used; and (3) the hydrophobic folding energy is incorporated into the potential. © 1993 John Wiley & Sons, Inc.