ACS Medicinal Chemistry Letters
Letter
BLT2 and the molecular mechanism responsible for ligand-
dependent efficacy for BLT2 remains poorly understood.
In this study, we aim to elucidate the mechanism of ligand
efficacies at BLT2 using a series of molecular modeling
techniques and experimental verification. First, we predicted
the three-dimensional (3D) structures of BLT2 in apo and
agonist- or antagonist-bound forms at an atomic scale using
molecular modeling. Based on the predicted structures, we
identified residues that could determine ligand efficacy.
Subsequently, we verified our predictions through site-directed
mutagenesis and cell motility experiments. By chemically
modifying a known agonist, we further confirmed our
proposed ligand efficacy mechanism thus enabling the design
of partial agonists. Additionally, we discuss in detail our
molecular modeling procedures of class A GPCR structure
prediction and GPCR ligand docking used to identify the
mechanism for ligand efficacy at BLT2.
according to their efficacies.16 Our predicted structure of BLT2
is based on an antagonist-bound structure of BLT1; therefore,
to accommodate the possibility of conformational changes
when an agonist binds BLT2, we introduced a structure
optimization step for protein−ligand complex after the
docking. It is well-known that the overall conformational
changes of GPCRs in active and inactive forms are difficult to
simulate.17 On the other hand, the structure optimization was
performed only on the binding site residues within 5 Å of the
ligand using the Prime18 module of the Schrodinger suite.
̈
We first predicted a binding pose of a native agonist
leukotriene B4 (LTB4). In the binding pose, LTB4 adopts a U-
shaped conformation in the binding site of BLT2 (Figure 1).
First, we predicted the 3D structure of human BLT2 (Figure
S1) using the GPCR homology modeling method we have
guinea pig BLT1 (PBD ID: 5X33) recently revealed was used
as the template for homology modeling. BLT2, like BLT1,
belongs to the leukotriene receptor family and recognizes the
same native ligands including LTB4. Human BLT2 and BLT1
exhibit 39.5% sequence identity and have similar amino acid
compositions in the ligand binding sites (Figure S2).
Seven transmembrane (TM) helices of GPCRs show
considerable variations in the helix shape because of kink
structures.13 Because conformations of TMs determine the
shape and volume of the ligand-binding site, accurately
predicting the helix conformation is crucial. Therefore, we
added two additional modeling steps. First, we ensured that the
position of proline or glycine that can form a kink in the helix
structure is consistent between template-BLT1 and target-
BLT2 (Figure S2). Second, we performed molecular dynamics
(MD) simulation to optimize the predicted structure of BLT2
under explicit membrane and solvent environments. MD
results showed that the predicted BLT2 structure rapidly
entered the equilibrium state at 200 ns and maintained stable
conformations to 1 μs (Figure S3). For comparison, we
performed another MD calculation on the crystal structure of
BLT1. At the equilibrium state, average heavy atom RMSD
values of the crystal structure of BLT1 and the predicted
structure of BLT2 were 4.84 and 4.86 Å, respectively.
Considering the thermal fluctuation of the BLT1 crystal
structure, one can conclude that the predicted BLT2 structure
did not undergo significant structural changes during the MD
simulation. This result also implies that the initial predicted
BLT2 structure is as stable as the crystal structure of BLT1.
The only notable structural difference between the predicted-
BLT2 and template-BLT1 is the conformation of extra-cellular
loop 2 (ECL2). Among the six loop structures in GPCRs,
ECL2 was found to be in direct interactions with the bound
ligand.14 The predicted BLT2 structure has a distinct parallel β
sheet in ECL2 like other members of γ-branch GPCRs
including BLT1 and opioid receptors.15 However, the ECL2 of
BLT2 is four residues shorter and located further inside the
ligand binding site than that of BLT1 (Figure S1).
Figure 1. Structural features of agonist LTB4 binding in BLT2 (green
surface: hydrophobic, cyan surface: polar uncharged).
The position 1 carboxyl group located in the binding site
entrance participates in various interactions with residues in
ELC2; the position 4 hydroxyl group forms hydrogen bonds
with Gln2987.33; and the position 12 hydroxyl group forms a
hydrogen bond with residue Tyr2716.51, located deep in the
binding site pocket. Long hydrophobic chain position 13−20
forms van der Waals contacts with the hydrophobic region in
the inner binding site. The predicted binding pose of LTB4 is
similar to that obtained from the nuclear magnetic resonance
(NMR) experiment.19 This agreement strongly suggests that
our predicted BLT2 structure and modeling methods are
adequate for studying BLT2 and ligand interactions.
We conducted additional docking simulations using the
agonists 12(S)-HETE, 12-HHT, and CAY1058320 to find
distinct interactions of agonists. 12(S)-HETE and 12-HHT are
native BLT2 agonists derived from arachidonic acids such as
LTB4, whereas CAY10583 is a synthetic agonist. Prior to the
interaction analysis, we calculated binding free energies of the
four agonists based on the docking results using the molecular
mechanics/generalized Born surface area (MM/GBSA)
method, and compared the calculation results with the
experimentally measured binding affinities to ensure the
reliability of the docking results. Experimentally, 12-HHT
exhibited the highest binding affinity for BLT2, followed by
LTB4 and 12(S)-HETE, respectively.21 Our calculated binding
energies were in the order consistent with the experimental
results (Table S1). The two native BLT2 agonists commonly
adopt a U-shaped conformation similar to LTB4 and
CAY10583 was bound in the transverse form of the binding
Following the structure prediction of BLT2, we performed a
docking simulation on the BLT2 structure to predict the
binding pose of ligands. GPCRs have various conformational
ensembles between active and inactive states, and the ligands
form complexes with their respective states of receptors
B
ACS Med. Chem. Lett. XXXX, XXX, XXX−XXX