Production of reliable MALDI spectra with quality threshold clustering of replicates

BitQT: a graph-based approach to the quality threshold clustering of molecular dynamics

MOTIVATION

Classical Molecular Dynamics (MD) is a standard computational approach to model time-dependent processes at the atomic level. The inherent sparsity of increasingly huge generated trajectories demands clustering algorithms to reduce other post-simulation analysis complexity. The Quality Threshold (QT) variant is an appealing one from the vast number of available clustering methods. It guarantees that all members of a particular cluster will maintain a collective similarity established by a user-defined threshold. Unfortunately, its high computational cost for processing big data limits its application in the molecular simulation field.

RESULTS

In this work, we propose a methodological parallel between QT clustering and another well-known algorithm in the field of Graph Theory, the Maximum Clique Problem. Molecular trajectories are represented as graphs whose nodes designate conformations, while unweighted edges indicate mutual similarity between nodes. The use of a binary-encoded RMSD matrix coupled to the exploitation of bitwise operations to extract clusters significantly contributes to reaching a very affordable algorithm compared to the few implementations of QT for MD available in the literature. Our alternative provides results in good agreement with the exact one while strictly preserving the collective similarity of clusters.

AVAILABILITY AND IMPLEMENTATION

The source code and documentation of BitQT are free and publicly available on GitHub (https://github.com/LQCT/BitQT.git) and ReadTheDocs (https://bitqt.readthedocs.io/en/latest/), respectively.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Quality Threshold Clustering of Molecular Dynamics: A Word of Caution

Clustering Molecular Dynamics trajectories is a common analysis that allows grouping together similar conformations. Several algorithms have been designed and optimized to perform this routine task, and among them, Quality Threshold stands as a very attractive option. This algorithm guarantees that in retrieved clusters no pair of frames will have a similarity value greater than a specified threshold, and hence, a set of strongly correlated frames are obtained for each cluster. In this work, it is shown that various commonly used software implementations are flawed by confusing Quality Threshold with another simplistic well-known clustering algorithm published by Daura et al. (Daura, X.; van Gunsteren, W. F.; Jaun, B.; Mark, A. E.; Gademann, K.; Seebach, D. Peptide Folding: When Simulation Meets Experiment. Angew. Chemie Int. Ed. 1999, 38 (1/2), 236-240). Daura's algorithm does not impose any quality threshold for the frames contained in retrieved clusters, bringing unrelated structural configurations altogether. The advantages of using Quality Threshold whenever possible to explore Molecular Dynamic trajectories is exemplified. An in-house implementation of the original Quality Threshold algorithm has been developed in order to illustrate our comments, and its code is freely available for further use by the scientific community.

Mass spectrometry imaging of lipids: untargeted consensus spectra reveal spatial distributions in Niemann-Pick disease type C1

Mass spectrometry imaging (MSI) is a tool to rapidly map the spatial location of analytes without the need for tagging or a reporter system. Niemann-Pick disease type C1 (NPC1) is a neurodegenerative, lysosomal storage disorder characterized by accumulation of unesterified cholesterol and sphingolipids in the endo-lysosomal system. Here, we use MSI to visualize lipids including cholesterol in cerebellar brain tissue from the NPC1 symptomatic mouse model and unaffected controls. To complement the imaging studies, a data-processing pipeline was developed to generate consensus mass spectra, thereby using both technical and biological image replicates to assess differences. The consensus spectra are used to determine true differences in lipid relative abundance; lipid distributions can be determined in an unbiased fashion without prior knowledge of location. We show the cerebellar distribution of gangliosides GM1, GM2, and GM3, including variants of lipid chain length. We also performed MALDI-MSI of cholesterol. Further analysis of lobules IV/V and X of the cerebellum gangliosides indicates regional differences. The specificity achieved highlights the power of MSI, and this new workflow demonstrates a universal approach for addressing reproducibility in imaging experiments applied to NPC1.

The quest for improved reproducibility in MALDI mass spectrometry

Reproducibility has been one of the biggest hurdles faced when attempting to develop quantitative protocols for MALDI mass spectrometry. The heterogeneous nature of sample recrystallization has made automated sample acquisition somewhat "hit and miss" with manual intervention needed to ensure that all sample spots have been analyzed. In this review, we explore the last 30 years of literature and anecdotal evidence that has attempted to address and improve reproducibility in MALDI MS. Though many methods have been attempted, we have discovered a significant publication history surrounding the use of nitrocellulose as a substrate to improve homogeneity of crystal formation and therefore reproducibility. We therefore propose that this is the most promising avenue of research for developing a comprehensive and universal preparation protocol for quantitative MALDI MS analysis. © 2016 Wiley Periodicals, Inc. Mass Spec Rev 37:217-228, 2018.

Biomarker Characterization by MALDI-TOF/MS

Mass spectrometric techniques frequently used in clinical diagnosis, such as gas chromatography-mass spectrometry, liquid chromatography-mass spectrometry, ambient ionization mass spectrometry, and matrix-assisted laser desorption ionization/time-of-flight mass spectrometry (MALDI-TOF/MS), are discussed. Due to its ability to rapidly detect large biomolecules in trace amounts, MALDI-TOF/MS is an ideal tool for characterizing disease biomarkers in biologic samples. Clinical applications of MS for the identification and characterization of microorganisms, DNA fragments, tissues, and biofluids are introduced. Approaches for using MALDI-TOF/MS to detect various disease biomarkers including peptides, proteins, and lipids in biological fluids are further discussed. Finally, various sample pretreatment methods which improve the detection efficiency of disease biomarkers are introduced.

BRAF pyrosequencing analysis aided by a lookup table

OBJECTIVES

BRAF mutations have substantial therapeutic, diagnostic, and prognostic significance, so detecting and specifying them is an important part of the workload of molecular pathology laboratories. Pyrosequencing assays are well suited for this analysis but can produce complex results. Therefore, we introduce a pyrosequencing lookup table based on Pyromaker that assists the user in generating hypotheses for solving complex pyrosequencing results.

METHODS

The lookup table contains all known mutations in the sequenced region and the positions in the dispensation sequence at which changes would occur with those mutations. We demonstrate the lookup table using a homebrew dispensation sequence for BRAF codons 596 to 605 as well as a commercially available kit-based dispensation sequence for codons 599 to 600.

RESULTS

These results demonstrate that the homebrew dispensation sequence unambiguously identifies all known BRAF mutations in this region, whereas the kit-based dispensation sequence has one unresolvable degeneracy that could be solved with the addition of two injections.

CONCLUSIONS

Using the lookup table and confirmatory virtual pyrogram, we unambiguously solved clinical pyrograms of the complex mutations V600K (c.1798_1799delGTinsAA), V600R (c.1798_1799delGTinsAG), V600D (c.1799_1800delTGinsAT), V600E (c.1799_1800delTGinsAA), and V600_K601delinsE (c.1799_1801delTGA). In addition, we used the approach to hypothesize and confirm a new mutation in human melanoma, V600_K601delinsEI (c.1799_1802delTGAAinsAAAT).

Comprehensive analysis of LC/MS data using pseudocolor plots

We have developed new applications of the pseudocolor plot for the analysis of LC/MS data. These applications include spectral averaging, analysis of variance, differential comparison of spectra, and qualitative filtering by compound class. These applications have been motivated by the need to better understand LC/MS data generated from analysis of human biofluids. The examples presented use data generated to profile steroid hormones in urine extracts from a Cushing's disease patient relative to a healthy control, but are general to any discovery-based scanning mass spectrometry technique. In addition to new visualization techniques, we introduce a new metric of variance: the relative maximum difference from the mean. We also introduce the concept of substructure-dependent analysis of steroid hormones using precursor ion scans. These new analytical techniques provide an alternative approach to traditional untargeted metabolomics workflow. We present an approach to discovery using MS that essentially eliminates alignment or preprocessing of spectra. Moreover, we demonstrate the concept that untargeted metabolomics can be achieved using low mass resolution instrumentation.