Metabolic network model guided engineering ethylmalonyl-CoA pathway to improve ascomycin production in Streptomyces hygroscopicus var. ascomyceticus
© The Author(s) 2017
Received: 14 June 2017
Accepted: 26 September 2017
Published: 3 October 2017
Ascomycin is a 23-membered polyketide macrolide with high immunosuppressant and antifungal activity. As the lower production in bio-fermentation, global metabolic analysis is required to further explore its biosynthetic network and determine the key limiting steps for rationally engineering. To achieve this goal, an engineering approach guided by a metabolic network model was implemented to better understand ascomycin biosynthesis and improve its production.
The metabolic conservation of Streptomyces species was first investigated by comparing the metabolic enzymes of Streptomyces coelicolor A3(2) with those of 31 Streptomyces strains, the results showed that more than 72% of the examined proteins had high sequence similarity with counterparts in every surveyed strain. And it was found that metabolic reactions are more highly conserved than the enzymes themselves because of its lower diversity of metabolic functions than that of genes. The main source of the observed metabolic differences was from the diversity of secondary metabolism. According to the high conservation of primary metabolic reactions in Streptomyces species, the metabolic network model of Streptomyces hygroscopicus var. ascomyceticus was constructed based on the latest reported metabolic model of S. coelicolor A3(2) and validated experimentally. By coupling with flux balance analysis and using minimization of metabolic adjustment algorithm, potential targets for ascomycin overproduction were predicted. Since several of the preferred targets were highly associated with ethylmalonyl-CoA biosynthesis, two target genes hcd (encoding 3-hydroxybutyryl-CoA dehydrogenase) and ccr (encoding crotonyl-CoA carboxylase/reductase) were selected for overexpression in S. hygroscopicus var. ascomyceticus FS35. Both the mutants HA-Hcd and HA-Ccr showed higher ascomycin titer, which was consistent with the model predictions. Furthermore, the combined effects of the two genes were evaluated and the strain HA-Hcd-Ccr with hcd and ccr overexpression exhibited the highest ascomycin production (up to 438.95 mg/L), 1.43-folds improvement than that of the parent strain FS35 (305.56 mg/L).
The successful constructing and experimental validation of the metabolic model of S. hygroscopicus var. ascomyceticus showed that the general metabolic network model of Streptomyces species could be used to analyze the intracellular metabolism and predict the potential key limiting steps for target metabolites overproduction. The corresponding overexpression strains of the two identified genes (hcd and ccr) using the constructed model all displayed higher ascomycin titer. The strategy for yield improvement developed here could also be extended to the improvement of other secondary metabolites in Streptomyces species.
Ascomycin (FK520) is a 23-membered polyketide macrolide antibiotic with immunosuppressant activity, which has similar characterization to FK506 (tacrolimus), exhibiting various therapeutic applications, including antifungal , immunosuppressive , antimalarial , antispasmodics  and several other uses [4, 5]. Meanwhile, many of its derivatives also have important pharmacological value, such as pimecrolimus, which has been confirmed to be effective in the treatment of various inflammatory skin diseases [4, 6]. Due to its broad clinical applicability and promising pharmaceutical prospects, ascomycin and its derivatives have recently attracted much attention of many researchers.
Recently, considerable efforts have been made to enhance the ascomycin titer of the production strains, including random mutagenesis and selection , fermentation technology optimization , comparative metabolic profile analysis , and so on. However, the relatively low yield of ascomycin in bio-fermentations is still a major obstacle for its further industrialization and commercialization. This problem may be due to our poor understanding of the mechanism of ascomycin overproduction as well as a lack of methods to efficiently identify the key-limiting steps as targets for strain engineering. Recently, genome-scale metabolic network model (GSMM)  has been developed as an effective strategy to analyze the intracellular metabolic behavior and identify metabolic engineering targets for efficient strain improvement. It has been successfully applied to the design of metabolic engineering strategies for a number of valuable products, such as succinic acid , spinosad , tacrolimus , riboflavin , dicarboxylic acid , putrescine , and so on.
Streptomyces coelicolor A3(2), the model organism of the Streptomyces genus, is the first one to be completely sequenced species of the genus Streptomyces . The GSMM of S. coelicolor A3(2) was first constructed in 2005 (iIB711) , and was updated in 2010 (iMA789)  and 2014 (iMK1208) . During the construction of the model iIB711, the authors analyzed the universality of the model by comparing the genome of S. coelicolor A3(2) with that of Streptomyces avermitilis . It was found that 78% ORFs in the model were synteny conserved with S. avermitilis and some of the non-conserved ORFs might have synteny-conserved isoenzymes with S. avermitilis. The results illustrated that metabolic genes were highly conservative in the above two Streptomyces strains. In addition, it has been demonstrated that the size of panmetabolism (the set of all reactions in the investigated strains) in Escherichia coli was limited to the set of reactions already known and approached a constant . This phenomenon may be due to the conservation of metabolic genes in these organisms, as well as our limited knowledge of gene functions and protein biochemical characteristics [21, 22]. However, in Streptomyces species, the size of panmetabolism might not approach a plateau because they contain large number of genes in their secondary metabolism. Nevertheless, the primary metabolic genes in Streptomyces species are more conserved, which has been illustrated by comparing genome and pangenome analysis [23, 24]. Thus, one can expect that a newly sequenced Streptomyces strain would not significantly expand the size of the set of known primary metabolic reactions in Streptomyces species. That is, the construction of a general metabolic network model for Streptomyces strains only requires the addition of limited diversity to the current primary metabolic network of the model strain S. coelicolor A3(2).
Nevertheless, the metabolic conservation of Streptomyces species has so far not been fully analyzed, and the primary reactions in model iMK1208 of S. coelicolor A3(2) have not been investigated to establish if these reactions are conserved and could be applied to construct a general metabolic network model for Streptomyces strains. To resolve these issues, the metabolic enzymes in the updated model iMK1208 of S. coelicolor A3(2) were compared with those of 31 Streptomyces strains. And high conservation of the primary metabolic enzymes had been found in the all selected Streptomyces strains. Then, a metabolic network model of S. hygroscopicus var. ascomyceticus for ascomycin biosynthesis was attempted to construct by using the universal primary metabolic reactions from model iMK1208 and adding the reactions associated with the synthesis of ascomycin in S. hygroscopicus var. ascomyceticus. Using the thus constructed metabolic network model, some potential targets for ascomycin overproduction were predicted with the aid of FBA  and MOMA . Among all the potential targets, two target genes (i.e. hcd and ccr) located in the ethylmalonyl-CoA pathway showed high promoting effects on ascomycin overproduction. After gene overexpression manipulation, all the engineered strains had shown higher ascomycin production than control, and the highest ascomycin titer reached to 438.95 mg/L in the strain HA-Hcd-Ccr. These findings not only suggested a general metabolic network model could be construed for Streptomyces species, but also provided an effective method for the titer improvement of other products in different Streptomyces species.
Conservation analysis of the proteins in iMK1208 against Streptomyces species
A total of 31 Streptomyces strains with completely sequenced genomes (deposited in the NCBI by September 2015) were used in the conservation analysis. All the available protein coding sequences of the 31 Streptomyces strains were retrieved from the National Center for Biotechnology Information (NCBI), and plasmids were omitted from the analysis. The accession numbers of all 31 strains are listed in Additional file 1. The protein sequences of S. coelicolor A3(2) in metabolic model iMK1208 were used as the reference, and the sequence similarities on protein level were calculated from a pair-wise comparison of homologous sequences between the protein sequences in the model iMK1208 and the protein coding sequences of the 31 Streptomyces strains using BLAST v2.2.30+. The top hit for each protein sequence (i.e. lowest e value) was selected as the optimal comparison result. The protein sequences were considered as conserved at the threshold of 40% identity, 50% of query protein length aligned and e value < 1e − 05.
Metabolic network construction
The initial model of S. hygroscopicus var. ascomyceticus FS35 was constructed based on the model iMK1208 of S. coelicolor A3(2). Several secondary metabolism subsystems actinorhodin, undecylprodigiosin, calcium dependent antibiotic, and germicidin biosynthesis were removed from the model iMK1208. Some specific secondary metabolites biosynthesis reactions in other subsystems such as membrane transport, cofactor and prosthetic group biosynthesis were also removed from the model. The reactions related with ascomycin specific precursors biosynthesis and overall biosynthesis reactions were added to the new model.
As a kind of constraints-based metabolic model, the necessary constraints were used to obtain the optimal flux solution space. For example, 6-phosphogluconate dehydrogenase was shown to be NADP-dependent and no detectable activity with NAD , the upper bounds of reaction GND2 (NAD-dependent) and GND (NADP-dependent) were set to 0 and 1000, respectively. Based on the way for glucose entry and phosphorylation, the upper bound of reaction G1DH was set to 0 [27, 28]. According to the analysis of the carbon-flux distribution in S. coelicolor , S. hygroscopicus var. ascomyceticus  and S. lividans , the lower bound of reaction PGL (catalyzed by 6-phosphogluconolactonase) was set to 0.2, which means that the lowest flux through pentose phosphate pathway was 20% of the available glucose-6-phosphate. Additionally, the lower bound of the anaplerotic reaction PPC was set to 0.1 according to the carbon-flux distribution of S. coelicolor . Menaquinone biosynthesis pathway was also further improved according the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The details of the all reactions and metabolites in the models are provided in Additional file 2.
Identification of the overexpression targets
Bacterial strains, plasmids, and cultural conditions
Bacterial strains and plasmids used in this study
S. hygroscopicus var. ascomyceticus FS35
Mutant, improved ascomycin producer
FS35 transformed with pIB139
FS35 transformed with pIB/ccr
FS35 transformed with pIB/hcd
FS35 transformed with pIB/hcd/ccr
E. coli DH5α
Plasmid construction and general cloning
E. coli ET12567
pUZ8002, nonmethylating plasmid donor, CmR, KanR
E. coli general cloning vector, AmpR
Integrative E. coli-Streptomyces shuttle plasmid containing oriT, attP, int, aac(3)IV and ermEp*
pIB139 based integrative plasmid containing ccr gene, AprR
pIB139 based integrative plasmid containing hcd gene, AprR
pUC18 based plasmid containing hcd gene
pUC18 based plasmid containing hcd and ccr gene
pIB139 based integrative plasmid containing hcd-ccr gene, AprR
Cloning, plasmid construction and transformation
General DNA manipulation and intergeneric conjugation of E. coli with Streptomyces were performed according to standard protocols . All the strains and primers used in this work are listed in Additional file 3: Table S1. A derivative of pSET152 , pIB139, was used for hcd and ccr gene overexpression. The plasmids pIB/ccr and pIB/hcd were constructed as follows. The coding sequence of hcd and ccr were amplified from the genomic DNA of S. coelicolor A3(2) by PCR using the primer pairs hcdF1/hcdR1 and ccrF1/ccrR1 respectively. The hcd and ccr PCR products were digested using the restriction enzymes NdeI and XbaI, and the resulting fragments produced by digestion were cloned into the pIB139 vector between NdeI and XbaI respectively to obtain pIB139R and pIB139D, respectively. The constructed plasmids pIB/hcd and pIB/ccr were both confirmed by DNA sequencing. The conformed plasmids pIB/hcd and pIB/ccr were introduced into the parent strain S. hygroscopicus var. ascomyceticus FS35 to obtain the hcd and ccr overexpression strains HA-Hcd and HA-Ccr, respectively. To construct the vector pIB/hcd/ccr, the primer pairs HcdF2 and HcdR2 were used for amplifying the gene hcd and the primer pairs CcrF2 and CcrR2 were used for amplifying the gene ccr. The hcd PCR product was digested using NdeI and PstI, and then the fragment was cloned into the pUC18 vector to form pUC/hcd. The resulted plasmid pUC/hcd was cut using PstI and XbaI, and linked with the ccr PCR product, which was digested with the same restriction enzymes, to form the plasmid pUC/hcd/ccr. The plasmid pUC/ccr/hcd was excised by NdeI/XbaI and the fragment including gene hcd and ccr was cloned into pIB139 to generate the plasmid pIB/hcd/ccr. The resulting plasmid was introduced into the parental strain S. hygroscopicus var. ascomyceticus FS35 via conjugal transfer, generating the overexpression strain HA-Hcd-Ccr. The parental strain S. hygroscopicus var. ascomyceticus FS35 with pIB139 (FS-PIB) was used as the negative control. All the positive mutants were confirmed by PCR amplification and DNA sequence analysis.
The culture samples for enzyme activity analysis in vitro were harvested at the end of exponential phase (72 h) by centrifugation at 8000×g for 10 min at 4 °C. Cell pellets were washed 3 times by 100 mM Tris·HCl (pH 7.0) containing 5 mM MnSO4, 20 mM KCl, 2 mM DTT and 0.1 mM EDTA, and then resuspended in the same buffer . Cell suspension was treated by ultrasonication for 10 min at 250 W for 5 cycles of 1 min each on ice bath. Supernatant obtained after centrifugation at 12,000×g for 20 min at 4 °C was stored at − 80 °C for assays of enzyme activity and total protein concentration. The concentration of total protein was detected by the Bradford method . One unit (1 U) of enzyme activity was defined as the transformation of 1 μM/min.
3-Hydroxybutyryl-CoA dehydrogenase (Hcd, EC:220.127.116.11) activity was assayed by monitoring the decrease in NADH concentration at 340 nm using acetoacetyl-CoA as the substrate, the reaction mixture contained 100 mM Tris–HCl buffer (pH 7.5), 100 μM acetoacetyl-CoA, 150 μM NADH, and extract. The reaction was started by the addition of acetoacetyl-CoA [46, 47]. Activity of the crotonyl-CoA carboxylase/reductase (Ccr, EC:18.104.22.168) was assayed by the method reported previously .
Fermentations of S. hygroscopicus var. ascomyceticus FS35 and its derivatives were carried out as described previously . Biomass concentration was determined by the measurement of dry cell weight (DCW). For the determination, 10 mL fermentation broth were sampled and then centrifuged at 8000×g for 10 min, washed once with 0.1 M HCl and twice with distilled water separately, and dried at 80 °C to constant weight. The concentration of ascomycin was quantitated by high-performance liquid chromatography (HPLC) analysis as described previously . The residual total sugar in fermentation broth was quantified according to the phenol–sulfuric acid method using glucose as the standard .
Metabolic conservation analysis in Streptomyces species
Metabolic network reconstruction and in silico identification of potential targets
According to the previous described methods (seen in the “Methods” section: Metabolic network construction), the metabolic network model of S. hygroscopicus var. ascomyceticus FS35 was constructed based on the model of S. coelicolor A3(2), and the details are listed in the Additional file 2. Using the maximum growth rate as an objective function, the constructed model was verified by comparing the specific growth rate from the simulation with experimental data. In the model simulation, the specific glucose uptake rate was set to 0.90 mM/g DCW/h and the specific ascomycin synthetic rate was set to 5 × 10−4 mM/g DCW/h respectively. Based on the above constraints, the maximal specific growth rate calculated by the model was 0.0653/h, which was very close to the experimental result of 0.0616/h (Additional file 3: Figure S3).
Improving ascomycin production by engineering the ethylmalonyl-CoA pathway
Specific activity of enzymes by parent strain S. hygroscopicus var. ascomyceticus FS35 and recombinants in batch cultures
Enzyme activites (U/mg protein)
0.25 ± 0.03
0.08 ± 0.01
0.49 ± 0.05
0.13 ± 0.01
0.27 ± 0.03
0.28 ± 0.03
0.55 ± 0.05
0.35 ± 0.04
In previous reports, GSMM was confirmed as an effective method for better understanding of intracellular metabolic behavior and identifying the key limiting steps for target metabolites biosynthesis. In Streptomyces species, the GSMM of S. coelicolor A3(2) has been constructed and updated for 2 times and successfully used in the guide of metabolic engineering targets identity [18–20]. It has also been used as a model for guiding the constructing other metabolic network models for some Streptomyces species [13, 59]. However, though several GSMM had been constructed in some Streptomyces species, the conservation of these metabolic reactions existed in the reported GSMM had not been investigated completely. Previous studies showed that primary metabolism might be highly conserved between Streptomyces species [23, 24]. To investigate whether the primary metabolic reactions in the metabolic network model of S. coelicolor A3(2) are conserved in other Streptomyces species, a more comprehensive, deeper metabolic gene analysis of Streptomyces species was performed in this study.
The metabolic conservation of Streptomyces species was investigated by comparing the metabolic enzymes in the latest reported model of S. coelicolor A3(2) with those of 31 completely sequenced Streptomyces strains. The results showed that the metabolic enzymes were highly conserved within the 31 Streptomyces strains, especially for primary metabolism. This was consistent with the previous comparative genome and pangenome analyses of Streptomyces species [18, 23, 24]. Furthermore, only minor differences were founded in the primary metabolism of Streptomyces species, and metabolic reactions showed higher conservation than the data obtained through protein homology analysis. Therefore, majority of the primary metabolic reactions in the model of S. coelicolor A3(2) are likely universal across Streptomyces species. In addition, due to the conservation of metabolic genes and our limited knowledge of gene functions and biochemical characteristics, the set of primary metabolic reactions in Streptomyces strains is expected to converge to a constant value with the addition of analyzed genomes, as was demonstrated for E. coli strains [21, 22]. Considering the universal nature of reactions in the primary metabolism and the fact that strain-specific portions of the primary metabolism remain largely uncharacterized in Streptomyces species, it is reasonable to reference the primary metabolic reactions in the model of S. coelicolor A3(2) to construct general metabolic network models for other Streptomyces strains.
As an important application of the general metabolic model, the GSMM of S. hygroscopicus var. ascomyceticus was firstly attempted to be constructed and applied to identify the genetic engineering targets for improving ascomycin production. The constructed metabolic model of S. hygroscopicus var. ascomyceticus facilitated the understanding of the ascomycin biosynthesis network. Using this model in conjunction with the FBA and MOMA algorithms, we identified 29 potential gene overexpression targets according to a ranking of their respective f PH values. Among all the predicted targets, several targets with high f PH values (ACACT1r, HACD1, ECOAH1 and CCCR) were found in the ethylmalonyl-CoA pathway, which had also been confirmed to be an limiting pathway for other products synthesis in some Streptomyces species . In fact, we also found that the ascomycin producing strain of S. hygroscopicus var. ascomyceticus FS35 could utilize acetate and ethylmalonic acid as a sole carbon source (data not shown), indicating the presence of the ethylmalonyl-CoA pathway in this strain. However, the ethylmalonyl-CoA pathway was incomplete in the most reported model of Streptomyces species. According to the conservative analysis of primary metabolism, the related metabolic reactions had been supplied into the general metabolic model by adding additional constraints and the missing reactions in the pathway (Additional file 4). Using the updated model, the targets related to the ethylmalonyl-CoA pathway were still predicted to have great relevance for ascomycin biosynthesis (Additional file 4). Therefore, the targets related to the ethylmalonyl-CoA pathway were successfully predicted to have great relevance for ascomycin biosynthesis (Additional file 4).
The two target genes hcd and ccr were overexpressed and co-overexpressed in S. hygroscopicus var. ascomyceticus. The experimental results suggested that all the engineered strains had a higher ascomycin production than control, and the highest ascomycin production could reach to 438.95 mg/L. The genetic manipulation about the target genes further verified the feasibility the strategy about model construction. However, the increase of production titer was different from the predicted improvements. This may be due to some limitations of the constraint-based metabolic network simulation, such as the lack of detailed knowledge on metabolic regulation, as well as thermodynamic and kinetic information. The strategy for metabolic model construction in Streptomyces species will make the construction process of metabolic network easier and make it possible for the strain without genomic information.
In addition, except for S. hygroscopicus var. ascomyceticus, the ethylmalonyl-CoA pathway is also found in many Streptomyces species, and the corresponding enzymes can be found in 30 of the 31 investigated Streptomyces strains (Additional file 3: Table S3). Enhancing the flux in ethylmalonyl-CoA pathway would improve the pools of many common precursor of in polyketide biosynthesis (Fig. 6), including ethylmalonyl-CoA, methylmalonyl-CoA, propionyl-CoA, crotonyl-CoA, butyryl-CoA, and so on . In addition, the ethylmalonyl-CoA pathway overlaps with several important metabolic pathways, such as the PHB pathway, TCA cycle and serine cycle, and many PKS gene clusters associated with butyryl-containing metabolites also contain another copy of the ccr gene to supply ethylmalonyl-CoA [60, 62]. The complicated relationship between the ethylmalonyl-CoA pathway and other metabolic pathways might open new avenues for the improvement of polyketide production and rational pathway design.
In this study, using the latest reported model of S. coelicolor A3(2) as model, a general metabolic model of Streptomyces species was attempted to be constructed based on the conservative analysis of the primary metabolic reactions of 31 Streptomyces strains. Using this general metabolic model, the GSMM of S. hygroscopicus var. ascomyceticus was reconstructed and applied to identify the potential targets for ascomycin overproduction. According to the simulated results, two targets hcd and ccr from the ethylmalonyl-CoA pathway were selected for overexpression, the experimental results indicated ascomycin production showed higher improvement in both of single- and double-overexpression strains. These findings could not only indicate that improving the supply of ethylmalonyl-CoA could improve ascomycin yield effectively, but also suggest that a general metabolic model could be existed and applied to guide for the yield improvement of other products in Streptomyces species.
JW carried out the experimental work, analyzed the data and wrote the manuscript. CW performed the partial of metabolic network model reconstruction and experimental analysis, helped in editing and revising the manuscript. KS helped in analyzing the data and writing the manuscript. JPW designed the experiments and supervised the work. All authors read and approved the final manuscript.
We want to thank Dr. Huanhuan Liu and Jiao Liu for the comments on our manuscript.
The authors declare that they have no competing interests.
Availability of data and supporting materials
Consent for publication
All authors consent to publish this manuscript.
Ethics approval and consent to participate
This manuscript does not involve the data collected from animals or humans.
This work was supported by the National 973 Project of China (No. 2013CB733600), the National Natural Science Foundation of China (No. 21676189), the Key Technologies R&D Program of Tianjin (No. 16YFZCSY00780) and the Key Program of National Natural Science Foundation of China (No. 21236005).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis 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.
- Arai T, Kouama Y, Suenaga T, Honda H. Ascomycin, an antifungal antibiotic. J Antibiot. 1962;15:231–2.Google Scholar
- Dumont FJ, Staruch MJ, Koprak SL, Siekierka JJ, Lin CS, Harrison R, et al. The immunosuppressive and toxic effects of FK-506 are mechanistically related: pharmacology of a novel antagonist of FK-506 and rapamycin. J Exp Med. 1992;176:751–60.View ArticleGoogle Scholar
- Monaghan P, Fardis M, Revill WP, Bell A. Antimalarial effects of macrolactones related to FK520 (Ascomycin) are independent of the immunosuppressive properties of the cmpounds. J Infect Dis. 2005;191:1342–9.View ArticleGoogle Scholar
- Jacobi A, Braeutigam M, Mahler V, Schultz E, Hertl M. Pimecrolimus 1% cream in the treatment of facial psoriasis: a 16-week open-label study. Dermatology. 2008;216:133–6.View ArticleGoogle Scholar
- Sierra-Paredes G, Sierra-Marcuno G. Ascomycin and FK506: pharmacology and therapeutic potential as anticonvulsants and neuroprotectants. CNS Neurosci Ther. 2008;14:36–46.View ArticleGoogle Scholar
- Boone B, Ongenae K, Van Geel N, Vernijns S, De Keyser S, Naeyaert JM. Topical pimecrolimus in the treatment of vitiligo. Eur J Dermatol. 2007;17:55–61.Google Scholar
- Qi H, Xin X, Li S, Wen J, Chen Y, Jia X. Higher-level production of ascomycin (FK520) by Streptomyces hygroscopicus var. ascomyceticus irradiated by femtosecond laser. Biotechnol Bioprocess E. 2012;17:770–9.View ArticleGoogle Scholar
- Parveen K, Kumar S, Kumar M. An improved fermentation process for preparing ascomycin. Patent WO 2007/29082 A2. 2007.Google Scholar
- Qi H, Zhao S, Fu H, Wen J, Jia X. Enhancement of ascomycin production in Streptomyces hygroscopicus var. ascomyceticus by combining resin HP20 addition and metabolic profiling analysis. J Ind Microbiol Biotechnol. 2014;41:1365–74.View ArticleGoogle Scholar
- Thiele I, Palsson BO. A protocol for generating a high-quality genome-scale metabolic reconstruction. Nat Protoc. 2010;5:93–121.View ArticleGoogle Scholar
- Lee KY, Park JM, Kim TY, Yun H, Lee SY. The genome-scale metabolic network analysis of Zymomonas mobilis ZM4 explains physiological features and suggests ethanol and succinic acid production strategies. Microb Cell Fact. 2010;9:94.View ArticleGoogle Scholar
- Wang X, Zhang C, Wang M, Lu W. Genome-scale metabolic network reconstruction of Saccharopolyspora spinosa for spinosad production improvement. Microb Cell Fact. 2014;13:41.View ArticleGoogle Scholar
- Huang D, Li S, Xia M, Wen J, Jia X. Genome-scale metabolic network guided engineering of Streptomyces tsukubaensis for FK506 production improvement. Microb Cell Fact. 2013;12:52.View ArticleGoogle Scholar
- Hao T, Han BB, Ma HW, Fu J, Wang H, Wang ZW, et al. In silico metabolic engineering of Bacillus subtilis for improved production of riboflavin, Egl-237, (R,R)-2,3-butanediol and isobutanol. Mol BioSyst. 2013;9:2034–44.View ArticleGoogle Scholar
- Mishra P, Park GY, Lakshmanan M, Lee HS, Lee H, Chang MW, et al. Genome-scale metabolic modeling and in silico analysis of lipid accumulating yeast Candida tropicalis for dicarboxylic acid production. Biotechnol Bioeng. 2016;113:1993–2004.View ArticleGoogle Scholar
- Park JM, Park HM, Kim WJ, Kim HU, Kim TY, Lee SY. Flux variability scanning based on enforced objective flux for identifying gene amplification targets. BMC Syst Biol. 2012;6:106.View ArticleGoogle Scholar
- Bentley SD, Chater KF, Cerdeno-Tarraga AM, Challis GL, Thomson NR, James KD, et al. Complete genome sequence of the model actinomycete Streptomyces coelicolor A3(2). Nature. 2002;417:141–7.View ArticleGoogle Scholar
- Borodina I, Krabben P, Nielsen J. Genome-scale analysis of Streptomyces coelicolor A3(2) metabolism. Genome Res. 2005;15:820–9.View ArticleGoogle Scholar
- Alam MT, Merlo ME, Hodgson DA, Wellington EM, Takano E, Breitling R. Metabolic modeling and analysis of the metabolic switch in Streptomyces coelicolor. BMC Genom. 2010;11:202.View ArticleGoogle Scholar
- Kim M, Sang Yi J, Kim J, Kim JN, Kim MW, Kim BG. Reconstruction of a high-quality metabolic model enables the identification of gene overexpression targets for enhanced antibiotic production in Streptomyces coelicolor A3 (2). Biotech J. 2014;9:1185–94.View ArticleGoogle Scholar
- Vieira G, Sabarly V, Bourguignon P-Y, Durot M, Le Fèvre F, Mornico D, et al. Core and panmetabolism in Escherichia coli. J Bacteriol. 2011;193:1461–72.View ArticleGoogle Scholar
- Baumler DJ, Peplinski RG, Reed JL, Glasner JD, Perna NT. The evolution of metabolic networks of E. coli. BMC Syst Biol. 2011;5:182.View ArticleGoogle Scholar
- Nett M, Ikeda H, Moore BS. Genomic basis for natural product biosynthetic diversity in the actinomycetes. Nat Prod Rep. 2009;26:1362–84.View ArticleGoogle Scholar
- Kim JN, Kim Y, Jeong Y, Roe JH, Kim BG, Cho BK. Comparative genomics reveals the core and accessory genomes of Streptomyces species. J Microbiol Biotechnol. 2015;25:1599–605.View ArticleGoogle Scholar
- Kauffman KJ, Prakash P, Edwards JS. Advances in flux balance analysis. Curr Opin Biotech. 2003;14:491–6.View ArticleGoogle Scholar
- Segre D, Vitkup D, Church GM. Analysis of optimality in natural and perturbed metabolic networks. Proc Natl Acad Sci USA. 2002;99:15112–7.View ArticleGoogle Scholar
- Coze F, Gilard F, Tcherkez G, Virolle M-J, Guyonvarch A. Carbon-flux distribution within Streptomyces coelicolor metabolism: a comparison between the actinorhodin-producing strain M145 and Its non-producing derivative M1146. PLoS ONE. 2013;8:e84151.View ArticleGoogle Scholar
- van Wezel GP, Konig M, Mahr K, Nothaft H, Thomae AW, Bibb M, et al. A new piece of an old jigsaw: glucose kinase is activated posttranslationally in a glucose transport-dependent manner in Streptomyces coelicolor A3(2). J Mol Microbiol Biotechnol. 2007;12:67–74.View ArticleGoogle Scholar
- Qi H, Lv M, Song K, Wen J. Integration of parallel 13C-labeling experiments and in silico pathway analysis for enhanced production of ascomycin. Biotechnol Bioeng. 2017;114:1036–44.View ArticleGoogle Scholar
- Avignone Rossa C, White J, Kuiper A, Postma PW, Bibb M, Teixeira de Mattos MJ. Carbon flux distribution in antibiotic-producing chemostat cultures of Streptomyces lividans. Metab Eng. 2002;4:138–50.View ArticleGoogle Scholar
- Schellenberger J, Que R, Fleming RMT, Thiele I, Orth JD, Feist AM, et al. Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0. Nat Protoc. 2011;6:1290–307.View ArticleGoogle Scholar
- Gilat A. MATLAB: an introduction with applications. New York: Wiley; 2009.Google Scholar
- Valipour M, Sefidkouhi MAG, Raeini M. Selecting the best model to estimate potential evapotranspiration with respect to climate change and magnitudes of extreme events. Agric Water Manag. 2017;180:50–60.View ArticleGoogle Scholar
- Valipour M, Montazar AA. An evaluation of SWDC and WinSRFR models to optimize of infiltration parameters in furrow irrigation. Am J Sci Res. 2012;69:128–42.Google Scholar
- Valipour M. Increasing irrigation efficiency by management strategies: cutback and surge irrigation. ARPN J Agric Biol Sci. 2013;8:35–43.Google Scholar
- Heavner BD, Smallbone K, Barker B, Mendes P, Walker LP. Yeast 5—an expanded reconstruction of the Saccharomyces cerevisiae metabolic network. BMC Syst Biol. 2012;6:55.View ArticleGoogle Scholar
- Boghigian BA, Armando J, Salas D, Pfeifer BA. Computational identification of gene over-expression targets for metabolic engineering of taxadiene production. Appl Microbiol Biotechnol. 2012;93:2063–73.View ArticleGoogle Scholar
- Qi H, Zhao S, Wen J, Chen Y, Jia X. Analysis of ascomycin production enhanced by shikimic acid resistance and addition in Streptomyces hygroscopicus var. ascomyceticus. Biochem Eng J. 2014;82:124–33.View ArticleGoogle Scholar
- Kieser T. Practical Streptomyces genetics. London: John Innes Foundation; 2000.Google Scholar
- MacNeil DJ, Gewain KM, Ruby CL, Dezeny G, Gibbons PH, MacNeil T. Analysis of Streptomyces avermitilis genes required for avermectin biosynthesis utilizing a novel integration vector. Gene. 1992;111:61–8.View ArticleGoogle Scholar
- Martinez E, Bartolomé B, de la Cruz F. pACYC184-derived cloning vectors containing the multiple cloning site and lacZα reporter gene of pUC8/9 and pUC18/19 plasmids. Gene. 1988;68:159–62.View ArticleGoogle Scholar
- Wilkinson CJ, Hughes-Thomas ZA, Martin CJ, Bohm I, Mironenko T, Deacon M, et al. Increasing the efficiency of heterologous promoters in actinomycetes. J Mol Microbiol Biotechnol. 2002;4:417–26.Google Scholar
- Bierman M, Logan R, O’Brien K, Seno ET, Rao RN, Schoner BE. Plasmid cloning vectors for the conjugal transfer of DNA from Escherichia coli to Streptomyces spp. Gene. 1992;116:43–9.View ArticleGoogle Scholar
- Wang B, Liu J, Liu H, Huang D, Wen J. Comparative metabolic profiling reveals the key role of amino acids metabolism in the rapamycin overproduction by Streptomyces hygroscopicus. J Ind Microbiol Biotechnol. 2015;42:949–63.View ArticleGoogle Scholar
- Bradford MM. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem. 1976;72:248–54.View ArticleGoogle Scholar
- Youngleson JS, Jones DT, Woods DR. Homology between hydroxybutyryl and hydroxyacyl coenzyme A dehydrogenase enzymes from Clostridium acetobutylicum fermentation and vertebrate fatty acid beta-oxidation pathways. J Bacteriol. 1989;171:6800–7.View ArticleGoogle Scholar
- Bond-Watts BB, Bellerose RJ, Chang MC. Enzyme mechanism as a kinetic control element for designing synthetic biofuel pathways. Nat Chem Biol. 2011;7:222–7.View ArticleGoogle Scholar
- Erb TJ, Berg IA, Brecht V, Müller M, Fuchs G, Alber BE. Synthesis of C5-dicarboxylic acids from C2-units involving crotonyl-CoA carboxylase/reductase: the ethylmalonyl-CoA pathway. Proc Natl Acad Sci USA. 2007;104:10631–6.View ArticleGoogle Scholar
- DuBois M, Gilles KA, Hamilton JK, Rebers P, Smith F. Colorimetric method for determination of sugars and related substances. Anal Chem. 1956;28:350–6.View ArticleGoogle Scholar
- Erb TJ, Rétey J, Fuchs G, Alber BE. Ethylmalonyl-CoA mutase from Rhodobacter sphaeroides defines a new subclade of coenzyme B12-dependent acyl-CoA mutases. J Biol Chem. 2008;283:32283–93.View ArticleGoogle Scholar
- Kupke T, Hernandez-Acosta P, Steinbacher S, Culianez-Macia FA. Arabidopsis thaliana flavoprotein AtHAL3a catalyzes the decarboxylation of 4′-phosphopantothenoylcysteine to 4′-phosphopantetheine, a key step in coenzyme A biosynthesis. J Biol Chem. 2001;276:19190–6.View ArticleGoogle Scholar
- Johnson MNR, Londergan CH, Charkoudian LK. Probing the phosphopantetheine arm conformations of acyl carrier proteins using vibrational spectroscopy. J Am Chem Soc. 2014;136:11240–3.View ArticleGoogle Scholar
- Volk K, Breunig SD, Rid R, Herzog J, Braeuer M, Hundsberger H, et al. Structural analysis and interaction studies of acyl-carrier protein (acpP) of Staphylococcus aureus, an extraordinarily thermally stable protein. Biol Chem. 2017;398:125–33.View ArticleGoogle Scholar
- Jung WS, Kim E, Yoo YJ, Ban YH, Kim EJ, Yoon YJ. Characterization and engineering of the ethylmalonyl-CoA pathway towards the improved heterologous production of polyketides in Streptomyces venezuelae. Appl Microbiol Biotechnol. 2014;98:3701–13.View ArticleGoogle Scholar
- Mo S, Ban Y-H, Park JW, Yoo YJ, Yoon YJ. Enhanced FK506 production in Streptomyces clavuligerus CKD1119 by engineering the supply of methylmalonyl-CoA precursor. J Ind Microbiol Biotechnol. 2009;36:1473–82.View ArticleGoogle Scholar
- Huang D, Xia M, Li S, Wen J, Jia X. Enhancement of FK506 production by engineering secondary pathways of Streptomyces tsukubaensis and exogenous feeding strategies. J Ind Microbiol Biotechnol. 2013;40:1023–37.View ArticleGoogle Scholar
- Chen D, Zhang Q, Zhang Q, Cen P, Xu Z, Liu W. Improvement of FK506 production in Streptomyces tsukubaensis by genetic enhancement of the supply of unusual polyketide extender units via utilization of two distinct site-specific recombination systems. Appl Environ Microbiol. 2012;78:5093–103.View ArticleGoogle Scholar
- Lu C, Zhang X, Jiang M, Bai L. Enhanced salinomycin production by adjusting the supply of polyketide extender units in Streptomyces albus. Metab Eng. 2016;35:129–37.View ArticleGoogle Scholar
- D’Huys P-J, Lule I, Vercammen D, Anne J, Van Impe JF, Bernaerts K. Genome-scale metabolic flux analysis of Streptomyces lividans growing on a complex medium. J Biotechnol. 2012;161:1–13.View ArticleGoogle Scholar
- Wu K, Chung L, Revill WP, Katz L, Reeves CD. The FK520 gene cluster of Streptomyces hygroscopicus var. ascomyceticus (ATCC 14891) contains genes for biosynthesis of unusual polyketide extender units. Gene. 2000;251:81–90.View ArticleGoogle Scholar
- Li CX, Florova G, Akopiants K, Reynolds KA. Crotonyl-coenzyme A reductase provides methylmalonyl-CoA precursors for monensin biosynthesis by Streptomyces cinnamonensis in an oil-based extended fermentation. Microbiology. 2004;150:3463–72.View ArticleGoogle Scholar
- Wilson MC, Moore BS. Beyond ethylmalonyl-CoA: the functional role of crotonyl-CoA carboxylase/reductase homologs in expanding polyketide diversity. Nat Prod Rep. 2012;29:72–86.View ArticleGoogle Scholar