Iván
Rivilla
a,
Abel
de Cózar
ab,
Thomas
Schäfer
bc,
Frank J.
Hernandez
c,
Alexander M.
Bittner
cd,
Aitziber
Eleta-Lopez
d,
Ali
Aboudzadeh
c,
José I.
Santos
e,
José I.
Miranda
e and
Fernando P.
Cossío
*a
aDepartment of Organic Chemistry I, Centro de Innovación en Química Avanzada (ORFEO-CINQA), Universidad del País Vasco/Euskal Herriko Unibertsitatea (UPV/EHU), Donostia International Physics Center (DIPC), Po Manuel Lardizabal 3, E-20018 Donostia/San Sebastián, Spain. E-mail: fp.cossio@ehu.es
bIkerbasque, Basque Foundation for Science, Ma Díaz de Haro 3, E-48013 Bilbao, Spain
cNanoBioSeparations Group, POLYMAT University of the Basque Country (UPV/EHU), Avda. Tolosa 72, E-20018 Donostia/San Sebastián, Spain
dCIC NanoGUNE, Consolider. Tolosa Hiribidea, 76, E-200018 Donostia/San Sebastián, Spain
eSGIker NMR Facility, Universidad del País Vasco/Euskal Herriko Unibertsitatea (UPV/EHU), Avda. Tolosa 72, E-20018 Donostia/San Sebastián, Spain
First published on 22nd August 2017
A novel catalytic system based on covalently modified DNA is described. This catalyst promotes 1,3-dipolar reactions between azomethine ylides and maleimides. The catalytic system is based on the distortion of the double helix of DNA by means of the formation of Pt(II) adducts with guanine units. This distortion, similar to that generated in the interaction of DNA with platinum chemotherapeutic drugs, generates active sites that can accommodate N-metallated azomethine ylides. The proposed reaction mechanism, based on QM(DFT)/MM calculations, is compatible with thermally allowed concerted (but asynchronous) [π4s + π2s] mechanisms leading to the exclusive formation of racemic endo-cycloadducts.
Despite the relevance of (3 + 2) cycloadditions in the chemical synthesis of five-membered rings,15 the enzymatic version of this reaction has not been identified in living systems.16,17 Only a very recent example of a possible enzymatic 1,3-dipolar reaction has been reported to date.18 In addition, nonenzymatic 1,3-dipolar reactions have been postulated in the biosynthesis of several alkaloids19 and natural products such as furanocembranoids20,21 and santolin Y.22 Therefore, to the best of our knowledge a biomolecule-assisted bona fide (3 + 2) cycloaddition between azomethine ylides and alkenes to produce unnatural proline derivatives has not been reported to date. Within this context, we herein describe the first example of a DNA-assisted 1,3-dipolar reaction in water.
Scheme 1 Basic design of DNA-based catalysts for 1,3-dipolar reactions. The distortion of DNA by coordination of G bases with Pt(II) is shown. |
First, we tested the feasibility of this basic design using computational methods. As starting point we took the crystal structure reported by Takahara et al.27 for a cisplatin-double-stranded oligodeoxynucleotide complex (pdb code: 1aio), where G*G* denotes a 1,2-intrastrand cis-[(H3N)2Pt-d(GpG)] crosslink. To this structure 2,2′-bipyrimidin (bipym) was added and the azomethine ylide derived from methyl (E)-2-(benzylidene-amino)acetate 4a was coordinated to a copper(II) metallic centre. The resulting structure was stabilized by incorporating 509 water molecules and 22 sodium cations (Fig. 1a). The whole ensemble was optimized using a hybrid QM/MM ONIOM28–30 scheme (ESI†). The full structure thus optimized was found to keep the folded geometry of the distorted double helix, where the Cu(II)–bipym–Pt(II)–G*G* system generated a cavity to which the azomethine ylide was coordinated.
Fig. 1 (a) Fully optimized (B3LYP/LanL2DZ:UFF level of theory) structure of a DNA double strand containing the d(CCTCTG*G*TCTCC)-d(GGAGACCAGAGG) sequence, in which the complex 1a is bound to the G*pG* unit and to the N–Cu(II) azomethine ylide derived from imine 4a (R2 = Ph, Scheme 1), surrounded by 509 water molecules and 22 sodium cations. The QM part of the optimization is within the circle highlighted in yellow. (b) The same optimized structure but showing the solvent accessible surface of the DNA fragment. The water molecules have been removed for clarity. Pt(II) and Cu(II) atoms are represented in light blue and green, respectively. |
As it can be seen by inspection of Fig. 1b, the resulting ensemble closely resembles the active site of a metalloenzyme. As a proof-of-concept experiment, we next examined the ability of (bipym)PtCl2 complex 1a to bind two equivalents of guanosine. Since 1a and its derivatives posed solubility problems, we monitored the different species in the solid state using Cross Polarization-Magic Angle Spinning (CP-MAS) spectroscopy.31 The NMR spectrum of bipym showed only one 15N-NMR signal, as expected from its D2h symmetry (Fig. 2). In contrast, the two nitrogen atoms of 1a coordinated to Pt(II) could be readily distinguished. Then, we used 5′-GMP as a suitable equivalent of G units in DNA and analysed the 1H–15N CP-MAS spectra at different mixing times to assign the five nitrogen atoms of the G-unit. Combination of 1a and 5′-GMP resulted in the formation of an adduct whose 15N chemical shift associated with the N7 atom of the purine system was considerably deshielded with respect to the signal recorded for 5′-GMP (Fig. 2).
We concluded that Pt(II) complex interacts with guanine units yielding a square planar complex, in which two equivalent G units bind the Pt(II) centre by means of the respective N7 atoms of the purine bicycle. We next studied the interaction between 1a and oligodeoxynucleotides in order to assess the binding abilities of different DNA sequences.
To this end, we used a quartz-crystal microbalance with dissipation monitoring (QCM-D) device.32 These QCM-D experiments measured changes in frequency, which correlate with mass adsorption on the sensor surface, and the dissipation of energy of the adsorbed film, which in turn correlates with its viscoelasticity.33,34 The dissipation vs. frequency points were collected for various times and ordered graphically from left to right and from bottom to top (Fig. 3). In a reference experiment denoted as “exp. 1” in Fig. 3a, the Pt(II) complex 1a was injected on bare gold (no DNA immobilized) in PBS buffer. Under these conditions, dissipation increased linearly with frequency, thus resulting in an incremental dissipation entirely correlated with the mass adsorbed on the surface. To carry out the experiments in the presence of DNA, we selected two different model oligodeoxynucleotides on the basis of the well-known binding ability of GpG pairs with platinum drugs.23,24 We tested in duplicate the oligomer containing 5′-thiol-AAAAATTAAATTAAA-3′ binding sequence (Fig. 3a, experiments 2 and 3). An overlay with the reference experiment 1 revealed that this response was indeed non-specific as both experiments showed practically identical profiles with respect to the reference. We interpreted these results as negative experiments in which there was no significant interaction between Pt(II) complex 1a and G-free deoxyoligonucleotide (Fig. 3b). We next examined the behavior of the oligomer containing 5′-thiol-AAAAAGGAAAGGAAA-3′ binding sequence (Fig. 3c, experiments 4 and 5). In this case, upon addition of 1a both plots showed an identical increase in dissipation while simultaneously no frequency change was observed, in sharp contrast with the previous blank and negative experiments 1–3. This indicated that first, there was no measurable non-specific binding occurring as it would be detected as a change in frequency (see the dotted lines in Fig. 3c); second, injection of 1a induced a strong change in dissipation, which can be only attributed to conformational changes taking place in the GpG pairs-containing immobilized and hybridized deoxyoligo-nucleotide. From these results, we concluded that GpG-containing DNA interacts specifically with Pt(II) complex 1a resulting in an increasingly more dissipative film, which in turn correlates with additive conformational changes of DNA upon addition of 1a (Fig. 3d).
Once we verified that complex 1a binds selectively G-containing oligonucleotides, we tested its binding ability to DNA strands. As a control experiment, we mixed 1a with 5′-GMP and the resulting complex was analyzed by MALDI-TOF mass spectrometry.35,36 We observed an ensemble of high mass-to-charge (m/z) signals distributed around a value of m/z ≈ 1042 a.u., as shown in Fig. 4a. These signals correlate satisfactorily with complex 7 (Fig. 4b), closely related to 6 (Fig. 2), in which two ketone units have been generated via dehydration-tautomerization of one hydroxy group of each ribose unit. After this control experiment, the same protocol was followed, but instead of 5′-GMP we used DNA sodium salt from salmon sperm (salmon sperm DNA in Fig. 4) with a % G-C content of 41.2% and a molecular mass of 1.3 × 106 Da (ca. 2000 bp). In this case, the same profile was observed in the corresponding MALDI-TOF mass spectrum (Fig. 4c), with a high m/z ensemble centered at ca. 1042 a.u. This response was interpreted in terms of ion 7 or, for two consecutive GpG units in the starting salmon sperm DNA, as the keto-enol phosphoric ester depicted as ion 8 in Fig. 4d. On the basis of these experiments, we concluded that the 1a–G2 complexes observed in monomeric and oligomeric G-containing species can be extended to double strain DNA.
We also studied the effect of Pt(II) complex 1a on the structure of DNA strands. Thus, samples of DNA (ca. 48 kb), from λ phage-infected E. coli were analyzed by atomic force microscopy (AFM).37 The corresponding DNA strands were unambiguously identified on oxidized silicon38 by the corresponding AFM images (Fig. 5a and b). When 1a was added, the AFM scans showed large morphological changes consisting of local kinks and crosslinks (Fig. 5c and d). On the basis of the previously presented experiments, these changes were attributed to the formation of intra- and interstrand cis-{(bipym)Pt(d[GpG + GpXpG + GpA])} adducts.39
These advanced analytical studies permitted us to establish the main features of the Pt(II)-mediated binding between molecule 1a and DNA. We concluded that the thus generated hybrid system could catalyse 1,3-dipolar reactions provided that (i) the heterobimetallic Cu(II)–Pt(II) centre is able to generate in situ the required N-metallated azomethine ylide derived from the corresponding imine 4; and (ii) the active site can accommodate the dipolarophile 3.
In order to test the catalytic ability in (3 + 2) cycloadditions of the DNA–heterobimetallic complex, we prepared catalyst 2 (Scheme 1) by using salmon sperm DNA in a buffered solution of (N-morpholino)propane-sulfonic acid (MOPS), to which 1a was added, followed by copper(II) triflate, triethylamine and the corresponding maleimide 3 and imine 4 (Scheme 2). After six days of reaction at room temperature, nine distinct endo-(3 + 2)-cycloadducts (5aa–5bf, Table 1) could be generated. This endo stereochemistry was secured on the basis of the spectroscopic data and, in two cases, by X-ray diffraction analysis (ESI†).
Scheme 2 Synthesis of (3 + 2) endo-cycloadducts 5aa–5bf catalysed by the heterobimetallic complex formed by salmon sperm DNA and 1a in the presence of copper(II) triflate. |
Entry | R1 (3) | R2 (4) | 5 | Yielda (mmol × 10−3) [%] |
---|---|---|---|---|
a Yields of isolated pure products correspond to averaged values obtained after three experiments. The limiting reactant was the corresponding maleimide 3 (0.05 mmol). b Structure confirmed by X-ray diffraction analysis (ESI).40 | ||||
1 | Me (3a) | Ph (4a) | 5aa | 10.5 ± 0.4 [21] |
2 | Ph (3b) | Ph (4a) | 5ba | 18.0 ± 0.5 [36] |
3 | 4-MeO–C6H4 (3c) | Ph (4a) | 5ca | 10.0 ± 0.8 [20] |
4 | Ph (3b) | 4-Me–C6H4 (4b) | 5bb | 17.9 ± 0.4 [23] |
5 | Ph (3b) | 4-MeO–C6H4 (4c) | 5bc | 17.5 ± 0.3 [35] |
6 | Ph (3b) | 4-F-C6H4 (4d) | 5bd | 12.5 ± 0.2 [25] |
7 | Ph (3b) | Cyclohexyl (4e) | 5be | 13.5 ± 0.4 [27] |
8b | Ph (3b) | 2-Thyenyl (4f) | 5bf | 16.0 ± 0.9 [32] |
In order to assess the relevance of each component of the catalytic system we tested 17 possibilities resulting from the combination of all the reagents except at least one (see Table S1 of the ESI†). In all these control experiments no 1,3-dipolar reaction was observed. In particular, when different combinations of triethylamine and the Cu(II) and Pt(II) salts were tested in the absence of salmon sperm DNA, the reaction did not proceed. Similarly, different combinations in the presence of salmon sperm DNA but in the absence of base, 2,2′-bipyrimidine and/or one of the metals did not produce any (3 + 2) cycloadduct. It is interesting to note that the combination of DNA, 2,2′-bipyrimidine and Cu(OTf)2 in the absence of Pt(II) was also unproductive (see Table S1 of the ESI,† entry 13), thus indicating that the powerful Roelfes–Feringa method consisting of a metallated intercalating heterocycle cannot catalyse this challenging 1,3-dipolar cycloaddition. In addition, these experiments demonstrate that there is no background reaction (see in particular entry 1 of Table S1 of the ESI†). In summary, these control experiments demonstrated that combination of 2,2′-bipyrimidine, Pt(II), Cu(II) and DNA, most likely by bonding to consecutive GG units, is required to achieve moderate yields of (3 + 2) racemic endo-cycloadducts 5 (see Table S1 of the ESI,† entry 10).
We interpreted our results as follows: Previously formed DNA-1a adducts bound Cu(OTf)2 and the resulting heterobimetallic complex 2a (Fig. 5a) coordinated imine 4 to form intermediate species INT1, from which the corresponding N-metallated azomethine ylide INT2 was formed via triethylamine-assisted deprotonation. This 1,3-dipole interacted with dipolarophile 3 to form the corresponding (3 + 2) cycloadduct and regenerating INT1via interaction with another equivalent of imine, thus completing the catalytic cycle. In order to understand the origins of the endo selectivity we optimized the possible endo- and exo-transition structures under the same computational framework used to optimize the structure of INT2 (Fig. 5b and c). Both saddle points were found to be quite asynchronous but still associated with a concerted [π4s + π2s] symmetry allowed mechanism, as indicated by the bond distances corresponding to the formation of the two new σ-bonds. We also found that endo-TS was stabilized by a strong electrostatic interaction between the nitrogen atom of the imide moiety and the metallic centre (Fig. 5b). As a result, exo-TS was calculated to be ca. 13 kcal mol−1 higher in energy than its endo congener, thus predicting the preferential formation of cycloadduct endo-5aa under kinetic control, in nice agreement with the experimental data.
We performed similar calculations in the absence of DNA. In these studies, the features of the computational model (including surrounding water molecules, see the ESI†) remained identical. The results are gathered in Fig. 6d and e. We observed that in the presence of Cu(II) and bipyrimidine, the chief features of the transition structure leading to endo-5aa are similar to those found for the 3a + 4a → 5aa reaction in the presence of DNA. However, the activation energy was found to be ca. 12 kcal mol−1 higher, which corresponds to a k(DNA–Cu–Pt)/k(Cu) ratio of ca. 4.8 × 108 (Fig. 6d). When the square planar diaqua–Pt(II) moiety was incorporated to the reaction coordinate (Fig. 6e), the shape of the corresponding saddle point did not change significantly and the activation energy was slightly lower than in the previous case, with a calculated k(DNA–Cu–Pt)/k(Cu–Pt) ratio of ca. 2.7 × 107. It is interesting to note that the DNA-free simulations also predict the preferential formation of the endo-cycloadduct (see Fig. S3 and Table S3 of the ESI†). These results are in agreement with the absence of reactivity observed when DNA was not present. We interpret the lower activation energies in the presence of DNA in terms of the destabilization of the substrate by restriction of the conformational freedom and by the electrostatic repulsion between the anionic polyphosphate environment and the anionic part of the starting 1,3-dipole, which is alleviated along the reaction coordinate leading to the non-zwitterionic cycloadduct.
Fig. 6 (a) Proposed catalytic cycle for (3 + 2) cycloadditions between azomethine ylides INT1 and maleimides 3 catalyzed by in situ formed DNA–bipym–Pt(II)–Cu(II) complexes. (b, c) Fully optimized geometries and relative energies of endo- and exo- transition structures (TS) associated with the 3a + 4a → 5aa reaction. Numbers in parentheses are the relative energies, in kcal mol−1. Bond distances are given in Å. Guanine units bound to Pt(II) are highlighted in blue. Calculations were performed at the B3LYP/LanL2DZ:UFF level of theory as indicated in Fig. 1b. Water molecules and sodium cations were included in the calculations (MM level, UFF force field) but are not shown for clarity. (d, e) Same transition structures associated with the formation of cycloadduct endo-5aa in the absence of DNA and Pt(II) (d), and in the absence of DNA (e). In the latter case, two binding water ligands have been included. Also in (d) and (e), water molecules were included in the calculations (MM level, UFF force field, not shown for clarity). |
Footnote |
† Electronic supplementary information (ESI) available: Experimental details, characterization of cycloadducts and intermediates, computational data. CCDC 1411372 and 1411373. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c7sc02311a |
This journal is © The Royal Society of Chemistry 2017 |