Open Access Article
Nicholas C.
Rose
ab,
Anaïs V.
Sanchez
ab,
Eve F.
Tipple
ab,
Jason M.
Lynam
a and
Christopher D.
Spicer
*ab
aDepartment of Chemistry, University of York, Heslington, YO10 5DD, UK. E-mail: chris.spicer@york.ac.uk
bYork Biomedical Research Institute, University of York, Heslington, YO10 5DD, UK
First published on 14th October 2022
Ortho-boronoaldehydes react with amine-based nucleophiles with dramatically increased rates and product stabilities, relative to unfunctionalised benzaldehydes, leading to exciting applications across biological and material chemistry. We have developed a novel Förster resonance energy transfer (FRET)-based assay to provide key new insights into the reactivity of these boronoaldehydes, allowing us to track conjugation with unprecedented sensitivity and accuracy under standardised conditions. Our results highlight the key role played by reaction pH, buffer additives, and boronoaldehyde structure in controlling conjugation speed and stability, providing design criteria for further innovations and applications in the field.
The ortho-boronic acid accelerates the rate of oBID formation (k1), relative to unfunctionalised benzaldehyde, by both activating the aldehyde to nucleophilic attack and accelerating the rate-determining dehydration step. Though the rate of hydrolysis (k−1) is also accelerated, stabilising B–N interactions lead to a significant overall shift in equilibrium towards product formation (decrease in dissociation constant, Kd). For example, while benzaldehyde reacts with alkyl amines with k1 ∼ 0.1 M−1 s−1 and Kd ∼ 300 M,4 the ortho-borono analogue ortho-formylphenylboronic acid (FBPA) forms analogous iminoboronates with a k1 ∼ 1000 M−1 s−1 and Kd ∼ 10 mM (Fig. 1a).5 This 5-order of magnitude increase in reaction rate, and 4-order of magnitude decrease in dissociation constant, typifies the unique reactivity of oBAs.
Through variation of the oBA and nucleophile coupling partners, Kds spanning 10−2–10−9 M have been reported.6 When coupled with high biocompatibility and potential for stimuli-responsive behaviour, this tunability greatly enhances the potential for biomedical applications of oBIDs. Despite this potential, much of our understanding of oBID chemistry comes from simple, unsubstituted model compounds such as FPBA, which poorly reflect the stereo-electronic characteristics of the functionalised oBAs necessary for translational applications.
Moreover, much of this understanding is pieced together from a large number of independent studies, performed and analysed under varying conditions, making comparison challenging. This is particularly important given the highly dynamic nature of oBID linkages, leading to a high sensitivity to environmental conditions and choice of analytical technique. Ultimately, this leads to significant discrepancies in reported rates of formation and stabilities of oBIDs, due to subtle differences in substrate choice, reaction conditions, or analysis method.7,8
To address this, we herein describe a versatile and sensitive FRET-assay for studying oBID formation, stability, and cleavage. This assay is highly tolerant of environmental conditions, allowing us to provide critical new insights into the effects of pH, additives, media, and oBA structure on conjugation. We therefore anticipate that this assay will find future use for the high-throughput screening of novel oBIDs and their stimuli-responsive behaviour.
(i) The platform should be operational at low concentrations. Many oBIDs form with k1 > 103 M−1 s−1, and at the substrate concentrations required for UV-vis/NMR analysis reactions are complete prior to the first measurements, preventing accurate quantification.9–11 Sensitive detection of oBID formation at low concentrations is therefore essential to allow a suitable number of time points to be recorded before equilibrium is reached;
(ii) The platform should function across a broad range of reaction conditions representative of the end applications, including the presence of complex additives, allowing us to study oBID formation under conditions relevant to biological applications;
And (iii) a readout of conjugation, rather than structure, should be provided. Following conjugation of an oBA and nucleophile, the oBID formed remains highly dynamic. Often, complex mixtures of rapidly interchanging species exist in equilibrium, in a manner that is highly sensitive to environmental conditions (e.g. pH dependent equilibrium between ortho-boronohydrazones and diazoborines;12 diastereomeric mixtures formed during thiazolidino boronate formation8). By focussing on conjugation, we would be able to gain an overall view of all species that contribute positively to oBID formation, which would not be provided by focussing on the formation of specific structures. UV-vis and NMR studies are again poorly suited to providing such an overview.
We reasoned that a FRET-based reporter system would fulfil these criteria (Fig. 1b).13 Upon conjugation of complementary dyes bearing either oBA or nucleophile reactive handles, the proximity of the dyes would allow energy transfer to take place under donor excitation. The donor
:
acceptor emission ratio (FRET ratio) would then directly correlate to reaction conversion, and would be largely independent of oBID structure. Notably, Schmidt et al. previously exploited a fluorescence quenching assay to monitor borono-oxime formation, but such experiments can be prone to interference from environmental conditions. In contrast, the ratiometric nature of FRET minimises such complications, and provides an extra level of robustness to kinetic measurements.14
The cyanine dyes Cy3 and Cy5 were chosen as a suitable FRET pair due to their chemical-, thermal-, and photo-stability, emission across a wide pH range, and known ability to undergo energy transfer.15 Due to the strong distance dependence of FRET (1/r6), we aimed to minimise the linker length from the dyes to the oBA/nucleophile functional handles, leading us to synthesise novel amino-cyanines 1 and 2 (Scheme 1). These derivatives, bearing a pendant sulfonate group that enhanced water solubility without affecting fluorescence,15 could be easily derivatised via amide coupling. Phthalimide protection of the amine (3 and 4) was required due to the harsh conditions needed for synthesis of the dye cores. Though traditional hydrazinolysis of the phthalimide was found to be problematic, due to the reaction of hydrazine with the conjugated backbone of Cy5, deprotection could be successfully achieved with high concentrations of methylamine to deliver 1 and 2 on a gram scale.
With the parent dyes in hand, they were functionalised with reactive handles for oBID formation. A 2-borono-4-alkoxy-benzaldehyde handle (5) was chosen as a representative oBA for attachment to Cy3 1. This choice was based on the ease with which 5 could be synthesised, relative to other oBA analogues, and therefore its position as a widely used functional oBA within the community.16–18 However, as discussed later, this choice of oBA plays a role in dictating conjugation rate and stability, and therefore has significant implications.
As nucleophiles, Cy5 analogues 6–10 were synthesised to represent the most common reagents for oBID formation. These analogues encompass the archetypal amine nucleophile (6) associated with iminoboronate formation (11),19 1,2-diamine (7)18 and 1,2-aminothiol (8)8 based nucleophiles which form cyclic imidazolidino-(IzB, 12) and thiazolidino-(TzB, 13) boronates respectively, and the ‘α-effect’ hydroxylamine (9)14 and hydrazine (10)20 handles that react to form stabilised oximes (14) and hydrazones/diazoborines (DAB, 15). A consistent linker structure was used to minimise any distance-dependent changes in FRET efficiency.
:
Cy5 emission ratios to conversions, reactions using an acetyl-capped Cy5-amine 16 (negative ‘0% conjugation’ control) and a covalently linked Cy3-Cy5 pair 17 (positive ‘100% conjugation’ control) were used as references (Fig. 2b and Table 1, see ESI Section 9† for details of data processing).
| pH | k 1/M−1 s−1 | k −1/s−1 | K d | |
|---|---|---|---|---|
| 11 | 6 | — | — | — |
| 12 | 6 | — | — | — |
| 13 | 6 | 250 ± 3 | 1.13 ± 0.04 × 10−3 | 4.5 μM |
| 14 | 6 | 3034 ± 27 | 6.04 ± 1.08 × 10−5 | 19 nM |
| 15 | 6 | 22 920 ± 839 |
8.07 ± 0.29 × 10−6 | 352 pM |
| 11 | 7.4 | — | — | — |
| 12 | 7.4 | 9882 ± 978 | 4.02 ± 0.45 × 10−2 | 4.3 μM |
| 13 | 7.4 | 613 ± 16 | 3.01 ± 0.11 × 10−3 | 4.9 μM |
| 14 | 7.4 | 4370 ± 164 | 5.61 ± 1.28 × 10−5 | 12 nM |
| 15 | 7.4 | 169 030 ± 6962 |
2.60 ± 0.55 × 10−6 | 15 pM |
| 11 | 8 | — | — | — |
| 12 | 8 | 17 521 ± 1211 |
2.72 × 10−2 | 1.5 μM |
| 13 | 8 | 21 261 ± 591 |
1.22 × 10−3 | 57 nM |
| 14 | 8 | 1528 ± 10 | 3.56 ± 0.11 × 10−4 | 233 nM |
| 15 | 8 | 37 471 ± 1238 |
8.05 ± 0.26 × 10−6 | 215 pM |
| 14 | 4 | 1517 ± 81 | 5.23 ± 0.56 × 10−5 | 34 nM |
| 15 | 4 | 876 ± 163 | 3.36 ± 0.70 × 10−2 | 38 μM |
As expected given reports of Kd ∼ 10 mM, iminoboronate 11 was not found to form to an appreciable extent at the low concentrations of these experiments.6 However, for the formation of IzB 12 and TzB 13 the data was found to fit a reversible second-order kinetic model, allowing us to extrapolate k1 and k−1, and thus calculate Kd.23 Interestingly, though the overall position of equilibrium for each reaction was similar (Kd ∼ 4–5 μM), IzB was found to form at an order of magnitude higher rate (∼9900 M−1 s−1vs. 610 M−1 s−1), highlighting thiol cyclisation as rate-limiting for the formation of a stable conjugate. Though our results are in line with previous reports on TzB formation,5 the rate of IzB formation was found to be accelerated by a factor of 10 relative to literature precedent.18 This difference can be attributed to the use of a functionalised 2-borono-4-alkoxy-benzaldehyde handle within our studies, in contrast to the use of FPBA models in the previous reports. This result therefore highlights the critical importance of studying derivatised oBIDs which more accurately reflect the substrates used in real-world applications. Moreover, our measurements were insensitive to the mixtures of diastereomers which complicate NMR analysis of IzB and TzB formation, providing further benefit over previous analyses.18,22
The formations of oxime 14 and hydrazone/DAB 15 were found to fit a second-order irreversible model, and so k1 was calculated from plots of inverse reagent concentration against time. Due to the fast kinetics with which 15 was formed, reactions were performed at high dilutions (750 nM) to provide sufficient data points within the linear region of the plot (See ESI Fig. S3†). Attempts to cleave 14 or 15 using an excess of a competitive nucleophile were unsuccessful, due to slow side-reactions with the conjugated framework of the dyes which became significant over the long timeframes of the equilibration. We therefore monitored the equilibration between propyl-amine functionalised oxime 18 or DAB 19 and an excess of a competitive nucleophile (20 or 21) via liquid chromatography analysis, allowing us to determine the rate of hydroxylamine/hydrazine exchange and thus calculate k−1 and Kd (Scheme 2, see ESI Fig. S4† for data).
The rate of oxime formation and dissociation constant (k1 ∼ 5000 M−1 s−1; Kd ∼ 10 nM) is in line with previous reports by Schmidt et al. (k1 ∼ 11
000 M−1 s−1; Kd ∼ 25 nM). Interestingly, this previous work studied the reactivity of 2-borono-5-alkoxy benzaldehydes, in contrast to the 4-alkoxy analogues used in this work, again highlighting the contributions of stereo-electronic arrangements to oBA reactivity.14
The reactivity of hydrazines with FPBA has been studied extensively by the Bane and Gillingham groups. They have shown that following the formation of an intermediate hydrazone, a second intramolecular cyclisation and dehydration step results in the eventual formation of an aromatic DAB at neutral pH. Though the rates of DAB formation have been determined (∼0.01 s−1), previous studies have had insufficient kinetic resolution to provide insight into the initial hydrazone formation.10,12,20 Our results indicate that this step is remarkably rapid, with k1 ∼ 105 M−1 s−1. Within the timeframe of our FRET experiments, it is unlikely that DAB formation makes a significant contribution, but would not be expected to significantly affect our FRET results, which report on conjugation not structure, regardless. This result again highlights the benefit of analysing all conjugated products collectively, rather than individual species which may otherwise be undergoing exchange or serving as intermediates towards an eventual thermodynamic product. Equilibration studies over the course of 1 week indicated the reaction remained dynamic, with 19 undergoing hydrolysis with a rate of ∼10−6 s−1 leading to a Kd for hydrazone/DAB formation of ∼10 pM (see ESI Fig. S4†). To the best of our knowledge, this represents the lowest Kd recorded for a dynamic oBID to date, in contrast to the nopoldiol-stabilised thiosemicarbazone reagents developed by the Hall group which show no reversibility.24
500 M−1 s−1), with a corresponding small increase in stability, but at pH 6 no conversion was observed. Accounting for error and noise within our calculations, this means Kd > 45 μM for IzB formation under such conditions. Notably, Li et al. previously reported that the carboxylic acid of diaminopropionic acid was able to cyclise to form a mixed boron-carbon anhydride structure which was favoured at pH 6 and stabilised IzB formation. However, when the 1,2-diamine is conjugated to a functional cargo, as in our system, no such interactions are possible demonstrating the importance of studying models relevant to functional applications.18
A similar increase in reaction was observed for the formation of TzB 13 at pH 8. In contrast to the results obtained at pH 7.4, TzB formation was actually faster than IzB formation at this pH (k1 ∼ 37
000 M−1 s−1vs. 21
000 M−1 s−1), presumably due to the accelerated cyclisation of the deprotonated thiolate that can form at higher pH. The stability of the TzB was also increased by 2-orders of magnitude (Kd ∼ 50 nM). At pH 6, the rate of formation and stability of 13 was comparable to the results obtained at pH 7.4, lending further support to the hypothesis that the protonation state of the thiol plays a key role in TzB formation.
For the formation of oxime 14, pH 7.4 was found to provide optimal stability within the conditions tested (Kd ∼ 10 nM). At pH 6 or 8, k1 was reduced slightly, but contributions from accelerated oxime hydrolysis to Kd were far more significant leading to Kds ∼ 20 nM and 200 nM at pH 6 and 8 respectively. These results are in line with recent work by Han and Domaille,25 as well as a prior study by Schmidt et al.,14 showing that borono-oxime stability is highly pH dependent. Interestingly, at pH 4 Kd ∼ 34 nM, indicating that acidic pH did not have a large effect on borono-oxime stability.
At both pH 6 and 8, the rate of reaction between oBA 5 and hydrazine 10 was reduced relative to the reaction at pH 7.4, with a corresponding increase in k−1 leading to an overall increase in Kd from ∼10 pM at pH 7.4 to ∼300 pM and ∼200 pM at pH 6 and 8 respectively. Overall, the reactions were still found to proceed very rapidly at both pHs though, with k1 22
000–37000 M−1 s−1 necessitating measurements at lowered concentrations (750 nM) to allow sufficient data points for analysis. At pH 4, the rate of formation of 15 was greatly decreased (∼850 M−1 s−1), with conjugate stability showing an even more dramatic 6-orders of magnitude drop (Kd ∼ 40 μM). It is noteworthy that at pH 4 Gu et al. reported that for the reaction of oBAs with related hydrazides, the open chain hydrazone oBID is dominant over the cyclic DAB.12 However, we found that the products formed from the hydrazine nucleophiles used in this study existed solely as the ring-closed DAB, with no evidence of the ring-open hydrazone by NMR, even at pH 4 (see ESI Section 11†). A switch between ring-closed and ring-open structure is therefore not responsible for the altered stability at lower pH.
| Additive | k 1/M−1 s−1 | k −1/s−1 | K d | |
|---|---|---|---|---|
| 14 | — | 4370 ± 164 | — | — |
| 14 | Glucose | 1340 ± 10 | 2.04 ± 0.12 × 10−4 | 152 nM |
| 14 | Fructose | 446 ± 4 | 2.54 ± 0.20 × 10−4 | 571 nM |
| 14 | 10% Serum | 64.0 ± 4.4 | 9.03 ± 2.22 × 10−4 | 14.1 μM |
| 15 | — | 169 030 ± 6962 |
2.60 ± 0.55 × 10−6 | 15 pM |
| 15 | Glucose | 45 750 ± 6614 |
2.36 ± 0.50 × 10−6 | 51 pM |
| 15 | Fructose | 16 549 ± 1338 |
2.35 ± 0.50 × 10−6 | 142 pM |
| 15 | 10% Serum | 153 ± 6 | 2.93 ± 0.20 × 10−3 | 19.2 μM |
In the case of hydrazones/DABs, the decrease in k1 led to an overall small increase in Kd (50–150 pM, vs. 15 pM in the absence of additives), despite there being no effect of either sugar on k−1. This may suggest that monosaccharides do not complex hydrazones/DABs once formed, nor does pre-complexation lead to a stable tri-component complex, as such scenarios would lead to differing k−1 values across the three samples. In contrast, the stability of oxime 14 was significantly decreased (Kd ∼ 100–500 nM), in contrast to previous reports that suggest minimal influence of sugar additives on the stability of pre-formed oximes.14 There are two possible origins for this effect: (i) The increased sensitivity of our platform may have given insight into reversibility that was not previously observed due to the use of high reagent concentrations; or (ii) The presence of glucose/fructose during the reaction altered the structure of the product formed. For example, it is possible that sugar complexation after oxime formation is no longer possible/favoured, and so tests of stability could give contrasting results to reactions that study formation in the presence of the additive. To investigate this, oxime 18 or DAB 19 were incubated at pH 7.4 with either glucose or fructose at a concentration of 9 mM and binding was assessed via NMR. In all cases, no binding was observed, suggesting sugar complexation is not a major contributing factor to oxime or DAB stability at concentrations relevant to this study (see ESI Section 12†). Further investigations are currently ongoing in our lab to study the interplay between oBA-sugar binding and reaction kinetics, and to elucidate the effects of other relevant additives on oBID formation, structure, and stability, exploiting the platform reported here.
To further demonstrate the benefits of our approach, we next studied the formation and stability of oxime and hydrazone/DAB constructs in serum (Table 2, see ESI Fig. S6†). Understanding oBID chemistry in this environment is essential for future in vivo applications, particularly for potential within drug delivery to be realised. However, studying conjugation chemistry within serum is highly challenging due to its complex composition. The limited number of previous studies into oBID formation have relied on high concentrations of reagents that may not truly reflect the end application.11,14 In contrast, as background fluorescence within the 520–700 nm window is minimal and can be accounted for, FRET provides us with an ideal means to monitor oBID formation in serum at low μM concentrations.
The rate of oxime formation was slowed significantly in the presence of 10% bovine serum (k1 ∼ 60 M−1 s−1vs. 4000 M−1 s−1 in PBS), and the conjugates were also found to be highly dynamic, with a Kd ∼ 15 μM being observed. The decrease in reaction rate was even more significant for hydrazine 15, with a three order of magnitude reduction in k1 being recorded (∼150 M−1 s−1). Furthermore, conjugation was also found to be dynamic with a Kd ∼ 20 μM. Since the rate of DAB formation is significantly slower than initial hydrazone formation, this data may indicate that this intermediate species is dynamic in the presence of serum, with implications for oBID formation in complex medias. In contrast, when DAB 15 was pre-formed in PBS prior to dilution in 10% serum, the FRET ratio was found to stay unchanged over a period of 25 min (see ESI Fig. S7†). This supports previous reports that DABs are stable in serum, with exciting opportunities for future in vivo applications. Interestingly, oximes pre-formed in PBS were also observed to undergo minimal change in FRET ratio upon addition of 10% serum, pointing towards a potential role for serum complexation to the oBA precursors, leading to the observed slow formation of 14.
Cumulatively, these results highlight the potential of our platform to give far greater insight into oBID chemistry under relevant conditions to their translational applications, and efforts to explore such factors in high throughput are currently underway in our group.
Reactions were run at a 1
:
1
:
1 ratio of 4-alkoxy oBA (25), oBA analogue, and the relevant nucleophile, with product ratios used to determine reaction parameters (Fig. 4). When 1,2-diamine 26 was used as a nucleophile, the rapid equilibration of the system (i.e. high k1 AND k−1) allowed us to calculate relative Kds, as described in the SI. The results obtained show a clear trend of increasing oBID stability as the substituent at the 4-position moved from electron donating (25) through to electron withdrawing (24). Unfortunately, efforts to perform analogous experiments using a hydrazine nucleophile as a kinetic trap proved unsuccessful, with the rapid rates of reaction leading insufficient resolution to distinguish differences in product distribution.
The intriguing insights provided by these studies of relative oBID stability open the door to the design of new oBA derivatives with tuneable properties. For example, rational design of highly electron-deficient oBAs may lead to the production of oBIDs with greatly increased stability relative to the current state of the art. Further efforts to explore these possibilities are currently being undertaken.
Using this platform, we were able to calculate k1,k−1, and Kd for the most widely used oBA-nucleophile pairs used in the literature across a range of conditions, allowing us to identify optimal pH ranges and interesting behaviour that provides new insight into oBID chemistry. The ability of our platform to function in the presence of buffer additives is particularly powerful, allowing us to study conjugation chemistry in complex environments such as serum. We therefore provide insight into the rate of formation of ortho-borono oximes and hydrazones/DABs for the first time. Our results also suggest that further tuning and optimization of oBID chemistry is possible through even small changes to the stereo-electronic configurations of the oBA and nucleophile precursors, and our platform provides an ideal means to test such structures. Efforts to do so, and deliver new functional methods relevant to bioconjugation and biomaterials chemistry, are currently underway in our lab, and we hope that our study will serve to stimulate further exploration and optimisation of the underlying chemistry of oBIDs in the future.
Footnote |
| † Electronic supplementary information (ESI) available: The supporting information contains all experimental details, including reagent synthesis, FRET measurements, and studies of reversibility. See DOI: https://doi.org/10.1039/d2sc04574e |
| This journal is © The Royal Society of Chemistry 2022 |