Quantitative two-dimensional NMR study of dermenkephalin, a highly potent and selective delta-opioid peptide

Toxins in anti-nociception and anti-inflammation

The use of toxins as novel molecular probes to study the structure-function relationship of ion-channels and receptors as well as potential therapeutics in the treatment of wide variety of diseases is well documented. The high specificity and selectivity of these toxins have attracted a great deal of interest as candidates for drug development. This review highlights the involvement of the proteins and peptide toxins as well as non-proteinaceous compounds derived from both venomous and non-venomous animals, in anti-nociception and anti-inflammation. The possible mechanisms of these potential therapeutic agents and possible clinical applications in the treatment of pain and inflammation are also summarized.

The mu-selective heptapeptide opioid dermorphin has two conformations around Phe3 psi with no head-to-tail interaction. A quantitative 2-D NMR and molecular modeling analysis

The mu opioid heptapeptide Dermorphin (DRM) is under 70 % of trans forms for the Tyr(5)-Pro(6) peptide bond in solution (CDCl(3)/DMSO-d(6) 1/1 vol/vol). Variations of NOE integrals at 5 temperatures show apparent correlation times of 0.8 to 0.9 ns (at 280 K) in that mixed solvent. Four NOE between non-adjacent residues reveal a large population of folded structures. However, in trans DRM, 4 adjacent NOE Phe(3)/Gly(4) can only be explained by an equilibrium between folded (psi(3) > 0) and extended (psi(3) > 0) conformations. Simulated annealing modeling gave about 60% (psi(3) > 0) and 40% (psi(3) > 0) of these conformer populations. Trans DRM study and previous studies on the heptapeptide opioids, dermenkephalin (DREK) and deltorphin-I (delta selective), and DREK(1-4)-DRM(5-7) hybrid (mu selective), show in folded structures more backbone bending of the first 4 residues in the mu opioids than in the delta peptides. Also, the main difference between mu- and delta-opioid peptides is a large fraction of extended conformations in mu heptapeptides. Either bending of the N-terminus, or extension of the C-terminal part in mu-opioid heptapeptides prevent the head-to-tail interactions which allow delta-opioid peptides to bind selectively to the delta-opioid receptor.

The delta-selective opioid peptide dermenkephalin and the mu-selective hybrid peptide dermenkephalin-[1-4]-dermophin-[5-7] display strikingly different conformations despite identical tetrapeptide N-termini. A quantitative 2-D NMR and molecular modeling analysis

The selective recognition of the aminoterminal binding pharmacophore Tyr-D-Xaa-Phe of the opioid heptapeptide dermorphin, Tyr-D-Ala-Phe-Gly-Tyr-Pro-Ser-NH2 (DRM)1, and of dermenkephalin, Tyr-D-Met-Phe-His-Leu-Met-Asp-NH2 (DREK), by the mu-opioid receptor and delta-opioid receptor, respectively, depends upon the constitution / conformation of the C-terminal tripeptide. The hybrid peptide DREK-[1-4]-DRM-[5-7] is very potent at, and exquisitely selective for the mu-opioid receptor, and differs only from dermenkephalin by its C-terminal tripeptide. Comparison of the structural features of DREK-[1-4]-DRM-[5-7] and dermenkephalin by nmr analysis and molecular modeling revealed striking differences, as well in the trans (Tyr5 - Pro6) isomer (population 75%) than in the cis isomer.. Whereas the folded C-terminal tail of dermenkephalin influenced the tertiary structure of the N-terminal tetrapeptide and placed the Tyr1 and Phe3 aromatic rings in definite orientations that are best suited for the delta-receptor, there were only weak contacts, as shown by NOE data, between the aminoterminal and carboxyterminal parts of the hybrid peptide. This promoted increased flexibility of the whole backbone and relaxed orientations for the side-chains of Tyr1 and Phe3 that are compatible with the mu-receptor but unsuitable for the delta-receptor. The steric hindrance introduced by Pro6 in DREK-[1-4]-DRM-[5-7], plus the absence of large hydrophobic side-chains in positions 5 and 6 may prevent close contacts between the N-terminal and C-terminal domains and reorientation of the main pharmacophoric elements Tyr1 and Phe3.

Opioid peptides from frog skin

The skin of the South American frogs Phyllomedusa secretes, in addition to numerous mammalian-like hormones and neuropeptides, several gene-encoded opioid peptides that contain a D-amino acid in position 2 of their sequence. Dermorphin, Tyr-D-Ala-Phe-Gly-Tyr-Pro-Ser-NH2, dermenkephalin/deltorphin A, Tyr-D-Met-Phe-His-Leu-Met-Asp-NH2 and the deltorphins, Tyr-D-Ala-Phe-Xaa-Val-Val-Gly-NH2 (where Xaa is either Asp or Glu) are highly potent at, and exquisitely selective, for the mu- and delta-opioid receptors. D-Ala and D-Met present in dermorphin and related peptides are coded for by the usual codons in the corresponding messenger RNAs. Prepro-dermorphin/dermenkephalin and prepro-deltorphins have considerable sequence identities to precursors encoding 10-46-residue-long antimicrobial peptides--dermaseptins, brevinins, temporins, esculentins and gaegurins--originating from various amphibian species. The similarity between the prepro-regions of precursors encoding end products with strikingly different structures and biological activities supports the suggestion that the genes encoding these peptides are all members of the same family.

What peptides these deltorphins be

The deltorphins are a class of highly selective delta-opioid heptapeptides from the skin of the Amazonian frogs Phyllomedusa sauvagei and P. bicolor. The first of these fascinating peptides came to light in 1987 by cloning of the cDNA of from frog skins, while the other members of this family were identified either by cDNA or isolation of the peptides. The distinctive feature of deltorphins is the presence of a naturally occurring D-enantiomer at the second position in their common N-terminal sequence, Tyr-D-Xaa-Phe, comparable to dermorphin, which is the prototype of a group of mu-selective opioids from the same source. The D-amino acid and the anionic residues, either Glu or Asp, as well as their unique amino acid compositions are responsible for the remarkable biostability, high delta-receptor affinity, bioactivity and peptide conformation. This review summarizes a decade of research from many laboratories that defined which residues and substituents in the deltorphins interact with the delta-receptor and characterized pharmacological and physiological activities in vitro and in vivo. It begins with a historical description of the topic and presents general schema for the synthesis of peptide analogues of deltorphins A, B and C as a means to document the methods employed in producing a myriad of analogues. Structure activity studies of the peptides and their pharmacological activities in vitro are detailed in abundantly tabulated data. A brief compendium of the current level of knowledge of the delta-receptor assists the reader to appreciate the rationale for the design of these analogues. Discussion of the conformation of these peptides addresses how structure leads to further hypotheses regarding ligand receptor interaction. The review ends with a broad discussion of the potential applications of these peptides in clinical and therapeutic settings.

Solution conformations of deltorphin-I obtained from combined use of quantitative 2D-NMR and energy calculations: a comparison with dermenkephalin

Deltorphin-I, Tyr-D-Ala-Phe-Asp-Val-Val-Gly-NH2 and dermenkephalin, Tyr-D-Met-Phe-His-Leu-Met-Asp-NH2, two highly related opioid peptides from frog skin, display very similar N-termini but strikingly different C-terminal tails. Nevertheless, both peptides are highly potent at, and exquisitely selective for the delta-opioid receptor. To identify common determinants concuring to the remarkably efficient targeting of deltorphin-I and dermenkephalin, combined use of quantitative two-dimensional nuclear magnetic resonance (53 dipolar interactions studied at four temperatures) and energy calculations using simulated annealing generated five groups of deltorphin-I conformers. These groups were pooled into two families whose overall conformation could be described either by a left-handed helix (Family I) or by a big loop (Family II), both stabilized by H-bonds. Proximity of D-Ala2-Phe3-Asp4 and Val5-Val6-Gly7 triads is an obvious structural similarity between almost all groups in both families of structures. Whereas differences between the two families originated mostly from a transition at psi Asp4 backbone dihedral angle, the backbone structures at segment 1-4 are similar and spatial arrangements of Tyr1 (t) and Phe3 (g-) are identical in one group of each family. Moreover, these two groups have a N-terminal tetrapeptide whose conformation most closely resembles that of a well-defined group of structures for dermenkephalin. Altogether, these results suggest that conformational attributes that are common to dermenkephalin and deltorphin-I, i.e., the backbone conformation of the N-terminal tetrapeptide and preferential orientations in the side-chain of Tyr1 (t) and Phe3 (g-) underlie their ability to bind with high selectivity to the delta-opioid receptor.

Folded conformations of the delta-selective opioid dermenkephalin with head-to-tail interactions. A simulated annealing study through NMR restraints

Despite similar tripeptide N-termini, dermorphin (Tyr-D-Ala-Phe-Gly-Tyr-Pro-Ser-NH2) and dermenkephalin (Tyr-D-Met-Phe-His-Leu-Met-Asp-NH2), naturally occuring opioid peptides from frog skin, exhibit high affinity but contrasting selectivity for the mu- and delta-opioid receptors, respectively. Structure-activity relationship studies have shown that the N-terminal tripeptide, Tyr-D-Xaa-Phe (where Xaa is either Ala or Met), is necessary for binding with both the mu- and delta-receptors while the nature and/or the conformation of the C-terminus His-Leu-Met-Asp-NH2 of dermenkephalin are responsible for addressing the peptide to the delta-receptor. In order to examine the conformational characteristics that are related to the selectivity of dermenkephalin towards the delta-receptor, 50 NOE restraints (10 between non-adjacent residues), and 7 dihedral angles, derived from a two-dimensional 1H-NMR study of dermenkephalin in dimethyl sulfoxide, were used in simulated annealing and energy minimization procedures. Twenty-four resulting conformers (60% of the generated structures) with no severe distance restraint violation were pooled into seven groups and three related families. These 24 conformers show close proximity between the two methionine residues, S-shaped structures, mean planes of N-terminal and C-terminal moieties almost at right angles to each other, a C-terminus region above the plane of the N-terminal region and g- as preferential orientation in the side chain of Phe. Aside these similarities, families of conformers differ by the preferential orientation in the side chain of Tyr (t or g-) and proximity between Tyr and Asp, or Tyr and the C-terminus. In contrast to previous models, practically no beta-turn structures exist for dermenkephalin, most of the NH hydrogen bonds participating to gamma-turns. The possible relationship between the conformational characteristics of dermenkephalin and the delta-opioid receptor selectivity is discussed.

Development of a model for the delta-opioid receptor pharmacophore: 3. Comparison of the cyclic tetrapeptide, Tyr-c[D-Cys-Phe-D-Pen]OH with other conformationally constrained delta-receptor selective ligands

We have previously proposed a model of the delta-opioid receptor bound conformation for the cyclic tetrapeptide, Tyr-c[D-Cys-Phe-D-Pen]OH (JOM-13) based on its conformational analysis and from conformation-affinity relationships observed for its analogues with modified first and third residues. To further verify the model, it is compared here with results of conformational and structure-activity studies for other known conformationally constrained delta-selective ligands: the cyclic pentapeptide agonist, Tyr-c[D-Pen-Gly-Phe-D-Phe]OH (DPDPE): the peptide antagonist, Tyr-Tic-Phe-PheOH (TIPP); the alkaloid agonist, 7-spiroindanyloxymorphone (SIOM); and the related alkaloid antagonist, oxymorphindole (OMI). A candidate delta-bound conformer is identified for DPDPE that provides spatial overlap of the functionally important N-terminal NH3+ and C-terminal COO- groups and the aromatic rings of the Tyr and Phe residues in both cyclic peptides. It is shown that all delta-selective ligands considered have similar arrangements of their pharmacophoric elements, i.e., the tyramine moiety and a second aromatic ring (i.e., the rings of Phe3, Phe4, and Tic2 residues in JOM-13, DPDPE, and TIPP, respectively; the indole ring system in OMI, and the indanyl ring system in SIOM). The second aromatic rings, while occupying similar regions of space throughout the analogues considered, have different orientations in agonists and antagonists, but identical orientations in peptide and alkaloid ligands with the same agonistic or antagonistic properties. These results agree with the previously proposed binding model for JOM-13, are consistent with the view that delta-opioid agonists and antagonists share the same binding site, and support the hypothesis of a similar mode of binding for opioid peptides and alkaloids.

New theoretical methodology for elucidating the solution structure of peptides from NMR data. II. Free energy of dominant microstates of Leu-enkephalin and population-weighted average nuclear Overhauser effects intensities

A small linear peptide in solution may populate several stable states (called here microstates) in thermodynamic equilibrium; elucidating its dynamic three dimensional structure by multi- dimensional nmr is complex since the experimentally measured nuclear Overhauser effect intensities (NOEs) represent averages over the individual contributions. We propose a new methodology based on statistical mechanical considerations for analyzing nmr data of such peptides. In a previous paper (called paper I, H. Meirovitch et al. (1995) Journal of Physical Chemistry, 99, 4847-4854] we have developed theoretical methods for determining the contribution to the partition function Z of the most stable microstates, i.e. those that pertain to a given energy range above the global energy minimum (GEM). This relatively small set of dominant microstates provides the main contribution to medium- and long-range NOE intensities. In this work the individual populations and NOEs of the dominant microstates are determined, and then weighted averages are calculated and compared with experiment. Our methodology is applied to the pentapeptide Leu-enkephalin H-Tyr-Gly-Gly-Phe-Leu-OH, described by the potential energy function ECEPP. Twenty one significantly different energy minimized structures are first identified within the range of 2 kcal/mol above the GEM by an extensive conformational search; this range has been found in paper I to contribute 0.6 of Z. These structures then become "seeds" for Monte Carlo (MC) simulations designed to keep the molecule relatively close to its seed. Indeed, the MC samples (called MC microstates) illustrate what we define as intermediate chain flexibility; some dihedral angles remain in the vicinity of their seed value, while others visit the full range of [-180 degrees, 180 degrees]. The free energies of the MC microstates (which lead to the populations) are calculated by the local states method, which (unlike other techniques) can handle any chain flexibility. The NOE of MC microstate i is calculated as the average <1/r(3)>i(2), and an effective interatomic distance ri(eff) is defined as ri(eff) = <l/r(3)>i(-1/3), where r is the distance between two protons. Under "initial rate approximation," and neglecting angular modulations, the overall I is the average over ri(eff-6), weighted by the populations of the MC microstates. This treatment is justified under the assumption that the rates at which conformations interconvert within, and among, microstates are faster and slower, respectively, than the rotational reorientation of the molecule. I(-6) leads to the virtual theoretical distances, compared to the corresponding virtual experimental distances, which were obtained previously from a cryoprotective solution of Leu-enkephalin at 280 K. A reasonable fit is found between theory and experiment. Future research directions are outlined.

NMR conformational study of a model tetradecapeptide mimicking the RXVRG consensus cleavage site of a Xenopus laevis skin endoprotease

A model tetradecapeptide, used for the purification of the RXVRG-endoprotease from Xenopus laevis skin exudate, has been studied by two-dimensional NMR, correlation (COSY) and NOE (NOESY) spectroscopy. This peptide has the 5-9 consensus sequence (RXVRG), along with an acidic moiety (1-4) and a hydrophobic domain (10-14). Variations with temperature of NH chemical shifts in a dimethyl sulfoxide solution (low thermal coefficients at residues 6, 7 and 8) and quantified NOE values from four spectra at different mixing times clearly showed a structural organization in the consensus domain with psi-angles around [-40, -10 degrees] for residues 7 and 8, and two NOE correlations of alpha HiNHi + 2 type (5-7 and 6-8). Moreover, a privileged rotamer in the side chain is established for three residues (Val2, Asp3 and Val7) and limited possibilities are discussed for seven others. Most of the folding trends were not observed in the [Ser7] derivative, underlying the relationship between the conformations and a full consensus sequence. In the model tetradecapeptide an equilibrium between two beta-turns of type I, fragments 4-7 and 5-8, seems the most probable. Comparison between this tetradecapeptide and its 4-14 fragment, also a substrate for RXVRG-endoprotease, shows that the 1-3 moiety (DVD) influences the consensus domain structure(s) and clearly stabilizes the folded one(s). Finally, two analytical methods are developed in order to determine: (1) the trifluoroacetic acid content of the peptide samples, on the basis of 19F NMR spectroscopy; (2) the mean phi- and psi-angles of each residue, from the whole set of NH/alpha H coupling constants (3JN alpha) and NOE data at a local level.

Two-dimensional 1H NMR study of a tetradecapeptide with the consensus sequence Arg5-Asp-Val-Arg-Gly9: structural effects of the outside substitution Ser12 by Ala12

Conformation of a tetradecapeptide with a RXVRG consensus sequence, Arg5-Asp-Val-Arg-Gly9, found in several precursors of antibacterian peptides, was investigated in dimethylsulfoxide solution by proton NMR spectroscopy. Complete resonance assignments and conformational parameters were obtained through correlated (COSY) and nuclear Overhauser (NOESY) techniques. The 3J(alpha H, beta H) coupling constants and the intramolecular NOE, NH...beta H, were used to analyse the conformers around the C alpha-C beta bond and, in four cases, to obtain stereospecific assignments. Use of restraints derived from NOE connectivities and 3J(NH, alpha H) coupling constants allows the determination of a range of phi and psi dihedral angles for all the residues in the sequence. The present NMR results provide favourable evidence for the formation of two bends in the consensus sequence of the tetradecapeptide. The first one has most of the features of a Glu4-Val7 beta-turn (low temperature coefficient of the Val7NH chemical shift, Arg5 alpha H...Val7NH and Asp6NH...Val7NH NOE correlations). The second one exhibits only the Asp6 alpha H...Arg7NH and Val7NH...Arg8NH NOE interactions. These consensus sequence organizations proposed were confirmed by molecular modeling based on low potential energy structure on the [4-9] fragment with high agreement of NOE data. Overall, the substitution of Ser12 by Ala12 shifts the conformation of the hydrophobic moiety [10-14] towards a quite random coil structure in this fragment and strongly destabilizes the folded structures of the consensus domain where only one NH (Val7) is solvent-shielded opposed to three (Asp6 to Arg8) in the [Ser12] tetradecapeptide. These conformational changes could be related to the processing enzyme activities on these model oligopeptides.

Conformational studies of an undecapeptide reproducing the consensus sequence around the cleavage site of the RXVRG endoprotease from Xenopus laevis skin

Two synthetic fragments, corresponding to the 4-9 and 4-14 sequences of a tetradecapeptide used as a model to test the RXVRG-endoprotease activity from Xenopus laevis skin, have been studied by two-dimensional nmr spectroscopies, correlated spectroscopy, and nuclear Overhauser effect (NOE) spectroscopy. Both peptides wore the 5-9 consensus sequence found in several hormonal precursors. The nmr data for the 4-9 hexapeptide did not indicate any particular organization, either in water or in dimethylsulfoxide (DMSO), whereas, the 4-14 undecapeptide, a substrate for the RXVRG endoprotease, showed, in DMSO solution, significant trends of structural organization involving the amino acids pertaining to the consensus domain. From variations of integrated NOE peaks with temperature, the apparent interproton correlation times tau c were estimated and the maxima observed with Val7, the central residue in the consensus sequence. A defined tertiary structure in that domain was also supported by medium- and long-range NOEs between Asp6 and Arg8, Glu4 and Gly9, and by the likely involvement of Arg8 and Gly9 NHs in intramolecular hydrogen bonds. Most of these observations could be rationalized by an equilibrium between a 5-8 beta-turn and a 9 > 4 H-bonded loop. The predominance of one rotamer for the C alpha-C beta bond was established in four residues. Finally, the average phi and psi angles were derived from two models taking, or not, into account variations in the correlation times along the sequence. This allowed us to discuss the artefacts generated by using an average correlation time through the whole molecule.

Comparative conformational analyses of mu-selective dermorphin and delta-selective deltorphin-II in aqueous solution by 1H-NMR spectroscopy

Two-dimensional 1H-NMR methods have been used to obtain complete proton resonance assignments and possible solution conformations of dermorphin (H-Tyr-D-Ala-Phe-Gly-Tyr-Pro-Ser-NH2) and deltorphin-II (H-Tyr-D-Ala-Phe-Glu-Val-Val-Gly-NH2), naturally occurring mu- and delta-selective opioids, respectively, in order to examine the conformational characteristics that are closely related to the selectivities towards mu/delta-opioid receptors. With the use of the proton-proton distances derived from ROESY measurements in aqueous solution, 50 possible 3D structures are generated by means of distance geometry calculations. The conformers which satisfy the distance constraints and the torsion angles estimated from JNHC alpha H vicinal coupling constants within the allowable range are then subjected to molecular dynamics simulations for 10 ps after equilibration. Although dermorphin and deltorphin-II are both in equilibrium among many flexible conformers, some conformational differences are observed between these peptides: many conformers of dermorphin show a structure rounded at the N-terminal Tyr-D-Ala-Phe-Gly-Tyr and C-terminal Gly-Tyr-Pro-Ser-NH2 moieties, which are almost at right angles to each other, while those of deltorphin-II are characterized by a 'hook'-shaped backbone structure in which the nearly extended conformation of the Val-Val-Gly-NH2 sequence is located under the folded conformation of the N-terminal Tyr-D-Ala-Phe-Glu sequence. The possible relationship between these conformational characteristics and the mu/delta-opioid receptor selectivities is discussed.