Bioactive Conformation in Bimolecular Recognition
A R T I C L E S
(Alwyn Jones, Upsala University, Sweden),16 and model refinement
was performed with the CNX software (Accelrys). The new model
was improved by a cycling procedure of model rebuilding and
refinement steps. The final model include two molecules of NS5B
(residues A1-A149, A154 to A563, B1 to B147, and B153 to B563).
The tip of the Λ2 finger loop (150-152) appeared to be flexible
and was not modeled. See the Table S1 in the Supporting
Information for the refinement statistics.
The structure of the complex involving compound 8 was solved
by MR using the apo structure as a model. On the second cycle,
the residues B18 to B35 were removed as the initial density map
was indicating a clear conformation change of that loop. The
subsequent difference e-density maps revealed a clear density for
the inhibitor molecule on the newly exposed surface of the thumb
domain. The final model includes two molecules of NS5B (residue
A1 to A149, A154 to A563, B1 to B17, B36 to B148, and B153 to
B563) and one molecule of compound 8 associated with NS5B
chain B. See Table S1 in the Supporting Information for refinement
statistics.
All data sets were processed using the TOPSPIN or WinNMR
software (Bruker Canada, Milton, Ontario) and software. Data sets
were typically zero-filled to yield 2048 (f2) × 1024 (f1) points
after Fourier transformation using a shifted sine-bell window
function. 1H and 13C spectra were chemical shift calibrated relative
to the standard values attributed to the DMSO peaks (2.49 and 39.51
ppm, respectively).
Molecular Dynamics Simulations and Docking. Molecular
dynamics (MD) simulations were performed for compounds 5 and
6 to generate free-state conformational histograms for angles S1
and S2 as shown in Figure 4. Each compound was first submitted
to a conformation sampling to generate 10 diverse conformations
that were then, each, submitted to 1 ns. MD run. The MD
simulations were run with MOE (versions 2008.10) under constant
volume and constant temperature conditions using MMFF94 (as
implemented in MOE version 2008.10). Each MD simulation
consisted of a heating phase that gradually increased the system
temperature from 0 K to 298 K in 10 ps, followed by a equilibration
phase at the target temperature of 298 K for 10 ps. The production
phase (data accumulation phase) consisted of 1 ns and snapshots
of the system were recorded every 0.5 ps. The nonbonded cutoff
radius was set to infinity and an implicit solvation model was used
to mimic an aqueous environment. The production phase trajectories
of the 10 MD runs were then combined to analyze the S1 and S2
angles with the Conformation Geometries module of MOE.
Compounds were docked to the allosteric “thumb pocket I” of
HCV polymerase using Glide (version 2007) in XP mode and with
the use of experimentally derived structures from X-ray crystal-
lography and bound-state information from NMR binding and
competition studies.
Despite the medium resolution of the data, most of the inhibitor
molecule was well-defined in the e-density maps with the exception
of the (N-acetyl substitution) which should be considered as flexible.
Also, the exact orientation of the furan group is not clearly defined
by the data as no H-bond interaction was found to differentiate the
oxygen atoms from the carbon of the ring.
X-ray data collection and refinement statistics are provided in
the Supporting Information. Also, both structures have been
deposited in the Protein Data Bank (accession code 3MWV and
3MWW for the apo and the inhibitor complex respectively).
NMR Binding and Competition Experiments. NMR experiments
were performed on Bruker Avance spectrometers (400 and/or 600
1
MHz H frequency) at 27 °C, unless otherwise stated.
Conclusions
Samples to be used for NMR binding and competition experi-
ments were prepared by adding concentrated inhibitor in DMSO-
d6 (typically 0.3 mgs. of compound) to an aqueous buffer composed
of 20 mM Tris-d11, 2 mM DTT-d10, 1 mM EDTA-d12, 300 mM
NaCl, and 10% (v/v) D2O spiked with TSP-2,2,3,3-d4. Buffer was
added to a final volume of 600 µL, and the pH was adjusted to 6.0
with diluted HCl. Spectra of free compound at 0.2 mM were then
acquired. A concentrated stock solution of HCV polymerase was
then added at less than stochiometric amounts to the NMR tube,
and spectra were again acquired at each increment such that the
inhibitor to polymerase ratio ranged from 40:1 to 10:1. The
polymerase stock solution contained 0.5 mM HCV polymerase
(NS5B∆21C-His6) in buffer consisting of 20 mM Tris-d11, 2 mM
DTT-d10, 1 mM EDTA-d12, 300 mM NaCl, and 10% (v/v) glycerol-
d8. NOESY spectra were also collected on these samples. For
competition experiments, concentrated stock solutions of other
compounds were added to the NMR tube followed by the
acquisition of spectra. TROSY NMR experiments were run using
26 µM NS5B∆21C-His6 (1b/J4) labeled with 15N in Tris-d11 25
mM, pH 7.5, NaCl 150 mM, glycerol-d8 5%, TCEP-d16 1 mM.
For spectra acquired on samples in buffer, suppression of the
solvent signal was achieved by the use of presaturation or by
inserting a 3-9-19 WATERGATE module prior to data acquisi-
tion. NOESY experiments on sample tubes containing no poly-
merase resulted in the observation of no significant cross-peaks.
However, NOESY spectra on sample tubes containing polymerase
resulted in the observation of many cross-peaks which contained
the valuable interhydrogen distance information of the compound
when bound to HCV polymerase. Thus, a series of NOESY spectra
were typically acquired on these latter samples which included the
following mixing times 50, 100, 200, and 300 ms (600 MHz).
NOESY crosspeak volumes were scaled and converted to apparent
interproton distances. The derived distances were normalized via
comparisons with the volumes from an internal fixed interhydrogen
distance.
In the context of deriving novel inhibitors of HCV poly-
merase, the bimolecular recognition we described here may be
represented by the collision of two flexible objects. Recognition
involves a disengaged “finger loop” from an allosteric site of
the thumb domain of the polymerase, and binding of ligands in
the bioactive “L” conformation (series A and B). From a ligand
design perspective, the work described here presents a good
example of how a ligand’s “structural hinge” influences potency,
and due to its solvent exposure in the complex, reasonable
structure-activity-relationships can be interpreted. Although the
other segments of the inhibitors in series A, B, and C are also
important for potency, such as the cyclohexyl and the lipophilic
benzimidazole/indole cores, the structural hinge can impart
significant potency that can be “fine tuned”. Moreover, variants
of the hinge also significantly influence cell culture activity and
other properties, which ultimately allowed us to successfully
transition from compounds that had poor physicochemical and
pharmacodynamic properties (series A) to a more promising
series (series B). As expected, the improvements in cell culture
potency observed in series A upon replacement of the benz-
imidazole by an indole core (compounds 3 to 4) were also
translated in the new series B (compounds 6 to 7), and likely
arise from an increase in lipophilic character (and cell membrane
permeability).
Overall, compound 7 maintains good enzymatic potency and
has an impressive cell-culture potency (IC50 ) 0.095 µM, EC50
) 0.072 µM,). Unlike compound 4, it also has beneficial
physicochemical and pharmacokinetic properties that promise
new opportunities for progress toward the development of novel
HCV therapeutics in the clinic. Indeed, a close analogue of
compound 7 has shown antiviral effects in clinic studies.3f
Acknowledgment. We are grateful to M.-A. Poupart for his
contributions, and to S. Lefebvre and G. McKercher for materials,
(16) Jones, T. A.; Zou, J.-Y.; Cowan, S. W.; Kjeldgaarrd, M. Acta
Crystallogr. 1991, A47, 110–119.
9
J. AM. CHEM. SOC. VOL. 132, NO. 43, 2010 15211