SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
This paper highlights how recent advances in computational identification of biosynthetic gene clusters (BGCs), and in their experimental manipulation, are opening new avenues to access novel secondary metabolites.
Diffusible small molecule microbial hormones drastically alter the expression profiles of antibiotics and other drugs in actinobacteria. For example, avenolide (a butenolide) regulates the production of avermectin, derivatives of which are used in the treatment of river blindness and other parasitic diseases. Butenolides and γ-butyrolactones control the production of pharmaceutically important secondary metabolites by binding to TetR family transcriptional repressors. Here, we describe a concise, 22-step synthetic strategy for the production of avenolide. We present crystal structures of the butenolide receptor AvaR1 in isolation and in complex with avenolide, as well as those of AvaR1 bound to an oligonucleotide derived from its operator. Biochemical studies guided by the co-crystal structures enable the identification of 90 new actinobacteria that may be regulated by butenolides, two of which are experimentally verified. These studies provide a foundation for understanding the regulation of microbial secondary metabolite production, which may be exploited for the discovery and production of novel medicines.
Highly pathogenic human coronavirus infection can cause a severe atypical, rapid onset pneumonia with a mortality rate of up to 10% for severe acute respiratory syndrome coronavirus 1 (SARS-CoV 1), 35% for Middle East respiratory syndrome coronavirus (MERS-CoV), or 1% for severe acute respiratory syndrome coronavirus 2 (SARS-CoV 2 causative agent of COVID 19). While medical countermeasures successfully controlled the worldwide SARS-CoV epidemic, the MERS-CoV epidemic is still ongoing and continues to be a concern for travelers in the Middle East and the multi-year SARS-CoV 2 pandemic underscore the importance of defining biomarkers that are diagnostic and/or predictive of severe disease outcomes for respiratory viruses. Systems biology approaches provide global snapshots of infection induced changes in host cells/tissues and provide extremely rich datasets for understanding host pathogen interactions. Metabolites, especially lipids, are critical for viral replication but less is understood about infection induced changes to lipids due to limits in lipid species detection and identification. To characterize how individual lipid species and proteins with lipid associated functions contribute to highly pathogenic human coronavirus replication and disease severity, existing datasets were probed to determine cell type specific lipid responses to MERS-CoV infection and verification studies were performed to determine if modification of the host lipid signature (by inhibiting enzymatic functions that produce specific lipid species) can perturb CoV replication in human lung cells. MERS-CoV infects human lung epithelial, endothelial and fibroblast cells. All three cell types were infected with MERS-CoV and samples collected to analyze lipids, proteins, metabolites, and transcripts from 12 to 48 hours post infection. Time matched mock-infected cells were collected in parallel for each cell type. Following sample and statistical analysis, functional enrichment and bioinformatic analysis was performed to identify differentially expressed lipids and proteins. Two lipid species were found to be significantly upregulated following MERS-CoV infection, ceramides, and triacylglycerol both of which are indicative of cells undergoing apoptotic cell death. In contrast, sphingomyelins (lipid molecules that can serve as a precursor for one pathway for ceramide synthesis) had significantly decreased expression in MERS-CoV infected cells. An inhibitor of acid sphingomyelinase (that converts sphingomyelin to ceramides and phosphorylcholine) reduced MERS-CoV replication suggesting that production of ceramides is key for successful viral replication and transmission. Acyl-CoA-synthetase 3 (ACSL3), the only protein whose function is lipid associated and had increased differential expression in our dataset, regulates the synthesis of triacylglycerol (increased expression). Inhibitors that directly block ACSL3 (Triacsin C) but not steps later in the triacylglycerol synthesis pathway (Etomoxir) inhibit MERS-CoV replication suggesting that triacylglycerol production and/or ACSL3 activity is also key for viral replication. ACSL3 expression is also upregulated in lung cancer cells and is predicted to promote continued cell viability which would also enhance viral replication. The differentially expressed lipid and lipid-associated protein species identified in our studies suggest that MERS-CoV infection is activating cellular death pathways to limit the number of cells that become infected but also stimulating the production of lipid-associated enzymes that can prolong host cell viability and the amount of time progeny virions can be produced and released. As the inhibitors that worked against MERS-CoV infection were also efficacious against SARS-CoV 2 infection, countermeasures that target these host pathways may provide novel ways to block highly pathogenic human coronavirus infection and/or prevent severe disease outcomes.
Aliphatic carboxylic acids, aldehydes, and ketones play diverse roles in microbial adaptation to their microenvironment, from excretion as toxins to adaptive metabolites for membrane fluidity. However, the spatial distribution of these molecules throughout biofilms and how microbes in these environments exchange these molecules remain elusive for many of these bioactive species due to inefficient molecular imaging strategies. Herein, we apply on-tissue chemical derivatization (OTCD) using 4-(2-((4-bromophenethyl)dimethylammonio)ethoxy)benzenaminium dibromide (4-APEBA) on a co-culture of a soil bacterium (Bacillus subtilis NCIB 3610) and fungus (Fusarium sp. DS 682) grown on agar as our model system. Using matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI), we spatially resolved more than 300 different metabolites containing carbonyl groups within this model system. Various spatial patterns are observable in these species, which indicate possible extracellular or intercellular processes of the metabolites and their up- or downregulation during microbial interaction. The unique chemistry of our approach allowed us to bring additional confidence in accurate carbonyl identification, especially when multiple isomeric candidates were possible, and this provided the ability to generate hypotheses about the potential role of some aliphatic carbonyls in this B. subtilis/Fusarium sp. interaction. The results shown here demonstrate the utility of 4-ABEBA-based OTCD MALDI-MSI in probing interkingdom interactions directly from microbial co-cultures, and these methods will enable future microbial interaction studies with expanded metabolic coverage.
Mass spectrometry (MS)-based, integrated proteomics, metabolomics, and lipidomics (collectively, multi-omics) studies provide a very detailed snap shot of virus-induced changes to the host following infection and can lead to the identification of novel prophylactic and therapeutic targets for preventing or lessening disease severity. Multi-omics studies with Middle East respiratory syndrome coronavirus (MERS-CoV) are challenging as the requirements of Biosafety Level 3 containment limit the numbers of samples that can be safely managed. To address these issues, the multi-omics sample preparation technique MPLEx (metabolite, protein, and lipid extraction) was developed to partition a single sample into three distinct parts (metabolites, proteins and lipids) for multi-omics analysis, while simultaneously inactivating MERS-CoV by solubilizing and disrupting the viral envelope and denaturing viral proteins. Here we describe the MPLEx protocol, highlight the step of inactivation, and describe the details of downstream processing, instrumental analysis of studies with Middle East respiratory syndrome coronavirus (MERS-CoV) are challenging as the requirements of Biosafety Level 3 containment limit the numbers of samples that can be safely managed. To address these issues, the multi-omics sample preparation technique MPLEx (metabolite, protein, and lipid extraction) was developed to partition a single sample into three distinct parts (metabolites, proteins and lipids) for multi-omics analysis, while simultaneously inactivating MERS-CoV by solubilizing and disrupting the viral envelope and denaturing viral proteins. Here we describe the MPLEx protocol, highlight the step of inactivation, and describe the details of downstream processing, instrumental analysis of the three separate analytes and their subsequent informatics pipelines.
Distinguishing the most important features (e.g. proteins, metabolites, etc.) per group (e.g. control and treatment) is a critical challenge in feature-rich multi-omics experiments, especially in soil data. Traditional feature identification and ranking approaches, such as differential expression, are based on single omics and thus not directly translatable to multi-omics experiments. Here, 5 multi-omics integration models (DIABLO, JACA, MOFA, MultiMLP, and SLIDE) that were not explicitly built for soil data applications were tested using a soil-based multi-omics experiment. The data were obtained from an experimental setup of an autoclaved soil system inoculated with 8 bacteria and using chitin as the carbon source and including samples collected at 0- (control), 4-, 8-, and 12-weeks post-inoculation. The omics data included metaproteomics, 16S rRNA sequencing, and LC-MS/MS metabolomics (in positive and negative mode). Each multi-omics integration model was implemented, and top features were compared to differential univariate statistics per omic type, demonstrating that integration approaches cut the potential number of top features from 2957 identified by differential statistics to 13-224 (a 99.6% to 92.4% reduction). Interestingly, most top features across integration models were not shared; though, scaling and averaging ranks across models shared similar patterns. This work highlights the usefulness of multi-omics integration models in soil-based microbial studies and the power of using multiple integration models together to interpret results.
Metabolomics provides a unique snapshot into the world of small molecules and the complex biological processes that govern the human, animal, plant, and environmental ecosystems encapsulated by the One Health modeling framework. However, this “molecular snapshot” is only as informative as the number of metabolites confidently identified within it. The spectral similarity (SS) score is traditionally used to identify compound(s) in mass spectrometry approaches to metabolomics, where spectra are matched to reference libraries of candidate spectra. Unfortunately, there is little consensus on which of the dozens of available SS metrics should be used. This lack of standard SS score creates analytic uncertainty and potentially leads to issues in reproducibility, especially as these data are integrated across other domains. In this work, we use metabolomic spectral similarity as a case study to showcase the challenges in consistency within just one piece of the One Health framework that must be addressed to enable data science approaches for One Health problems. Here, using a large cohort of datasets comprising both standard and complex datasets with expert-verified truth annotations, we evaluated the effectiveness of 66 similarity metrics to delineate between correct matches (true positives) and incorrect matches (true negatives). We additionally characterize the families of these metrics to make informed recommendations for their use. Our results indicate that specific families of metrics (the Inner Product, Correlative, and Intersection families of scores) tend to perform better than others, with no single similarity metric performing optimally for all queried spectra. This work and its findings provide an empirically-based resource for researchers to use in their selection of similarity metrics for GC-MS identification, increasing scientific reproducibility through taking steps towards standardizing identification workflows.
Lignocellulosic conversion residue (LCR) is the material remaining after deconstructed lignocellulosic biomass is subjected to microbial fermentation and treated to remove the biofuel. Technoeconomic analyses of biofuel refineries have shown that further microbial processing of this LCR into other bioproducts may help offset the costs of biofuel generation. Identifying organisms able to metabolize LCR is an important first step for harnessing the full chemical and economic potential of this material. In this study, we investigated the aerobic LCR utilization capabilities of 71 Streptomyces and 163 yeast species that could be engineered to produce valuable bioproducts. The LCR utilization by these individual microbes was compared to that of an aerobic mixed microbial consortium derived from a wastewater treatment plant as representative of a consortium with the highest potential for degrading the LCR components and a source of genetic material for future engineering efforts. We analyzed several batches of a model LCR by chemical oxygen demand (COD) and chromatography-based assays and determined that the major components of LCR were oligomeric and monomeric sugars and other organic compounds. Many of the Streptomyces and yeast species tested were able to grow in LCR, with some individual microbes capable of utilizing over 40% of the soluble COD. For comparison, the maximum total soluble COD utilized by the mixed microbial consortium was about 70%. This represents an upper limit on how much of the LCR could be valorized by engineered Streptomyces or yeasts into bioproducts. To investigate the utilization of specific components in LCR and have a defined media for future experiments, we developed a synthetic conversion residue (SynCR) to mimic our model LCR and used it to show lignocellulose-derived inhibitors (LDIs) had little effect on the ability of the Streptomyces species to metabolize SynCR. We found that LCR is rich in carbon sources for microbial utilization and has vitamins, minerals, amino acids and other trace metabolites necessary to support growth. Testing diverse collections of Streptomyces and yeast species confirmed that these microorganisms were capable of growth on LCR and revealed a phylogenetic correlation between those able to best utilize LCR. Identification and quantification of the components of LCR enabled us to develop a synthetic LCR (SynCR) that will be a useful tool for examining how individual components of LCR contribute to microbial growth and as a substrate for future engineering efforts to use these microorganisms to generate valuable bioproducts.
Large-scale kinetic models provide the computational means to dynamically link metabolic reaction fluxes to metabolite concentrations and enzyme levels while also conforming to substrate level regulation. However, the development of broadly applicable frameworks for efficiently and robustly parameterizing models remains a challenge. Challenges arise due to both the heterogeneity, paucity, and difficulty in obtaining flux and/or concentration data but also due to the computational difficulties of the underlying parameter identification problem. Both the computational demands for parameterization, degeneracy of obtained parameter solutions and interpretability of results has so far limited widespread adoption of large-scale kinetic models despite their potential. Herein, we introduce the Kinetic Estimation Tool Capturing Heterogeneous Datasets Using Pyomo (KETCHUP), a flexible parameter estimation tool that leverages a primal-dual interior-point algorithm to solve a nonlinear programming (NLP) problem that identifies a set of parameters capable of recapitulating the (non)steady-state fluxes and concentrations in wild-type and perturbed metabolic networks. KETCHUP is benchmarked against previously parameterized large-scale kinetic models demonstrating an at least an order of magnitude faster convergence than the tool K-FIT while at the same time attaining better data fits. Finally, this versatile toolbox accepts different kinetic descriptions, metabolic fluxes, enzyme levels and metabolite concentrations, under either steady-state or instationary conditions to enable robust kinetic model construction and parameterization. KETCHUP supports the SBML format and can be accessed at https://github.com/maranasgroup/KETCHUP.
Even as the field of microbiome research has made huge strides in mapping microbial community composition in a variety of environments and organisms, explaining the phenotypic influences on the host by microbial taxa – both known and unknown – and their specific functions still remain major challenges. A pressing need is the ability to assign specific functions in terms of enzymes and small molecules to specific taxa or groups of taxa in the community. This knowledge will be crucial for advancing personalized therapies based on the targeted modulation of microbes or metabolites that have predictable outcomes to benefit the human host. This perspective article advocates for the combined use of standards-free metabolomics and activity-based protein profiling strategies to address this gap in functional knowledge in microbiome research via the identification of novel biomolecules and the attribution of their production to specific microbial taxa.
Current approaches to analyzing metabolomic data often rely on matching MS/MS fragmentation data to sparse libraries or databases. This approach results in limited identification of features, often with less than 10% of the dataset being annotated. A complementary approach is to assign molecular formula to features based on accurate mass measurements, but the platforms commonly used for metabolomics do not have the needed accuracy or resolving power to do this robustly, particularly for larger molecules. Using our newly modified analysis tool, CoreMS, we generated a library of molecular formula from pooled samples analyzed with LC-21T FT-ICR MS. This library successfully annotated approximately 53.2% of features identified from the exometabolome of marine diatom Phaeodactylum tricornutum – a nearly ten-fold increase over the 5.9% annotation rate achieved using a conventional MS/MS library matching approach. Using this FT-ICR MS library approach, we were able to differentiate differences in the exometabolome of P. tricornutum in iron replete and iron limited conditions, with 668 metabolites being differentially expressed (p < 0.05, 2 x intensity difference) under these conditions. The traditional MS/MS fragmentation-based annotation approach only annotated 61 of these metabolites, while our novel pipeline annotated 450 metabolites and revealed 12 metabolites that were significantly more abundant under low iron conditions. Our results demonstrate the utility of ultrahigh resolution mass spectrometry for generating more comprehensive and confident molecular annotations.
Sesquiterpene cyclases (STC) catalyse the cyclization of the C15 molecule farnesyl diphosphate into a vast variety of mono- or polycyclic hydrocarbons and, for a few enzymes, oxygenated structures, with diverse stereogenic centres. The huge diversity in sesquiterpene skeleton structures in nature is primarily the result of the type of cyclization driven by the STC. Despite the phenomenal impact of fungal sesquiterpenes on the ecology of fungi and their potentials for applications, the fungal sesquiterpenome is largely untapped. The identification of fungal STC is generally based on protein sequence similarity with characterized enzymes. This approach has improved our knowledge on STC in a few fungal species, but it has limited success for the discovery of distant sequences. Besides, the tools based on secondary metabolite biosynthesis gene clusters have shown poor performance for terpene cyclases. Here, we used four sets of sequences of fungal STC that catalyse four types of cyclization, and specific amino acid motives to identify phylogenetically related sequences in the genomes of basidiomycetes fungi from the order Polyporales. We validated that four STC genes newly identified from the genome sequence of Leiotrametes menziesii, each classified in a different phylogenetic clade, catalysed a predicted cyclization of farnesyl diphosphate. We built HMM models and searched STC genes in 656 fungal genomes genomes. We identified 5605 STC genes, which were classified in one of the four clades and had a predicted cyclization mechanism. We noticed that the HMM models were more accurate for the prediction of the type of cyclization catalysed by basidiomycete STC than for ascomycete STC.
Complex interactions exist between microbiomes and their hosts. Increasingly, defensive metabolites that have been attributed to host biosynthetic capability are now being recognized as products of host-associated microbes. These unique metabolites often have bioactivity targets in human disease and can be purposed as pharmaceuticals. Polyketides are a complex family of natural products that often serve as defensive metabolites for competitive or pro-survival purposes for the producing organism, while demonstrating bioactivity in human diseases as cholesterol lowering agents, anti-infectives, and anti-tumor agents. Marine invertebrates and microbes are a rich source of polyketides. Palmerolide A, a polyketide isolated from the Antarctic ascidian Synoicum adareanum, is a vacuolar-ATPase inhibitor with potent bioactivity against melanoma cell lines. The biosynthetic gene clusters (BGCs) responsible for production of secondary metabolites are encoded in the genomes of the producers as discrete genomic elements. A candidate palmerolide BGC was identified from a S. adareanum microbiome-metagenome based on a high degree of congruence with a chemical structure-based retrobiosynthetic prediction. Protein family homology analysis, conserved domain searches, active site and motif identification were used to identify and propose the function of the ~75 kbp trans-acyltransferase (AT) polyketide synthase-non-ribosomal synthase (PKS-NRPS) domains responsible for the stepwise synthesis of palmerolide A. Though PKS systems often act in a predictable co-linear sequence, this BGC includes multiple trans-acting enzymatic domains, a non-canonical condensation termination domain, a bacterial luciferase-like monooxygenase (LLM), and is found in multiple copies within the metagenome-assembled genome (MAG). Detailed inspection of the five highly similar pal BGC copies suggests the potential for biosynthesis of other members of the palmerolide chemical family. This is the first delineation of a biosynthetic gene cluster from an Antarctic microbial species, recently proposed as Candidatus Synoicihabitans palmerolidicus. These findings have relevance for fundamental knowledge of PKS combinatorial biosynthesis and could enhance drug development efforts of palmerolide A through heterologous gene expression.
Ochratoxin A (OTA) is a well-known mycotoxin with wide distribution in food and feed. Fungal genome sequencing has great utility for identifying secondary metabolites gene clusters for known and novel compounds. A comparative analysis of the OTA-biosynthetic cluster in A. steynii, A. westerdijkiae, A. niger, A. carbonarius, and P. nordicum has revealed a high synteny in OTA cluster organization in five structural genes (otaA, otaB, ota, otaR1, and otaD). Moreover, a recent detailed comparative genome analysis of Aspergilli OTA producers led to the identification of a cyclase gene, otaY, located in the OTA cluster between the otaA and otaB genes, encoding for a predicted protein with high similarity to SnoaLs domain. These proteins have been shown to catalyze ring closure steps in the biosynthesis of polyketide antibiotics produced in Streptomyces. In the present study, we demonstrated an upregulation of the cyclase gene in A. carbonarius under OTA permissive conditions, consistent with the expression trends of the other OTA cluster genes and their role in OTA biosynthesis by complete gene deletion. Our results pointed out the involvement of a cyclase gene in OTA biosynthetic pathway for the first time. They represent a step forward in the understanding of the molecular basis of OTA biosynthesis in A. carbonarius.
Cognitive computing is revolutionizing the way big data are processed and integrated, with artificial intelligence (AI) natural language processing (NLP) platforms helping researchers to efficiently search and digest the vast scientific literature. Most available platforms have been developed for biomedical researchers, but new NLP tools are emerging for biologists in other fields and an important example is metabolomics. NLP provides literature-based contextualization of metabolic features that decreases the time and expert-level subject knowledge required during the prioritization, identification and interpretation steps in the metabolomics data analysis pipeline. Here, we describe and demonstrate four workflows that combine metabolomics data with NLP-based literature searches of scientific databases to aid in the analysis of metabolomics data and their biological interpretation. Additionally, the four procedures can be used in isolation or consecutively, depending on the research questions. The first, used for initial metabolite annotation and prioritization, creates a list of metabolites that would be interesting for follow-up. The second workflow finds literature evidence of the activity of metabolites and metabolic pathways in governing the biological condition on a systems biology level. The third is used to identify candidate biomarkers, and the fourth looks for metabolic conditions or drug-repurposing targets that the two diseases have in common. The protocol can take 1–4 h or more to complete, depending on the processing time of the various software used.
This feasibility study project aimed to lay further groundwork for studying and understanding how microorganisms interact with each other at the molecular level using model algal-bacterial co-cultures. A better understanding of such interactions is critical to optimizing and operating biotic platforms for a number of applications, including bioenergy, algal bioproducts and agriculture, and carbon capture. Here we explored the feasibility of characterizing and following these interactions by identifying and tracking volatile metabolites and biomarkers in the complex gaseous headspace of such cultures. We set up an experimental platform using a co-culture of Phaeodactylum tricornutum and Marinobacter subspecies 3-2 as a model system, used solid-phase microextraction fibers for volatile collection at various time points throughout culture growth, and employed gas chromatography mass spectrometry-based instrumentation available at LLNL for compound identification and quantification. We found that it is feasible to determine volatile profiles from the growth of Phaeodactylum tricornutum and Marinobacter subspecies 3-2 using our methodology and monitor profile changes over time. We also observed that volatiles distinguishing between cultures of the individual species and the co-cultures could be detected. We expect that results of this project will feed into the ongoing LLNL’s DOE BER funded Biofuels Scientific Focus Area (SFA) on algal-bacterial interactions and pave the way for future expanded volatile studies supported by DOE. The approach and technologies developed here would also be applicable and transferable to other microbial communities under study, such as cyanobacteria, rhizosphere communities and biofilms.
Phosphorus (P) deficiency in plants creates a variety of metabolic perturbations that decrease photosynthesis and growth. Phosphorus deficiency is especially challenging for the production of bioenergy feedstock plantation species, such as poplars (Populus spp.), where fertilization may not be practically or economically feasible. While the phenotypic effects of P deficiency are well known, the molecular mechanisms underlying whole-plant and tissue-specific responses to P deficiency, and in particular the responses of commercially valuable hardwoods, are less studied. We used a multi-tissue and multi-omics approach using transcriptomic, proteomic, and metabolomic analyses of the leaves and roots of black cottonwood (Populus trichocarpa) seedlings grown under P-deficient (5 µM P) and replete (100 µM P) conditions to assess this knowledge gap and to identify potential gene targets for selection for P efficiency. In comparison to seedlings grown at 100 µM P, P-deficient seedlings exhibited reduced dry biomass, altered chlorophyll fluorescence, and reduced tissue P concentrations. In line with these observations, growth, C metabolism, and photosynthesis pathways were downregulated in the transcriptome of the P-deficient plants. Additionally, we found evidence of strong lipid remodeling in the leaves. Metabolomic data showed that the roots of P-deficient plants had a greater relative abundance of phosphate ion, which may reflect extensive degradation of P-rich metabolites in plants exposed to long-term P-deficiency. With the notable exception of the KEGG pathway for Starch and Sucrose Metabolism (map00500), the responses of the transcriptome and the metabolome to P deficiency were consistent with one another. No significant changes in the proteome were detected in response to P deficiency. Collectively, our multi-omic and multi-tissue approach enabled the identification of important metabolic and regulatory pathways regulated across tissues at the molecular level that will be important avenues to further evaluate for P efficiency. These included stress-mediating systems associated with reactive oxygen species maintenance, lipid remodeling within tissues, and systems involved in P scavenging from the rhizosphere.
Leveraging predicted origin sequences from a previously characterized groundwater plasmidome, we constructed a barcoded plasmid library to screen for previously unknown origins. Testing this library against a panel of representative bacterial strains led to the identification of 3 previously unknown origins that replicate in gram-negative bacteria not previously associated with these origin sequences. Experimental validation confirmed that a plasmid bearing origin 6911 as the sole origin could replicate with a copy number of 9 (±2) in Pantoea sp. MT58, a fast growing and metal tolerant, environmentally important bacterium. Plasmids based on this new origin were used to express the reporter protein GFP, and non-native metabolite pathways for the natural product indigoidine and the terpenoid compound isoprenol. Functional previously unknown origins of replication in such non-model organisms can expand the toolkit for genetic manipulations of both model and less-studied bacteria.