Highly accelerated inverse electron-demand cycloaddition of electron-deficient azides with aliphatic cyclooctynes

Abstract

Strain-promoted azide-alkyne cycloaddition (SPAAC) as a conjugation tool has found broad application in material sciences, chemical biology and even in vivo use. However, despite tremendous effort, SPAAC remains fairly slow (0.2-0.5 M^-1 s^-1 ) and efforts to increase reaction rates by tailoring of cyclooctyne structure have suffered from a poor trade-off between cyclooctyne reactivity and stability. We here wish to report tremendous acceleration of strain-promoted cycloaddition of an aliphatic cyclooctyne (bicyclo[6.1.0]non-4-yne, BCN) with electron-deficient aryl azides, with reaction rate constants reaching 2.0-2.9 M^-1 s^-1 . A remarkable difference in rate constants of aliphatic cyclooctynes versus benzoannulated cyclooctynes is noted, enabling a next level of orthogonality by a judicious choice of azide-cyclooctyne combinations, which is inter alia applied in one-pot three-component protein labelling. The pivotal role of azide electronegativity is explained by density-functional theory calculations and electronic-structure analyses, which indicates an inverse electron-demand mechanism is operative with an aliphatic cyclooctyne.

ARTICLE

Received 22 Jul 2014 | Accepted 25 Sep 2014 | Published 10 Nov 2014

Strain-promoted azidealkyne cycloaddition (SPAAC) as a conjugation tool has found broad application in material sciences, chemical biology and even in vivo use. However, despite tremendous effort, SPAAC remains fairly slow (0.20.5 M 1 s 1) and efforts to increase reaction rates by tailoring of cyclooctyne structure have suffered from a poor trade-off between cyclooctyne reactivity and stability. We here wish to report tremendous acceleration of strain-promoted cycloaddition of an aliphatic cyclooctyne (bicyclo[6.1.0]non-4-yne, BCN) with electron-decient aryl azides, with reaction rate constants reaching 2.02.9 M 1 s 1.

A remarkable difference in rate constants of aliphatic cyclooctynes versus benzoannulated cyclooctynes is noted, enabling a next level of orthogonality by a judicious choice of azide cyclooctyne combinations, which is inter alia applied in one-pot three-component protein labelling. The pivotal role of azide electronegativity is explained by density-functional theory calculations and electronic-structure analyses, which indicates an inverse electron-demand mechanism is operative with an aliphatic cyclooctyne.

DOI: 10.1038/ncomms6378

Highly accelerated inverse electron-demand cycloaddition of electron-decient azides with aliphatic cyclooctynes

Jan Dommerholt1, Olivia van Rooijen2, Annika Borrmann1, Clia Fonseca Guerra2, F. Matthias Bickelhaupt1,2 & Floris L. van Delft1

1 Institute for Molecules and Materials, Radboud University Nijmegen, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands. 2 Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, VU University Amsterdam, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands. Correspondence and requests for materials should be addressed to F.M.B. (email: mailto:[email protected]

Web End [email protected] ) or to F.L.v.D. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 5:5378 | DOI: 10.1038/ncomms6378 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6378

Recent years have seen a tremendous interest in a class of chemistry known as strain-promoted cycloadditions1. The high mutual reactivity of reaction components, yet

(relative) inertness for a large repertoire of other molecular functionalities, explains the fast-growing popularity of strain-promoted cycloadditions for application in material sciences2, surface functionalization3, bioconjugation4, chemical biology and even in vivo use (known as bioorthogonal chemistry)5,6.

The major focus in the eld of strain-promoted cycloadditions generally dened1 as a bimolecular reaction of a strained system (typically an unsaturated carbocycle) with a complementary component (typically a 1,3-dipole or a diene)has centred around cyclooctyne, the smallest stable cycloalkyne. Its reactivity with 1,3-dipoles, 1,3-dienes, carbenes and other functional groups has been reported throughout the years1, and inspired Bertozzi et al.7 to develop strain-promoted alkyneazide cycloaddition (SPAAC) as a copper-free version of the popular click reaction for application in cellular assays and bioconjugation. Subsequently, a whole family of cyclooctynes has been developed (Fig. 1), mostly with the aim to enhance the relatively slow reaction rate constant of plain cyclooctyne, for example, by introduction of electron-poor uoride substituents as in diuorocyclooctyne (DIFO, 1) (ref. 8), or by dibenzoannulation as in dibenzocyclooctyne (DIBO, 2) (ref. 9). Dibenzo-aza-cyclooctyne (DIBAC, 3) rst developed by us10, also called aza-dibenzocyclooctyne (ADIBO)11 or dibenzocyclooctyne (DBCO)12 displayed signicantly enhanced reactivity and is currently the most broadly applied cyclooctyne for strain-promoted cycloadditions. Yet more reactive probes such as biarylazacyclooctynones (BARAC)13 (4) and tetramethylthiacycloheptyne (TMTH)14 (5) exist, but suffer from poor stability. Bicyclo[6.1.0]non-4yne (BCN, 6) is a stable, non-benzoannulated cyclooctyne that is synthetically readily accessible15, but generally regarded to be less reactive than DIBAC. Finally, Sondheimer diyne (7) is unique in the sense that it is able to react sequentially with two azides with distinctly different reaction rates16. In general it may be concluded that, despite tremendous synthetic effort, optimization of cycloaddition of azides with (stable) cyclooctynes has, depending on specic conditions, leveled off around0.20.5 M 1 s 1 (refs 1,4).

We here wish to report tremendous acceleration of strain-

promoted cycloaddition with cyclooctynes, only by modulation of the azide structure. In particular, it is found that BCN reacts with electron-poor aryl azides up to 29 times faster than with aliphatic azides, with reaction rate ratios reaching 2.02.9 M 1 s 1.

Moreover, a judicious choice of azide enables near-complete tailoring of BCN:DIBAC reaction rate ratios from 0.24 to 425, which is applied in the rst orthogonal SPAAC reported to date. The remarkable difference in rate constants of aliphatic versus benzoannulated cyclooctynes, as well as the pivotal role of

substituent electronegativity are explored and explained by density-functional theory (DFT) calculations17 and electronic-structure analyses18.

ResultsAzidecyclooctyne structurereactivity relationship. Our investigations on the structurereactivity relationship of azides in strain-promoted cycloadditions were stimulated by the fact that much effort has been devoted to the development of more reactive cyclooctynes8,1015, but only scant information is available on the inuence of the azide structure on SPAAC rates19,20. In addition, it struck us that SPAAC reaction rate constants are nearly always determined with an aliphatic azide (typically benzyl azide), but rarely with an aromatic azide. The reason to generally avoid aromatic azides for SPAAC presumably lies in the reported sevenfold reduced reactivity of phenyl azide in comparison with an aliphatic azide, at least in conjunction with DIBAC19. The observation that reaction rates of aromatic azides are hardly inuenced by changing the electronic nature of substituents (as determined for p-methoxy and p-CF3-phenyl azide)20, has provided further ground to avoid aromatic azides for SPAAC. However, because an in-depth structurereactivity relationship of azides with cyclooctynes is lacking to date, we designed a comparative study of aliphatic and aromatic azides in their reactivity with an aliphatic cyclooctyne (BCN) and a dibenzoannulated cyclooctyne (DIBAC).

IR-based reaction rate determination. To facilitate the evaluation of a large number of azides in reaction with cyclooctynes, the availability of a fast and straightforward method to follow conversion rates was desirable. However, the common nuclear magnetic resonance-based method for reaction rate determination is laborious and unsuitable for fast reactions (k41 M 1 s 1) as a consequence of the time-consuming data acquisition.

An alternative method based on ultraviolet21,22 is fast, but is also limited to cases where a difference exists in absorption between substrate(s) and product(s) at a specic wavelength. We reasoned that infrared (IR) detection should enable direct monitoring of substrate-to-product conversion, simply by integration of the distinct azide stretch vibration (around 2,100 cm 1). Indeed, it was found that IR-monitoring of reactions performed in a 9:1 mixture of THF/H2O enables fast and accurate monitoring of azide disappearance, as exemplied for reaction of benzyl azide and BCN alcohol 6a (Fig. 2a,b; Supplementary Fig. 1). From this plot, a reaction rate constant of 0.07 M 1 s 1 could be readily derived (Fig. 2c; Supplementary Methods). Other solvents were also explored, for example, 20% aqueous THF or pure MeOH, but in these cases signal-to-noise ratios were too small to measure azide stretch vibration.

F

F

S

O

DIFO (1)

N N

O

O

DIBO (2)

TMTH (5)

H H

O

DIBAC/DBCO (3) BARAC (4) BCN (6)

Sondheimer diyne (7)

Figure 1 | Known cycloalkynes for bioconjugation. Cyclooctynes employed in strain-promoted alkyneazide cycloaddition (SPAAC) most typically include DIBO (2), DIBAC (3) and BCN (6).

2 NATURE COMMUNICATIONS | 5:5378 | DOI: 10.1038/ncomms6378 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6378 ARTICLE

N

N

N

N3

H H

OH

H H

OH

6a

IR-integral (azide)

kt

10

8

6

4

2

0 0 2,000 3,000 4,000

1,000

300

250

200

150

100

50

0

0 1,000 2,000 4,000

3,000

t (s)

t (s)

kt = 0.075 t

r2 = 0.99475

Figure 2 | Reaction rate determination based on IR. (a) Model SPAAC reaction between benzyl azide and non-functionalized BCN alcohol 6a. (b) Data points collected by integration of N3 stretch vibration of benzyl azide during reaction with BCN alcohol 6a. (c) Determination of the second-order rate constant k for the reaction of benzyl azide and BCN alcohol 6a by linear regression.

Cycloadditions of azides with BCN or DIBAC. Having established a practical protocol for the determination of reaction rate constants, a small range of aliphatic azides (AC, for synthesis see Supplementary Methods 2) was evaluated rst (Fig. 3; Table 1, entries 13, and Supplementary Table 1). Not surprisingly, near-identical rate constants were determined for azides AC in reaction with BCN alcohol (6a) or with a DIBAC derivative (3a). Similar to earlier reported data10,15, an approximately three- to fourfold higher reactivity was observed for DIBAC versus BCN (kBCN:kDIBAC). It must be noted, however, that the absolute rate constants as measured by this protocol are about twofold lower

than those reported earlier in other solvent systems1. One sensible explanation lies in the fact that 1,3-dipolar cycloadditions generally proceed faster in a more polar environment (for example, 50% aqueous CH3CN)23, as for example was also demonstrated by us for strain-promoted nitrone cycloadditions24. Indeed, competition experiments in a 2:1 mixture of CD3CN/D2O (Supplementary

Methods 5) indicate k values (numbers in brackets) that are a factor 1.52.5 higher than those obtained by IR in 9:1 THF/H2O.

The rst unexpected result was obtained while determining the rate constants for phenyl azide (D, entry 4). While the sevenfold drop in reactivity with DIBAC was in line with earlier observations19,20, reaction with BCN unexpectedly proceeded2.9 faster in comparison with benzyl azide (A, entry 1). Thus,

in a side-by-side comparison, reaction of phenyl azide with BCN is in fact 6 faster than with DIBAC, quite opposite to aliphatic

azides. Stimulated by this nding, we were curious to investigate whether cycloaddition of BCN with 2,6-diisopropylphenyl azide (E) would be further accelerated, as reported for Sondheimer diyne 7 (ref. 20). However, in this case, only a modest effect was observed(1.5 faster than unsubstituted phenyl azide, D), while DIBAC, a

dibenzoannulated cyclooctyne, like Sondheimer diyne, showed a large increase in reaction rate (23 faster). Puzzled by these

contradictory results, a range of substituted aryl azides (FN) was subsequently synthesized (Supplementary Methods 2) and evaluated in reaction with BCN and DIBAC (entries 614). It was found that replacing the o-isopropyl groups with a less sterically hindered but more electronegative halogen group (Cl orF) was favourable for the reaction with BCN, but not for DIBAC. The positive effect of an electronegative substituent on reactivity with BCN was also established for p-nitrophenyl azide (H, entry 8), but quite opposite for reaction with DIBAC, thereby reaching a BCN:DIBAC reactivity ratio of a factor 27. Entries 9 and 10 further underline the positive correlation between electronegativity and cyclooaddition reactivity for BCN, but not for DIBAC.

Finally, four particularly electron-poor aryl azides (KN) were evaluated in a reaction with BCN and DIBAC (entries 1114). Much to our delight, N-propyl 4-azido-2,3,5,6-tetrauorobenzamide (L, entry 12) and 4-azido-1-methylpyridinium iodide (N, entry 14) reached reaction rate constants of 1.23 and2.0 M 1 s 1 (or 2.9 M 1 s 1 by extrapolation from the competition experiment, Supplementary Methods 5), respectively, the highest absolute number for SPAAC (with a stable cycloalkyne), as well as relative BCN-to-DIBAC ratio (40

faster) reported to date.

Frontier molecular orbital considerations. The marked increase in reaction rate for electron-poor azides is in sharp contrast to earlier ndings of studies focusing on the inuence of electron-withdrawing substituents on SPAAC, all of which have been suggestive of a HOMOazideLUMOcyclooctyne interaction. For

example, Hosoya et al. noted a 3.9 increase in reaction rate for

p-MeOC6H4N3 (I) versus PhN3 in reaction with 7 (ref. 20). A positive effect on reactivity (1.6 faster) upon cyclooctyne

uorination, as in DIFO8 or BARAC25, or upon DIBAC halogenation (up to 2.2 faster)26 further supports this notion.

Interestingly, the fast reaction of electron-poor azides with an aliphatic cyclooctyne (BCN) as reported here is highly indicative of an alternative frontier orbital interaction. Further experimental support for the latter hypothesis is found in the much faster reaction of cyclooctyne 8 (Supplementary Methods 3) with p-nitrophenyl azide than with benzyl azide (6 faster), while

only a negligible reaction rate difference was noted for the more electron-decient cyclooctyne 9 (Fig. 4).

Our relativistic DFT calculations17 (at ZORA-OLYP/TZ2P, Supplementary Methods 8)2729 conrm the above reactivity

NATURE COMMUNICATIONS | 5:5378 | DOI: 10.1038/ncomms6378 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6378

Aliphatic azides Aromatic azides

N3

N3

Cl

F

N3 N3 N3 N3 N3

N3

A

Cl

F

R

H

I

R = NO2 R = MeO

D

E

F

G

O

O O O M

J

F

n-PrNH

N3 N3

N3

NO2 N N

R

B

C

N3

n-PrNH n-PrNH

HO

+

F

R

I

K

L

R = H

N

R = F

Figure 3 | Azide structure determines reaction rate. Structures of aliphatic (AC) and (substituted) aromatic (DN) azides.

Table 1 | Experimentally determined absolute (M 1 s 1) and relative rate constants (BCN 6a versus DIBAC 3a) for SPAAC with azides A-N.

Entry Azide H

HOH

kBCN:kDIBAC

kBCN k(rel) kDIBAC k(rel)1 A 0.07 1 0.24 1 0.292 B 0.05 0.7 0.20 0.83 0.243 C 0.05 0.7 0.21 0.88 0.244 D 0.20 (0.5)* 2.9 0.033 0.14 65 E 0.30 4.3 0.77 3.2 0.46 F 0.41 6.0 0.28 1.2 1.57 G 0.63 9.1 0.14 0.58 4.58 H 0.53 (1.1)* 7.6 0.02 0.08 279 I 0.18 (0.6)* 2.6 0.06 0.25 3.110 J 0.38 (0.7)* 5.4 0.018 0.08 2111 K 0.73 (1.9)* 11 0.11 0.46 6.612 L 1.23 18 0.16 0.67 7.713 M 0.68 10 0.05 0.21 1414 N 2.0 (2.9)* 29 0.05 0.21 40

SPAAC, Strain-promoted azidealkyne cycloaddition.*Reaction rate constants in brackets are derived from competition experiments in CH CN/H O (see Supplementary Information).

O

N 3a

6a O NH2

trend of the reactions involving BCN (see Supplementary Table 24 for full structural and energy details). Thus, the barrier Ea for the reaction of BCN with A, D, E, F and G decreases monotonically from 16.9 to 14.4 kcal mol 1 as the rate constant increases from 0.07 to 0.63 M 1 s 1 (see Table 2). The formation of the new CN bonds occurs in an asynchronous manner: in the transition state, the bond length d2 involving the aryl- or alkyl-substituted azide nitrogen runs somewhat behind that (d1) involving the other nitrogen atom. All reactions are pronouncedly exothermic with reaction energies Er between 55

and 60 kcal mol 1.

Interestingly, the DFT computations reveal that these fast cycloadditions with electron-poor azides proceed via an inverse electron-demand SPAAC mechanism (IED SPAAC). This is in line with the fact that BCN undergoes extremely fast IED Diels Alder cycloaddition with 1,2,4,5-tetrazines20,30,33, while DIBAC is inert31. Fig. 5 is a schematic representation of the frontier molecular orbital interactions between the BCN and phenyl azide reactants in the transition state that emerge from our quantitative KohnSham electronic-structure analyses. Activation-strain analyses32 conrm that the overall orbital interactions indeed contribute ca 35 kcal mol 1 of stabilization, but that net

interactions between the reactants in the TS are small, at most

0.8 kcal mol 1, due to substantial Pauli repulsion between closed-shell orbitals on each of the two fragments. The dominant

donoracceptor orbital interactions are between the p electrons of BCN and the relatively low-energy LUMO of the azide. The associated inverse HOMOBCNLUMOazide gap of 1.82.4 eV is in all cases signicantly smaller than the normal HOMOazide

LUMOBCN gap of 4.54.9 eV (see Table 2). Besides the in-plane HOMO of BCN and azides, also the corresponding perpendicular

O

NHn-Pr

8

O

F H

N

O

9

NHn-Pr

Figure 4 | Cyclooctyne structures determines reaction rate. Structures of a regular cyclooctyne 8 and an electron-poor cyclooctyne 9 as used in a competition experiment between p-NO2C6H4N3 (H) and benzyl azide (A).

4 NATURE COMMUNICATIONS | 5:5378 | DOI: 10.1038/ncomms6378 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6378 ARTICLE

Table 2 | Quantum-chemically computed key parameters of selected SPAAC reactions for BCN*.

Entry Azide k Ea Er HOMOazideLUMOBCN HOMOBCNLUMOazide d1 d2

1 A 0.07 16.9 60.1 4.90 2.37 2.184 2.397

2 D 0.20 15.9 55.4 4.55 1.93 2.136 2.458

3 E 0.30 15.8 58.8 4.51 2.43 2.179 2.411

4 F 0.41 14.6 57.3 4.96 1.83 2.111 2.509

5 G 0.63 14.4 56.9 4.91 1.76 2.132 2.481

SPAAC, Strain-promoted azidealkyne cycloaddition.*Computed at ZORA-OLYP/TZ2P. Experimental rate constant (k in M 1 s 1), activation and reaction energies (E , E in kcal mol 1), HOMOLUMO gaps (in eV) and newly forming CN bond lengths (d , d in ; d is the CN bond of the aryl or alkyl-substituent azide nitrogen atom).

HOMO

HOMO-1

LUMO

LUMO

HOMO

HOMO-2

HO H

H

BCN TS

R

Phenyl azide

N

N N +

R

Figure 5 | Frontier orbital interactions in SPAAC or IED SPAAC. Schematic representation of the BCNazide frontier molecular orbital interactions emerging from our quantitative KohnSham electronic-structure analyses.

p orbitals located on the BCN triple bond and azide N3 moiety, respectively, contribute to the donoracceptor interactions (see Fig. 5). The inverse electron donation from BCN HOMO and HOMO1 to the azide is signicantly larger, 0.220.26 electrons, than the normal electron donation from azide HOMO and HOMO2 (or HOMO3) to BCN, 0.070.10 electrons (not shown in Table 2). This agrees with the smaller inverse HOMO LUMO gap and conrms the predominant IED character of azide with aliphatic cyclooctynes.

One-pot three-component cycloadditions by orthogonal SPAAC. Finally, the high reactivity of DIBAC and BCN with aliphatic or (electron-poor) aromatic azides, respectively, opens up the unique possibility of orthogonal SPAAC with two azides. First, we subjected 2,6-diuorophenyl azide derivative K to an excess of BCN alcohol 6a and a DIBAC derivative 3b (ve equiv. each, Supplementary Fig. 2) in CD3CN/D2O (3:1). As expected, analysis of the crude reaction mixture with 1H nuclear magnetic

resonance led to the sole detection of the (3 2) cycloadduct

formed by reaction with BCN, but no DIBAC-derived triazole (Supplementary Fig. 10). A further level of selectivity was attained by mixing compound 10 (Fig. 6), bearing one aliphatic and one electron-decient aryl azide, with a twofold excess of BCN alcohol 6a and DIBAC derivative 3c (Supplementary Fig. 3; Supplementary Methods 6), which led to the formation of a mixed BCN,DIBAC-bistriazole (12) as a pure compound in 83% yield after silica gel column chromatography, with only minor amounts (o6%) of double-addition products of 10 with BCN or

DIBAC (Supplementary Fig. 23). The concept of a judicious choice of azidecyclooctyne combination for orthogonal SPAAC was nally extended to a biological context by simultaneous incubation of green uorescent protein (GFP) bearing a genetically encoded BCN33 (Fig. 6) with stoichiometric 10 and DIBACTAMRA (3d, Supplementary Methods 7). The reaction mixture was analysed by protein gel electrophoresis and in-gel uorescence detection (Supplementary Fig. 4). Incubation with 10 led to efcient labelling in only 2 h, while only minor labelling

NATURE COMMUNICATIONS | 5:5378 | DOI: 10.1038/ncomms6378 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6378

N3

H N

O

O2N

O

O

10

11

N

O

TAMRA

N3 N3

N3

O

O

NH

O

O

4

HN

3d

O

GFPBCN + 11 + 3d

GFPBCN + 10 + 3d

t (min)

15

30 60 120 15 30 60 120

M

GFPBCN

R =

3d (1 equiv.) +

10 or 11 (1 equiv.)

25 kDa

50 kDa

25 kDa

N N

N

H N

N

N

R =

R

O2N

O O

N

O

N

O

O

NH 4

TAMRA

O

HN

O

Figure 6 | Judicious choice of cyclooctyne-azide combination enables selective SPAAC. Labeling of GFPBCN (20 mM, 37 C) with stoichiometric DIBACTAMRA (3d) in the presence of a bivalent linker containing either aliphatic and aromatic azide 10 or two aliphatic azides 11 shows signicantly more labeling with linker 10 than with 11.

was observed upon similar incubation with bis-azido-triethylene glycol reagent 11. In this respect, the here presented one-pot sequential SPAAC for labelling of a BCN-bearing protein bears resemblance, yet is nicely complementary, to the procedure reported by Hosoya et al.16 for three-component labelling of an azido-bearing protein in the presence of Sondheimer diyne and an azido-uorophore.

DiscussionIn conclusion, it was established that aliphatic and aromatic azides display a striking difference in reactivity, with strong preference for benzoannulated and aliphatic cyclooctynes, respectively. The latter experimental nding is strongly supported by relativistic DFT calculation, which indicate that aromatic azides and aliphatic cyclooctynes (in particular BCN) react via an IED mechanism, with the dominant donoracceptor orbital interactions involving the p electrons of BCN (HOMO) and the

LUMO of the azide, opposite to what currently is assumed to be the operative mechanism during SPAAC. The latter mode of reactivity enabled the development of highly accelerated SPAAC by means of introduction of electron-withdrawing substituents on the aromatic azide, thereby leading to unprecedented reaction rates (42 M 1 s 1), up to 30-fold faster than traditional

SPAAC. In this respect, use of IED SPAAC also compares favourably with the so-called supersensitive copper-catalysed click chemistry, which can be accelerated only 46 times by including picolyl azides instead of regular benzyl azides34,35. It must be noted that several research groups reported on a sizable effect of electronegative substituents on aryl azides in

cycloaddition with alkenes already 50 years ago23,36,37. For example, an approximate fourfold reaction rate enhancement was reported by Nowack et al. for nitro-substituted versus unsubstituted phenyl azide in reaction with norbornene36. Even more markedly, picryl azide was found to react with norbornene almost 1,000 times faster than phenyl azide, which is even further enhanced in reaction with electron-rich vinyl ethers and enamines as reported by Bailey and White37. The latter phenomenon was investigated in-depth by Huisgen et al.23 and it was proposed that electronegative aryl substituents are able to lower the energy of the (concerted) transition state in cycloaddition with electron-rich olen-like enamines by partial (negative) charge delocalization. However, why such a substituent effect is also operative in cycloaddition with neutral alkenes (for example, norbornene) and why the effect is so much smaller for alkenes with an electron-withdrawing substituent, remains unexplained. The current study rst of all reports absolute reaction rate constants for cycloaddition with BCN that are signicantly faster than those reported for alkenes, for example, a factor 1,000 times higher than for norbornene23. Moreover, a similar Hammett effect in reaction with alkynes was to the best of our knowledge unexplored to date. Furthermore, we have elucidated the bonding mechanism between the reactants in terms of quantitative KohnSham molecular orbital theory based on DFT computations. These analyses link the reactivity trend of our SPAAC reactions involving BCN to the inverse HOMOBCN

LUMOazide gap. We are condent that the latter insight will

provide solid ground for the design of yet faster strain-promoted cycloadditions.

6 NATURE COMMUNICATIONS | 5:5378 | DOI: 10.1038/ncomms6378 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6378 ARTICLE

A practical convenience of the IED reactions reported herein furthermore lies in the fact that BCN and several of the aromatic azides are readily synthesized or commercially available, such as, for example, photoafnity labelling reagents based on azidotetrauorobenzoic acid38. It is clear that the reaction rate constants of the here presented accelerated IED SPAAC still lag far behind those involved with IED cycloadditions of tetrazines.1,4,22,30 However, besides specic in vivo applications that proceed at extremely small concentrations, in the vast majority of circumstances such high reaction rates as pertaining to tetrazine ligation are not a prerequisite. Hence, also given the high popularity of (copper-free) click chemistry, broad application of the here reported accelerated SPAAC may be foreseen, in scientic areas spanning from polymer and materials science, chemical and cell biology and beyond. Furthermore, it is to be expected that with growing understanding of IED SPAAC, even faster probes may be readily designed in forth-coming years. Efforts in this direction, as well as further computational analyses with the aim to more rationally design IED SPAAC, are currently under investigation in our laboratories.

References

1. Debets, M. F. et al. Bioconjugation with strained alkenes and alkynes. Acc. Chem. Res. 44, 805815 (2011).

2. Wilson, J. T., Krishnamurthy, V. R., Cui, W., Qu, Z. & Chaikof, E. L. Non-covalent cell surface engineering with cationic graft copolymers. J. Am. Chem. Soc. 131, 1822818229 (2009).

3. Manova, R., van Beek, T. A. & Zuilhof, H. Surface functionalization by strain-promoted alkyne-azide click reactions. Angew. Chem. Int. Ed. 50, 54285430 (2011).

4. Debets, M. F., van Hest, J. C. C. & Rutjes, F. P. J. T. Bioorthogonal labelling of biomolecules: new functional handles and ligation methods. Org. Biomol. Chem. 11, 64396455 (2013).

5. Sletten, E. M. & Bertozzi, C. R. Bioorthogonal chemistry: shing for selectivity in a sea of functionality. Angew. Chem. Int. Ed. 48, 69746998 (2009).

6. Borrmann, A. & van Hest, J. C. M. Bioorthogonal chemistry in living organisms. Chem. Sci. 5, 21232134 (2014).

7. Agard, N. J., Prescher, J. A. & Bertozzi, C. R. A strain-promoted [3 2]

azide-alkyne cycloaddition for covalent modication of biomolecules in living systems. J. Am. Chem. Soc. 126, 1504615047 (2004).8. Baskin, J. M. et al. Copper-free click chemistry for dynamic in vivo imaging. Proc. Natl Acad. Sci. USA 104, 1679316797 (2007).

9. Ning, X., Guo, J., Wolfert, M. A. & Boons, G.-J. Visualizing metabolically labeled glycoconjugates of living cells by copper-free and fast Huisgen cycloadditions. Angew. Chem. Int. Ed. 47, 22532255 (2008).

10. Debets, M. F. et al. Aza-dibenzocyclooctynes for fast and efcient enzyme PEGylation via copper-free (3 2) cycloaddition. Chem. Commun. 46, 9799

(2010).11. Kuzmin, A., Poloukhtine, A., Wolfert, M. A. & Popik, V. V. Surface functionalization using catalyst-free azide-alkyne cycloaddition. Bioconj. Chem. 21, 20762085 (2010).

12. Campbell-Verduyn, L. S. et al. Strain-promoted copper-free click chemistry for 18F radiolabeling of bombesin. Angew. Chem. Int. Ed. 123, 1131311316 (2011).

13. Jewett, J. C., Sletten, E. M. & Bertozzi, C. R. Rapid Cu-free click chemistry with readily synthesized biarylazacyclooctynones. J. Am. Chem. Soc. 132, 36883690 (2010).

14. de Almeida, G., Sletten, E. M., Nakamura, H., Palaniappan, K. K. & Bertozzi, C.R. Thiacycloalkynes for copper-free click chemistry. Angew. Chem. Int. Ed. 51, 24432447 (2012).15. Dommerholt, J. et al. Readily accessible bicyclononynes for bioorthogonal labeling and three-dimensional imaging of living cells. Angew. Chem. Int. Ed. 49, 94229425 (2010).

16. Kii, I. et al. Strain-promoted double-click reaction for chemical modication of azido-biomolecules. Org. Biomol. Chem. 8, 40514055 (2010).

17. te Velde, G. et al. Chemistry with ADF. J. Comput. Chem. 22, 931967 (2001).18. Bickelhaupt, F. M. & Baerends, E. J. in Reviews in Computational Chemistry Vol. 15 (eds Lipkowitz, K. B. & Boyd, D. B.) 186 (Wiley-VCH, 2000).

19. Zimmerman, E. S. et al. Production of site-specic antibody drug conjugates

using optimized non-natural amino acids in a cell-free expression system. Bioconj. Chem. 25, 351361 (2014).

20. Yoshida, S. et al. Enhanced clickability of doubly sterically-hindered aryl azides. Sci. Rep. 82, 14 (2011).

21. Poloukhtine, A. A., Mbua, N. E., Wolfert, M. A., Boons, G.-J. & Popik, V. P. Selective labeling of living cells by a photo-triggered click reaction. J. Am. Chem. Soc. 131, 1576915776 (2009).

22. Lang, K. et al. Genetic encoding of bicyclononynes and trans-cyclooctenes for site-specic protein labeling in vitro and in live mammalian cells via rapid uorogenic Diels Alder reactions. J. Am. Chem. Soc. 134, 1031710320

(2012).23. Huisgen, H., Szeimies, G. & Mbius, L. Kinetik der additionen organischer azide an CC-Mehrfach-bindungen. Chem. Ber. 100, 24942507 (1967).

24. Ning, X. et al. Protein modication by strain-promoted alkynenitrone cycloaddition. Angew. Chem. Int. Ed. 49, 30653068 (2010).

25. Gordon, C. G. et al. Reactivity of biarylazacyclooctynones in copper-free click chemistry. J. Am. Chem. Soc. 134, 91999208 (2012).

26. Debets, M. F. et al. Synthesis of DIBAC analogues with excellent SPAAC rate constants. Org. Biomol. Chem. 12, 50315037 (2014).

27. Handy, N. C. & Cohen, A. J. A dynamical correlation functional. J. Chem. Phys. 116, 5411 (2002).

28. Lee, C., Yang, W. & Parr, R. G. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev.B. 37, 785789 (1988).29. van Lenthe, E., van Leeuwen, R., Baerends, E. J. & Snijders, J. G. Relativistic regular two-component Hamiltonians. Int. J. Quantum Chem. 57, 281293 (1996).

30. Chen, W., Wang, D., Dai, C., Hamelberg, D. & Wang, B. Clicking 1,2,4,5-tetrazine and cyclooctynes with tunable reaction rates. Chem. Commun. 48, 17361738 (2012).

31. Karver, M. R., Weissleder, R. & Hilderbrand, S. A. Bioorthogonal reaction pairs enable simultaneous, selective, multi-target imaging. Angew. Chem. Int. Ed. 51, 920922 (2012).

32. Fernandez, I. & Bickelhaupt, F. M. The activation strain model and molecular orbital theory: understanding and designing chemical reactions. Chem. Soc. Rev. 43, 49534967 (2014).

33. Borrmann, A. et al. Genetic encoding of a bicyclo[6.1.0]nonyne-charged amino acid enables fast cellular protein imaging by metal-free ligation. Chembiochem 13, 20942099 (2012).

34. Uttamapinant, C. et al. Fast, cell-compatible click chemistry with copper-chelating azides for biomolecular labeling. Angew. Chem. Int. Ed. 51, 58525856 (2012).

35. Jiang, H. et al. Monitoring dynamic glycosylation in vivo using supersensitive click chemistry. Bioconj. Chem. 25, 698706 (2014).

36. Scheiner, P., Schomaker, J. H., Deming, S., Libbey, W. J. & Nowack, G. P. The addition of aryl azides to norbornene. A kinetic investigation. J. Am. Chem. Soc. 87, 306311 (1965).

37. Bailey, A. S. & White, J. E. Reaction of picryl azide with olens. J. Chem. Soc. B 819822 (1966).

38. Liu, L. & Yang, M. Peruorophenyl azides: new applications in surface functionalization and nanomaterial synthesis. Acc. Chem. Res. 43, 14341443 (2010).

Acknowledgements

This work was supported by the National Research School CombinationCatalysis and the Netherlands Organization for Scientic Research (NWO-EW and NWO-CW).

Author contributions

J.D., A.B. and F.L.v.D. conceived and designed the experiments and analysed the data. O.v.R., C.F.G. and F.M.B. conceived and designed the DFT calculations and analysed the data. J.D. performed the kinetic experiments and A.B. performed the protein experiments. O.v.R. performed the DFT calculations. F.L.v.D. and F.M.B. wrote the paper, and all the authors edited and commented on the manuscript.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: F.L.v.D. is the CSO of SynAfx BV. The remaining authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Dommerholt, J. et al. Highly accelerated inverse electron-demand cycloaddition of electron-decient azides with aliphatic cyclooctynes.

Nat. Commun. 5:5378 doi: 10.1038/ncomms6378 (2014).

NATURE COMMUNICATIONS | 5:5378 | DOI: 10.1038/ncomms6378 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

& 2014 Macmillan Publishers Limited. All rights reserved.