Mégane
Debiais
a,
Alejandro
Gimenez Molina
a,
Sabine
Müller
b,
Jean-Jacques
Vasseur
a,
Ivan
Barvik
c,
Carine
Baraguey
*a and
Michael
Smietana
*a
aInstitut des Biomolécules Max Mousseron, Université de Montpellier, CNRS, ENSCM, 1919 route de Mende, 34095 Montpellier, France. E-mail: carine.baraguey@umontpellier.fr; michael.smietana@umontpellier.fr
bUniversity Greifswald, Institute for Biochemistry, Greifswald, Germany
cInstitute of Physics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 5, Prague 2, 121 16, Czech Republic
First published on 16th March 2022
Inspired by the ability of boronic acids to bind with compounds containing diol moieties, we envisioned the formation in solution of boronate ester-based macrocycles by the head-to-tail assembly of a nucleosidic precursor that contains both a boronic acid and the natural 2′,3′-diol of ribose. DOSY NMR spectroscopy experiments in water and anhydrous DMF revealed the dynamic assembly of this precursor into dimeric and trimeric macrocycles in a concentration-dependent fashion as well as the reversibility of the self-assembly process. NMR experimental values and quantum mechanics calculations provided further insight into the sugar pucker conformation profile of these macrocycles.
Based on the formation of reversible covalent bonds, the synthesis of macrocycles using dynamic covalent chemistry is characterized by its ability to be controlled by changing experimental conditions and external stimuli. The arsenal of reversible reactions is vast and has been described in detail in dedicated recent reviews.23–27 Among these reactions, boronic acid–diol interactions leading to the formation of stable cyclic boronates have received increasing attention for the elaboration of macrocycles and cages.28–34
Traditionally, approaches to form macrocycles via boronic ester formation involve to remove water from a mixture composed of a functionalized mono- or di-boronic acid and a polyol.29,35,36 By contrast, the construction of boronate-based macrocycles using head-to-tail assembly of a molecule carrying both a boronic acid and a diol moiety are rather scarce.37,38
Our laboratory has long been interested in the study of DNA-templated formation of reversible boronate internucleoside linkages generated from two strands, one having the natural cis-diol functions of ribonucleotides at its 3′-end, while the second one was substituted by a boronic acid at its 5′-extremity.39–42 Whereas earlier studies relied on the incorporation of 2′-deoxyborononucleotide phosphoramidites into DNA sequences,43–45 we recently extended the formation of boronate internucleoside linkages to RNA sequences.46 This was achieved through the synthesis of a 6′-boronoribonucleotidic 2′-O-pivaloyloxymethyl (O-PivOM) phosphoramidite building block 2 prepared in 8 steps starting from commercially available 2′-O-pivaloyloxymethyl-5′-O-dimethoxytrityluridine 1 (Fig. 1A). These results led us to consider the synthesis of a 6′-boronic acid uridine analogue 3 functionalized by both a cis-diol moiety and a boronic acid and to evaluate its ability to self-assemble in a head-to-tail (htt) fashion to generate boronate-linked macrocyclic nucleotides (Fig. 1B).
![]() | ||
Fig. 1 Schematic illustration of boronic acid-mediated assembly of nucleic acids (DMTr: dimethoxytrityl, U = uracyl). |
Decreasing the pD back to 6.4 restored the initial spectrum corresponding to the monomer, thus excluding the observed changes to result from any degradation and demonstrating the reversibility of the process. While these results suggest the reversible self-assembly of 3 in basic aqueous media, detailed insights into the structure were difficult to determine due to the broad spectral features observed.
Compound 3 was thus studied in anhydrous DMF-d7, a solvent expected to promote the ligation between the boronic acid and the 2′-3′ diol moieties. The 1H NMR spectrum of 3 at 37 mM exhibited two datasets in a 2:
1 ratio indicating the presence of two species (bottom spectra Fig. 3 and ESI Fig. S13†). Selective saturation experiments (data not shown) confirmed the existence of a slow exchange at the NMR time-scale between the two species since saturation of a signal belonging to the major species affected the corresponding resonance in the minor species.
![]() | ||
Fig. 3 1H NMR spectra of 3 in anhydrous DMF-d7 with increasing quantities of water (for clarity, parts of the spectra free of signal were removed). |
1H NMR experiments revealed that neither the 2′ and 3′ hydroxyl protons, nor those from the boronic acid function could be observed. Moreover, the boron-bound methylene group at the 6′ position presented a multiplicity more complex than the triplet resulting from the free rotation of the boronic extremity as observed in D2O. Few aliquots of H2O were then gradually added to the sample. As expected, this resulted in the disappearance of the two signals around 5.20 and 4.80 ppm assigned to H2′ and H3′ respectively, that were progressively replaced by two signals at higher field corresponding to the same protons in the free form.
At the same time, three new signals characteristic of the exchangeable hydroxyl protons merged at 5.21, 5.52 and 7.70 ppm (Fig. 3 upper part and Table 1). After addition of 20 μL of water (4%v/v) a unique data set corresponding to the monomer was observed.
Atom | Anhydrous DMF-d7 | DMF-d7/H2Oa | Δδb | |
---|---|---|---|---|
Major speciesc | Minor species | (3) | ||
a H2O: 4%. b Difference (in ppm) between the mean value of the chemical shifts of the 2 species existing in anhydrous DMF-d7 (water < 0.05%) and the chemical shift of 3 observed in wet DMF-d7 (water 4%). c Major species population estimated to 73% from the H6 peaks integration and taking into account the di- and trimeric nature of the 2 structures (vide infra). | ||||
H1′ | 5.88 | 5.86 | 5.89 | −0.02 |
H2′ | 5.19 | 5.18 | 4.20 | +1.02 |
H3′ | 4.81 | 4.80 | 3.92 | +0.89 |
H4′ | 4.05 | 3.90 | 3.82 | +0.16 |
H5′/5′′ | 1.93 | 1.91 | 1.81/1.74 | |
H6′/6′′ | 1.10–0.95 | 1.10–0.95 | 0.84/0.78 | +0.20 |
OH2′ | — | — | 5.52 | — |
OH3′ | — | — | 5.21 | — |
H5 | 5.71 | 5.72 | 5.71 | +0.01 |
H6 | 7.85 | 7.81 | 7.71 | +0.12 |
NH | 11.40 | 11.40 | 11.33 | +0.07 |
B(OH)2 | — | — | 7.70 | — |
Additional 1H NMR spectra were recorded in anhydrous DMF-d7 for diluted samples, showing the concentration dependence of the two species distribution. Lowering the concentration led to a gradual decrease of the minor species signals. The ratio of the integrations of H6 peaks for the major and minor species dropped from 1.7:
1 at 37 mM to 6.7
:
1 at 1.4 mm. At 0.14 mM, only the major species could be detected. No significant chemical shift change was noticed in this concentration range (0.14 to 37 mM), excluding any self-association due to intermolecular π-stacking interactions (ESI, Fig. S13†).
11B NMR was also investigated and two spectra were acquired under strictly identical conditions for the same sample before and after addition of water (Fig. S12†). The spectrum acquired in DMF-d7 containing 20 μL of water showed a single and relatively narrow peak at 31.9 ppm characteristic of the free boronic acid with a trigonal geometry around the boron atom. In anhydrous DMF-d7, a very broad band ranging between 20 and 50 ppm was observed. No significant chemical shift variation could be reported, the signal linewidth in the bound state preventing to detect the chemical shift with accuracy, a characteristic that is well established for 11B atom.48
The symmetry loss around the boron center as well as the slower molecular tumbling, both consecutive to the boron-ribose ligation might explain the flattened signal observed in anhydrous DMF. 1H and 13C resonances of the two species coexisting in anhydrous DMF-d7 were individually assigned and compared to the ones of the monomer 3 generated by addition of water (Tables 1 and 2 and ESI Fig. S14†). Under anhydrous conditions, H2′ and H3′ protons were both down-field shifted by about 1 ppm (Table 1) compared to the monomer, while the corresponding carbons were shifted by more than 10 ppm (Table 2). Numerous studies on boronic acids complexation with sugars reported that boronate formation leads to a significant down-field shift of both the 1H and 13C resonances of the boron-bound diol moiety in both water and organic solvent.49,50 Consistent with the literature, these observations showed that the 2′ and 3′ oxygen atoms are bound to boron in both of the two structures coexisting in anhydrous DMF-d7. The comparison of the signals in anhydrous DMF-d7 indicated that the two species were structurally close (Tables 1 and 2).
Atom | Anhydrous DMF-d7 | DMF-d7/H2Oa | Δδb | |
---|---|---|---|---|
Major species | Minor species | (3) | ||
a H2O: 4%. b Difference (in ppm) between the mean value of the chemical shifts of the 2 species existing in anhydrous DMF-d7 (water < 0.05%) and the chemical shift of 3 observed in wet DMF-d7 (water 4%). | ||||
C1′ | 92.4 | 92.9 | 88.8 | +3.9 |
C2′ | 84.9 | 85.3 | 73.9 | +11.2 |
C3′ | 84.6 | 83.3 | 73.5 | +10.5 |
C4′ | 87.4 | 88.1 | 86.1 | +1.7 |
C5′ | 27.7 | 27.3 | 28.7 | −1.2 |
C6′ | 6.3 | 6.2 | 11.3 | −5.1 |
C2 | 151.0 | 150.9 | 151.2 | −0.3 |
C4 | 163.6 | 163.6 | 163.6 | 0 |
C5 | 102.6 | 102.6 | 102.2 | +0.4 |
C6 | 142.9 | 143.2 | 141.2 | +1.9 |
In 1H NMR, the strongest variation between the 2 assemblies was found for H4′ (4.05 vs. 3.90 ppm, Table 1) while on the 13C spectrum moderate deviations were noticed for the sugar ring atoms, with a maximum variation of 1.3 ppm for C4′. The evidence of boronate ester formation with the 2′-3′ cis-diol moiety and the absence of hydroxyl protons in anhydrous DMF-d7 suggested the presence of cyclic assemblies rather than linear oligomers.
To determine the size of these structures, diffusion ordered spectroscopy experiments (DOSY) were carried out and confirmed that the 2 sets of signals belonged to different species (Fig. 4). The diffusion coefficient measurements showed that both assemblies diffused at different but still relatively close rates (2.81 × 10−10 m2 s−1 and 3.26 × 10−10 m2 s−1 for the minor and major species respectively, Fig. S20†). The translational diffusion coefficient of the monomeric form was also determined from the sample in wet DMF-d7 to be 3.84 × 10−10 m2 s−1 (Fig. S21†). Considering that the addition of 4% of water in the sample did not significantly modify the medium viscosity, the diffusion coefficients were used to estimate the molecular weight of each structure using the Stokes–Einstein equation after correction of the f-factor according to Gierer and Wirtz (Table 3 and ESI†).51 From this calculation the molecular weights of the major and minor species were estimated to be 467 and 730 g mol−1 respectively. These values were in agreement with calculated molecular weights for dimeric (htt2) and trimeric (htt3) head-to-tail macrocycles (500 and 750 g mol−1 respectively).
![]() | ||
Fig. 4 DOSY spectrum of 3 in anhydrous DMF-d7 (for clarity purpose, parts of the spectra free of signal were removed). Red and blue dotted lines refer to minor and major species respectively. |
Compound | 3 | ||||||
---|---|---|---|---|---|---|---|
a Maximum standard deviation: 0.08 × 10−10 m2 s−1. | |||||||
D (×10−10 m2 s−1) | 3.84 | 3.19 | 2.96 | 3.29 | 2.86 | 3.26 | 2.81 |
To assess the efficiency of the DOSY experiment for the characterization of these dynamic assemblies, we synthesized linear (dinucleotide U2 and trinucleotide U3) and cyclic (c(T2) and c(T3)) reference compounds and measured their diffusion coefficients in anhydrous DMF-d7. The linear compounds were prepared using standard automatized phosphoramidite chemistry while the cyclic derivatives were obtained according to the procedure described by Smietana and Kool.52 The reference compounds were selected for practical reasons assuming that their diffusion coefficients would not differ significantly from c(U2) and c(U3) which would require more steps to achieve their synthesis. Indeed, an excellent correlation of the diffusion coefficients between the cyclic c(T2) and c(T3) reference compounds and the corresponding boronate assemblies was obtained (Table 3). Together, these results clearly indicate the head-to-tail formation of boronate-linked macrocyclic nucleotides analogues. Given the rapid opening of these macrocycles in the presence of water, their separation and isolation by HPLC techniques is illusive. We have already demonstrated in the past during DNA-templated ligation experiments that the generated internucleoside boronate linkages were disrupted upon denaturation i.e. in the absence of the template.41,46 The very fact that the boronate-linked macrocycle is formed and observed by NMR in the absence of a matrix is already remarkable.
We then investigated the impact of the boronate linkage on the intrinsic conformational preference of the sugar moiety. First structural insights about the macrocyclic dimer htt2 were obtained from molecular modeling. The geometry of the structure was optimized by means of quantum mechanics (QM) calculations using the Gaussian 03 software package with the MP2/6-31G(d) method (Fig. 5). Analysis of the sugar ring pucker for this model indicated a flattened C4′-exo conformation (puckering amplitude Φm = 26°). In spite of our efforts we did not succeed in obtaining crystals of htt2 or htt3. Therefore, to experimentally support the molecular modeling results, vicinal coupling constants were analyzed. In free nucleosides, the ribose adopts a conformation among the most representative forms, namely C2′-endo (south) and C3′-endo (north) usually subjected to a rapid interconversion equilibrium.
For the C2′-endo conformation, typical values of the vicinal coupling constants are found about 8–10 Hz for 3J1′-2′ and 1–2 Hz for 3J3′-4′ while these orders of magnitude are reversed for the C3′-endo conformation. The three ribose coupling constants for native uridine and its boronic analogue 3 are given in Table S4.† In water, these two compounds presented comparable values, characteristic of mononucleotides operating a dynamic north–south equilibrium. In wet DMF-d7, the higher 3J1′-2′ value indicated a slight preference of the south conformation. In anhydrous DMF on the contrary, markedly lower 3J1′-2′ constants were observed (3.2 and 2.5 Hz for htt2 and htt3 respectively, Table S4†) suggesting a significant prevalence of the north conformation. The population of the north conformers for 3 were estimated to be 53% in water, 47% in wet DMF-d7 and 68% and 75% in anhydrous DMF-d7 for htt2 and htt3 respectively (calculated from % north = 100 − 10 × 3J1′-2′).54
On the other hand, 3J2′-3′ coupling constant (∼5 Hz in free mononucleotides) is almost unaffected by the north–south equilibrium but strongly depends on the ring pucker.55,56 For htt2, both the rather high 3J2′-3′ value (7.8 Hz) and the lower value of 3J1′-2′ + 3J3′-4′ (8.3 Hz compared to 9.5 Hz observed in the classical north–south equilibrium), indicated a flattening of the ribose ring and a possible atypical conformation.
The coupling constants were then calculated from a generalized Karplus equation and empirical relationships between torsion angles and pseudorotation parameters (P and Φ) defined by Altona and coworkers for β-D-ribose,53 and compared to the experimental values measured for htt2 (Table 4). The results for ideal C3′-endo (P = 18°, Φm = 37°) and C2′-endo conformations (P = 162°, Φm = 37°) led to predictable markedly high RMS deviations (Table 4, entries 1 and 2). Assuming a single conformation, the phase angle and the ring pucker were then optimized to fit the experimental data (Table 4, entry 3) and P and Φm values corresponding to a flattened C4′-exo conformer, together with a reasonable RMSD (0.21 Hz) were obtained. To evaluate if the coupling constants could result from an equilibrium between flattened C2′-endo and C3′-endo conformations, the calculations were then repeated assuming a symmetrical two states model (e.g. ΦN = ΦS) and optimizing both the molar fraction of the C3′-endo conformer (XN) and the Φ values. In this case, calculated values were found to fit the experimental data with a reasonable deviation (0.14 Hz) when reducing the pucker amplitude to 15° and for a population of the north conformer of 70% (Table 4, entry 4). Finally, a last calculation was run assuming an equilibrium between the C4′-exo conformation and any other south conformer, optimizing PN, Φm and XN values (Table 4, entry 5). In that case, the best RMSD deviations was reached (0.03 Hz) for a C4′-exo ↔ C4′-endo type equilibrium, with a population of 69% of the C4′-exo form.
Entry | Model |
P
N![]() |
Φ
N![]() |
P
S![]() |
Φ
S![]() |
X
N![]() |
J
1′-2′![]() |
J
2′-3′![]() |
J
3′-4′![]() |
RMSDd (Hz) |
---|---|---|---|---|---|---|---|---|---|---|
a P and Φ are the amplitude puckering and the phase angle of pseudorotation respectively, as defined by Haasnoot et al.53 b X N represents the mole fraction of the north conformer. c Experimental values obtained for htt2: 3J1′-2′ = 3.2 Hz; 3J2′-3′ = 7.8 Hz; 3J3′-4′ = 5.1 Hz. d RMSD = root-mean-square deviation. | ||||||||||
1 | C3′-endo | 18 | 37 | 1.20 | 5.14 | 8.79 | 2.87 | |||
2 | C2′-endo | 162 | 37 | 7.85 | 5.20 | 1.24 | 3.80 | |||
3 | Single conformer | 52 | 9 | 3.13 | 8.15 | 5.12 | 0.21 | |||
4 | C2′-endo ↔ C3′-endo | 18 | 15 | 162 | 15 | 0.70 | 3.06 | 7.66 | 4.97 | 0.14 |
5 | C4′-exo ↔ C4′-endo | 54 | 20 | 234 | 20 | 0.69 | 3.17 | 7.85 | 5.11 | 0.03 |
Both the single conformation (Table 3, entry 3) and the north ↔ south equilibrium models (Table 4, entries 4 and 5) gave reliable fitting results, although in the former case the calculations led to an overestimation of the vicinal 3J2′-3′ constant (0.3 Hz).
To determinate whether the sugar maintains a certain flexibility or is conformationally locked, the coupling constants temperature-dependence was assessed. 1H NMR spectra were recorded at 273 K and 333 K and almost unchanged J values were measured (ΔJmax = 0.2 Hz). This confirmed that the sugar was frozen in a rigid C4′-exo conformation as highlighted by the molecular modeling study. Nucleosides bridged at the 2′ and 3′ position either by an alkyl chain57 or by a phosphate group have been already reported in the literature.58,59 In solution, 2′-3′-O-methoxymethyluridine was reported as a flexible system engaged in a C2′-endo ↔ C3′-endo equilibrium with a slightly flattened pucker (Φm = 34°)57 and 2′-3′-phosphate cyclic mononucleosides were considered as similar flexible systems.58 By contrast 2′-3′-isopropylidene nucleosides were described as rigid structures.60
Whenever the structures were studied in liquid or in solid state, the formation of a fused five-membered ring bridging the 2′ and 3′ oxygen atoms clearly constrains the ribose, thus severely reducing the pucker amplitude (Φm between 17 and 32°) compared to that commonly found in unmodified nucleosides and nucleotides (Φm about 39°).61 The macrocycle described here presented an additional constraint due to the 5′-2′/3′ head-to-tail cyclization that probably contribute to the structure rigidity.
Comparison of the dihedral exocyclic angles from QM geometry optimization and for the 2 models resulting from the fitting of the NMR data are given in Table S5.† Both NMR and QM gave comparable angles values, characteristic of a favored north C4′-exo conformation with a reduced pucker amplitude thus corroborating the formation of the cyclic homodimer htt2.
Several nucleoside derivatives adopting a favored C4′-exo conformations have already been reported in the literature. In all cases, these derivatives are constrained either by chemical modifications of the sugar,62–64 or because of an association of a nucleotide with a protein target.65
The 6′-boronic acid uridine analogue 3 described here is an example of a highly constrained structure, due to both the formation of a five-membered-ring at the 2′-3′ extremities of the sugar and the head-to-tail cyclization.
Footnote |
† Electronic supplementary information (ESI) available. See https://doi.org/10.1039/d2ob00232a |
This journal is © The Royal Society of Chemistry 2022 |