ACS Medicinal Chemistry Letters
Letter
with multiple structural modifications in the Phe−Phe portion,
comprising amino acid homologation and variation upon the
configuration of the incorporated β2hPhe-NH2 at the fourth
position.
and 4 seems to go hand in hand with an essentially extended
conformation, such as that adopted as a result of β3-hPhe-β2-
hPhe insertion in place of the third and fourth EM-2 residues.
To evaluate the pharmacophoric ability of our analogues 3
and 4, we decided to explore their conformational behavior by
means of computational tools, focusing our attention on the
intramolecular distances outlined by the pharmacophore model,
assembled for bioactivity and MOR selectivity. It suggests as
ideal the key distance of 10−13 Å (C distance) between the
aromatic rings of residues 1 and 3. In addition, two other
distances concerning the message domain have been
considered conclusive for activity, namely, the Tyr1 N−Tyr1
O distance (about 8 Å) and the Tyr1 N−Phe3 ring distance
(about 7 Å).25,26
First a conformational analysis was performed, then the
lowest energy conformers were submitted to a 10 ns molecular
dynamics at 300 K to test their conformational stability.
Calculations were executed on the Tyr1 protonated form of 3
and 4 in water and in DMSO. The conformations were also
clustered according to folded or extended conformations by
measuring the interatomic Cα distance between the first and
the fourth residue. Results showed a preference for folded
conformations in water for the two ligands, both after Monte
Carlo conformational analysis (MC) and molecular dynamics
(MD) simulations. In DMSO, a substantial equilibrium
between folded and extended conformations is established,
with a slight preference for the latter. The pharmacophore A, B,
and C distances were measured and averaged on all found
conformations for both compounds. A good agreement of all
values with the proposed model could be observed for water
environment, while lower C distances have been measured in
DMSO. In both environments, comparison between com-
pounds 3 and 4 highlights a greater value of the C distance for
compound 4, accordingly with its slightly higher activity.
In order to better evaluate the binding orientations and to
study the key ligand−receptor interactions involved in the
molecular recognition process with the μ opioid receptor
(MOR), most active compounds 3 and 4 were submitted to
docking studies. Recently the crystal structure of the mouse
MOR with the irreversible morphinan antagonist β-funaltrex-
amine (β-FNA) has been reported.27 Since it can be supposed
that the reported structure was referred to an inactive form of
the receptor, we rather decided to use an active form of the
MOR model as described by Mosberg.28 The ligands were
docked in to the MOR model flexibly with the Molegro Virtual
Docker software (Figure 3).
The conformationally constrained mimics 1, 2, and 5 are not
as active as 3 and 4, disclosing how deep modifications in the
key address fragment (C-terminal Phe4-NH2), leading to an
appreciable reduction of the conformational freedom, would
result in a decrease in affinity. Finally, the low activity of
peptidomimetic 7 clearly highlights the detrimental effect of β2-
homologation at Phe3, which alters significantly the pharmaco-
phore distances21 of the aromatic ring in the message fragment
(N-terminal tripeptide unit). This date confirms once again that
the proper aromatic ring distance and spatial orientation of the
third residue is highly influencing the binding process,
discriminating a potent lead from a much less active derivative.
Structural studies, in particular 1D and 2D proton NMR,
molecular dynamics, and docking analysis, have been performed
in order to fully investigate the three-dimensional properties of
synthesized analogues (see Supporting Information). At
present, a large number of data have been reported on the
conformational analysis of endomorphins, but conclusions are
rather contradictory as both extended and folded structures
have been suggested as the bioactive conformation. For
instance, the importance of the Tyr1−Pro2 amide bond
conformation of EM-2 in its bioactive form is a controversial
and maybe an overstressed problem. In fact, if the opioid
receptor protein selects its ligand by conformational selection
in a dynamic environment, the conformation of a certain amide
bond alone cannot be a strict determinant for a stable ligand−
receptor interaction. However, relying on total energy measure-
ments, it was recently proposed that the trans isomer of EM-2
and analogues could be the mainly bioactive form and that the
cis isomer is mainly an artifact under the solution conditions.22
We performed NMR studies in DMSO-d6 solution, as this
solvent is considered a good physical approximation of
transport fluids environments23 and is claimed to be a better
approximation for the mechanical and electrostatic environ-
ment of binding to the MOR than D2O is.24 Moreover, being a
good hydrogen bond acceptor, it may reveal intrinsic
conformational preferences, mimicking at the same time the
physical circumstances of receptor−ligand interactions. The
spectra show the presence of two conformers (cis and trans
with respect to the Tyr−Pro peptide bond) for all compounds.
A higher prevalence of the trans conformers (trans/cis ratio
ranging from 1 for most compounds to 4 for 7) has been
observed according to the intensities of OH peaks in 1H NMR
spectra.
Both compound 3 and 4 displayed high negative MolDock
score (−260 and −339, respectively), indicating a strong
favorable interaction of the ligands with the receptor, with a
lower energy for 4, thus in accordance with biological activity
results. Most of the interactions found to be relevant for the
EM-2 activity are observed in the binding mode of 3 and 4. The
Tyr1 residue is deeply placed in the binding site providing
strong interactions with the receptor. An H-bond/electrostatic
interaction between the protonated amino group and the
carboxylate moiety of the Asp147 residue is present, and the
Tyr-OH is involved in an H-bond with the nitrogen of His297.
A π−π interaction with Trp293 is established. Further H-bonds
are formed with Tyr148 and Glu229. Also a stabilizing
interaction with Lys 233 is established. Other important
interactions are found between lipophilic regions of the ligand
and Met151, Phe152, Ile296, Val300, and Ile322. The binding
mode of 4 showed additional π−π interactions between Tyr1
From relevant inter-residue ROE interactions, some common
conformational behaviors can be highlighted within all
compounds. In particular, the Tyr1 aromatic ring seems to be
always folded over the Pro2 segment (ROE between Pro H-δ
and various Tyr protons), and in most cases, the NH of the β-
hPhe3 interacts with the Pro’s H-α. Constrained compounds 1
and 2 present additional ROE contacts, involving the β2-hTbac4
residue from one side, and the benzylic or H-α protons of β3-
hPhe3, on the other side, and clearly underlining a folding
between the third and fourth residue. No significant inter-
residue ROE interaction could be observed for linear
derivatives, revealing preferentially extended conformations.
The only exception regards compound 6 where the presence of
multiple ROE contacts would suggest a prominent folding of
the overall structure. Back to biological results, high activity of 3
797
dx.doi.org/10.1021/ml400189r | ACS Med. Chem. Lett. 2013, 4, 795−799