Mass spectrometry based metabolomics for in vitro systems pharmacology: pitfalls, challenges, and computational solutions.

INTRODUCTION: Mass spectrometry based metabolomics has become a promising complement and alternative to transcriptomics and proteomics in many fields including in vitro systems pharmacology. Despite several merits, metabolomics based on liquid chromatography mass spectrometry (LC-MS) is a developing area that is yet attached to several pitfalls and challenges. To reach a level of high reliability and robustness, these issues need to be tackled by implementation of refined experimental and computational protocols. OBJECTIVES: This study illustrates some key pitfalls in LC-MS based metabolomics and introduces an automated computational procedure to compensate for them. METHOD: Non-cancerous mammary gland derived cells were exposed to 27 chemicals from four pharmacological classes plus a set of six pesticides. Changes in the metabolome of cell lysates were assessed after 24 h using LC-MS. A data processing pipeline was established and evaluated to handle issues including contaminants, carry over effects, intensity decay and inherent methodology variability and biases. A key component in this pipeline is a latent variable method called OOS-DA (optimal orthonormal system for discriminant analysis), being theoretically more easily motivated than PLS-DA in this context, as it is rooted in pattern classification rather than regression modeling. RESULT: The pipeline is shown to reduce experimental variability/biases and is used to confirm that LC-MS spectra hold drug class specific information. CONCLUSION: LC-MS based metabolomics is a promising methodology, but comes with pitfalls and challenges. Key difficulties can be largely overcome by means of a computational procedure of the kind introduced and demonstrated here. The pipeline is freely available on www.github.com/stephanieherman/MS-data-processing.

Detection of cell aggregation and altered cell viability by automated label-free video microscopy: a promising alternative to endpoint viability assays in high-throughput screening.

Automated phase-contrast video microscopy now makes it feasible to monitor a high-throughput (HT) screening experiment in a 384-well microtiter plate format by collecting one time-lapse video per well. Being a very cost-effective and label-free monitoring method, its potential as an alternative to cell viability assays was evaluated. Three simple morphology feature extraction and comparison algorithms were developed and implemented for analysis of differentially time-evolving morphologies (DTEMs) monitored in phase-contrast microscopy videos. The most promising layout, pixel histogram hierarchy comparison (PHHC), was able to detect several compounds that did not induce any significant change in cell viability, but made the cell population appear as spheroidal cell aggregates. According to recent reports, all these compounds seem to be involved in inhibition of platelet-derived growth factor receptor (PDGFR) signaling. Thus, automated quantification of DTEM (AQDTEM) holds strong promise as an alternative or complement to viability assays in HT in vitro screening of chemical compounds.

NMR spectroscopy-based metabolic profiling of drug-induced changes in vitro can discriminate between pharmacological classes.

Drug-induced changes in mammalian cell line models have already been extensively profiled at the systemic mRNA level and subsequently used to suggest mechanisms of action for new substances, as well as to support drug repurposing, i.e., identifying new potential indications for drugs already licensed for other pharmacotherapy settings. The seminal work in this field, which includes a large database and computational algorithms for pattern matching, is known as the "Connectivity Map" (CMap). However, the potential of similar exercises at the metabolite level is still largely unexplored. Only recently, the first high-throughput metabolomic assay pilot study was published, which involved screening the metabolic response to a set of 56 kinase inhibitors in a 96-well format. Here, we report results from a separately developed metabolic profiling assay, which leverages (1)H NMR spectroscopy to the quantification of metabolic changes in the HCT116 colorectal cancer cell line, in response to each of 26 compounds. These agents are distributed across 12 different pharmacological classes covering a broad spectrum of bioactivity. Differential metabolic profiles, inferred from multivariate spectral analysis of 18 spectral bins, allowed clustering of the most-tested drugs, according to their respective pharmacological class. A more-advanced supervised analysis, involving one multivariate scattering matrix per pharmacological class and using only 3 spectral bins (3 metabolites), showed even more distinct pharmacology-related cluster formations. In conclusion, this type of relatively fast and inexpensive profiling seems to provide a promising alternative to that afforded by mRNA expression analysis, which is relatively slow and costly. As also indicated by the present pilot study, the resulting metabolic profiles do not seem to provide as information-rich signatures as those obtained using systemic mRNA profiling, but the methodology holds strong promise for significant refinement.

Label free high throughput screening for apoptosis inducing chemicals using time-lapse microscopy signal processing.

Label free time-lapse microscopy has opened a new avenue to the study of time evolving events in living cells. When combined with automated image analysis it provides a powerful tool that enables automated large-scale spatiotemporal quantification at the cell population level. Very few attempts, however, have been reported regarding the design of image analysis algorithms dedicated to the detection of apoptotic cells in such time-lapse microscopy images. In particular, none of the reported attempts is based on sufficiently fast signal processing algorithms to enable large-scale detection of apoptosis within hours/days without access to high-end computers. Here we show that it is indeed possible to successfully detect chemically induced apoptosis by applying a two-dimensional linear matched filter tailored to the detection of objects with the typical features of an apoptotic cell in phase-contrast images. First a set of recorded computational detections of apoptosis was validated by comparison with apoptosis specific caspase activity readouts obtained via a fluorescence based assay. Then a large screen encompassing 2,866 drug like compounds was performed using the human colorectal carcinoma cell line HCT116. In addition to many well known inducers (positive controls) the screening resulted in the detection of two compounds here reported for the first time to induce apoptosis.

Label-free detection and dynamic monitoring of drug-induced intracellular vesicle formation enabled using a 2-dimensional matched filter.

Analysis of vesicle formation and degradation is a central issue in autophagy research and microscopy imaging is revolutionizing the study of such dynamic events inside living cells. A limiting factor is the need for labeling techniques that are labor intensive, expensive, and not always completely reliable. To enable label-free analyses we introduced a generic computational algorithm, the label-free vesicle detector (LFVD), which relies on a matched filter designed to identify circular vesicles within cells using only phase-contrast microscopy images. First, the usefulness of the LFVD is illustrated by presenting successful detections of autophagy modulating drugs found by analyzing the human colorectal carcinoma cell line HCT116 exposed to each substance among 1266 pharmacologically active compounds. Some top hits were characterized with respect to their activity as autophagy modulators using independent in vitro labeling of acidic organelles, detection of LC3-II protein, and analysis of the autophagic flux. Selected detection results for 2 additional cell lines (DLD1 and RKO) demonstrate the generality of the method. In a second experiment, label-free monitoring of dose-dependent vesicle formation kinetics is demonstrated by recorded detection of vesicles over time at different drug concentrations. In conclusion, label-free detection and dynamic monitoring of vesicle formation during autophagy is enabled using the LFVD approach introduced.