Interpretation of alkyl diselenide and selenosulfenate mass spectra

Masserstein: Linear regression of mass spectra by optimal transport

RATIONALE

The linear regression of mass spectra is a computational problem defined as fitting a linear combination of reference spectra to an experimental one. It is typically used to estimate the relative quantities of selected ions. In this work, we study this problem in an abstract setting to develop new approaches applicable to a diverse range of experiments.

METHODS

To overcome the sensitivity of the ordinary least-squares regression to measurement inaccuracies, we base our methods on a non-conventional spectral dissimilarity measure, known as the Wasserstein or the Earth Mover's distance. This distance is based on the notion of the cost of transporting signal between mass spectra, which renders it naturally robust to measurement inaccuracies in the mass domain.

RESULTS

Using a data set of 200 mass spectra, we show that our approach is capable of estimating ion proportions accurately without extensive preprocessing of spectra required by other methods. The conclusions are further substantiated using data sets simulated in a way that mimics most of the measurement inaccuracies occurring in real experiments.

CONCLUSIONS

We have developed a linear regression algorithm based on the notion of the cost of transporting signal between spectra. Our implementation is available in a Python 3 package called masserstein, which is freely available at https://github.com/mciach/masserstein.

Topological and quantum molecular descriptors as effective tools for analyzing cytotoxic activity achieved by a series of (diselanediyldibenzene-4,1-diylnide)biscarbamate derivatives

A molecular modeling study has been carried out on a previously reported series of (diselanediyldibenzene-4,1-diylnide)biscarbamate derivatives that show cytotoxic and antiproliferative in vitro activity against MCF-7 human cell line; radical scavenging properties were also confirmed when these compounds were tested for their ability to scavenge DPPH and ABTS radicals. The data obtained allowed us to classify the compounds into two different groups: (a) aliphatic carbamates for which the activity could be related with a first nucleophilic attack (mediated by H₂O, for example) on the selenium atoms of the central scaffold, followed by the release of the alkyl N-(4-selanylphenyl) and N-(4-selenenophenyl)carbamate moieties. Then, a second nucleophilic attack on the carbamate moiety, to yield 4-aminobenzeneselenol and 4-selenenoaniline respectively, which can ultimately be responsible for the activity of the compounds; (b) aromatic carbamates, for which we propose a preferred nucleophilic attack on the carbamate moiety, yielding 4-[(4-aminophenyl)diselanyl]aniline, the common structural fragment for this series, for which we have previously demonstrated its cytotoxic profile. Then, selenium atoms of the central fragment may later undergo a new nucleophilic attack, to yield 4-selenenoaniline and 4-aminobenzeneselenol. The phenolic moieties released in this process may also have a synergistic cytotoxic and redox activity. The data that support this connection include the conformational behavior and the molecular topography of the derivatives which can influence the accessibility of the hydrolysis points, and some quantum descriptors (bond order, atomic charges, total valences, ionization potential, electron affinity, HOMO 0 and LUMO 0 location, etc.) that have been related to the biological activity of the compounds.

Electron capture dissociation of disulfide, sulfur-selenium, and diselenide bound peptides

To examine the electron capture dissociation (ECD) behavior of disulfide (S-S), sulfur-selenium (S-Se), and diselenide (Se-Se) bonds-containing peptides, a series of free cysteine (Cys) and selenocysteine (Sec) containing peptides were reacted to form interchain S-S, S-Se, and Se-Se bonds, and then studied using ECD with Fourier transform ion cyclotron mass spectrometry (FTICR MS). These results demonstrate that the radical has higher tendency to stay at selenium rather than sulfur after the cleavage of Se-S bonds by ECD. In addition, -SH (-33), -S (-32), and -S + H (-31) small neutral losses were all observed from the cleavage of C-S bonds of a disulfide bound peptide. Similar, but minor, fragments were also detected in S-Se bound peptides. In contrast, the cleavage of C-Se bonds of the Se-Se species mainly forms fragments with neutral loss of -Se + H (-78.90868), and the radical tends to stay on the selenium of its corresponding complementary pair. Although the electron affinities of S atom (2.07 eV) and Se atom (2.02 eV) are very close; they have very different reactivity towards electrons. The replacement of sulfur with selenium greatly increases the electron affinities of S-Se and Se-Se bonds comparing to S-S bonds (with an increase of electron affinity by about 0.20 eV by replacing a sulfur with a selenium) (Int J Quantum Chem 110:513-523, 2010), which in turn leads to different ECD fragmentation behavior and mechanisms. Our results are in good agreement with previously published ab initio calculations on Se-Se compounds by other groups.

Mathematical tools in analytical mass spectrometry

Over the last few decades, mass spectrometry has become a powerful tool for exploring various aspects of molecular processes occurring in biological systems. Such exploration is leading to a greater understanding of various complex life processes; unraveling these processes poses the greatest challenge to contemporary bioscience. With due respect to sample preparation, data analysis is rapidly becoming a major obstacle to the conversion of experimental knowledge into valid conclusions. It is interesting to note that many problems related to mass spectrometry can be solved using techniques from computer science, graph theory and discrete mathematics. The aim of this manuscript is to recollect several essays that demonstrate the power and the need to apply such skills to mass spectrometry data interpretation. Special attention is paid to situations where traditional chemical analysis reaches its limits but mathematical reasoning can still allow us to reach valid conclusions.

Interpretation of butyltin mass spectra using isotope pattern reconstruction for the accurate measurement of isotope ratios from molecular clusters

The fragmentation patterns of butyltin compounds (mono-, di-, and tributyltin) in an electron impact ion source were studied using an isotope pattern reconstruction algorithm with emphasis on isotope ratio measurements from molecular clusters. For this purpose, standards of natural tin isotope abundance and a (119)Sn-enriched mixture of the three compounds were both ethylated and propylated using sodium tetraalkylborates. The corresponding mass spectra of the various tetraalkyltin compounds prepared were obtained by GC/MS after their extraction with hexane. The results showed that pure interference-free molecular clusters were obtained only for certain R(3)Sn(+) ions where no isobaric overlap with R(2)SnH(+) ions occurred (e.g. BuEt(2)Sn(+) overlaps with Bu(2)SnH(+)). These ions are ideal candidates for accurate Sn isotope ratio measurements, while isotope pattern perturbing interferences are observed for other molecular fragments down to Sn(.)(+). Isotope pattern reconstruction algorithm thus can be used as an analytical tool to ensure the absence of molecular interferences--a requirement for accurate isotope ratio measurements from molecular clusters. The relevance of these studies for the determination of butyltin compounds in environmental samples by isotope dilution GC/MS is also discussed.