Skip to main content
Download PDF

Elucidation of the glycosylation steps during biosynthesis of antitumor macrolides PM100117 and PM100118 and engineering for novel derivatives

Microbial Cell Factories volume 15, Article number: 187 (2016) Cite this article

Abstract

Background

Antitumor compounds PM100117 and PM100118 are glycosylated polyketides derived from the marine actinobacteria Streptomyces caniferus GUA-06-05-006A. The organization and characterization of the PM100117/18 biosynthesis gene cluster has been recently reported.

Results

Based on the preceding information and new genetic engineering data, we have outlined the pathway by which PM100117/18 are glycosylated. Furthermore, these genetic engineering experiments have allowed the generation of novel PM100117/18 analogues. Deletion of putative glycosyltranferase genes and additional genes presumably involved in late biosynthesis steps of the three 2,6-dideoxysugars appended to the PM100117/18 polyketide skeleton, resulted in the generation of a series of intermediates and novel derivatives.

Conclusions

Isolation and identification of the novel compounds constitutes an important contribution to our knowledge on PM100117/18 glycosylation, and set the basis for further characterization of specific enzymatic reactions, additional genetic engineering and combinatorial biosynthesis approaches.

Background

Natural products (NPs) are an invaluable source of biologically active drugs with therapeutic significance [1, 2]. A considerable portion of these products present as part of their chemical structure deoxysugar (DOH) moieties, mainly 6-deoxyhexoses (6DOH) and 2,6-dideoxyhexoses (2,6DOH) [3, 4]. Many studies aimed to unravel the structure–activity relationship of numerous NPs have shown the relevant and diverse functional roles played by DOHs on a wide range of biological activities, including antibiotic, antifungal, antitumor and antiparasitic [5–8]. In accordance with these important roles, removal of DOH moieties might lead to various degrees of reduction in activity, ranging from a complete activity loss to the retention of certain activity level as the glycosidic fraction is eliminated [4].

Glycobiology research has led to discover and delineate the biosynthesis pathway of many natural sugars [9]. As a result, a core set of enzymatic activities has been found to synthesize the diverse sugar structures observed in nature. Various common enzyme reactions, which include deoxygenations, transaminations, ketoreductions, C-, N- or O-methylations, epimerizations, and less frequently oxygenations, are required to generate the surprisingly high structural diversity of natural DOHs [3, 9]. The complete enzyme set required for DOH biosynthesis is, in most cases, encoded in the biosynthesis gene cluster of the respective NP. Therefore, because the great majority of sugars present in bioactive NPs are DOHs, identification of clusters containing genes coding for these activities may constitute a useful strategy for the discovery of clinical attractive drugs [10, 11]. Some of these enzymes have been thoroughly characterized in vitro [12–14] and in vivo [15] in multiple attempts to modify sugar structures and alter the glycosylation pattern of NP. Upon DOHs biosynthesis, transfer onto specific positions of the corresponding aglycone is controlled by glycosyltransferase (GT) enzymes. These enzymes require an activated DOH donor, typically a NDP or NMP-DOH, for transfer of the glycosyl residue to the acceptor molecule [16]. In this way, a staggering diversity of sugars can be transferred to a no less structurally diverse array of acceptor molecules, including polyketide, nonribosomal peptide, indolocarbazole and aminocoumarin aglyca among many others, by the various hundreds of putative GTs described in databases as involved in NPs biosynthesis [16] (http://www.cazy.org/GlycosylTransferases.html). Most GTs catalyze O-glycosidic bond formation between the sugar and its natural aglycone. Other linkages such as C–C [17] or C–N [18] glycosidic bonds are relatively unusual. Substrate flexibility is an interesting feature of some sugar biosynthesis and GT enzymes. This is especially remarkable in GTs which can show different degrees of promiscuity towards the acceptor and/or the donor molecules [3, 8, 19, 20]. Similarly, other GTs possesses the unusual ability to transfer sugars to different positions of the same aglycone sometimes using different type of linkages [21, 22]. GT flexibility is a key aspect for the success of many glycosylation pattern engineering strategies aimed to afford and develop new NPs. Deletion and heterologous expression of DOH biosynthesis and/or GTs genes have allowed the generation of a handful of interesting derivatives from structurally diverse NPs, such as erythromycin [23], methymicin/pikromycin [24], daunorubicin/doxorubicin [25], rebeccamycin [26] or elloramycin [27–29], among others. In some cases, novel derivatives resulting from these strategies show an enhanced bioactivity and/or improved pharmacological properties. Relevant examples of glycosylation-pattern modifications leading to these improvements are the aureolic acids demycarosyl-3D-β-d-digitoxosyl-mithramycin SK, demycarosyl-mithramycin SDK and demycarosyl-3D-β-d-digitoxosyl-mithramycin SDK, which derived from the clinically used anticancer agent mithramycin by combinatorial biosynthesis in Streptomyces argillaceus [30].

PM100117 and PM100118 are glycosylated macrolides produced by the marine actinobacteria Streptomyces caniferus GUA-06-05-006A [31, 32]. The glycosylation pattern of PM100117/18 consists of three 2,6DOHs (l-axenose, 2-deoxy-l-fucose and l-rhodinose), which along with a naphthoquinone chromophore form a side chain attached to the PM100117/18 polyketide skeleton (Fig. 1a). The remarkable antitumor activity of PM100117/18, and the generation of structural analogues with improved bioactivity recently reported [32], prompted us to undertake the engineering of these compounds looking for more derivatives. Here, we describe the generation of seven novel analogues by a glycoengineering approach consisting in the deletion of PM100117/18 biosynthetic genes previously characterized in silico as coding for putative GTs and key sugar biosynthetic enzymes (Fig. 1b). In addition, the outcome of these genetic engineering manipulations provided valuable information on the PM100117/18 glycosylation steps.

Fig. 1
figure 1

PM100117 and PM100118 2,6-DOH moieties and biosynthesis gene cluster. a PM100117/18 chemical structure. DOH moieties are indicated: a l-axenose b l-rhodinose (R = H) or 2-deoxi-l-fucose (R = OH) c l-rhodinose. b Genetic organization of PM100117/18 biosynthesis gene cluster showing genes involved in 6DOH biosynthesis (black arrows) and transfer (grey arrows)

Full size image

Results

Deletion of 2,6-dideoxysugar biosynthesis genes

The PM100117/18 gene cluster encodes most of the putative functions required to catalyze the biosynthesis of the three 2,6DOH moieties appended to the PM100117/18 macrolactone ring (Fig. 1). From these activities, as deduced in silico (Additional file 1), biosynthesis of l-axenose, 2-deoxy-l-fucose and l-rhodinose can be envisaged to occur as depicted in Fig. 2. Early biosynthesis steps of 6DOHs involves intermediacy of nucleoside diphosphate (NDP)-4-keto-6DOH, which is synthesized from a hexose-1-phospate, most possibly d-glucose, in two reaction steps catalyzed by a NDP-d-hexose synthase (GonGS) and a NDP-d-hexose-4,6-dehydratase (GonD2), respectively. To achieve, 2-deoxy-l-fucose and l-rhodinose, NDP-4-keto-6DOH could be converted into a NDP-4-keto-2,6DOH intermediate with an equatorial hydroxyl group at carbon C3, by the sequential activity of NDP-d-hexose-2,3-dehydratase (GonD1) and NDP-d-hexose-3-ketoreductase (GonR3), in accordance with the high resemblance of GonR3 to other 3-ketoreductases that render this configuration (Additional file 1). Similarly, l-axenose biosynthesis would involve a NDP-4-keto-2,6DOH intermediate with a hydroxyl group at carbon C3 in axial stereo configuration, however the corresponding NDP-d-hexose-3-ketoreductase is not present in the PM100117/18 gene cluster. Later steps towards PM100117/18 DOHs biosynthesis from this key 2,6DOH intermediate might involve four reaction types outlined in Fig. 2, which include C-methylation at C3 (GonCM), 3-deoxygenation (GonD3), 5- or 3,5-epimerization (GonE) and 4-ketoreduction (GonR1/GonR2).

Fig. 2
figure 2

Proposed biosynthesis pathway of PM100117 and PM100118 2,6-DOH moieties. 3KR: ketoreductase

Full size image

Herein, we sought to generate PM100117/18 structural analogues by manipulating the biosynthesis of their pendant 2,6DOHs. For that purpose we generated the mutant strains ΔgonCM, ΔgonD3, ΔgonE, and ΔgonR1, in which genes involved in late steps of 2,6DOH biosynthesis, gonCM, gonD3, gonE, and gonR1, respectively, were deleted. To detect the biosynthesis of PM100117/18 derivatives, whole-culture extracts from these strains were analyzed by UPLC (Fig. 3) and HPLC/MS. UPLC chromatograms from strains ΔgonCM, ΔgonD3 and ΔgonE showed suppression of PM100117 (1, UPLC R t  = 5.182 min, m/z 1601.9 [M + Na]+) and PM100118 (2, UPLC R t  = 5.536 min, m/z 1585.9 [M + Na]+) biosynthesis, and the appearance of several peaks corresponding to potential analogues (Fig. 3a). Mass assessments by HPLC/MS suggested that the novel compounds detected in ΔgonCM are compatible with PM100117/18 analogues in which the first sugar moiety, l-axenose, has been substituted for either 2-deoxy-l-fucose or l-rhodinose. This was confirmed by elucidating the chemical structure of products 3 [UPLC R t  = 6.035 min, m/z 1571.9 (M + Na)+] and 4 [UPLC R t  = 6.102 min, m/z 1555.9 (M + Na)+] by NMR, showing that both are PM100118 derivatives with an l-rhodinose unit in place of l-axenose as first sugar moiety (Fig. 3b and Additional file 2: Fig. S1 and Table S1). The presence of a hydroxyl group at C18 of the polyketide ring in 3 (absent in 4 and also in 1 and 2) suggests that compound 3 might be a shunt product of PM100117/18 biosynthesis after undergoing a hydroxylation reaction at this position. These results are consistent with the involvement of GonCM in l-axenose biosynthesis (Fig. 2). Likewise, compound 5 [UPLC R t  = 4.941 min, m/z 1487.8 (M + Na)+], produced by strain ΔgonD3 (Fig. 3a), was also purified and subjected to NMR analysis (Additional file 2: Fig. S1 and Table S2), confirming that, as expected from the hypothetical 2,6DOH biosynthetic pathways outlined in Fig. 2, this product structurally derives from PM100117 by lacking the third sugar moiety, l-rhodinose. The estimated molecular mass for compounds 6 [UPLC R t  = 5.721 min, m/z 297 (M + Na)+] and 7 [UPLC R t  = 5.796 min, m/z 239 (M + Na)+] from strain ΔgonE was strikingly low in comparison to that of the parental products, suggesting that these compounds might not be related to PM100117/18 biosynthesis. However, structure elucidation (Additional file 2: Fig. S1 and Table S3) revealed that 6 corresponds to the naphthoquinone chromophore, while 7 is the precursor of 6 to which a propionate side chain is incorporated (Fig. 3b). Masses displayed by compounds 8 [UPLC R t  = 6.012 min, m/z 1457.8 (M + Na)+] and 9 [UPLC R t  = 6.098 min, m/z 1441.8 (M + Na)+], fitted with PM100117 lacking both l-rhodinose as third sugar moiety and oxygenations at C16 and/or C17. Nonetheless, neither these structural features nor the configuration of the sugar moieties presumably present in 8 and 9 have been further confirmed because the low production level and high instability of these products, which made impossible the purification of quantities suitable for NMR analysis. On the other hand, mutant ΔgonR1 produced PM100117/18 at similar levels to those observed in wild type S. caniferus GUA-06-05-006A (Fig. 3a), indicating that this gene is not essential for PM100117/18 2,6DOHs biosynthesis. It is important to note the presence in the PM100117/18 gene cluster of an additional putative NDP-4-keto-6-deoxyhexose reductase gene, gonR2 (Fig. 1b). We set to determine whether this gene codes for the 4-ketoreductase activity presumably involved in the final biosynthesis step of the three PM100117/18 2,6DOHs (Fig. 2). Nevertheless, after recurrent attempts, we fail to produce a gonR2 deletion mutant. As an alternative approach, gonR2 was cloned under the control of the ermE*p promoter to produce pCgonR2, which was introduced in S. caniferus GUA-06-05-006A. The metabolic profile of the resultant strain, OEgonR2, was then analyzed by UPLC, showing no differences with respect to wild type strain, GUA-pSETH, carrying the respective empty plasmid (Additional file 3: Fig. S2).

Fig. 3
figure 3

Production of PM100117 (1) and PM100118 (2) derivatives 3–9 in 2,6-DOH biosynthesis gene mutant strains. a UPLC analysis (254 nm) of metabolite production in S. caniferus GUA-06-05-006A wild-type and ΔgonCM, ΔgonD3, ΔgonE and ΔgonR1 deletion mutants. NR: peak non-related to PM100117/18. b Chemical structure of various PM100117/18 analogues

Full size image

Deletion of glycosyltransferase genes

The PM100117/18 biosynthesis gene cluster contains four ORFs (gonG1, gonG2, gonG3 and gonG4) coding for putative GT proteins [32], each of which might be involved in transfer of at least one DOH moiety to a defined position of PM100117/18. As additional route to produce PM100117/18 structural analogues, we pursued to manipulate 2,6DOH transfer to the PM100117/18 macrolactone aglycone by independently deleting these putative GT genes. Thus, we produced mutant strains ΔgonG1, ΔgonG2, ΔgonG3 and ΔgonG4. Metabolite production of these strains was examined as described above. It is known that loss of the naphthoquinone chromophore unit causes a shift in the PM100117/18 maximum absorption wavelength from 254 to 216 nm [32]. Thereby, since deletion of GT genes could lead to glycosylation pattern changes involving loss of this structural element, screening for novel derivatives entailed the analysis of chromatograms at both wavelengths. Deletion of gonG1 or gonG2 abolished PM100117/18 production (Fig. 4a) and induced the accumulation of various potential analogues with a maximum absorption wavelength at 216 nm, hence lacking the naphthoquinone moiety (Fig. 4b). HPLC/MS analysis showed that novel compounds produced by strain ΔgonG1 possess masses corresponding to different variations of the PM100117/18 non-glycosylated aglycone, while masses of those produced by mutant strain ΔgonG2 matched putative PM100117/18 analogues lacking the second 2,6DOH moiety. Structural elucidation by NMR (Additional file 2: Fig. S1 and Table S4 and S5) revealed that compound 10 (UPLC R t  = 3.808 min, m/z 957.6 [M + Na]+) from ΔgonG1 indeed corresponded to the PM100117/18 aglycone, while 11 (UPLC R t  = 3.842 min, m/z 1102.4 [M + Na]+) and 12 (UPLC R t  = 3.905 min, m/z 1117.7 [M + Na]+) from strain ΔgonG2 were PM100117/18 derivatives retaining only the first sugar moiety, l-axenose (Fig. 4c). These results implicate the GTs GonG1 and GonG2 in transfer of the first and second 2,6DOH moiety, respectively, during PM100117/18 biosynthesis. UPLC and HPLC/MS analysis of ΔgonG3 metabolite biosynthesis showed the production of derivative 5 (Fig. 4b) as well as an additional compound, 13 (UPLC R t  = 5.143 min, m/z 1471.8 [M + Na]+), with absorption (254 nm) and mass spectra compatible with a PM100118 derivative lacking the third sugar moiety. Confirmation of this structure by NMR analysis (Fig. 4c; Additional file 2) endorsed the involvement of the GT GonG3 in transfer of the third 2,6DOH, l-rhodinose. Lastly, deletion of gonG4 did not show any effect either on PM100117/18 biosynthesis, which were produced at comparable levels to those observed in the wild type strain S. caniferus GUA-06-05-006A, or on the accumulation of novel PM100117/18 derivatives (Fig. 4a, b). The study of the possible role of this gene in PM100117/18 biosynthesis is addressed in section below.

Fig. 4
figure 4

Production of PM100117 (1) and PM100118 (2) biosynthesis intermediates 10–13 in glycosyltransferase gene mutant strains. UPLC analysis at 254 (a) and 216 nm (b) wavelength of metabolite production in S. caniferus GUA-06-05-006A wild-type and ΔgonG1, ΔgonG2, ΔgonG3and ΔgonG4 deletion mutants. NR: peak non-related to PM100117/18. c Chemical structure of various PM100117/18 analogues

Full size image

We conducted a series of experiments to investigate the possible functional role of the putative GT GonG4. Firstly, a precise quantification of PM100117/18 production in strain ΔgonG4 and wild type S. caniferus GUA-06-05-006A confirmed that gonG4 deletion has no impact on the production level of these compounds, suggesting that putative GT GonG4 might not be involved in their biosynthesis (Fig. 5a). To assess the potential of GonG4 to transfer 2,6DOHs to PM100117/18, we increased its presence in S. caniferus GUA-06-05-006A by over-expressing gonG4 using the integrative plasmid pSGonG4. Precise UPLC-based quantification of PM100117/18 production in the resultant strain, OEgonG4, did not show any change with respect to the wild type strain (Fig. 5b; Additional file 3: Fig. S3a). Next, we reasoned that, under certain conditions, transfer of 2-deoxy-l-fucose and l-rhodinose as second 2,6DOH units to form PM100117 and PM100118, respectively, might involve two different GTs. However no potential derivative containing the second sugar moiety was detected when gonG4 was over-expressed in ΔgonG2 (Additional file 3: Fig. S3b). On the other hand, a defensive function consisting of transferring of glucose units to antibiotic agents, which as a result become inactivated, has been proposed for several GTs [33–35]. Notwithstanding, the involvement of GonG4 in PM100117/18 inactivation by glycosylation has not been proven by experiments consisting of gonG4 over-expression and subsequent feeding with PM100118. Therefore, the role of gonG4 on PM100117/18, if any, remains elusive.

Fig. 5
figure 5

Analysis of GonG4 functional role on PM100117 and PM100118 biosynthesis. a Determination of PM100117/18 UPLC peak areas in strains S. caniferus GUA-06-05-006A (white bars) and ΔgonG4 (grey bars). b Peak area quantification of PM100117/18 in strains GUA-pSETH (white bars) and OEgonG4 (grey bars). In all cases, data are means ± standard deviations from three independent experiments

Full size image

In vitro cytotoxicity analysis of PM100117 and PM100118 derivatives

We assessed the in vitro cytotoxic activity of the novel analogues generated against cancer cell lines A549 (human lung carcinoma cells), PSN1 (pancreas carcinoma), MDA-MB-231 (human breast adenocarcinoma) and HT29 (human colorectal carcinoma). The GI50 (50% inhibition on cell growth) and TGI (total growth inhibition) concentrations of compounds 3 and 4 increased slightly with respect to PM100117/18, indicating a drop in the antitumor activity of the derivatives relative to the parental products (Table 1). The remaining novel products, 5-13, showed a more remarkable antitumor activity reduction, showing GI50 concentrations up to ten times higher than the natural compounds.

Table 1 In vitro cytotoxic activity of compounds 1–13
Full size table

Discussion

Up to thirteen ORFs from the PM100117 and PM100118 biosynthesis gene cluster can reliably be predicted to code for enzymes potentially involved in glycosylation reactions. Although bioinformatics analysis of the entire PM100117/18 gene cluster was previously addressed [32], here we provide additional in silico data supporting function assignment to putative DOH biosynthesis and GT enzymes (Additional file 1). The biosynthesis pathways of the PM100117/18 2,6DOHs moieties (Fig. 2) that can be deduced from these functions are congruent with 2,6DOHs biosynthesis reaction schemes previously reported [36]. Only a gene coding for an EryBII-type 3-ketoreductase (3KR, Fig. 2) presumably involved in the formation of a NDP-4-keto-2,6DOH intermediate with an axial hydroxyl group at C3 has not been located within the gene cluster. The absence of functions required for 2,6DOHs biosynthesis has been also described in the sequence analysis of the gene cluster of other glycosylated NPs, such as spinosyn [37], brasilicardin A [38] or mithramycin [39]. However, in the later example, gene mtmC has been shown to code for a ketoreductase/methyltrasnferase bifunctional enzyme [40, 41]. Thus, 3-C-methyltransferase GonCM might be also responsible for the 3-ketoreduction reaction leading to l-axenose. Alternatively, methylation of the NDP-4-keto-2,6-dideoxyl-d-Glucose intermediate by GonCM could be sequential to GonR3 activity (Fig. 2). Furthermore, the involvement of 3-C-methyltransferase GonCM and 3-dehydratase GonD3 in reaction steps of l-axenose and l-rhodinose biosynthesis, respectively, is also supported by genetic engineering experiments (Fig. 3). Likewise, implication of 3,5-epimerase GonE in PM100117/18 biosynthesis has been also confirmed. Nevertheless, putative isomerase activity of this enzyme has not been further verified because of the impossibility to conduct NMR analysis of compounds 8 and 9. It is noteworthy the presence in the PM100117/18 gene cluster of two genes, gonR1 and gonR2, coding for putative 4-ketoreductase enzymes. However, the involvement of these enzymes in the biosynthesis of the three PM100117/18 DOH moieties has not been fully demonstrated by genetic experiments. The inactivation of gonR1 has no apparent effect on PM100117/18 biosynthesis, and its lack could be complemented by the presence of an additional 4-ketoreductase encoded in S. caniferus chromosome, being gonR2 a good candidate. On the other hand, the impossibility of producing a gonR2 deletion mutant suggests the involvement of this gene in important primary processes such as, for example, supply of essential sugars for cell-wall formation. A similar connection of primary and secondary metabolism via DOH biosynthesis genes has been previously reported for the biosynthesis of l-rhamnose required for spinosyn A production in Saccharopolyspora spinosa [42].

Gene deletion experiments have also enabled the identification of GTs responsible for transfer of each 2,6DOH sugar moiety to its corresponding position of the PM100117/18 glycosylation pattern. Putative GT GonG1 is involved in transfer of the first sugar unit, l-axenose, since upon deletion of gonG1 no glycosylated PM100117/18 derivatives were produced (Fig. 4). Following l-axenose transfer, GonG2 is responsible for the introduction of the second sugar moiety, 2-deoxy-l-fucose and l-rhodinose, to PM100117 and PM100118 biosynthesis, respectively. Since two different 2,6DOH moieties can be transferred as second sugar unit, we previously predicted the involvement in this step of two GTs, each dedicated to transfer only one 2,6DOH. Nevertheless, UPLC and MS analysis did not detect production in ΔgonG2 strain of PM100117/18 analogues carrying the second DOH unit (Fig. 4), thus pointing at GonG2 as the only GT responsible for the introduction of both 2,6DOHs at second position of the PM100117/18 glycosylation pattern. Finally, the structure of analogues derived from ΔgonG3 (Fig. 4c), which lack the third sugar moiety, suggests that GonG3 is the GT in charge of transferring l-rhodinose to the naphthoquinone moiety of the PM100117/18 glycosylation profile. Surprisingly, the functional role of the putative GT GonG4 remains elusive. As a primary hypothesis, we speculated a role for GonG4 on transferring 2-deoxy-l-fucose or l-rhodinose as second sugar moiety. However, neither PM100117 nor PM100118 production was detected upon gonG4 over-expression in ΔgonG2 mutant strain (Fig. 5d). Alternatively, GonG4 might be inactive, but sequence alignment with functional glycosyltransferases did not reveal obvious differences supporting that possibility, or require post-translational regulatory signals for its activation. Moreover, it seems plausible also to contemplate that under certain environmental conditions, GonG4 might have a role on transferring 2,6DOH moieties other than l-axenose, 2-deoxy-l-fucose or l-rhodinose, leading to the production of natural PM100117/18 derivatives. This is consistent with the presence in the PM100117/18 gene cluster of two ORFs, gonR1 and gonR2, coding for putative 4-ketoreductase enzymes. Their presence in the cluster allows us to speculate with the biosynthesis of additional 2,6DOHs. For instance, last step of 2-deoxy-l-fucose an l-rhodinose biosynthesis could be branched towards the formation of l-olivose and l-amicetose, respectively, if either GonR1 or GonR2 were able to produce hydroxyl groups in equatorial configuration at C4. However, based on sequence analysis data, it is not possible to forecast the resulting stereo configuration after GonR1 and/or GonR2 activity. Otherwise, we cannot rule out the biosynthesis of other 2,6DOHs by the activity of enzymes encoded in other locations of S. caniferus chromosome.

Standing by the chemical structure of PM100117/18 analogues produced by the mutant strains analyzed in this work, we can make interesting conclusions on the substrate flexibility of GonG1, GonG2 and GonG3. When genes coding for these putative GTs are individually deleted, none of the two remaining GTs is flexible to catalyze the transfer of its own natural sugar to the acceptor aglycone. This apparently low promiscuity of the three GTs towards the acceptor substrate contrasts with the flexibility showed by GonG1 to transfer l-rhodinose (compounds 3 and 4) as first sugar moiety to the PM100117/18 polyketide ring (Fig. 3) in the ΔgonCM mutant, in which l-axenose biosynthesis is presumably suppressed. Another interesting observation from the structure of the compounds detected in mutants ΔgonCM and ΔgonG2 is the production of a derivative with a hydroxyl functional group at C18 of the PM100117/18 macrolactone ring (compounds 3 and 11), a feature also observed in the polyketide skeleton of the structurally related antibiotic brasilinolide A [43]. This seems to reinforce the structural and biosynthetic relationships between both compounds. On the other hand, accumulation of the naphthoquinone moiety (compound 6) is exclusively detected in strain ΔgonE, and not in other mutant strains, such as ΔgonG1 or ΔgonG2, in which this chromophore is neither incorporated into the PM100117/18 structure (Figs. 3, 4). Based on our current knowledge on PM100117/18 biosynthesis, no suitable explanation to this observation can be exposed. Moreover, the naphthoquinone chromophore and its derivative compound 7 do not show in vitro cytotoxic activity (Table 1), in spite of the structural and biosynthetic connection of these products to menaquinone, whose antitumor properties have been repeatedly noted [44, 45]. Furthermore, as it was predictable from a previous report [32], removal of large portions from PM100117/18 structure, as displayed in different degrees by compounds 5, 10, 11, 12 and 13, resulted in a remarkable decrease of in vitro cytotoxicity. Only compounds 3 and 4, which preserve chemical structures containing three 2,6DOHs, retain activity levels comparable to those showed by the parental products (Table 1).

Given their biological activity, PM100117/18 constitute appealing potential targets of biochemical, genetic and medical research. Findings on PM100117/18 biosynthesis reported in this work will enable the manipulation of the described glycosylation machinery by genetic methods in order to generate additional structural diversity of these compounds.

Conclusions

PM100117 and PM100118 exhibit a composition of structural elements with biosynthetic interest, such as a naphthoquinone chromophore and three 2,6DOHs. Previously, we showed the vital role of the naphthoquinone unit on the PM100117/18 biological activity. Herein, the 2,6DOHs moieties are proven to be equally essential to PM100117/18 cytotoxicity. Therefore, manipulations aimed to biosynthesize novel PM100117/18 antitumor derivatives must contemplate the preservation of these structural components. In addition, we have demonstrated that PM100117/18 glycosylation can be genetically manipulated to direct the biosynthesis of novel derivatives. At least one GT coded within the PM100117/18 gene cluster has shown certain flexibility towards its natural 2,6DOH donor substrate, widening future perspectives to successfully fulfill the engineering of novel analogues by combinatorial biosynthesis methods. On the other hand, the uncertain functional role of the putative GT GonG4, and the presence in the PM100117/18 gene cluster of two ORFs coding for putative 4-ketoreductase enzymes, invites to speculate with the biosynthesis, under certain environmental conditions, of additional natural products structurally related to PM100117/18. The role of these enzymes deserves further investigations, which might require performance of biochemical and enzymatic analyses.

Methods

Bacterial strains, media and culture conditions

The PM100117/18 producer S. caniferus GUA-06-05-006A [31] and its isogenic strains were routinely cultivated on solid medium A [46] or in TSB [47] medium. Cultures for PM100117/18 production experiments were carried out in supplemented MS medium (SMS) [32] at 30 °C and rotary shaking at 200 rpm for 7 days. A more detailed description of these fermentations can be found in a previous report [32]. The Escherichia coli strains used for gene cloning (DH10B) [48] and intergeneric conjugation (ET12567/pUB307) [47] were grown at 37 °C in 2xTY medium. All media were supplemented, when required, with the appropriate antibiotic for plasmid selection following standard directions [47].

Analysis of PM100117 and PM100118 production

Whole-culture samples (2 ml) from S. caniferus GUA-06-05-006A wild type, or the corresponding mutant strains, were mixed with the same volume of ethyl acetate and shacked vigorously at room temperature for 1 h. The organic phase was collected and evaporated to dryness with a Speed-Vac. The residue was re-dissolved in 60 µl of methanol:DMSO (1:1), clarified by centrifugation and then subjected to UPLC and LC/MS analysis as described elsewhere [49]. To assess the PM100117/18 production level, SMS medium was vacuum-filtrated with nitrocellulose filters of 0.45 µm pore size (Millipore). After 7-days cultivation of the respective strains in clarified SMS medium, 2-ml whole-culture samples were collected by triplicate and centrifuged. The cell pellet was washed twice with milliQ water and then incubated for 72 h at 100 °C to determine cell dry-weight of the culture. Metabolite extraction and UPLC analysis were performed as mentioned above. The final sample solution was appropriately diluted to ensure reliable linear-range PM100117/18 quantifications. Peak areas were determined by triplicate and normalized with respect to the corresponding dry-weight values.

Purification of PM100117 and PM100118 derivatives

For purification of compounds 3–13, cultures of the corresponding producing strains were centrifuged and the supernatants were filtered and applied to a solid-phase extraction cartridge (Sep-Pak Vac C18, 10 g, Waters). The retained material in each case was eluted with a mixture of methanol and 0.05% TFA in water. A linear gradient from 0 to 100% methanol in 60 min, at 10 ml/min, was used. Fractions were taken every 5 min and analyzed by UPLC. The fractions containing the desired compounds were evaporated in vacuo and subsequently re-dissolved in a small volume of a mixture of DMSO and methanol (50:50). For purification of compounds 6 and 7 the extracted broth was acidified by adding formic acid up to 1% and re-extracted with one volume of ethyl acetate. Additionally, purification of compounds 3 and 4, required extraction from cell pellets with ethyl acetate. The resulting pellet extracts were similarly evaporated and re-dissolved. The compounds of interest were purified by preparative HPLC using a SunFire C18 column (10 µm, 10 × 250 mm, Waters). Compounds were chromatographed with mixtures of acetonitrile or methanol and 0.05% TFA in water in isocratic conditions optimized for each peak, at 7 ml/min, and when needed (compounds 3, 4, 6 and 7) collected on 0.1 M phosphate buffer, pH 7.0. After every purification step, the collected compounds were diluted fourfold with water and then applied to a solid-phase extraction cartridge (Sep-Pak C18, Waters). The cartridge was washed with water, the retained compound was eluted with methanol and dried in vacuo. Once the purification was finished, the compounds were dissolved in a mixture of tert-butanol and water (1:1) and lyophilized.

In vitro cytotoxicity assay

The procedure to determine the in vitro antitumor activity of the PM100117/18 analogues obtained in this work was previously described [31, 32].

DNA manipulation and plasmids construction

DNA manipulation experiments in Escherichia coli and S. caniferus GUA-06-05-006A were carried out according to standard protocols [47, 50]. PCR amplifications were performed by using Herculase II Fusion polymerase (Agilent Technologies) with the touchdown PCR procedure previously reported [32]. Primers and plasmids used in this work are described in (Additional file 4: Table S6). To accomplish single deletion of genes gonCM, gonD3, gonE, gonR1, gonR2, gonG1, gonG2, gonG3 and gonG4, the downstream DNA sequence of the corresponding target ORFs were amplified with the primer pairs EcoRV-CM/NdeI-CM, EcoRV-D3/BamHI-D3, EcoRV-E/BamHI-E, BamHI-R1/EcoRV-R1, EcoRV-R2/BamHI-R2, BamHI-G1/EcoRV-G1, BglII-G2/EcoRV-G2, NdeI-G3/EcoRV-G3 and BamHI-G4/EcoRV-G4, respectively, and cloned in the designated restriction sites (bold) of plasmid pEFBA-oriT [51], downstream to the aac(3)IV gene, which confers resistance to apramycin (ApmR). Then, the upstream sequence of the cited target genes were likewise amplified with the primer pairs NsiI-CM/SpeI-CM, NsiI-D3/SpeI-D3, NsiI-E/SpeI-E, SpeI-R1/NsiI-R1, NsiI-R2/SpeI-R2, SpeI-G1/NsiI-G1, SpeI-G2/NsiI-G2, XbaI-G3/NsiI-G3 and SpeI-G4/NsiI-G4, respectively, and cloned upstream to the aac(3)IV gene. Finally, the hygromycin B resistance (HygR) gene marker, hyg, was extracted by XbaI/NheI digestion from plasmid pLHyg [52] and introduced in the XbaI site of the deletion plasmids. Gene hyg allows recognizing clones in which a complete gene replacement by a double cross-over has taken place (HygS ApmR) from those in which a single cross-over event has integrated the deletion plasmid into the chromosome (HygR ApmR). To construct the pSETH plasmid backbone, used to achieve gene complementation and over-expression experiments, a 1.6-kb DNA fragment was amplified from plasmid pDR4 [53] with the primer pair NcoI-dHGR/NcoI-rvHGR and cloned in the NcoI site of plasmid pSETec [49], which harbors the constitutive ermE*p promoter. Plasmids intended to accomplish gene complementation and over-expression were constructed by inserting the respective target ORF in the BamHI/EcoRV sites of pSETH (Additional file 4: Table S6). Correct fragment insertion in deletion and complementation plasmids was verified by PCR and sequencing.

Gene mutation and complementation by intergeneric conjugation

Deletion plasmids pDgonCM, pDgonD3, pDgonE, pDgonR1, pDgonG1, pDgonG2, pDgonG3 and pDgonG4 (Table S6) were introduced in S. caniferus GUA-06-05-006A by intergeneric conjugation [32] to generate deletion strains ΔgonCM, ΔgonD3, ΔgonE, ΔgonR1, ΔgonG1, ΔgonG2, ΔgonG3 and ΔgonG4, respectively. ORF replacement in the resulting mutant strains was confirmed by PCR (Additional file 4: Fig. S5) with the primer pairs indicated in Table S6 (Additional file 4: Table S6). Complementation plasmids pCgonCM, pCgonD3, pCgonE, pCgonG1, pCgonG2 and pCgonG3 were then introduced by the same procedure in the corresponding mutant strains to afford strains CgonCM, CgonD3, CgonE, CgonG1, CgonG2, CgonG3, in which PM100117/18 biosynthesis was partially restored (Additional file 5: Fig S6). Plasmids, pSETH, pCgonR2 and pCgonG4 were also transferred to S. caniferus GUA-06-05-006A to produce strains GUA-pSETH, OEgonR2 and OEgonG4, respectively. On the other hand, transfer of pSETH and pCgonG4 to mutant ΔgonG2 yielded strains ΔG2-pSETH and ΔG2-OEG4, respectively. The exconjugant selection procedure followed in this work was thoroughly described in a previous report [31].

Abbreviations

ApmR :

apramycin resistant phenotype

DMSO:

dimethyl sulfoxide

DOH:

deoxysugar

2,6DOH:

2,6-dideoxyhexoses

6DOH:

6-deoxyhexoses

GT:

glycosyltransferase

HPLC:

 high-performance liquid chromatography

HygR :

 hygromycin B resistant phenotype

Hygs :

 hygromycin B sensitive phenotype

LC:

 liquid chromatography

MS:

 mass spectrometry

NDP:

nucleoside diphosphate

NMP:

nucleoside monophosphate

HPLC:

high-performance liquid chromatography

NMR:

 nuclear magnetic resonance

NP:

natural product

ORF:

 open reading frame

Rt:

 retention time

TFA:

 trifluoroacetic acid

UPLC:

 ultra-performance liquid chromatography

References

  1. Cragg G. Natural products in drug discovery and development. J Nat Prod. 1997;60:52–60.

    Article  CAS  Google Scholar 

  2. Schmitt EK, Hoepfner D, Krastel P. Natural products as probes in pharmaceutical research. J Ind Microbiol Biotechnol. 2016;43:249–60.

    Article  CAS  Google Scholar 

  3. Salas JA, Méndez C. Biosynthesis pathways for deoxysugars in antibiotic-producing actinomycetes: isolation, characterization and generation of novel glycosylated derivatives. J Mol Microbiol Biotechnol. 2005;9:77–85.

    Article  CAS  Google Scholar 

  4. Mendez C, Luzhetskyy A, Bechthold A, Salas JA. Deoxysugars in bioactive natural products: development of novel derivatives by altering the sugar pattern. Curr Top Med Chem. 2008;8:710–24.

    Article  CAS  Google Scholar 

  5. Weymouth-Wilson AC. The role of carbohydrates in biologically active natural products. Nat Prod Rep. 1997;14:99–110.

    Article  CAS  Google Scholar 

  6. Kren V, Martínková L. Glycosides in medicine: “The role of glycosidic residue in biological activity”. Curr Med Chem. 2001;8:1303–28.

    Article  CAS  Google Scholar 

  7. Gantt RW, Peltier-Pain P, Thorson JS. Enzymatic methods for glyco(diversification/randomization) of drugs and small molecules. Nat Prod Rep. 2011;28:1811–53.

    Article  CAS  Google Scholar 

  8. Méndez C, Salas JA. Altering the glycosylation pattern of bioactive compounds. Trends Biotechnol. 2001;19:449–56.

    Article  Google Scholar 

  9. Thibodeaux CJ, Melançon CE, Liu H. Natural-product sugar biosynthesis and enzymatic glycodiversification. Angew Chem Int Ed Engl. 2008;47:9814–59.

    Article  CAS  Google Scholar 

  10. Sohng JK, Oh TJ, Kim CG. Method for cloning biosynthetic genes of secondary metabolites including deoxysugar from actinomycetes. J Biochem MolBiol. 1998;31:475–83.

    CAS  Google Scholar 

  11. Luzhetskyy A, Weiss H, Charge A, Welle E, Linnenbrink A, Vente A, Bechthold A. A strategy for cloning glycosyltransferase genes involved in natural product biosynthesis. Appl Microbiol Biotechnol. 2007;75:1367–75.

    Article  CAS  Google Scholar 

  12. Amann S, Dräger G, Rupprath C, Kirschning A, Elling L. (Chemo)enzymatic synthesis of dTDP-activated 2,6-dideoxysugars as building blocks of polyketide antibiotics. Carbohydr Res. 2001;335:23–32.

    Article  CAS  Google Scholar 

  13. Chen H, Yamase H, Murakami K, Chang C-W, Zhao L, Zhao Z, Liu H-W. Expression, purification, and characterization of two N, N-dimethyltransferases, tylM1 and desVI, involved in the biosynthesis of mycaminose and desosamine. Biochemistry. 2002;41:9165–83.

    Article  CAS  Google Scholar 

  14. Hong L, Zhao Z, Melançon CE, Zhang H, Liu H. In Vitro Characterization of the Enzymes Involved in TDP-d-Forosamine Biosynthesis in the Spinosyn Pathway of Saccharopolyspora spinosa. J Am Chem Soc. 2008;130:4954.

    Article  CAS  Google Scholar 

  15. Blanchard S, Thorson JS. Enzymatic tools for engineering natural product glycosylation. Curr Opin Chem Biol. 2006;10:263–71.

    Article  CAS  Google Scholar 

  16. Lairson LL, Henrissat B, Davies GJ, Withers SG. Glycosyltransferases: structures, functions, and mechanisms. Annu Rev Biochem. 2008;77:521–55.

    Article  CAS  Google Scholar 

  17. Mittler M, Bechthold A, Schulz GE. Structure and action of the C–C Bond-forming glycosyltransferase UrdGT2 involved in the biosynthesis of the antibiotic urdamycin. J Mol Biol. 2007;372:67–76.

    Article  CAS  Google Scholar 

  18. Salas AP, Zhu L, Sánchez C, Braña AF, Rohr J, Méndez C, Salas JA. Deciphering the late steps in the biosynthesis of the anti-tumour indolocarbazole staurosporine: sugar donor substrate flexibility of the StaG glycosyltransferase. Mol Microbiol. 2005;58:17–27.

    Article  CAS  Google Scholar 

  19. Rupprath C, Schumacher T, Elling L. Nucleotide deoxysugars: essential tools for the glycosylation engineering of novel bioactive compounds. Curr Med Chem. 2005;12:1637–75.

    Article  CAS  Google Scholar 

  20. Luzhetskyy A, Méndez C, Salas JA, Bechthold A. Glycosyltransferases, important tools for drug design. Curr Top Med Chem. 2008;8:680–709.

    Article  CAS  Google Scholar 

  21. Trefzer A, Blanco G, Remsing L, Künzel E, Rix U, Lipata F, Braña AF, Méndez C, Rohr J, Bechthold A, Salas JA. Rationally designed glycosylated premithramycins: hybrid aromatic polyketides using genes from three different biosynthetic pathways. J Am Chem Soc. 2002;124:6056–62.

    Article  CAS  Google Scholar 

  22. Dürr C, Hoffmeister D, Wohlert S-E, Ichinose K, Weber M, Von Mulert U, Thorson JS, Bechthold A. The glycosyltransferase UrdGT2 catalyzes both C- and O-glycosidic sugar transfers. Angew Chem Int Ed Engl. 2004;43:2962–5.

    Article  Google Scholar 

  23. Summers RG, Donadio S, Staver MJ, Wendt-Pienkowski E, Hutchinson CR, Katz L. Sequencing and mutagenesis of genes from the erythromycin biosynthetic gene cluster of Saccharapolyspora erythraea that are involved in l-mycarose and d-desosamine production. Microbiology. 1997;143:3251–62.

    Article  CAS  Google Scholar 

  24. Hong JSJ, Park SH, Choi CY, Sohng JK, Yoon YJ. New olivosyl derivatives of methymycin/pikromycin from an engineered strain of Streptomyces venezuelae. FEMS Microbiol Lett. 2004;238:391–9.

    CAS  Google Scholar 

  25. Madduri K, Kennedy J, Rivola G, Inventi-Solari A, Filippini S, Zanuso G, Colombo A, Gewain K, Occi J, MacNeil D, Hutchinson C. Production of the antitumor drug epirubicin (4′-epidoxorubicin) and its precursor by a genetically engineered strain of Streptomyces peucetius. Nat Biotechnol. 1998;16:291–4.

    Article  Google Scholar 

  26. Sánchez C, Zhu L, Braña AF, Salas AP, Rohr J, Méndez C, Salas JA. Combinatorial biosynthesis of antitumor indolocarbazole compounds. Proc Natl Acad Sci USA. 2005;102:461–6.

    Article  Google Scholar 

  27. Rodriguez L, Oelkers C, Aguirrezabalaga I, Braña AF, Rohr J, Méndez C, Salas JA. Generation of hybrid elloramycin analogs by combinatorial biosynthesis using genes from anthracycline-type and macrolide biosynthetic pathways. J Mol Microbiol Biotechnol. 2000;2(3):271–6.

    CAS  Google Scholar 

  28. Lombó F, Gibson M, Greenwell L, Braña AF, Rohr J, Salas JA, Méndez C. Engineering biosynthetic pathways for deoxysugars: branched-chain sugar pathways and derivatives from the antitumor tetracenomycin. Chem Biol. 2004;11:1709–18.

    Article  Google Scholar 

  29. Pérez M, Lombó F, Zhu L, Gibson M, Braña AF, Rohr J, Salas JA, Méndez C. Combining sugar biosynthesis genes for the generation of l-and d-amicetose and formation of two novel antitumor tetracenomycins. Chem Commun. 2005;12:1604–6.

    Article  Google Scholar 

  30. Núñez LE, Nybo SE, González-Sabín J, Pérez M, Menéndez N, Braña AF, Shaaban KA, He M, Morí-S F, Salas JA, Rohr J, Méndez C. A novel mithramycin analogue with high antitumor activity and less toxicity generated by combinatorial biosynthesis. J Med Chem. 2012;55:5813–25.

    Article  Google Scholar 

  31. Pérez M, Schleissner C, Fernández R, Rodríguez P, Reyes F, Zuñiga P, de la Calle F, Cuevas C. PM100117 and PM100118, new antitumor macrolides produced by a marine Streptomyces caniferus GUA-06-05-006A. J Antibiot. 2016;69(5):388–94.

    Article  Google Scholar 

  32. Salcedo RG, Olano C, Gómez C, Fernández R, Braña AF, Méndez C, de la Calle F, Salas JA. Characterization and engineering of the biosynthesis gene cluster for antitumor macrolides PM100117 and PM100118 from a marine actinobacteria: generation of a novel improved derivative. Microb Cell Fact. 2016;15:44.

    Article  Google Scholar 

  33. Vilches C, Hernandez C, Mendez C, Salas JA. Role of glycosylation and deglycosylation in biosynthesis of and resistance to oleandomycin in the producer organism, Streptomyces antibioticus. J Bacteriol. 1992;174:161–5.

    CAS  Google Scholar 

  34. Quirós LM, Aguirrezabalaga I, Olano C, Méndez C, Salas JA. Two glycosyltransferases and a glycosidase are involved in oleandomycin modification during its biosynthesis by Streptomyces antibioticus. Mol Microbiol. 1998;28:1177–85.

    Article  Google Scholar 

  35. Zhao L, Beyer NJ, Borisova SA, Liu HW. Beta-glucosylation as a part of self-resistance mechanism in methymycin/pikromycin producing strain Streptomyces venezuelae. Biochemistry. 2003;42:14794–804.

    Article  CAS  Google Scholar 

  36. Lombó F, Olano C, Salas JA, Méndez C. Sugar biosynthesis and modification. Methods Enzymol. 2009;458:277–308.

    Article  Google Scholar 

  37. Waldron C, Madduri K, Crawford K, Merlo DJ, Treadway P, Broughton MC, Baltz RH. A cluster of genes for the biosynthesis of spinosyns, novel macrolide insect control agents produced by Saccharopolyspora spinosa. Antonie van Leeuwenhoek, Int J Gen Mol Microbiol. 2000;78:385–90.

    Article  CAS  Google Scholar 

  38. Komaki H, Nemoto A, Tanaka Y, Takagi H, Yazawa K, Mikami Y, Shigemori H, Kobayashi J, Ando A, Nagata Y. Brasilicardin A, a new terpenoid antibiotic from pathogenic Nocardia brasiliensis: fermentation, isolation and biological activity. J Antibiot. 1999;52:13–9.

    Article  CAS  Google Scholar 

  39. González A, Remsing LL, Lombó F, Fernández MJ, Prado L, Braña AF, Künzel E, Rohr J, Méndez C, Salas JA. The mtmVUC genes of the mithramycin gene cluster in Streptomyces argillaceus are involved in the biosynthesis of the sugar moieties. Mol Gen Genet. 2001;264:827–35.

    Article  Google Scholar 

  40. Wang G, Pahari P, Kharel MK, Chen J, Zhu H, Vanlanen SG, Rohr J. Cooperation of two bifunctional enzymes in the biosynthesis and attachment of deoxysugars of the antitumor antibiotic mithramycin. Angew Chemie—Int Ed. 2012;51:10638–42.

    Article  CAS  Google Scholar 

  41. Chen JM, Hou C, Wang G, Tsodikov OV, Rohr J. Structural insight into MtmC, a bifunctional ketoreductase-methyltransferase involved in the assembly of the mithramycin trisaccharide chain. Biochemistry. 2015;54:2481–9.

    Article  CAS  Google Scholar 

  42. Madduri K, Waldron C, Merlo DJ. Rhamnose biosynthesis pathway supplies precursors for primary and secondary metabolism in Saccharopolyspora spinosa. J Bacteriol. 2001;183:5632–8.

    Article  CAS  Google Scholar 

  43. Tanaka Y, Komaki H, Yazawa K, Mikami Y, Nemoto A, Tojyo T, Kadowaki K, Shigemori H, Kobayashi J. Brasilinolide A, a new macrolide antibiotic produced by Nocardia brasiliensis: producing strain, isolation and biological activity. J Antibiot (Tokyo). 1997;50:1036–41.

    Article  CAS  Google Scholar 

  44. Hitomi M, Yokoyama F, Kita Y, Nonomura T, Masaki T, Yoshiji H, Inoue H, Kinekawa F, Kurokohchi K, Uchida N, Watanabe S, Kuriyama S. Antitumor effects of vitamins K1, K2 and K3 on hepatocellular carcinoma in vitro and in vivo. Int J Oncol. 2005;26:713–20.

    CAS  Google Scholar 

  45. Samykutty A, Shetty AV, Dakshinamoorthy G, Kalyanasundaram R, Zheng G, Chen A, Bosland MC, Kajdacsy-Balla A, Gnanasekar M. Vitamin K2, a naturally occurring menaquinone, exerts therapeutic effects on both hormone-dependent and hormone-independent prostate cancer cells. Evid Based Complement Altern Med. 2013;2013:287358.

    Article  Google Scholar 

  46. Fernández E, Weißbach U, Reillo CS, Braña AF, Méndez C, Rohr J, Salas JA. Identification of two genes from Streptomyces argillaceus encoding glycosyltransferases involved in transfer of a disaccharide during biosynthesis of the antitumor drug mithramycin. J Bacteriol. 1998;180:4929–37.

    Google Scholar 

  47. Kieser T, Bibb MJ, Buttner MJ, Chater KF, Hopwood DA. Practical Streptomyces genetics. Norwich: The John Innes Foundation; 2000.

    Google Scholar 

  48. Grant SGN, Jessee J, Bloom FR, Hanahan D. Differential plasmid rescue from transgenic mouse DNAs into Escherichia coli methylation-restriction mutants. Proc Natl Acad Sci USA. 1990;87:4645–9.

    Article  CAS  Google Scholar 

  49. Cano-Prieto C, García-Salcedo R, Sánchez-Hidalgo M, Braña AF, Fiedler H-P, Méndez C, Salas JA, Olano C. Genome Mining of Streptomyces sp. Tü 6176: characterization of the Nataxazole Biosynthesis Pathway. Chem BioChem. 2015;16:1461–73.

    CAS  Google Scholar 

  50. Sambrook J. W Russell D: Molecular cloning: A laboratory manual. 3rd ed. New York: Cold Spring Harbor Laboratory Press; 2001.

    Google Scholar 

  51. Horna DH, Gómez C, Olano C, Palomino-Schätzlein M, Pineda-Lucena A, Carbajo RJ, Braña AF, Méndez C, Salas JA. Biosynthesis of the RNA polymerase inhibitor streptolydigin in Streptomyces lydicus: tailoring modification of 3-methyl-aspartate. J Bacteriol. 2011;193:2647–51.

    Article  CAS  Google Scholar 

  52. Olano C, Moss SJ, Braña AF, Sheridan RM, Math V, Weston AJ, Méndez C, Leadlay PF, Wilkinson B, Salas JA. Biosynthesis of the angiogenesis inhibitor borrelidin by Streptomyces parvulus Tü4055: insights into nitrile formation. Mol Microbiol. 2004;52:1745–56.

    Article  CAS  Google Scholar 

  53. Wang W, Li X, Wang J, Xiang S, Feng X, Yang K. An engineered strong promoter for streptomycetes. Appl Environ Microbiol. 2013;79:4484–92.

    Article  CAS  Google Scholar 

Download references

Authors’ contributions

JAS, FDC, CM and CO conceived and designed the project; RGS conducted experiments and analyzed the data; RGS and CO performed sequence in silico analysis; AFB carried out compound purifications; RF performed structural elucidation of novel derivatives; RGS drafted the manuscript and CO, JAS and CM contributed to preparing the final version of the paper. All authors read and approved the final manuscript.

Acknowledgements

We thank Fundación Bancaria Caja de Ahorros de Asturias for financial support to CO.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The datasets on which the findings and conclusions of this article is based upon are all included in this manuscript and the Additional files associated with it.

Funding

This work was supported by projects INNPACTO IPT-2011-0752-900000 and BIO2015-64161-R (to J.A.S) of the Spanish Ministry of Economy and Competitiveness (MINECO).

Author information

Authors and Affiliations

  1. Departamento de Biología Funcional e Instituto Universitario de Oncología del Principado de Asturias (I.U.O.P.A), Universidad de Oviedo, Asturias, 33006, Oviedo, Spain

    Raúl García Salcedo, Carlos Olano, Alfredo F. Braña, Carmen Méndez & José A. Salas

  2. Drug Discovery Area, PharmaMar S.A., Avda. de los Reyes, Colmenar Viejo, 128770, Madrid, Spain

    Rogelio Fernández & Fernando de la Calle

Authors
  1. Raúl García Salcedo
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Carlos Olano
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Rogelio Fernández
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Alfredo F. Braña
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Carmen Méndez
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Fernando de la Calle
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. José A. Salas
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to José A. Salas.

Additional files

Additional file 1. In silico analysis of sugar biosynthesis and glycosyltransferase genes.

Additional file 2. Fig. S1. Mass and NMR spectra compounds 3-13. Tables S1–S5. NMR data of compounds 3-13.

12934_2016_591_MOESM3_ESM.pdf

Additional file 3. Fig. S2. Effect of gonR2 over-expression on PM100117/18 biosynthesis. Fig. S3. Effects of gonG4 over-expression on S. caniferus GUA-06-05-006A and ΔgonG2. Methods: Feeding experiment. Fig. S4. Heterologous expression of gonG4 in Streptomyces albus J1074 and feeding with PM100118.

12934_2016_591_MOESM4_ESM.pdf

Additional file 4. Table S6. Primers and plasmids used in this work. Fig. S5. Gene deletion confirmation by PCR in mutant strains ∆gonG1, ∆gonG2, ∆gonG3, ∆gonG4, ∆gonR1, ∆gonE, ∆gonCM and ∆gonD3. 

Additional file 5. Fig. S6. Genetic complementation of mutant strains ΔgonCM, ΔgonD3, ΔgonE, ΔgonG1, ΔgonG2 and ΔgonG3.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Salcedo, R.G., Olano, C., Fernández, R. et al. Elucidation of the glycosylation steps during biosynthesis of antitumor macrolides PM100117 and PM100118 and engineering for novel derivatives. Microb Cell Fact 15, 187 (2016). https://doi.org/10.1186/s12934-016-0591-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12934-016-0591-7

Keywords

Microbial Cell Factories

ISSN: 1475-2859

Contact us