Despite decades of accumulated knowledge about proteins and their post-translational modifications (PTMs), numerous questions remain regarding their molecular composition and biological function. One of the most fundamental queries is the extent to which the combinations of DNA-, RNA- and PTM-level variations explode the complexity of the human proteome. Here, we outline what we know from current databases and measurement strategies including mass spectrometry-based proteomics. In doing so, we examine prevailing notions about the number of modifications displayed on human proteins and how they combine to generate the protein diversity underlying health and disease. We frame central issues regarding determination of protein-level variation and PTMs, including some paradoxes present in the field today. We use this framework to assess existing data and to ask the question, "How many distinct primary structures of proteins (proteoforms) are created from the 20,300 human genes?" We also explore prospects for improving measurements to better regularize protein-level biology and efficiently associate PTMs to function and phenotype.

The nature of the “toxic gain of function” that results from amyotrophic lateral sclerosis (ALS)-, Parkinson-, and Alzheimer-related mutations is a matter of debate. As a result no adequate model of any neurodegenerative disease etiology exists. We demonstrate that two synergistic properties, namely, increased protein aggregation propensity (increased likelihood that an unfolded protein will aggregate) and decreased protein stability (increased likelihood that a protein will unfold), are central to ALS etiology. Taken together these properties account for 69% of the variability in mutant Cu/Zn-superoxide-dismutase-linked familial ALS patient survival times. Aggregation is a concentration-dependent process, and spinal cord motor neurons have higher concentrations of Cu/Zn-superoxide dismutase than the surrounding cells. Protein aggregation therefore is expected to contribute to the selective vulnerability of motor neurons in familial ALS.

Single-cell analysis is gaining popularity in the field of mass spectrometry as a method for analyzing protein and peptide content in cells. The spatial resolution of MALDI mass spectrometry (MS) imaging is by a large extent limited by the laser focal diameter and the displacement of analytes during matrix deposition. Owing to recent advancements in both laser optics and matrix deposition methods, spatial resolution on the order of a single eukaryotic cell is now achievable by MALDI MS imaging. Provided adequate instrument sensitivity, a lateral resolution of ~10 µm is currently attainable with commercial instruments. As a result of these advances, MALDI MS imaging is poised to become a transformative clinical technology. In this article, the crucial steps needed to obtain single-cell resolution are discussed, as well as potential applications to disease research.

Drug transit through the blood-brain barrier (BBB) is essential for therapeutic responses in malignant glioma. Conventional methods for assessment of BBB penetrance require synthesis of isotopically labeled drug derivatives. Here, we report a new methodology using matrix assisted laser desorption ionization mass spectrometry imaging (MALDI MSI) to visualize drug penetration in brain tissue without molecular labeling. In studies summarized here, we first validate heme as a simple and robust MALDI MSI marker for the lumen of blood vessels in the brain. We go on to provide three examples of how MALDI MSI can provide chemical and biological insights into BBB penetrance and metabolism of small molecule signal transduction inhibitors in the brain – insights that would be difficult or impossible to extract by use of radiolabeled compounds.

Diagnosis of prostate cancer is based on histologic evaluation of tumor architecture using a system known as the "Gleason score." This diagnostic paradigm, while the standard of care, is time-consuming, shows intraobserver variability, and provides no information about the altered metabolic pathways, which result in altered tissue architecture. Characterization of the molecular composition of prostate cancer and how it changes with respect to the Gleason score (GS) could enable a more objective and faster diagnosis. It may also aid in our understanding of disease onset and progression. In this work, we present mass spectrometry imaging for identification and mapping of lipids and metabolites in prostate tissue from patients with known prostate cancer with GS from 6 to 9. A gradient of changes in the intensity of various lipids was observed, which correlated with increasing GS. Interestingly, these changes were identified in both regions of high tumor cell density, and in regions of tissue that appeared histologically benign, possibly suggestive of precancerous metabolomic changes. A total of 31 lipids, including several phosphatidylcholines, phosphatidic acids, phosphatidylserines, phosphatidylinositols, and cardiolipins were detected with higher intensity in GS (4þ3) compared with GS (3þ4), suggesting they may be markers of prostate cancer aggression. Results obtained through mass spectrometry imaging studies were subsequently correlated with a fast, ambient mass spectrometry method for potential use as a clinical tool to support imageguided prostate biopsy. Implications: In this study, we suggest that metabolomic differences between prostate cancers with different Gleason scores can be detected by mass spectrometry imaging.

Mass spectrometry imaging (MSI) is an emerging technology that holds potential for improving, biomarker discovery, metabolomics research, pharmaceutical applications and clinical diagnosis. Despite many solutions being developed, the large data size and high dimensional nature of MSI, especially 3D datasets, still pose computational and memory complexities that hinder accurate identification of biologically relevant molecular patterns. Moreover, the subjectivity in the selection of parameters for conventional pre-processing approaches can lead to bias. Therefore, we assess if a probabilistic generative model based on a fully connected variational autoencoder can be used for unsupervised analysis and peak learning of MSI data to uncover hidden structures. The resulting msiPL method learns and visualizes the underlying non-linear spectral manifold, revealing biologically relevant clusters of tissue anatomy in a mouse kidney and tumor heterogeneity in human prostatectomy tissue, colorectal carcinoma, and glioblastoma mouse model, with identification of underlying m/z peaks. The method is applied for the analysis of MSI datasets ranging from 3.3 to 78.9 GB, without prior pre-processing and peak picking, and acquired using different mass spectrometers at different centers.

Multimodal integration between mass spectrometry imaging (MSI) and radiology-established modalities such as magnetic resonance imaging (MRI) would allow the investigations of key questions in complex biological systems such as the central nervous system. Such integration would provide complementary multiscale data to bridge the gap between molecular and anatomical phenotypes, potentially revealing new insights into molecular mechanisms underlying anatomical pathologies presented on MRI. Automatic co-registration between 3D MSI/MRI is a computationally challenging process due to dimensional complexity, MSI data sparsity, lack of direct spatial-correspondences, and non-linear tissue deformation. Here, we present a new computational approach based on stochastic neighbor embedding to non-linearly align 3D MSI to MRI data, identify and reconstruct biologically-relevant molecular patterns in 3D, and fuse the MSI datacube to the MRI space. We demonstrate our method using multimodal high-spectral resolution MALDI 9.4 Tesla MSI and 7 Tesla in vivo MRI data, acquired from a patient-derived, xenograft mouse brain model of glioblastoma following administration of the EGFR inhibitor drug of Erlotinib. Results show the distribution of some identified molecular ions of the EGFR inhibitor

Therapeutic options for the treatment of glioblastoma remain inadequate despite concerted research efforts in drug development. Therapeutic failure can result from poor permeability of the blood-brain barrier, heterogeneous drug distribution, and development of resistance. Elucidation of relationships among such parameters could enable the development of predictive models of drug response in patients and inform drug development. Complementary analyses were applied to a glioblastoma patient-derived xenograft model in order to quantitatively map distribution and resulting cellular response to the EGFR inhibitor erlotinib. Mass spectrometry images of erlotinib were registered to histology and magnetic resonance images in order to correlate drug distribution with tumor characteristics. Phosphoproteomics and immunohistochemistry were used to assess protein signaling in response to drug, and integrated with transcriptional response using mRNA sequencing. This comprehensive dataset provides simultaneous insight into pharmacokinetics and pharmacodynamics and indicates that erlotinib delivery to intracranial tumors is insufficient to inhibit EGFR tyrosine kinase signaling.