Identification of molecules in complex natural matrices relies on matching the fragmentation spectra of ions under investigation and the spectra acquired for the corresponding analytical standards. Currently, there are many databases of experimentally measured tandem mass spectrometry spectra (such as NIST, MzCloud, and Metlin), and considerable progress has been made in the development of software for predicting tandem mass spectrometry fragments in silico using combinatorial, machine learning, and quantum chemistry approaches (such as MetFrag, CFM‐ID, and QCxMS). However, the electrospray ionization molecules can be ionized at different sites (protonated or deprotonated), and the fragmentation spectra of such ions are different. Here, we are using the combination of the in‐ESI source hydrogen/deuterium exchange reaction and MSn fragmentation for the investigation of the fragmentation pathways for different protomers of organic molecules. It is shown that the distribution of the deuterium in the fragment ions reflects the presence of different protomers. For several molecules, the distribution of deuterium was traced up to the MS5 level of fragmentation revealing many unusual and unexpected effects. For example, we investigated the loss of HF from the ciprofloxacin and norfloxacin ions and observed that for ions protonated at –COOH group, the eliminating hydrogen always comes from –NH group. When ions are protonated at another site, the elimination of hydrogen with a probability of 30% occurs from the –NH group, and with a probability of 70%, it originates from other sites on the molecule. Such effects were not described previously. Quantum chemical simulation was used for the verification of the protonated structures and simulation of the corresponding fragmentation spectra.