8
Y. Yamashita et al. / Bioorg. Med. Chem. xxx (2014) xxx–xxx
tiplet). HR-MS were measured on a Jeol JMS-700 MStation. The
progress of all reactions was monitored by thin-layer chromatogra-
phy with silica gel glass plates (60 F254) (Merck Ltd, Tokyo, Japan),
and spots were visualized with ultraviolet light (254 nm) and
stained with 5% ethanolic phosphomolybdic acid. Column chroma-
tography was performed using silica gel 60N (Kanto Chemical Co.,
Tokyo, Japan).
HPLC-UV chromatograms were acquired in a Waters Alliance
2695 chromatographer equipped with a Waters 2996 photodiode
array detector (Waters, Milford, MA). HPLC analysis was conduced
according to method A (see below) with the retention time
expressed in min detected at 220 nm. For HPLC method A, chroma-
(10À10 to 3 Â 10À5 M) of each mepenzolate were tested in triplicate
to generate competition curves. Non-specific binding was deter-
mined in the presence of atropine (2.5 lM). The samples were
passed through a GF/C filter (Filtermat A, PerkinElmer Life and
Analytical Sciences, Boston, MA) that was pre-incubated for 1 h
with wash buffer (50 mM Tris/HCl (pH 7.4), 100 mM NaCl) contain-
ing 1.0% polyethylenimine, and washed four times with ice-cold
wash buffer. Filters were then dried for 30 min before attachment
to MeltiLex A (melt-on scintillation sheet; PerkinElmer Life and
Analytical Sciences, Boston, MA). The radioactivity remaining on
the filter was monitored with MicroBeta Trilux microplate scintil-
lation counter (PerkinElmer Life and Analytical Sciences, Boston,
MA). Affinities at equilibrium were determined as equilibrium
antagonist dissociation constant (Ki) values after correcting the
experimentally determined IC50 values with the experimentally
determined Kd value of NMS for hM2R or hM3R and the concentra-
tion of NMS, as described previously.35 The Ki value was obtained
from three independent curves. All adjustments were performed
using Prism (GraphPad Software, Inc., San Diego, CA).
tography was performed on
a
Daicel Chiralpak IA
(250 mm  46 mm). The mobile phase, at a flow of 1.0 mL/min,
was a binary gradient of water (containing 0.1 M potassium hexa-
fluorophosphate) and acetonitrile, 30:70.
5.2.1. (R)-Mepenzolate ((R)-1)
To a solution of benzilic acid (830 mg, 3.6 mmol) in N,N-dimeth-
ylformamide (DMF) (8 mL), CDI (883 mg, 5.4 mmol) was added and
the mixture was stirred for 15 min at room temperature. To this
5.4. Homology modelling
mixture, a solution of (R)-3 (500 lL, 4.3 mmol) in DMF (4 mL)
was added dropwise at 80 °C and the resulting mixture was stirred
for 18 h at the same temperature. After cooling to room tempera-
ture, the reaction was quenched with water and organic materials
were extracted with ethyl acetate. The combined extract was
washed with brine, dried over sodium sulfate and concentrated
in vacuo. The residue was then applied to a short silica gel column,
eluted with ethyl acetate and concentrated in vacuo to give (R)-4 as
a yellow oil, which was used for the next step without further
purification.
To a solution of (R)-4 (860 mg, 2.6 mmol) in acetonitrile
(10 mL), methyl bromide (2.0 M in tetrahydrofuran, 6.2 mL,
12.4 mmol) was added and stirred for 5 h at room temperature.
The precipitates were filtered off and re-crystallized from dichloro-
methane and methanol to give (R)-1 as colorless fine needles
The structure of rat M3R bound to its antagonist, tiotropium,
was recently solved by X-ray crystallography (PDB code 4DAJ.22
Although a T4-lysozyme was fused to intracellular loop 3 (ICL3)
for enhancing crystallization, there were no residue gaps and few
residue mismatches between the crystal structure and hM3R
amino acid sequence except for both terminus regions and the
ICL3. Therefore, we constructed an hM3R model by fixing the miss-
ing atoms in the crystal and replacing the mismatched amino acid
residues with human ones. The fused T4-lysozyme is replaced with
a ten amino acid linker sequence (GGGGSGGGGS), because this
region has no reliable template for modelling. These modelling
procedures were performed using Modeller v9.4.36
5.5. Molecular dynamics
(980 mg, 64% 2 steps). IR
mmax: 3432, 3315, 1735, 1216, 1093,
1068, 923, 705 cmÀ1 1H NMR (DMSO-d6): d = 7.28–7.37 (m,
;
MD simulations were carried out in an explicit membrane and
water system by using myPresto/cosgene.37 The whole structure
consists of the hM3R model, (R)-1 or (S)-1, four cholesterols, 146
1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine and 10,212
water molecules (TIP3P model) including 21 Na+ and 35 ClÀ parti-
cles.38 NaCl was added to achieve the saline condition and exces-
sive ClÀ was added to neutralize the total charge of the system.
The AMBER99 force field parameters39 were used for hM3R. Force
field parameters for (R)- or (S)-mepenzolate were obtained accord-
ing to the generalized AMBER force field procedure40 with partial
charges derived from quantum chemical calculations by GAUSSIAN
03 at the HF/6-31G⁄ level of theory. Periodic boundary conditions
were applied and Berendsen’s method for temperature and pres-
sure coupling was adopted (300 K and 1 atm, respectively). After
performing 400 steps of energy minimization by the steepest des-
cent method and the conjugate gradient method, the NPT (constant
particle number, constant pressure, and constant temperature)41
for 200 ps was performed to obtain an equilibrated system with
the following conditions: periodic boundary, 300 K, 1 atm, cut-
10H), 6.82 (s, 1H), 5.28 (m, 1H), 3.63 (dd, J = 13.5, 3.5 Hz, 1H),
3.49 (dd, J = 13.5, 4.7 Hz, 1H), 3.37 (m, 1H), 3.33 (m, 1H), 3.10 (s,
3H), 2.82 (s, 3H), 1.63–1.85 (m, 4H); 13C NMR (DMSO-d6):
d = 171.9, 143.0, 142.7, 128.0, 127.9, 127.7, 127.6, 127.0, 80.7,
26
66.8, 61.6, 60.8, 25.2, 16.2; [
a]
À8.12 (c 1.00, methanol); mp
D
224.2–224.9 °C; HR-MS (FAB): calcd for
C
21H26O3N: [M+1]+:
340.1913; found: m/z = 340.1905. HPLC analysis was done accord-
ing to method A (retention time, 16.8 min; single peak).
5.2.2. (S)-Mepenzolate ((S)-1)
The title compound ((S)-1) was synthesized using the same pro-
cedure described for the preparation of (R)-mepenzolate, except
for the use of (S)-3-hydroxy-1-methylpiperidine instead of (R)-iso-
mer. The final sample was a colorless needle (934 mg, 61% 2 steps).
25
[a
]
+8.33 (c 1.00, methanol). Other physical and spectral data of
D
this enantiomer were in good accordance with those of (R)-1. HPLC
analysis was done according to method
15.0 min; single peak).
A (retention time,
off = 12 Å, and 1.0 fs/step. The ZD method42 = 0.0 ÅÀ1) was
(a
5.3. Filter-binding assay
applied for the calculation of electrostatic force terms. Further
200 ps equilibrium calculations were performed under NVT (con-
stant particle number, constant volume, and constant tempera-
ture) conditions with the cell size maintained (300 K, 1.0 fs/step
for 50 ps, 0.75 fs/step for 60 ps, 1.0 fs/step for 40 ps and 2.0 fs/step
for 50 ps). After these preliminary calculations, a production run at
300 K was performed with the same cell size and the time step was
set to 2 fs. The SHAKE algorithm was applied to fix all the
bonds involving hydrogen during the NVT simulations. Snapshot
The filter-binding assay was done as described previously34
with some modifications. Membrane fractions prepared from
CHO-K1 cells expressing hM2R or hM3R (Membrane Target Sys-
tems, Perkin–Elmer Life and Analytical Sciences, Boston, MA; pro-
tein concentration, 10
[3H]NMS (85.5 Ci/mmol) at room temperature for 2 h in 200
PBS in the presence of each mepenzolate. A range of concentrations
lg/well) were incubated with 2 nM
lL