Advancing Anticancer Drug Discovery: Leveraging Metabolomics and Machine Learning for Mode of Action Prediction by Pattern Recognition.

A bottleneck in the development of new anti-cancer drugs is the recognition of their mode of action (MoA). Metabolomics combined with machine learning allowed to predict MoAs of novel anti-proliferative drug candidates, focusing on human prostate cancer cells (PC-3). As proof of concept, 38 drugs are studied with known effects on 16 key processes of cancer metabolism, profiling low molecular weight intermediates of the central carbon and cellular energy metabolism (CCEM) by LC-MS/MS. These metabolic patterns unveiled distinct MoAs, enabling accurate MoA predictions for novel agents by machine learning. The transferability of MoA predictions based on PC-3 cell treatments is validated with two other cancer cell models, i.e., breast cancer and Ewing's sarcoma, and show that correct MoA predictions for alternative cancer cells are possible, but still at some expense of prediction quality. Furthermore, metabolic profiles of treated cells yield insights into intracellular processes, exemplified for drugs inducing different types of mitochondrial dysfunction. Specifically, it is predicted that pentacyclic triterpenes inhibit oxidative phosphorylation and affect phospholipid biosynthesis, as confirmed by respiration parameters, lipidomics, and molecular docking. Using biochemical insights from individual drug treatments, this approach offers new opportunities, including the optimization of combinatorial drug applications.

Hordedane diterpenoid phytoalexins restrict Fusarium graminearum infection but enhance Bipolaris sorokiniana colonization of barley roots.

Plant immunity is a multilayered process that includes recognition of patterns or effectors from pathogens to elicit defense responses. These include the induction of a cocktail of defense metabolites that typically restrict pathogen virulence. Here, we investigate the interaction between barley roots and the fungal pathogens Bipolaris sorokiniana (Bs) and Fusarium graminearum (Fg) at the metabolite level. We identify hordedanes, a previously undescribed set of labdane-related diterpenoids with antimicrobial properties, as critical players in these interactions. Infection of barley roots by Bs and Fg elicits hordedane synthesis from a 600-kb gene cluster. Heterologous reconstruction of the biosynthesis pathway in yeast and Nicotiana benthamiana produced several hordedanes, including one of the most functionally decorated products 19-ß-hydroxy-hordetrienoic acid (19-OH-HTA). Barley mutants in the diterpene synthase genes of this cluster are unable to produce hordedanes but, unexpectedly, show reduced Bs colonization. By contrast, colonization by Fusarium graminearum, another fungal pathogen of barley and wheat, is 4-fold higher in the mutants completely lacking hordedanes. Accordingly, 19-OH-HTA enhances both germination and growth of Bs, whereas it inhibits other pathogenic fungi, including Fg. Analysis of microscopy and transcriptomics data suggest that hordedanes delay the necrotrophic phase of Bs. Taken together, these results show that adapted pathogens such as Bs can subvert plant metabolic defenses to facilitate root colonization.

Combinatorial biosynthesis in yeast leads to over 200 diterpenoids.

Diterpenoids form a diverse group of natural products, many of which are or could become pharmaceuticals or industrial chemicals. The modular character of diterpene biosynthesis and the promiscuity of the enzymes involved make combinatorial biosynthesis a promising approach to generate libraries of diverse diterpenoids. Here, we report on the combinatorial assembly in yeast of ten diterpene synthases producing (+)-copalyl diphosphate-derived backbones and four cytochrome P450 oxygenases (CYPs) in diverse combinations. This resulted in the production of over 200 diterpenoids. Based on literature and chemical database searches, 162 of these compounds can be considered new-to-Nature. The CYPs accepted most substrates they were given but remained regioselective with few exceptions. Our results provide the basis for the systematic exploration of the diterpenoid chemical space in yeast using sequence databases.

Dynamics of Reactive Carbonyl Species in Pea Root Nodules in Response to Polyethylene Glycol (PEG)-Induced Osmotic Stress.

Drought dramatically affects crop productivity worldwide. For legumes this effect is especially pronounced, as their symbiotic association with rhizobia is highly-sensitive to dehydration. This might be attributed to the oxidative stress, which ultimately accompanies plants' response to water deficit. Indeed, enhanced formation of reactive oxygen species in root nodules might result in up-regulation of lipid peroxidation and overproduction of reactive carbonyl compounds (RCCs), which readily modify biomolecules and disrupt cell functions. Thus, the knowledge of the nodule carbonyl metabolome dynamics is critically important for understanding the drought-related losses of nitrogen fixation efficiency and plant productivity. Therefore, here we provide, to the best of our knowledge, for the first time a comprehensive overview of the pea root nodule carbonyl metabolome and address its alterations in response to polyethylene glycol-induced osmotic stress as the first step to examine the changes of RCC patterns in drought treated plants. RCCs were extracted from the nodules and derivatized with 7-(diethylamino)coumarin-3-carbohydrazide (CHH). The relative quantification of CHH-derivatives by liquid chromatography-high resolution mass spectrometry with a post-run correction for derivative stability revealed in total 194 features with intensities above 1 × 105 counts, 19 of which were down- and three were upregulated. The upregulation of glyceraldehyde could accompany non-enzymatic conversion of glyceraldehyde-3-phosphate to methylglyoxal. The accumulation of 4,5-dioxovaleric acid could be the reason for down-regulation of porphyrin metabolism, suppression of leghemoglobin synthesis, inhibition of nitrogenase and degradation of legume-rhizobial symbiosis in response to polyethylene glycol (PEG)-induced osmotic stress effect. This effect needs to be confirmed with soil-based drought models.

Does Protein Glycation Impact on the Drought-Related Changes in Metabolism and Nutritional Properties of Mature Pea (Pisum sativum L.) Seeds?

Protein glycation is usually referred to as an array of non-enzymatic post-translational modifications formed by reducing sugars and carbonyl products of their degradation. The resulting advanced glycation end products (AGEs) represent a heterogeneous group of covalent adducts, known for their pro-inflammatory effects in mammals, and impacting on pathogenesis of metabolic diseases and ageing. In plants, AGEs are the markers of tissue ageing and response to environmental stressors, the most prominent of which is drought. Although water deficit enhances protein glycation in leaves, its effect on seed glycation profiles is still unknown. Moreover, the effect of drought on biological activities of seed protein in mammalian systems is still unstudied with respect to glycation. Therefore, here we address the effects of a short-term drought on the patterns of seed protein-bound AGEs and accompanying alterations in pro-inflammatory properties of seed protein in the context of seed metabolome dynamics. A short-term drought, simulated as polyethylene glycol-induced osmotic stress and applied at the stage of seed filling, resulted in the dramatic suppression of primary seed metabolism, although the secondary metabolome was minimally affected. This was accompanied with significant suppression of NF-kB activation in human SH-SY5Y neuroblastoma cells after a treatment with protein hydrolyzates, isolated from the mature seeds of drought-treated plants. This effect could not be attributed to formation of known AGEs. Most likely, the prospective anti-inflammatory effect of short-term drought is related to antioxidant effect of unknown secondary metabolite protein adducts, or down-regulation of unknown plant-specific AGEs due to suppression of energy metabolism during seed filling.

The inconspicuous gatekeeper: endophytic Serendipita vermifera acts as extended plant protection barrier in the rhizosphere.

In nature, beneficial and pathogenic fungi often simultaneously colonise plants. Despite substantial efforts to understand the composition of natural plant-microbe communities, the mechanisms driving such multipartite interactions remain largely unknown. Here we address how the interaction between the beneficial root endophyte Serendipita vermifera and the pathogen Bipolaris sorokiniana affects fungal behaviour and determines barley host responses using a gnotobiotic soil-based split-root system. Fungal confrontation in soil resulted in induction of B. sorokiniana genes involved in secondary metabolism and a significant repression of genes encoding putative effectors. In S. vermifera, genes encoding hydrolytic enzymes were strongly induced. This antagonistic response was not activated during the tripartite interaction in barley roots. Instead, we observed a specific induction of S. vermifera genes involved in detoxification and redox homeostasis. Pathogen infection but not endophyte colonisation resulted in substantial host transcriptional reprogramming and activation of defence. In the presence of S. vermifera, pathogen infection and disease symptoms were significantly reduced despite no marked alterations of the plant transcriptional response. The activation of stress response genes and concomitant repression of putative effector gene expression in B. sorokiniana during confrontation with the endophyte suggest a reduction of the pathogen's virulence potential before host plant infection.

Strigolactone Levels in Dicot Roots Are Determined by an Ancestral Symbiosis-Regulated Clade of the PHYTOENE SYNTHASE Gene Family.

Strigolactones (SLs) are apocarotenoid phytohormones synthesized from carotenoid precursors. They are produced most abundantly in roots for exudation into the rhizosphere to cope with mineral nutrient starvation through support of root symbionts. Abscisic acid (ABA) is another apocarotenoid phytohormone synthesized in roots, which is involved in responses to abiotic stress. Typically low carotenoid levels in roots raise the issue of precursor supply for the biosynthesis of these two apocarotenoids in this organ. Increased ABA levels upon abiotic stress in Poaceae roots are known to be supported by a particular isoform of phytoene synthase (PSY), catalyzing the rate-limiting step in carotenogenesis. Here we report on novel PSY3 isogenes from Medicago truncatula (MtPSY3) and Solanum lycopersicum (SlPSY3) strongly expressed exclusively upon root interaction with symbiotic arbuscular mycorrhizal (AM) fungi and moderately in response to phosphate starvation. They belong to a widespread clade of conserved PSYs restricted to dicots (dPSY3) distinct from the Poaceae-PSY3s involved in ABA formation. An ancient origin of dPSY3s and a potential co-evolution with the AM symbiosis is discussed in the context of PSY evolution. Knockdown of MtPSY3 in hairy roots of M. truncatula strongly reduced SL and AM-induced C13 α-ionol/C14 mycorradicin apocarotenoids. Inhibition of the reaction subsequent to phytoene synthesis revealed strongly elevated levels of phytoene indicating induced flux through the carotenoid pathway in roots upon mycorrhization. dPSY3 isogenes are coregulated with upstream isogenes and downstream carotenoid cleavage steps toward SLs (D27, CCD7, CCD8) suggesting a combined carotenoid/apocarotenoid pathway, which provides "just in time"-delivery of precursors for apocarotenoid formation.

Multi-Omics of Tomato Glandular Trichomes Reveals Distinct Features of Central Carbon Metabolism Supporting High Productivity of Specialized Metabolites.

Glandular trichomes are metabolic cell factories with the capacity to produce large quantities of secondary metabolites. Little is known about the connection between central carbon metabolism and metabolic productivity for secondary metabolites in glandular trichomes. To address this gap in our knowledge, we performed comparative metabolomics, transcriptomics, proteomics, and 13C-labeling of type VI glandular trichomes and leaves from a cultivated (Solanum lycopersicum LA4024) and a wild (Solanum habrochaites LA1777) tomato accession. Specific features of glandular trichomes that drive the formation of secondary metabolites could be identified. Tomato type VI trichomes are photosynthetic but acquire their carbon essentially from leaf sucrose. The energy and reducing power from photosynthesis are used to support the biosynthesis of secondary metabolites, while the comparatively reduced Calvin-Benson-Bassham cycle activity may be involved in recycling metabolic CO2 Glandular trichomes cope with oxidative stress by producing high levels of polyunsaturated fatty acids, oxylipins, and glutathione. Finally, distinct mechanisms are present in glandular trichomes to increase the supply of precursors for the isoprenoid pathways. Particularly, the citrate-malate shuttle supplies cytosolic acetyl-CoA and plastidic glycolysis and malic enzyme support the formation of plastidic pyruvate. A model is proposed on how glandular trichomes achieve high metabolic productivity.

Early responses of mature Arabidopsis thaliana plants to reduced water potential in the agar-based polyethylene glycol infusion drought model.

Drought is one of the most important environmental stressors resulting in increasing losses of crop plant productivity all over the world. Therefore, development of new approaches to increase the stress tolerance of crop plants is strongly desired. This requires precise and adequate modeling of drought stress. As this type of stress manifests itself as a steady decrease in the substrate water potential (ψw), agar plates infused with polyethylene glycol (PEG) are the perfect experimental tool: they are easy in preparation and provide a constantly reduced ψw, which is not possible in soil models. However, currently, this model is applicable only to seedlings and cannot be used for evaluation of stress responses in mature plants, which are obviously the most appropriate objects for drought tolerance research. To overcome this limitation, here we introduce a PEG-based agar infusion model suitable for 6-8-week-old A. thaliana plants, and characterize, to the best of our knowledge for the first time, the early drought stress responses of adult plants grown on PEG-infused agar. We describe essential alterations in the primary metabolome (sugars and related compounds, amino acids and polyamines) accompanied by qualitative and quantitative changes in protein patterns: up to 87 unique stress-related proteins were annotated under drought stress conditions, whereas further 84 proteins showed a change in abundance. The obtained proteome patterns differed slightly from those reported for seedlings and soil-based models.

Autofluorescence as a Signal to Sort Developing Glandular Trichomes by Flow Cytometry.

The industrial relevance of a number of metabolites produced in plant glandular trichomes (GTs) has spurred research on these specialized organs for a number of years. Most of the research, however, has focused on the elucidation of secondary metabolite pathways and comparatively little has been undertaken on the development and differentiation of GTs. One way to gain insight into these developmental processes is to generate stage-specific transcriptome and metabolome data. The difficulty for this resides in the isolation of early stages of development of the GTs. Here we describe a method for the separation and isolation of intact young and mature type VI trichomes from the wild tomato species Solanum habrochaites. The final and key step of the method uses cell sorting based on distinct autofluorescence signals of the young and mature trichomes. We demonstrate that sorting by flow cytometry allows recovering pure fractions of young and mature trichomes. Furthermore, we show that the sorted trichomes can be used for transcript and metabolite analyses. Because many plant tissues or cells have distinct autofluorescence components, the principles of this method can be generally applicable for the isolation of specific cell types without prior labeling.

Isoprenoid and metabolite profiling of plant trichomes.

Plant glandular trichomes are specialized secretory structures located on the surface of the aerial parts of plants with large biosynthetic capacity, often with terpenoids as output molecules. The collection of plant trichomes requires a method to separate trichomes from leaf epidermal tissues. For metabolite profiling, trichome tissue needs to be rapidly quenched in order to maintain the indigenous state of intracellular intermediates. Appropriate extraction and chromatographic separation methods must be available, which address the wide-ranging polarity of metabolites. In this chapter, a protocol for trichome harvest using a frozen paint brush is presented. A work flow for broad-range metabolite profiling using LC-MS(2) analysis is described, which is applicable to assess very hydrophilic isoprenoid precursors as well as more hydrophobic metabolites from trichomes and other plant tissues.

Nitrogen as an indicator of mass transfer during in-situ gas sparging.

Aiming at the stimulation of intrinsic microbial activity, pulses of pure oxygen or pressurized air were recurrently injected into groundwater polluted with chlorobenzene. To achieve well-controlled conditions and intensive sampling, a large, vertical underground tank was filled with the local unconfined sandy aquifer material. In the course of two individual gas injections, one using pure oxygen and one using pressurized air, the mass transfer of individual gas species between trapped gas phase and groundwater was studied. Field data on the dissolved gas composition in the groundwater were combined with a kinetic model on gas dissolution and transport in porous media. Phase mass transfer of individual gas components caused a temporary enrichment of nitrogen, and to a lower degree of methane, in trapped gas leading to the formation of excess dissolved nitrogen levels downgradient from the dissolving gas phase. By applying a novel gas sampling method for dissolved gases in groundwater it was shown that dissolved nitrogen can be used as a partitioning tracer to indicate complete gas dissolution in porous media.

Interplay between oxygen demand reactions and kinetic gas-water transfer in porous media.

Gas-water phase transfer associated with the dissolution of trapped gas in porous media is a key process that occurs during pulsed gas sparging operations in contaminated aquifers. Recently, we applied a numerical model that was experimentally validated for abiotic situations, where multi-species kinetic inter-phase mass transfer and dissolved gas transport occurred during pulsed gas penetration-dissolution events [Balcke, G.U., Meenken, S., Hoefer, C. and Oswald, S.E., 2007. Kinetic gas-water transfer and gas accumulation in porous media during pulsed oxygen sparging. Environmental Science & Technology 41(12), 4428-4434]. Here we extend the model by using a reactive term to describe dissolved oxygen demand reactions via the formation of a reaction product, and to study the effects of such an aerobic degradation process on gas-water mass transfer and dissolution of trapped gas in porous media. As a surrogate for microbial oxygen reduction, first-order oxygen demand reactions were based on the measured oxidation of alkaline pyrogallol in column experiments. This reaction allows for adjusting the rate to values close to expected biodegradation rates and detection of the reaction product. The experiments and model consistently demonstrated accelerated oxygen gas-water mass transfer with increasing oxygen demand rates associated with an influence on the partitioning of other gases. Thus, as the oxygen demand accelerates, less gas phase residues, consisting mainly of nitrogen, are observed, which is in general beneficial to the performance of field biosparging operations. Model results additionally predict how oxygen demand influences oxygen mass transfer for a range of biodegradation rates. A typical field case scenario was simulated to illustrate the observed coupling of oxygen consumption and gas bubble dissolution. The model provides a tool to improve understanding of trapped gas behavior in porous media and contributes to a model-assisted biosparging.

Functional metaproteome analysis of protein extracts from contaminated soil and groundwater.

Using proteins from soil or groundwater as functional biomarkers requires efficient extraction. We developed an extraction method in which the separation of proteins from the inorganic and organic constituents of the soil matrix was achieved by a combination of 0.1 M NaOH treatment and phenol extraction. Incubation with NaOH released humic acids and proteins from soil minerals, and simultaneously, disrupted microorganisms. The subsequent phenol extraction separated the proteins from the humic organic matter. Protein extracts were applied to sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and 2D-electrophoresis (2-DE). Spots and bands were excised and individual proteins identified by liquid chromatography online linked to mass spectrometry (MS) via electrospray ionization source (LC-ESI-MS). To assess the suitability of the method for the functional analysis of environmental metaproteomes, it was applied to soil that had been enriched in chlorophenoxy acid-degrading bacteria by incubation with 2,4-dichlorophenoxy acetic acid (2,4-D) for 22 days. The method was also used to analyze groundwater from the aquifer of a chlorobenzene-contaminated site. The identification of enzymes such as chlorocatechol dioxygenases was consistent with bacterial metabolic pathways expected to be expressed in these samples. The protocol enabled thus the analysis of the metaproteome of soil and groundwater samples. It thereby provides a means to study the diversity of environmental microbial communities while addressing functional aspects more directly than metagenome or even metatranscriptome analysis.

Kinetic gas-water transfer and gas accumulation in porous media during pulsed oxygen sparging.

Gas-water mass transfer and the transport of dissolved gases in variably saturated porous media are key processes for in-situ remediation by pulsed gas sparging. In this context, gas dissolution tests were conducted during pulsed oxygen gas injection into sand columns. The columns were recharged with anoxic water, effluents were analyzed for dissolved O2, and tracer tests were performed to detect accumulation of trapped gas. In a second series oxygen gas was blended with sulfur hexafluoride (SF6), and O2 and SF6 breakthrough curves were recorded. To interpret experimental results, a numerical model was applied that simulates multi-species kinetic mass transfer during gas dissolution. The model predicted breakthrough curves of dissolved gas species and delivered spatially resolved values for gas phase accumulation and composition, which are not directly accessible experimentally. It was shown how dissolved nitrogen accumulates increasingly in trapped gas phase and inhibits its complete dissolution, in case the pulsed gas injections were operated based on O2 breakthrough only. Accumulation of nitrogen also retarded dissolved oxygen transport and thus oxygen breakthrough. Experiments plus modeling demonstrated that SF6 measurements are highly sensitive to the gas dissolution processes, and provide a more sensitive criterion for determining gas injection frequencies during pulsed biosparging.

Determination of naturally occurring MTBE biodegradation by analysing metabolites and biodegradation by-products.

Methyl tert-butyl ether (MTBE) is one of the main additives in gasoline. Its degradation is known to be difficult in natural environments. In this study, significant MTBE degradation is demonstrated at a contaminated site in Leuna (eastern Germany). Since the extent of the plume appeared to be constant over the last 5 years, an extended study was performed to elucidate the degradation processes. Special attention was paid to the production, accumulation and degradation of metabolites and by-products. Groundwater samples from 105 monitoring wells were used to measure 20 different substances. During the degradation process, several intermediates such as tert-butyl alcohol (TBA), tert-butyl formate, formate and lactate were produced. However, the potentially carcinogenic by-product methacrylate was not detected in several hundred samples. At the Leuna site, MTBE degradation occurred under microaerobic conditions. In contrast to hydrocarbons and BTEX, there was no evidence for anaerobic MTBE degradation. Among the degradation products, TBA was found to be a useful intermediate to identify MTBE degradation, at least under microaerobic conditions. TBA accumulation was strongly correlated to MTBE degradation according to the kinetic properties of both degradation processes. Since maximum degradation rates (v(max)) and k(m) values were higher for MTBE (v(max)=2.3 mg/l/d and k(m)=3.2 mg/l) than for TBA (v(max)=1.35 mg/l/d and k(m)=0.05 mg/l), TBA significantly accumulated as an intermediate by-product. The field results were supported by bench scale model aquifer experiments.

Chlorobenzene biodegradation under consecutive aerobic-anaerobic conditions.

The biodegradation of monochlorobenzene, the main contaminant in a quaternary aquifer at Bitterfeld, Central Germany, was studied in microcosm experiments employing either original groundwater or defined mineral media together with the indigenous microbial community from the polluted site. The impact of consecutive aerobic-anaerobic-aerobic incubations on monochlorobenzene biodegradation, microbial diversity, and pH development was examined. The related changes in microbial community composition were analyzed by 16S rRNA gene-based single-strand conformation polymorphism (SSCP) fingerprints and sequencing of dominant bands and by quantitative analysis of bacterial respiratory chain quinones as biomarkers. Under aerobic conditions, the indigenous microbial community of the groundwater degraded monochlorobenzene mainly via the modified ortho-pathway. Respiratory chain quinones and SSCP analysis suggested dominance of the genera Acidovorax and Pseudomonas. A shift to anoxic conditions resulted in monochlorobenzene biotransformation but no dechlorination. The ability to degrade monochlorobenzene aerobically remained preserved throughout a fortnightly anoxic period at sufficiently high buffer capacity. Acidification, caused by monochlorobenzene biodegradation, was alkalinity-controlled. At low initial alkalinity a substantial decrease in pH, monochlorobenzene degradation, and total counts of live cells, accompanied by a change of the microbial community composition, was observed.