Harikrishnan M.
Kurup
ab,
Maksim V.
Kvach
a,
Stefan
Harjes
a,
Geoffrey B.
Jameson
ab,
Elena
Harjes
*ab and
Vyacheslav V.
Filichev
*ab
aSchool of Natural Sciences, Massey University, Private Bag 11 222, Palmerston North 4442, New Zealand. E-mail: e.harjes@massey.ac.nz; v.filichev@massey.ac.nz
bMaurice Wilkins Centre for Molecular Biodiscovery, Auckland 1142, New Zealand
First published on 7th June 2023
The APOBEC3 (APOBEC3A-H) enzyme family as a part of the human innate immune system deaminates cytosine to uracil in single-stranded DNA (ssDNA) and thereby prevents the spread of pathogenic genetic information. However, APOBEC3-induced mutagenesis promotes viral and cancer evolution, thus enabling the progression of diseases and development of drug resistance. Therefore, APOBEC3 inhibition offers a possibility to complement existing antiviral and anticancer therapies and prevent the emergence of drug resistance, thus making such therapies effective for longer periods of time. Here, we synthesised nucleosides containing seven-membered nucleobases based on azepinone and compared their inhibitory potential against human cytidine deaminase (hCDA) and APOBEC3A with previously described 2′-deoxyzebularine (dZ) and 5-fluoro-2′-deoxyzebularine (FdZ). The nanomolar inhibitor of wild-type APOBEC3A was obtained by the incorporation of 1,3,4,7-tetrahydro-2H-1,3-diazepin-2-one in the TTC loop of a DNA hairpin instead of the target 2′-deoxycytidine providing a Ki of 290 ± 40 nM, which is only slightly weaker than the Ki of the FdZ-containing inhibitor (117 ± 15 nM). A less potent but notably different inhibition of human cytidine deaminase (CDA) and engineered C-terminal domain of APOBEC3B was observed for 2′-deoxyribosides of the S and R isomers of hexahydro-5-hydroxy-azepin-2-one: the S-isomer was more active than the R-isomer. The S-isomer shows resemblance in the position of the OH-group observed recently for the hydrated dZ and FdZ in the crystal structures with APOBEC3G and APOBEC3A, respectively. This shows that 7-membered ring analogues of pyrimidine nucleosides can serve as a platform for further development of modified ssDNAs as powerful A3 inhibitors.
![]() | ||
Fig. 1 (A) Deamination of dC in ssDNA by A3 enzymes. (B) Previously described powerful CDA inhibitors. |
Despite low sequence identity, cytidine deaminase (CDA) and A3 share a similar overall structural topology and close structural homology for the Zn2+-containing active site. Consequently, these enzymes are believed to have a similar mechanism of cytosine deamination. The major difference between them is that CDA accepts only individual nucleosides as substrates, whereas A3s act mainly on ssDNA having at least four nucleotides, one of which is 2'-deoxycytidine.20
CDA is a key enzyme of primary metabolism that catalyses the deamination of cytidine or 2′-deoxycytidine to uridine or 2′-deoxyuridine, respectively.21,22 The spontaneous deamination of cytidine is a very slow process; the accelerated rate of deamination by the enzyme is through the formation of a hydrated intermediate species.23 CDA forms a part of the pyrimidine salvage pathway and full inhibition of CDA leads to accumulation of toxic pyrimidine catabolism intermediates.24,25 However, several cytosine-based chemotherapeutics such as cytarabine,26 gemcitabine27 and decitabine28 have been shown to be deaminated by CDA in the liver and spleen, which reduces their potency. Local and temporary inhibition of CDA in these organs is therefore a useful strategy that can mitigate detrimental degradation of chemotherapeutics. In July 2020, a combination of the known CDA inhibitor cedazuridine, (4R)-2′-deoxy-2′,2′-difluoro-3,4,5,6-tetrahydrouridine,29 and decitabine as an oral pill (C-DEC or ASTX727) developed by Astex Pharmaceuticals was approved by the US Food and Drug Administration (FDA) for the treatment of patients with intermediate/high-risk myelodysplastic syndromes (MDS) and chronic myelomonocytic leukemia (CMML).30
The most powerful inhibitors of CDA based on nucleosides described so far are 1,3-diazepin-2-one riboside (rDiAzep, Fig. 1B),31–36 phosphapyrimidine riboside (unstable in water),37 tetrahydrouridine (THU)38 and zebularine (Z).39,40 We demonstrated that the use of 2′-deoxy analogues of zebularine (dZ) and 5-fluorozebularine (FdZ) instead of dC in the preferred ssDNA substrate motif resulted in competitive inhibitors of A3.17–19 Recently, we have shown that cross-linking of two distant nucleotides at +1 and −2 positions resulted in a more powerful dZ-containing A3A inhibitor than the linear ssDNA.41 DNA hairpins with short loops have also been shown to provide selective inhibitors of A3A when the target dC in the loop was changed to dZ-derivatives.42–44 We also demonstrated that dZ and THU as free nucleosides did not inhibit A3 enzymes, which indicates that ssDNA delivers dZ into the active site of A3.17
In the past, the diazepinone riboside was described as a more powerful inhibitor of CDA than dZ (Ki = 25 nM for 1,3-diazepin-2-one31–35 and 2.9 μM for dZ).45 Hence, we decided to expand a series of 2′-deoxy nucleoside analogues of CDA inhibitors31,34,46 used in ssDNA for A3 inhibition and investigate 1,3-diazepin-2-one riboside (dDiAzep) and its hydrated derivatives (dDiAzep-OH)31–35 in this context.
Out of all the reported potent CDA inhibitors, 1,3-diazepin-2-one riboside31–35 is notably different structurally as it has no functionality at C5 to coordinate with Zn2+ or to undergo hydration in the active site of CDA. In the crystal structure of 1,3-diazepin-2-one riboside with human CDA (hCDA, PDB ID: 1MQ0),35 the protein self-associates to form a dimer of dimers (tetramer) (Fig. S2 in the ESI†) which is known to be active.19 Each subunit is bound to one molecule of diazepinone riboside but each active site is made up of three of the four hCDA monomers and all of them contribute to the recognition of the inhibitor. It is interesting that diazepinone is not coordinated to Zn2+. Instead, a water molecule is bound to Zn2+. The nucleobase N3–H is engaged in interactions with Glu67 of hCDA and the oxygen of the carbonyl C2O is hydrogen bonded to the main chain NH of Ala66 of the first subunit. The seven-membered nucleobase is puckered and the C5–C6 double bond of the molecule is involved in a canonical π–π interaction with Phe137 residue of the third subunit of hCDA. Reduction of the double bond to the corresponding saturated nucleoside led to significant loss in binding affinity to CDA which highlights the importance of the double bond for the inhibition of hCDA.31 Interestingly, aromatic amino acids Tyr130 and Tyr313 are present at the same position in the active site of monomeric A3A and A3BCTD, respectively, instead of Phe137 in hCDA (Fig. 2). We hypothesised that a π–π interaction between the double bonds of diazepinone and Tyr130(A3A)/Tyr313(A3Bctd) will be preserved for A3A and A3B and potentially result in potent inhibition of these enzymes by diazepinone. These observations make this compound interesting to be synthesised as a 2′-deoxy analogue (dDiAzep, Fig. 1B) for testing as an A3 inhibitor when incorporated into ssDNA.
![]() | ||
Fig. 2 Active site residues (grey) and diazepinone riboside (green) of hCDA (1MQ0) overlaid with the active site of A3A (5KEG) (cyan). Tyr 130 (cyan) of A3A is on the same place as Phe 137 (orange) of hCDA interacting with diazepinone; zinc2+ atoms are shown as spheres. Structures were overlaid in pymol. The CDA tetramer is shown in ESI Fig. S2.† |
To mimic the respective tetrahedral substrate intermediate in the context of 7-membered ring nucleobase, diazepinone derivatives carrying an OH-group on carbon C5 have been prepared as R and S isomers in the past (rDiAzep-OH, Fig. 1B).31 One of the isomers was identified as a slow-binding inhibitor of CDA, as potent as THU. Unfortunately, the authors were unable to determine which isomer was more potent. The lack of proper stereochemical information prompted us to attempt the synthesis, isolation, and characterisation of the stereochemistry of both dDiAzep-OH isomers before their incorporation into ssDNA as A3 inhibitors.
As reported in the literature, acylation of the amide 3 was required to lock the allyl groups in a cis orientation for the successful RCM.47 The amide 3 was protected with benzoyl chloride to provide compound 4. The subsequent RCM with the GreenCat™ catalyst provided the required seven-membered ring 5 in 41% isolated yield. Resuspending and washing the precipitate 5 in methanol gave a >99% pure β anomer. The toluoyl groups on compound 5 were selectively deprotected with a solution of aqueous ammonia in methanol followed by selective protection of the 5′-hydroxyl group using standard DMT protection conditions to provide compound 6 in 82% yield. Phosphoramidite 7 was prepared using standard phosphitylation conditions with an isolated yield of 82%.
The fully deprotected nucleoside, dDiAzep (8), was obtained by the deprotection of compound 5 in concentrated aqueous ammonia and isolated by preparative thin-layer chromatography (TLC) in 58% yield.
The above failed reaction on substrate 5 prompted us to investigate the possibility of hydroboration/oxidation on a protected nucleobase 12. We started the synthesis with the previously obtained bis-allyl urea 2. Benzoyl protection of the compound 2 followed by tert-butyloxy carbonyl (boc) protection allowed the facile RCM of compound 11 producing the protected nucleobase 12 in 50% overall yield over three steps. The attempted hydroboration/oxidation protocol resulted only in benzoyl deprotection. Next, we tried the same hydroboration/oxidation protocol on a free nucleobase 14 which was reported previously for the synthesis of compound 15.32 Unfortunately, the reported procedure did not work in our study and the starting material was isolated (Scheme 2). In the past this protocol was not chosen for the synthesis of the target rDiAzep-OH. Instead, the formation of a diastereomeric mixture of ribosides of R- and S-5-hydroxy-diazepin-2-ones (Fig. 3A) relied on the use of toxic mercury salts (HgO/HgBr2) and perhydro-1,3-diazepine-2,5-dione that has to be synthesised.31
![]() | ||
Fig. 3 (A) Previously described ribosides of R- and S-5-hydroxy-1,3-diazepin-2-ones. (B) Proposed structure of 2′-deoxyribosides of 5-hydroxy-azepin-2-ones. |
We envisioned that a potent CDA and A3 inhibitor can be based on the nucleoside carrying 1,3-diazepine-2-one in which the N3 atom is replaced with CH2. One should note that the IUPAC nomenclature assignment of R and S isomers is opposite for 1,3-diazepin-2-one and azepin-2-one (Fig. 3).
The synthesis started from commercially available 1,4-cyclohexanedione cyclic ethylene monoketal (16) that was converted to oxime and then, using a Beckmann rearrangement, provided the required protected nucleobase as a 7-membered cyclic amide (17).50 This compound was coupled with Hoffer's chlorosugar51 in the presence of K-tBuO in 1,4-dioxane at room temperature to provide the required nucleoside 18 with a β:
α ratio of 99
:
1 as determined by NOESY NMR. However, deprotection of the ketal in compound 18 under various acidic conditions failed. The difficulties of deprotecting the cyclic ketal in 7-membered ring compounds have been reported in the past.32 To circumvent these difficulties, we decided to change the synthetic route (Scheme 3) and rely on the benzyl-protecting group on the alcohol at C5 of the hexahydro azepin-2-one (21).
A pathway reported in the current paper includes the synthesis of the protected azepan nucleobase 21 from the cyclic ketal 16 using the Beckmann rearrangement on oxime 20 (Scheme 3).50 The cyclic amide 21 and Hoffer's chlorosugar were coupled in the presence of K-tBuO in 1,4-dioxane at room temperature to provide a pair of diastereomers 22R and 22S in moderate yields (47 and 24%, respectively) and a β:
α ratio of 99
:
1. These diastereomers were purified using silica gel column chromatography and crystallised. The stereochemistry of each compound was assigned based on X-ray crystal data (Fig. 4). Here onwards, these diastereomers were subjected to the same sequence of reactions separately. Benzyl deprotection of 22 was performed by catalytic hydrogenation using 10% Pd/C at room temperature to give 23R and 23S with 93 and 88%, respectively. One part of the benzyl-deprotected compounds was treated with 30% aq. NH3 in MeOH at r.t. for 3 days to provide free nucleosides 24R and 24S with 72% and 79% yield, respectively, for evaluation as inhibitors of hCDA. The rest of each compound was protected with tert-butyl diphenylsilyl (TBDPS). After the removal of toluoyl groups, 5′-O-DMT-protected 3′-O-phosphoramidites 28R and 28S were prepared employing standard protocols.
![]() | ||
Fig. 4 X-ray structures of 22S (A) and 22R (B). Arrows indicate the stereocenters of these compounds where they differ. The structures have been deposited with the CCDC with deposition numbers 2207174 (22S) and 2207173 (22R).† |
We incorporated the modified nucleosides at the location of dC in the preferred A3 substrate motifs. Detailed procedures for oligodeoxynucleotide synthesis are provided in the ESI.† All oligodeoxynucleotides were purified by reverse-phase HPLC with detection at 260 nm. ESI-MS (Table 1) confirmed their compositions.
Name | DNA sequence, 5′ → 3′ | Retention timea (min) | ESI-MS [Da], found/calculated |
---|---|---|---|
a Reverse-phase HPLC on a 250/4.6 mm, 5 μm, 300 Å C18 column (Thermo Fisher Scientific) using a gradient of CH3CN (0 → 20% for 20 min, 1.3 mL min−1) in 0.1 M TEAA buffer (pH 7.0) with detection at 260 nm. | |||
7-mer dZ-linear | TTTT dZ AT | 14.2 | 2045.37/2045.37 |
9-mer dDiAzep | ATTT dDiAzep ATTT | 15.9 | 2679.48/2678.50 |
7-mer R-dAzep-OH | TTTT R-dAzep-OH AT | 14.9 | 2078.39/2078.41 |
7-mer S-dAzep-OH | TTTT S-dAzep-OH AT | 14.8 | 2078.38/2078.41 |
FdZ-hairpin | (GC)2TT FdZ (GC)2 | 14.3 | 3311.50/3310.60 |
dDiAzep-hairpin | T(GC)2TT dDiAzep (GC)2T | 14.3 | 3916.65/3916.69 |
S-dAzep-OH-hairpin | T(GC)2TT S-dAzep-OH (GC)2T | 14.0 | 3934.65/3933.71 |
By analysing the initial rate of deamination, we observed that the S-isomer of 24 was more potent than the R-isomer at all inhibitor concentrations with and without preincubation with hCDA. However, both analogues were not as powerful hCDA inhibitors as the previously described dZ. Therefore, no detailed analysis was performed. Instead, we focused our attention on the 1,3-diazepin-2-one 2′-deoxyriboside 8.
The 1,3-diazepin-2-one riboside was previously reported as a powerful inhibitor of CDA. In order to accurately compare the 2′-deoxy form of 1,3-diazepin-2-one riboside and dZ we monitored the deamination of dC at 286 nm and analysed the kinetic profiles at various inhibitor concentrations using a global regression analysis of the kinetic data using Lambert's W function.52 This method provides better estimates for Km and Vmax than non-linear regression analysis of the initial rate or any of the known linearised transformations of the Michaelis–Menten equation, such as Lineweaver–Burk, Hanes–Woolf and Eadie–Hofstee transformations.52 Then the Michaelis–Menten constant (KM) for the substrate and the inhibition constant (Ki) for each inhibitor were calculated (Table 2) assuming a competitive nature of inhibitors.
Indeed, dDiAzep nucleoside is an at least 20 times more powerful inhibitor of hCDA than dZ if Km/Ki values are compared. This warrants evaluation of dDiAzep incorporated in ssDNA as an inhibitor of A3 enzymes.
The results revealed that to achieve a similar level of inhibition induced by 7-mer-dZ we had to increase the concentration of oligodeoxynucleotides containing azepinone nucleotides from 8 to 36 μM. Despite the fact that dDiAzep was a more powerful inhibitor of hCDA in comparison with dZ, it did not exhibit a similar effect on engineered A3BCTD. Interestingly, the inhibition by R-dAzep-OH and S-dAzep-OH containing oligodeoxynucleotides was poor compared to the standard inhibitor (7-mer-dZ), but the difference in the inhibitory potential observed between R and S isomers was similar to that seen for hCDA (Fig. 5). This suggests the orientation of –OH on the nucleobase which mimics the intermediate formed after the attack of H2O at C4 of cytosine by A3BCTD-QM-ΔL3-AL1swap.
Inhibitor | K i (nM) |
---|---|
Conditions: 50 mM Na+/K+ phosphate buffer, pH 7.4 supplemented with 100 mM NaCl, 1 mM TCEP, 100 μM DSS and 10% D2O at 20 °C; enzyme concentration: 140 nM; substrate concentration (dC-hairpin): 500 μM. Km for the dC-hairpin against A3A is 21 ± 10 μM (ref. 43) (33 ± 7 μM in experiments with the S-dazep-OH-hairpin). Full structural characterisation of the FdZ-hairpin in the complex with A3A and the activity of the FdZ-hairpin in cells will be published elsewhere.43 | |
FdZ-linear | 2400 ± 940 |
FdZ-hairpin43 | 117 ± 15 |
dDiAzep-hairpin | 290 ± 40 |
S-dAzep-OH-hairpin | 6000 ± 1000 |
The Ki values show that FdZ and dDiAzep-containing hairpins are superior inhibitors of the wild-type A3A in comparison with the best linear oligodeoxynucleotide (FdZ-linear). It is also evident from these data that FdZ and dDiAzep are more potent inhibitors of A3 when used instead of dC in the structure of the hairpin than in the linear 9-mer sequence (Fig. 6 and Table 3). A hairpin carrying S-dAzep-OH instead of dC has a Ki of 6 μM which is worse than the Ki for FdZ-linear. These data support our assertion that nucleoside-like inhibitors of CDA can be converted into inhibitors of A3 when used instead of dC in a DNA hairpin.
Recently, we have shown that the incorporation of cytosine-like 2′-deoxyzebularine45,57 and 5-fluoro-2′-deoxyzebularine (dZ and FdZ on Fig. 1B) into short, linear ssDNA led to the first selective A3 inhibitors17–19 with low micromolar inhibition constants (Ki). We have proposed that the inhibitory potential of ssDNAs can be further improved through the incorporation of more potent inhibitors of CDA into ssDNA.19 Here we have considered several diazepin derivatives known to inhibit CDA and incorporated them into ssDNA as possible inhibitors of A3 enzymes. It was pleasing to note that the RCM methodology allowed us to obtain and evaluate diazepinone 2′-deoxyriboside (dDiAzep, Fig. 1B) as an inhibitor of hCDA, engineered A3BCTD and A3A. In comparison with dZ, dDiAzep with a Ki of 430 ± 80 nM was a 20 times more potent inhibitor of recombinantly expressed hCDA but it did not reach the Ki of 25 nM against human liver CDA reported forty years ago for the diazepinone riboside.36 It is plausible that different sugar puckering of ribose and 2′-deoxyribose had an effect here, as well as a different source of hCDA enzyme. The use of dDiAzep instead of dC in the loop of a DNA hairpin led to an A3A inhibitor with nM potency, exhibiting at least 70 times tighter binding to A3A than a dC-hairpin if their Ki and Km values are compared. However, the inhibitory effect of this nucleoside against hCDA was substantially higher (Km/Ki = 512, Table 2). This may be due to the differences in the π–π interactions of dDiAzep with Phe and Tyr because the binding energy of π–π interactions depends on interacting pairs.58
As a control for the evaluation of dDiAzep in the structure of a DNA hairpin, we synthesised and tested an FdZ-hairpin that showed powerful inhibition of wild-type A3A with a Ki of 117 nM.43 Full characterisation of this compound and its crystal structure with wild-type A3A will be published elsewhere.43
Unfortunately, due to chemical difficulties we had to abandon the synthesis of the 2′-deoxy form of dDiAzep-OH (Fig. 1B). Instead, we synthesised the R and S isomers of 1,3-diazepine-2-one 2′-deoxyriboside in which the N3 atom is replaced with CH2. The replacement of the nitrogen atom by CH2 might be a reason for a considerable decrease in the inhibitory potential of nucleosides against hCDA and engineered A3BCTD. However, the S-isomer was more potent than the R-isomer against both hCDA and engineered A3BCTD. Very recently, using X-ray crystallography it was revealed that FdZ and dZ form tetrahedral intermediates (4-(R)-hydroxy-3,4-dihydro-2′-deoxy-5-fluorozebularine and 4-(R)-hydroxy-3,4-dihydro-2′-deoxyzebularine) in the enzyme's active site with wild-type A3A42,43 and catalytically active C-terminal domain of A3G.59 The position of OH in hydrated FdZ and dZ is the same in the complex with CDA and A3 with carbon C4 having R-stereochemistry. Considering that the nomenclature changes when we replace N3 in 1,3-diazepin-2-one to CH2 in azepin-2-one (Fig. 3), our observation that the S-isomer is more active than the R-isomer for 5-hydroxy-azepin-2-ones is consistent with the same conformation seen for the S-isomer in the model (Fig. 7A) and for the hydrated FdZ in the crystal structure (Fig. 7B).43 This observation also implies that the 5-hydroxy group of S-dAzep-OH coordinates to Zn2+ and possibly displaces a water molecule that is usually coordinated to Zn2+ in A3's active site (Fig. 7).
![]() | ||
Fig. 7 (A) Chemical structure of S-dAzep-OH in ssDNA; (B) hydration of FdZ and dZ by A3 results in 4-(R)-hydroxy-3,4-dihydro-2′-deoxyzebularines. |
In the past, high-throughput screening efforts led to the identification of covalent small molecule A3G inhibitors with low micromolar potencies.60,61 During our study two articles have been published on drug discovery efforts to obtain small molecule inhibitors of A3.62,63 None approaches the nM potency exhibited by our FdZ43 and dDiAzep-containing hairpins. This highlights the interest and need for powerful A3 inhibitors. It is gratifying that the 7-membered ring nucleobase upon incorporation into DNA inhibits A3A, which provides another useful scaffold for further development of A3 inhibitors.
Footnote |
† Electronic supplementary information (ESI) available: Experimental details of the synthesis of nucleosides and modified oligodeoxynucleotides and enzymatic assays; and the sequence alignment of proteins used in this study, 1H, 13C, 31P NMR, IR and HRMS (ESI) spectra of new compounds synthesised, RP-HPLC profiles and HRMS (ESI) spectra of oligodeoxynucleotides. CCDC 2207173 and 2207174. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d3ob00392b |
This journal is © The Royal Society of Chemistry 2023 |