MedChemComm
Concise Article
Fig. 5 HINT H-bond component score for ring interactions with Cys241b. Closed
squares represent compounds possessing both H-bond acceptors and appropri-
ately placed hydrophobic groups. These compounds generally possess superior
pEC50s. Open triangles represent compounds with weak or no acceptors. Open
circles represent compounds with steric issues and/or lacking key hydrophobic
interactions.
Fig. 4 Specific hydrogen bonding (yellow) and hydrophobic (green) interactions
in subpocket A. Compounds 2 (solid white), 5h (translucent purple), 5i (translucent
green) and 5m (translucent red) are shown. Notes: (1) the key H-bond interaction is
with Cys241b, which is strongest with the O of methoxy in the ring's meta position;
(2) some analogues, e.g., with F, can weakly H-bond with the NH of Leu252b; (3)
the CH3 of p-methoxy has key hydrophobic interactions with Leu242b; (4) Ile378b
has hydrophobic interactions with m-methoxy or the rings of 5h or 5i.
two-fold; rst, the hydrogen bonding interaction with Cys241b is
the key predictor, absent of steric clashes, for the microtubule
inhibitory activity for this set of analogues; second, other inter-
actions in the pocket, i.e., hydrophobic, are also necessary, but
competitive with this weakly scored hydrogen bonding.
direction. The better hydrogen bonding for a meta-position
substituent explains the activity of 5f compared to 5e and other
similar cases. As for 5h and 5i, the distal (from the pyrrole core)
rings were located directly beneath the thiol hydrogen, thus
acting as acceptors for the weak but critical hydrogen bond, but
pocket steric issues cancelled this advantage. The uorine atom
of 5m is also located at the meta-position, but the docking study
suggested that an 180ꢄ ring ip shied its position in space
such that, although the uorine was anchored by the backbone
NH of Leu252b, it provided no additional bonding to Cys241b
compared to 5e. In the case of 5g, the detrimental effect of the
The importance of both hydrophobic interactions and
hydrogen bonding in subpocket A was seen in the SAR analysis
and modeling studies. The latter dictates whether the C-4
analogues of pyrrole-based antitubulin agents display micro-
tubule depolymerizing activity and the strength of that activity,
while the character of the pocket requires predominantly
hydrophobic moieties. Underestimation of the Cys241b inter-
action was one probable reason that the total HINT score was a
poor predictor of microtubule depolymerizing activity. This
thiol group acts as a hydrogen bond donor and while this type of
hydrogen bonding interaction is generally regarded as weak and
is thusly parameterized by HINT, it is not even considered by
many other scoring functions. For the downstream biological
effect, microtubule depolymerization, the interaction assumed
to be weak surprisingly stands out as a key factor. In fact, its
absence might produce a different mechanism of action even
when other portions of the structure are exactly the same, as
shown particularly by 5d with potent antiproliferative activity
(0.312 mM) but weaker microtubule depolymerization activity
(ꢂ94 mM). Cys241b has been previously identied as an
important target residue for colchicine site agents.16 In a study
of 15 structurally diverse colchicine site inhibitors, the docked
binding modes of all included hydrogen bonding to Cys241b
˚
third methoxy is visually apparent: the tight distance (3.58 A)
between the backbone of Leu252b and the phenyl ring of the
ligand can lead to signicant steric clashes with a large
substituent such as the 5-methoxy.
The total HINT scores of C-4 analogues fail to show a tight
relationship with pEC50 (ESI, Fig. S6†). However, isolating
the HINT score for hydrogen bonding interactions involving
Cys241b for a subset of analogues (2, 5e, 5f, 5j, 5k, 5m and 5q)
that place, as separate entities, appropriately positioned hydro-
phobic groups and a hydrogen bond acceptor in the subpocket
(while not inducing steric clashes),‡ reveals a linear relation with
respect to these compounds' pEC50s (Fig. 5). The implications are
‡ Compounds 5a and 5b do not possess H-bond acceptor substituents to the
phenyl ring. The halogens of 5c, 5d and 5l are considered as weak acceptors by (Cys239b in that study).17 Our combined SAR and modeling
HINT and, while their dominant property is hydrophobicity, these halogen
study conrms the importance of that cysteine. It should be
noted that there is potentially a systematic error in our proce-
dure. As GOLD optimizes ligand placement with a different
atoms are not at ideal hydrophobic contact distance. Similarly, the aromatic
rings of 5h and 5i can simultaneously be weak hydrogen bond acceptors and
strongly hydrophobic, but these analogues are inappropriately shaped to make
forceeld (set of rules) than used by HINT in scoring, subtle
all necessary hydrophobic contacts in the pocket (see Fig. 4). Compounds 5g,
5n–5p have a number of structural or conformational issues as they are docked structural effects, or in this case, the interplay of several of
in subpocket A: there is insufficient room for the trimethoxyphenyl of 5g or the
them, are not well scored post-docking as none of the models
generated by GOLD capture the set of features in a single model
that HINT would score highest. This is likely to be a general
ethoxyl-2-napthyl of 5n; and although the two ether oxygens of 5o and 5p can
participate in hydrogen bonding with Cys241b as in 2, the methylene and
ethylene, respectively, are forced to be close to the thiol sulfur and thus
observation in docking/rescoring studies, irrespective of
utilized scoring functions, when subtle effects are at play.
produce unfavorable hydrophobic–polar interactions. In contrast, the methyls
of 2 can jacknife away from the sulfur.
420 | Med. Chem. Commun., 2013, 4, 417–421
This journal is ª The Royal Society of Chemistry 2013