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

Conformational heterogeneity of a leucine enkephalin analogue in aqueous solution and sodium dodecyl sulfate micelles: comparison of time-resolved FRET and molecular dynamics simulations

We have undertaken time-resolved Förster resonance energy transfer (FRET) and molecular dynamics simulations to analyze conformations and conformational heterogeneity of an analogue of leucine enkephalin in solution and in the presence of sodium dodecyl sulfate (SDS) micelles. Enkephalins are opioid pentapeptides that interact with opioid receptors in the central nervous system. We used time-correlated single-photon counting to detect energy transfer between the N-terminal tyrosine and a tryptophan residue substituted for phenylalanine at the 4 position. FRET from Tyr to Trp was measured over a temperature range from 5 to 55 degrees C in aqueous solution. By taking into account Tyr rotamer interconversion rates measured previously, we determined average distances between Tyr and Trp for the two populated rotameric conformations of Tyr. Molecular dynamics simulations (100 ns) support this analysis and indicate extensive conformational heterogeneity. The simulations also predict that the FRET orientational factor is correlated with the Tyr-Trp separation. Failure to account for the correlation between orientation and distance results in errors that appear to be largely offset in the leucine enkephalin analogue (YGGWL) by a weighting bias inherent in the R(-6) dependence of the energy-transfer rate. The Tyr lifetimes decrease upon titration of the peptides with SDS, indicating formation of compact conformations of the peptide in the micelle environment. This result is consistent with the conjecture that the lipid environment may induce formation of bioactive conformations of the peptide.

Quantitative conformationally sampled pharmacophore for delta opioid ligands: reevaluation of hydrophobic moieties essential for biological activity

Recent studies have indicated several therapeutic applications for delta opioid agonists and antagonists. To exploit the therapeutic potential of delta opioids developing a structural basis for the activity of ligands at the delta opioid receptor is essential. The conformationally sampled pharmacophore (CSP) method (Bernard et al. J. Am. Chem. Soc. 2003, 125, 3103-3107) is extended here to obtain quantitative models of delta opioid ligand efficacy and affinity. Quantification is performed via overlap integrals of the conformational space sampled by ligands with respect to a reference compound. Iterative refinement of the CSP model identified hydrophobic groups other than the traditional phenylalanine residues as important for efficacy and affinity in DSLET and ICI 174 864. The obtained models for a structurally diverse set of peptidic and nonpeptidic delta opioid ligands offer good predictions with R2 values>0.9, and the predicted efficacy for a set of test compounds was consistent with the experimental values.

A multidimensional 1H NMR investigation of the conformation of methionine-enkephalin in fast-tumbling bicelles

Enkephalins are pentapeptides found in the central nervous system. It is believed that these neuropeptides interact with the nerve cell membrane to adopt a conformation suitable for their binding to an opiate receptor. In this work, we have determined the three-dimensional structure of methionine-enkephalin (Menk) in fast-tumbling bicelles using multidimensional (1)H NMR. Bicelles were selected as model membranes because both their bilayer organization and composition resemble those of natural biomembranes. The effect of the membrane composition on the peptide conformation was explored using both zwitterionic (PC bicelles) and negatively charged bicelles (Bic/PG). Pulsed field gradient experiments allowed the determination of the proportion of Menk bound to the model membranes. Approximately 60% of the water-soluble enkephalin was found to associate to the bicellar systems. Structure calculations from torsion angle and NOE-based distance constraints suggest the presence of both micro - and delta-selective conformers of Menk in each system and slightly different conformers in PC bicelles and Bic/PG. As opposed to previous studies of enkephalins in membrane mimetic systems, our results show that these opiate peptides could adopt several conformations in a membrane environment, which is consistent with the flexibility and poor selectivity of enkephalins.

Opioid receptor three-dimensional structures from distance geometry calculations with hydrogen bonding constraints

Three-dimensional structures of the transmembrane, seven alpha-helical domains and extracellular loops of delta, mu, and kappa opioid receptors, were calculated using the distance geometry algorithm, with hydrogen bonding constraints based on the previously developed general model of the transmembrane alpha-bundle for rhodopsin-like G-protein coupled receptors (Biophys. J. 1997. 70:1963). Each calculated opioid receptor structure has an extensive network of interhelical hydrogen bonds and a ligand-binding crevice that is partially covered by a beta-hairpin formed by the second extracellular loop. The binding cavities consist of an inner "conserved region" composed of 18 residues that are identical in delta, mu, and kappa opioid receptors, and a peripheral "variable region," composed of 19 residues that are different in delta, mu, and kappa subtypes and are responsible for the subtype specificity of various ligands. Sixteen delta-, mu-, or kappa-selective, conformationally constrained peptide and nonpeptide opioid agonists and antagonists and affinity labels were fit into the binding pockets of the opioid receptors. All ligands considered have a similar spatial arrangement in the receptors, with the tyramine moiety of alkaloids or Tyr1 of opioid peptides interacting with conserved residues in the bottom of the pocket and the tyramine N+ and OH groups forming ionic interactions or H-bonds with a conserved aspartate from helix III and a conserved histidine from helix VI, respectively. The central, conformationally constrained fragments of the opioids (the disulfide-bridged cycles of the peptides and various ring structures in the nonpeptide ligands) are oriented approximately perpendicular to the tyramine and directed toward the extracellular surface. The results obtained are qualitatively consistent with ligand affinities, cross-linking studies, and mutagenesis data.

Development of a common 3D pharmacophore for delta-opioid recognition from peptides and non-peptides using a novel computer program

A unified three-dimensional (3D) pharmacophore for recognition of the delta-opioid receptor by families of structurally diverse delta-opioid ligands, including peptides and non-peptides, has been determined. An additional structural feature required for delta-selectivity was also characterized using a subset of these ligands that are highly selective for the delta-opioid receptor. To obtain these pharmacophores, we have used a recently developed computer program that performs systematic and automated comparisons of molecules to determine whether any common 3D relationships exist among candidate recognition moieties in high-affinity analogs. All the low-energy conformations of each ligand are included in these comparisons. The program developed should be applicable in general to molecular superimposition problems in rational drug design and to develop both 3D recognition and activation pharmacophores for any receptor for which high- and low-affinity analogs and agonists and antagonists have been identified.