The integration of mathematical modeling and experimental testing is emerging as a powerful approach for improving our understanding of the regulation of metabolic pathways. In this study, we report on the development of a kinetic mathematical model that accurately simulates the developmental patterns of monoterpenoid essential oil accumulation in peppermint (Mentha×piperita). This model was then used to evaluate the biochemical processes underlying experimentally determined changes in the monoterpene pathway under low ambient-light intensities, which led to an accumulation of the branchpoint intermediate (+)-pulegone and the side product (+)-menthofuran. Our simulations indicated that the environmentally regulated changes in monoterpene profiles could only be explained when, in addition to effects on biosynthetic enzyme activities, as yet unidentified inhibitory effects of (+)-menthofuran on the branchpoint enzyme pulegone reductase (PR) were assumed. Subsequentin vitroanalyses with recombinant protein confirmed that (+)-menthofuran acts as a weak competitive inhibitor of PR (K_i= 300 μM). To evaluate whether the intracellular concentration of (+)-menthofuran was high enough for PR inhibitionin vivo, we isolated essential oil-synthesizing secretory cells from peppermint leaves and subjected them to steam distillations. When peppermint plants were grown under low-light conditions, (+)-menthofuran was selectively retained in secretory cells and accumulated to very high levels (up to 20 mM), whereas under regular growth conditions, (+)-menthofuran levels remained very low (<400 μM). These results illustrate the utility of iterative cycles of mathematical modeling and experimental testing to elucidate the mechanisms controlling flux through metabolic pathways.

We have previously reported the use of a combination of computational simulations and targeted experiments to build a first generation mathematical model of peppermint (Mentha × piperita) essential oil biosynthesis. Here, we report on the expansion of this approach to identify the key factors controlling monoterpenoid essential oil biosynthesis under adverse environmental conditions. We also investigated determinants of essential oil biosynthesis in transgenic peppermint lines with modulated essential oil profiles. A computational perturbation analysis, which was implemented to identify the variables that exert prominent control over the outputs of the model, indicated that the essential oil composition should be highly dependent on certain biosynthetic enzyme concentrations [(+)-pulegone reductase and (+)-menthofuran synthase], whereas oil yield should be particularly sensitive to the density and/or distribution of leaf glandular trichomes, the specialized anatomical structures responsible for the synthesis and storage of essential oils. A microscopic evaluation of leaf surfaces demonstrated that the final mature size of glandular trichomes was the same across all experiments. However, as predicted by the perturbation analysis, differences in the size distribution and the total number of glandular trichomes strongly correlated with differences in monoterpenoid essential oil yield. Building on various experimental data sets, appropriate mathematical functions were selected to approximate the dynamics of glandular trichome distribution/density and enzyme concentrations in our kinetic model. Based on a χ² statistical analysis, simulated and measured essential oil profiles were in very good agreement, indicating that modeling is a valuable tool for guiding metabolic engineering efforts aimed at improving essential oil quality and quantity.

Species of the genus Streptomyces are known for their ability to produce multiple secondary metabolites; their genomes have been extensively explored to discover new bioactive compounds. The richness of genomic data currently available allows filtering for high quality genomes, which in turn permits reliable comparative genomics studies and an improved prediction of biosynthetic gene clusters (BGCs) through genome mining approaches. In this work, we used 121 genome sequences of the genus Streptomyces in a comparative genomics study with the aim of estimating the genomic diversity by protein domains content, sequence similarity of proteins and conservation of Intergenic Regions (IGRs). We also searched for BGCs but prioritizing those with potential antibiotic activity. Our analysis revealed that the pan-genome of the genus Streptomyces is clearly open, with a high quantity of unique gene families across the different species and that the IGRs are rarely conserved. We also described the phylogenetic relationships of the analyzed genomes using multiple markers, obtaining a trustworthy tree whose relationships were further validated by Average Nucleotide Identity (ANI) calculations. Finally, 33 biosynthetic gene clusters were detected to have potential antibiotic activity and a predicted mode of action, which might serve up as a guide to formulation of related experimental studies.

Streptomyces clavuligerus is a gram-positive filamentous bacterium notable for producing clavulanic acid (CA), an inhibitor of β-lactamase enzymes, which confers resistance to bacteria against several antibiotics. Here we present a comparative analysis of the morphological and metabolic response of S. clavuligerus linked to the CA production under low and high shear stress conditions in a 2D rocking-motion single-use bioreactor (CELL-tainer ^®) and stirred tank bioreactor (STR), respectively. The CELL-tainer^® guarantees high turbulence and enhanced volumetric mass transfer at low shear stress, which (in contrast to bubble columns) allows the investigation of the impact of shear stress without oxygen limitation. The results indicate that high shear forces do not compromise the viability of S. clavuligerus cells; even higher specific growth rate, biomass, and specific CA production rate were observed in the STR. Under low shear forces in the CELL-tainer^® the mycelial diameter increased considerably (average diameter 2.27 in CELL-tainer^® vs. 1.44 µm in STR). This suggests that CA production may be affected by a lower surface-to-volume ratio which would lead to lower diffusion and transport of nutrients, oxygen, and product. The present study shows that there is a strong correlation between macromorphology and CA production, which should be an important aspect to consider in industrial production of CA.