Nearline acquisition and processing of liquid chromatography-tandem mass spectrometry data

Liquid chromatography–mass spectrometry (LC–MS) is a commonly used analytical platform for non-targeted metabolite profiling experiments. Although data acquisition, processing and statistical analyses are almost routine in such experiments, further annotation and subsequent identification of chemical compounds are not. For identification, tandem mass spectra provide valuable information towards the structure of chemical compounds. These are typically acquired online, in data-dependent mode, or offline, using handcrafted acquisition methods and manually extracted from raw data. Here, we present several methods to fast-track and improve both the acquisition and processing of LC–MS/MS data. Our nearly online (nearline) data-dependent tandem MS strategy creates a minimal set of LC–MS/MS acquisition methods for relevant features revealed by a preceding non-targeted profiling experiment. Using different filtering criteria, such as intensity or ion type, the acquisition of irrelevant spectra is minimized. Afterwards, LC–MS/MS raw data are processed with feature detection and grouping algorithms. The extracted tandem mass spectra can be used for both library search and de-novo identification methods. The algorithms are implemented in the R package MetShot and support the export to Bruker, Agilent or Waters QTOF instruments and the vendor-independent TraML standard. We evaluate the performance of our workflow on a Bruker micrOTOF-Q by comparison of automatically acquired and extracted tandem mass spectra obtained from a mixture of natural product standards against manually extracted reference spectra. Using Arabidopsis thaliana wild-type and biosynthetic gene knockout plants, we characterize the metabolic products of a biosynthetic pathway and demonstrate the integration of our approach into a typical non-targeted metabolite profiling workflow.     

Discovery of novel indole-based aroylhydrazones as anticonvulsants: Pharmacophore-based design

A number of novel melatonin derivatives, containing aroylhydrazone moieties, were synthesized and explored in vivo for anticonvulsant activity, neurotoxicity in ICR mice as well as in-vitro for cytoxicity and oxidative stress in rats. The structures and configurations were confirmed by NMR, FTIR, HRMS and crystal X-ray diffraction method. For selection of potent structures for synthesis a pharmacophore model was used. Two compounds 3e, with a 2-furyl moiety fragment and 3f with 2-thienyl fragment, showed a potency in maximal electroshock (MES) test (ED₅₀ = 50.98 mg kg⁻¹, PI > 5.88 and ED₅₀ = 108.7 mg kg⁻¹; PI > 2.76), respectively, higher than melatonin (ED₅₀ = 160.3 mg kg⁻¹, PI > 1.87). The compounds 3c, 3e, 3f and 3i suppressed psychomotor seizures in the 6 Hz test and 3c was the most potent with higher ED₅₀ = 13.98 mg kg⁻¹ and PI of > 21.46 compared to that of melatonin (ED₅₀ = 49.76 mg kg⁻¹ and PI of > 6.03) in mice. None of the compounds displayed neurotoxicity in the rota-rod test. The novel melatonin derivatives exerted weak cytotoxic effects while 3f showed the lowest hepatoxic effects comparable to that of the positive control melatonin in rats. The high affinities to the elucidated pharmacophore model of the novel melatonin compounds derived from the inclusion of aroylhydrazone moiety in the indole scaffold yielded suitable candidates with anticonvulsant activity in the MES and 6 Hz test of psychomotor seizures.     

Using MS-FINDER for identifying 19 natural products in the CASMI 2016 contest

In its fourth year, the CASMI 2016 contest was organized to evaluate current chemical structure identification strategies for 19 natural products using high-resolution LC–MS and LC–MS/MS challenge datasets using automated methods with or without the combination of other tools. These natural products originate from plants, fungi, marine sponges, algae, or micro-algae. Every compound annotation workflow must start with determination of elemental compositions. Of these 19 challenges, one was excluded by the organizers after submission. For the remaining 18 challenges, three software programs were used. MS-FINDER version 1.62 was able to correctly identify 89% of the molecular formulas using an internal database that comprised of 13 metabolomics repositories with 45,181 formulas. SIRIUS correctly identified 61% compositions using PubChem formulas and Seven Golden Rules correctly identified 83% by using the Dictionary of Natural Products as a targeted database. Next, we performed structural dereplication for which we used the consensus formula from the three software programs. We submitted two solution sets for these challenges. In the first solution set, avaniya001, we only used the internal MS-FINDER functions for predicting and ranking structures, correctly identifying 53% of the structures as top-hit, 72% within the top-3 structures, and 78% within the top-10 hits. For our second set, avaniya002, we used both MS-FINDER predictions as well as MS/MS queries against the commercial NIST 14, METLIN, and the public MassBank of North America libraries. Here we correctly identified 78% of the structures as top-hit and 83% within the top-3 hits. Three challenge spectra remained unidentified in either of our submissions within the top-10 hits.     

Mass spectrometric characterization of siderophores produced by <Emphasis Type="Italic">Pseudomonas taiwanensis</Emphasis> VLB120 assisted by stable isotope labeling of nitrogen source

The structures of three previously unknown siderophores produced by the fluorescent, biotechnologically relevant Pseudomonas taiwanensis (P. taiwanensis) VLB120 bacteria were elucidated by means of hydrophilic interaction liquid chromatography (HILIC) hyphenated to high-resolution tandem mass spectrometry (HRMS/MS). They could be verified as iron(III)- and aluminum(III) complexes as well as the protonated molecules of the siderophores formed by in-source fragmentation. The siderophores were separated according to their different acyl side chains and additionally according their central ions. One of the siderophores was identified as pyoverdine with a malic acid (hydroxy succinic acid) amide side chain and a peptide moiety consisting of Orn-Asp-OHAsn-Thr-AcOHOrn-Ser-cOHOrn. The other analytes were assigned to an azotobactin with the identical peptide chain linked to the characteristic chromophoric unit and a pyoverdine with a variation in the amino acid sequence. Proline is directly linked to the pyoverdine chromophore instead of ornithine. Acidic and enzymatic hydrolyses were carried out to analyze the individual amino acids. Beside OHAsn, each amino acid of the peptide part was identified by HILIC–HRMS and comparison to authentic standards. Additionally, ¹⁵N-labeled cellular supernatants were analyzed by means of HRMS/MS. The results of the MS/MS experiments complemented by accurate mass data facilitated elucidation of the structures studied in this work and provided further confirmation of the three recently described pyoverdines of P. taiwanensis VLB120 (Baune et al. in Biometals 30:589–597, 2017.  https://doi.org/10.1007/s10534-017-0029-7).     

Studies on Core-Shell Nanocapsules of Felodipine: <Emphasis Type="Italic">In Vitro</Emphasis>-<Emphasis Type="Italic">In Vivo</Emphasis> Evaluations

The present study aimed for in vitro-in vivo-in silico simulation studies of experimentally designed (3²-factorial) Capmul PG-8-cored, Eudragit RSPO-Lutrol F 127 nanocapsules to ferry felodipine using GastroPlus?. The in silico parameter sensitivity analysis for pharmacokinetic parameters was initially assessed to justify the preparation of felodipine-loaded nanocapsules (FLNs) with enhanced solubility to overcome the bioavailability issues of felodipine. The overall integrated desirability ranged between 0.8187 and 0.9488 for three optimized FLNs when analyzed for mean particle size, zeta potential, encapsulation efficiency, and in vitro dissolution parameters. The morphological evaluation (SEM, TEM, and AFM) demonstrated spherical nanoparticles (200–300 nm). Validated LC-MS/MS analysis demonstrated enhanced relative bioavailability (13.37-fold) of optimized FLN as compared to suspension. The simulated regional absorption of the FLN presented significant absorption from the cecum (26.3%) and ascending colon (20.1%) with overall absorption of 67.4% from the GIT tract. Furthermore, in vitro-in vivo correlation demonstrated the Wagner-Nelson method as the preferred model as compared to mechanistic and numerical deconvolution on the basis of least mean absolute prediction error, least standard error of prediction, least mean absolute error, and maximum correlation coefficient (r ² = 0.920). The study demonstrated enhanced oral absorption of felodipine-loaded nanocapsules, and GastroPlus? was found to be an efficient simulation tool for in vitro-in vivo-in silico simulations.     

Characterization of low-molecular-weight-glutenin subunit genes from the D-genome of Triticum aestivum,Aegilops crassa,Ae. cylindrica and Ae. tauschii

Twenty low-molecular-weight-glutenin subunit (LMW-GS) gene sequences from the D-genome from Aegilops crassa (2n = 4x = 28), Ae. cylindrica (2n = 4x = 28), Ae. tauschii (2n = 2x = 14) and Triticum aestivum (2n = 6x = 42) were obtained using five sets of specific allele primer pairs. Only the sequences of the first primer pair were complete coding sequences (cds) of LMW-GS, and had 305, 304, 306 and 305 LMW-m amino acid residues in Ae. crassa, Ae. cylindrica, Ae. tauschii and T. aestivum, respectively. The repetitive domain and repeat motif numbers of all LMW glutenin subunits showed eight conserved cysteine residues that lead to the same functional activity in different genome. Based on DNA and predicted protein sequences, phylogenetic trees for all sets of sequences were drawn. At the DNA level, the species closest to T. aestivum for the second, third, fourth and fifth set of sequences were Ae. cylindrica, Ae. tauschii and Ae. crassa, respectively. At the protein level, the species closest to T. aestivum based on the first, second and fifth set of sequences were Ae. cylindrica, Ae. crassa and Ae. crassa, respectively. For other sets of sequences, bread wheat proved to be a distinct species. The LMW-GS gene sequences have been recorded in the GenBank with accession numbers JQ726549–JQ726568JQ726549JQ726550JQ726551JQ726552JQ726553JQ726554JQ726555JQ726556JQ726557JQ726558JQ726559JQ726560JQ726561JQ726562JQ726563JQ726564JQ726565JQ726566JQ726567JQ726568.     

Strategies for dereplication of natural compounds using high-resolution tandem mass spectrometry

Complete structural elucidation of natural products is commonly performed by nuclear magnetic resonance spectroscopy (NMR), but annotating compounds to most likely structures using high-resolution tandem mass spectrometry is a faster and feasible first step. The CASMI contest 2016 (Critical Assessment of Small Molecule Identification) provided spectra of eighteen compounds for the best manual structure identification in the natural products category. High resolution precursor and tandem mass spectra (MS/MS) were available to characterize the compounds. We used the Seven Golden Rules, Sirius2 and MS-FINDER software for determination of molecular formulas, and then we queried the formulas in different natural product databases including DNP, UNPD, ChemSpider and REAXYS to obtain molecular structures. We used different in-silico fragmentation tools including CFM-ID, CSI:FingerID and MS-FINDER to rank these compounds. Additional neutral losses and product ion peaks were manually investigated. This manual and time consuming approach allowed for the correct dereplication of thirteen of the eighteen natural products.     

Annotation of the human cerebrospinal fluid lipidome using high resolution mass spectrometry and a dedicated data processing workflow

Introduction

Due to its proximity with the brain, cerebrospinal fluid (CSF) could be a medium of choice for the discovery of biomarkers of neurological and psychiatric diseases using untargeted analytical approaches.

Objectives

This study explored the CSF lipidome in order to generate a robust mass spectral database using an untargeted lipidomic approach.

Methods

Cerebrospinal fluid samples from 45 individuals were analyzed by liquid chromatography coupled to high-resolution mass spectrometry method (LC-HRMS). A dedicated data processing workflow was implemented using XCMS software and adapted filters to select reliable features. In addition, an automatic annotation using an in silico lipid database and several MS/MS experiments were performed to identify CSF lipid species.

Results

Using this complete workflow, 771 analytically relevant monoisotopic lipid species corresponding to 550 unique lipids which represent five major lipid families (i.e., free fatty acids, sphingolipids, glycerophospholipids, glycerolipids, and sterol lipids) were detected and annotated. In addition, MS/MS experiments enabled to improve the annotation of 304 lipid species. Thanks to LC-HRMS, it was possible to discriminate between isobaric and also isomeric lipid species; and interestingly, our study showed that isobaric ions represent about 50 % of the total annotated lipid species in the human CSF.

Conclusion

This work provides an extensive LC/HRMS database of the human CSF lipidome which constitutes a relevant foundation for future studies aimed at finding biomarkers of neurological disorders.

     

Purification and characterization of the extracellular lipase Lip2 from Yarrowia lipolytica

An extracellular lipase from Yarrowia lipolytica (YlLip2) has been purified by ion exchange chromatography on Q sepharose FF, followed by hydrophobic interaction chromatography on butyl sepharose FF. SDS-PAGE showed that the molecular weight of this lipase is about 38 kDa. N-terminal amino acid sequencing and MALDI-TOF mass spectral analysis showed that this lipase is encoded by gene LIP2 (GenBank accession no. AJ012632). Enzymatic deglycosylation showed that this lipase is a glycosylated protein which contains about 12% sugar. The corresponding deglycosylated lipase remained 88% specific activity of untreated lipase. There was a high amino acid sequence identity (91%) between YlLip2 and Candida deformans lipase CdLip1 (GenBank accession no. AJ428393). The optima temperature and pH for the purified lipase was 40 °C and 8.0, respectively. The lipase showed a preference for long chain fatty acid methyl esters (C12–C16), with the highest activity toward methyl myristate (C14). Lipase activity was stimulated by Ca²⁺ and Mg²⁺ and inhibited by Zn²⁺, Ni²⁺ and Cu²⁺, whereas EDTA had no effect on its activity. A 0.1% of Tween 80 and Span 65 increased slightly the enzyme activity and SDS inhibited it.     

Profiling and inhibiting reversible palmitoylation

1. Download : Download high-res image (278KB)
2. Download : Download full-size image

Highlights► New methods allow direct enrichment of palmitoylated proteins for LC–MS annotation. ► SILAC proteomics methods enhance the annotation of low abundance palmitoyl proteins. ► Alkynyl fatty acids facilitate pulse-chase labeling for measuring palmitoylation dynamics. ► Thioesterase inhibitors are promising new tools to dissect palmitoylation dynamics.     

An Integrated Mass-Spectrometry Pipeline Identifies Novel Protein Coding-Regions in the Human Genome

Background

Most protein mass spectrometry (MS) experiments rely on searches against a database of known or predicted proteins, limiting their ability as a gene discovery tool.

Results

Using a search against an in silico translation of the entire human genome, combined with a series of annotation filters, we identified 346 putative novel peptides [False Discovery Rate (FDR)<5%] in a MS dataset derived from two human breast epithelial cell lines. A subset of these were then successfully validated by a different MS technique. Two of these correspond to novel isoforms of Heterogeneous Ribonuclear Proteins, while the rest correspond to novel loci.

Conclusions

MS technology can be used for ab initio gene discovery in human data, which, since it is based on different underlying assumptions, identifies protein-coding genes not found by other techniques. As MS technology continues to evolve, such approaches will become increasingly powerful.     

A novel stable isotope labelling assisted workflow for improved untargeted LC–HRMS based metabolomics research

Many untargeted LC–ESI–HRMS based metabolomics studies are still hampered by the large proportion of non-biological sample derived signals included in the generated raw data. Here, a novel, powerful stable isotope labelling (SIL)-based metabolomics workflow is presented, which facilitates global metabolome extraction, improved metabolite annotation and metabolome wide internal standardisation (IS). The general concept is exemplified with two different cultivation variants, (1) co-cultivation of the plant pathogenic fungi Fusarium graminearum on non-labelled and highly ¹³C enriched culture medium and (2) experimental cultivation under native conditions and use of globally U-¹³C labelled biological reference samples as exemplified with maize and wheat. Subsequent to LC–HRMS analysis of mixtures of labelled and non-labelled samples, two-dimensional data filtering of SIL specific isotopic patterns is performed to better extract truly biological derived signals together with the corresponding number of carbon atoms of each metabolite ion. Finally, feature pairs are convoluted to feature groups each representing a single metabolite. Moreover, the correction of unequal matrix effects in different sample types and the improvement of relative metabolite quantification with metabolome wide IS are demonstrated for the F. graminearum experiment. Data processing employing the presented workflow revealed about 300 SIL derived feature pairs corresponding to 87–135 metabolites in F. graminearum samples and around 800 feature pairs corresponding to roughly 350 metabolites in wheat samples. SIL assisted IS, by the use of globally U-¹³C labelled biological samples, reduced the median CV value from 7.1 to 3.6 % for technical replicates and from 15.1 to 10.8 % for biological replicates in the respective F. graminearum samples.     

Endophytic <Emphasis Type="Italic">Streptomyces</Emphasis> sp. AC35, a producer of bioactive isoflavone aglycones and antimycins

In this research, a microbial endophytic strain obtained from the rhizosphere of the conifer Taxus baccata and designated as Streptomyces sp. AC35 (FJ001754.1 Streptomyces, GenBank) was investigated. High 16S rDNA gene sequence similarity suggests that this strain is closely related to S. odorifer. The major fatty acid profile of intracellular lipids was also carried out to further identify this strain. Atomic force microscopy and scanning acoustic microscopy were used to image our strain. Its major excreted substances were extracted, evaluated for antimicrobial activity, purified, and identified by ultraviolet–visible spectroscopy (UV–vis), liquid chromatography–mass spectrometry (LC–MS/MS) and nuclear magnetic resonance as the bioactive isoflavone aglycones—daidzein, glycitein and genistein. Batch cultivation, performed under different pH conditions, revealed enhanced production of antimycin components when the pH was stable at 7.0. Antimycins were detected by HPLC and identified by UV–vis and LC–MS/MS combined with the multiple reaction monitoring. Our results demonstrate that Streptomyces sp. AC35 might be used as a potential source of effective, pharmaceutically active compounds.     

Dereplication of natural products from complex extracts by regression analysis and molecular networking: case study of redox-active compounds from <Emphasis Type="Italic">Viola alba</Emphasis> subsp. <Emphasis Type="Italic">dehnhardtii</Emphasis>

Introduction

In natural product research, bioassay-guided fractionation was previously widely employed but is now judged to be inadequate in terms of time and cost, particularly if only known compounds are ultimately isolated. The development of metabolomics, along with improvements in analytical tools, allows comprehensive metabolite profiling. This enables dereplication to target unknown active compounds early in the purification workflow.

Objectives

Starting from an ethanolic extract of violet leaves, this study aims to predict redox active compounds within a complex matrix through an untargeted metabolomics approach and correlation analysis.

Methods

Rapid fractionation of crude extracts was carried out followed by multivariate data analysis (MVA) of liquid chromatography–high resolution mass spectrometry (LC–HRMS) profiles. In parallel, redox active properties were evaluated by the capacity of the molecules to reduce 2,2-diphenyl-1-picrylhydrazyl (DPPH^·) and superoxide (O₂ ^·?) radicals using UV–Vis and electron spin resonance spectroscopies (ESR), respectively. A spectral similarity network (molecular networking) was used to highlight clusters involved in the observed redox activities.

Results

Dereplication on Viola alba subsp. dehnhardtii highlighted a reproducible pool of redox active molecules. Polyphenols, particularly O-glycosylated coumarins and C-glycosylated flavonoids, were identified and de novo dereplicated through molecular networking. Confirmatory analyses were undertaken by thin layer chromatography (TLC)–DPPH–MS assays and nuclear magnetic resonance (NMR) spectra of the most active compounds.

Conclusion

Our dereplication strategy allowed the screening of leaf extracts to highlight new biologically active metabolites in few steps with a limited amount of crude material and reduced time-consuming manipulations. This approach could be applied to any kind of natural extract for the study of various biological activities.

     

Tandem mass spectrometry for measuring stable-isotope labeling

1. Download : Download high-res image (187KB)
2. Download : Download full-size image

Highlights► Tandem MS combines efficient separation, high sensitivity and selectivity. ► Tandem MS is used for measuring stable-isotopic labeling for ¹³C-flux analysis. ► Complete positional isotopomer distribution of aspartate was measured by tandem MS. ► Measurement accuracy is improved by 2-fold to10-fold compared to GC–MS and LC–MS. ► Validated libraries of parent–daughter fragments are needed for accurate ¹³C-analysis.     

<Emphasis Type="Italic">In silico</Emphasis> and <Emphasis Type="Italic">in situ</Emphasis> characterization of the zebrafish (<Emphasis Type="Italic">Danio rerio</Emphasis>) <Emphasis Type="Italic">gnrh3</Emphasis> (sGnRH) gene

Background

Gonadotropin releasing hormone (GnRH) is responsible for stimulation of gonadotropic hormone (GtH) in the hypothalamus-pituitary-gonadal axis (HPG). The regulatory mechanisms responsible for brain specificity make the promoter attractive for in silico analysis and reporter gene studies in zebrafish (Danio rerio).

Results

We have characterized a zebrafish [Trp⁷, Leu⁸] or salmon (s) GnRH variant, gnrh 3. The gene includes a 1.6 Kb upstream regulatory region and displays the conserved structure of 4 exons and 3 introns, as seen in other species. An in silico defined enhancer at -976 in the zebrafish promoter, containing adjacent binding sites for Oct-1, CREB and Sp1, was predicted in 2 mammalian and 5 teleost GnRH promoters. Reporter gene studies confirmed the importance of this enhancer for cell specific expression in zebrafish. Interestingly the promoter of human GnRH-I, known as mammalian GnRH (mGnRH), was shown capable of driving cell specific reporter gene expression in transgenic zebrafish.

Conclusions

The characterized zebrafish Gnrh3 decapeptide exhibits complete homology to the Atlantic salmon (Salmo salar) GnRH-III variant. In silico analysis of mammalian and teleost GnRH promoters revealed a conserved enhancer possessing binding sites for Oct-1, CREB and Sp1. Transgenic and transient reporter gene expression in zebrafish larvae, confirmed the importance of the in silico defined zebrafish enhancer at -976. The capability of the human GnRH-I promoter of directing cell specific reporter gene expression in zebrafish supports orthology between GnRH-I and GnRH-III.