RESEARCH
| REPORT
observed at eight atoms (catalyst F), providing
the desired phosphoramidate in 90% yield and
could address this shortcoming, closing the gap
to “evolutionary perfection.”
18. L.-C. Lee, Y.-L. Lee, R.-J. Leu, J.-F. Shaw, Biochem. J. 397,
9–76 (2006).
6
1
9. Merck Sharp & Dohme Corp., Efficacy and Safety of
grazoprevir (MK-5172) and MK-3682 with elbasvir (MK-8742)
or ruzasvir (MK-8408) for chronic hepatitis C virus (HCV)
genotype (GT) 3, GT4, GT5, and GT6 Infection (MK-3682-012).
in ClinicalTrials.gov. [Internet] Bethesda (MD): National Library
of Medicine (US). 2000- [cited 21 November 2016]. Available
from: https://clinicaltrials.gov/ct2/show/NCT02332720; NLM
Identifier NCT02332720.
97:3 d.r. and validating our computationally driven
Although this catalyst was optimized for use
in a single ProTide therapeutic (MK-3682), the
key principles that govern selectivity proved gen-
erally applicable to other nucleoside analogs. A
variety of ribofuranose-based nucleosides could
be phosphorylated with high P stereoselectivity
by means of this method (Fig. 4). Neither the nu-
cleobase nor the C3′ hydroxyl is an important
driver for selectivity, exemplified in the highly
stereoselective phosphoramidation of dihydro-
furanones 5 and 9. Deoxyribofuranose derivatives,
lacking a C2′ substituent, are tolerated; notably,
the anti-HIV nucleoside derivative azidothymidine
(AZT) (7) could be converted to the prodrug in
93:7 d.r. In a seminal paper, McGuigan et al. re-
ported this nucleoside to yield only a 2:1 mix-
ture of diastereomers when phosphoramidated
in the presence of NMI (33).
Through a distinctive activation mode, we have
designed a practical small-molecule catalyst that
mimics the complex function of an enzyme, culmi-
nating in phosphoramidation that is highly ste-
reoselective at the phosphorus center. The general
applicability of the catalyst to a broad range of
therapeutically relevant nucleophiles invites its
rapid adoption in the preparation of these clin-
ically important molecules.
design principle (Fig. 3A, entry 5). Maximum se-
lectivity was achieved at nine atoms (catalyst G,
98:2 d.r.); further lengthening of the linker resulted
in erosion of both stereoselectivity and chemo-
selectivity, likely on account of undesirable flex-
ibility. Given this observation, we anticipated that
further conformational restriction of the linker
would result in increased selectivity through the
reduction of rotational degrees of freedom. As
such, catalyst I, which contains a central 1,3-phenyl
group, was synthesized and displayed optimal re-
activity [relative rate krel = 10 versus (R)-B (fig. S16)]
and further improved stereoselectivity (99:1 d.r.).
To better understand the effect of catalyst link-
age on stereoselectivity, we reoptimized the transi-
tion states with linked catalyst I (Fig. 3B). Relative
transition state energies predict the desired P-(R)
stereochemistry to be favored by 2.6 kcal/mol, in
agreement with the higher experimentally ob-
served selectivity (31). The same structural prefer-
ences are conserved and further amplified in this
structure (table S15).
Experimental validation of intramolecular co-
operativity was garnered by the observation of
a change to first-order dependence on catalyst
concentration in this system (fig. S14). We were
initially concerned that linking catalysts might
result in deactivation through competitive bis-
phosphoramidation of the dimer. Surprisingly,
P–O bond formation proved considerably faster
than chlorophosphoramidate activation, exem-
plified by first- and zero-order kinetics of chloro-
phosphoramidate 1 and nucleoside 2, respectively
20. P. M. Weintraub, P. L. Tiernan, J. C. Huffman, J. Het. Chem. 24,
61–563 (1987).
1. Z. Zhang, F. Xie, J. Jia, W. Zhang, J. Am. Chem. Soc. 132,
5939–15941 (2010).
2. A. S. De Miranda et al., J. Mol. Catal. B 91, 77–80 (2013).
5
2
2
1
23. D. G. Blackmond, Angew. Chem. Int. Ed. 44, 4302–4320
2005).
(
2
4. The less-selective catalyst (R)-A, based on the same bicyclic
framework but lacking the carbamate functionality, exhibited
a first-order dependence (fig. S10).
25. S. Xu et al., Chemistry 11, 4751–4757 (2005).
2
2
2
6. A traditional nonlinear-effect study was not meaningful in this
case because of the diastereomeric relationship of the products
and the complex equilibria associated with the chiral reactants.
7. Owing to the complexity in sampling van der Waals
precomplexes, we report relative activation free energies in
relation to a common interconverting lowest-energy species.
a a
8. Estimates of the pK (where K is the acid dissociation
constant) of this C–H moiety suggest that it is more acidic
than an amide or imidazole N–H, frequently implicated in
oxyanion hole–type stabilization.
29. R. G. Konsler, J. Karl, E. N. Jacobsen, J. Am. Chem. Soc. 120,
1
0780–10781 (1998).
0. S. E. Denmark, S. K. Ghosh, Angew. Chem. Int. Ed. 40,
759–4762 (2001).
3
4
31. The experimentally determined selectivities for catalysts B
and I correspond to energy differences at 263 K of 1.66 and
2.40 kcal/mol, respectively.
REFERENCES AND NOTES
3
2. W. J. Albery, J. R. Knowles, Biochemistry 15, 5631–5640
1976).
1.
L. P. Jordheim, D. Durantel, F. Zoulim, C. Dumontet, Nat. Rev.
Drug Discov. 12, 447–464 (2013).
(
33. C. McGuigan, R. N. Pathirana, N. Mahmood, A. J. Hay, Bioorg.
Med. Chem. Lett. 2, 701–704 (1992).
2.
D. Cahard, C. McGuigan, J. Balzarini, Mini Rev. Med. Chem. 4,
371–381 (2004).
3
4
.
.
M. J. Sofia et al., J. Med. Chem. 53, 7202–7218 (2010).
U. Pradere, E. C. Garnier-Amblard, S. J. Coats, F. Amblard,
R. F. Schinazi, Chem. Rev. 114, 9154–9218 (2014).
B. S. Ross, P. G. Reddy, H.-R. Zhang, S. Rachakonda,
M. J. Sofia, J. Org. Chem. 76, 8311–8319 (2011).
F. Pertusati, C. McGuigan, Chem. Commun. 51, 8070–8073 (2015).
K. Tran et al., J. Org. Chem. 80, 4994–5003 (2015).
J. M. McCabe Dunn et al., Chem. Sci. 8, 2804–2810 (2017).
B. R. Sculimbrene, S. J. Miller, J. Am. Chem. Soc. 123,
10125–10126 (2001).
ACKNOWLEDGMENTS
We thank Y. Yu for synthesis of catalyst libraries, J. Liu and M. Biba
for separations support, Y. Liu for NMR support, J. Forbes and
H. Gurukar for high-performance computing support, and
E. N. Jacobsen for helpful discussions. Structural parameters for
compounds 3, (R)-B, (R)-J, and (R)-G are available free of charge
from the Cambridge Crystallographic Data Centre under accession
numbers CCDC-1525797, CCDC-1523526, CCDC-1523527, and
CCDC-1523528, respectively. Full details of the kinetic, spectroscopic,
and computational studies are available in the supplementary
materials. Stationary points for all computed structures can be
located in data S1.
(figs. S11 and S12). These data suggest a change in
5.
the turnover-limiting step to chlorophosphorami-
date activation. Given that P–O bond formation
is faster than activation of 1, it seemed plausible
that equilibration may be a stereo-determining fac-
6
.
7.
8
9.
.
31
tor in this system. However, P NMR studies re-
vealed a single observable species, consistent with
fast equilibration and a lower barrier to P inver-
sion through a cooperative displacement mech-
anism (figs. S32 and S33). As argued by Albery
and Knowles in their view of the “evolutionary
perfection” of an enzyme, the enzyme must con-
front the problem of finding a mechanism that
lowers the barrier to each elementary step (32).
We have demonstrated that this catalyst lowers
barriers to P–O bond formation and P(V) inver-
sion but is turnover-limited by the rate of catalyst
addition to 1. Arguably, further development
10. C. M. Longo, Y. Wei, M. F. Roberts, S. J. Miller, Angew. Chem.
Int. Ed. 48, 4158–4161 (2009).
11. J. W. Evans, M. B. Fierman, S. J. Miller, J. A. Ellman, J. Am.
Chem. Soc. 126, 8134–8135 (2004).
12. S. Liu et al., Tetrahedron Asymmetry 23, 329–332 (2012).
3. J. Steinreiber, K. Faber, H. Griengl, Chem. Eur. J. 14,
060–8072 (2008).
4. W. B. Farnham, R. A. Lewis, R. K. Murray Jr., K. Mislow, J. Am.
SUPPLEMENTARY MATERIALS
1
www.sciencemag.org/content/356/6336/426/suppl/DC1
8
Materials and Methods
Figs. S1 to S41
Tables S1 to S17
References (34–63)
NMR Spectra
1
Chem. Soc. 92, 5808–5809 (1970).
15. M. Dutartre, J. Bayardon, S. Jugé, Chem. Soc. Rev. 45,
5771–5794 (2016).
1
6. C. Brown, J. A. Boudreau, B. Hewitson, R. F. Hudson, J. Chem.
Soc. Chem. Commun. 1975, 504–505 (1975).
7. M. Ryan, T. Liu, F. W. Dahlquist, O. H. Griffith, Biochemistry 40,
9743–9750 (2001).
Data S1
1
17 January 2017; accepted 9 March 2017
10.1126/science.aam7936
DiRocco et al., Science 356, 426–430 (2017) 28 April 2017
4 of 4