Effects of sodium cationization versus protonation on the conformations and N-glycosidic bond stabilities of sodium cationized Urd and dUrd: solution conformation of [Urd+Na]+ is preserved upon ESI

Solvation effects on glyphosate protonation and deprotonation states evaluated by mass spectrometry and explicit solvation simulations

Glyphosate is a widely used herbicide, and its protonation and deprotonation sites are fundamental to understanding its properties. In this work, the sodiated, protonated, and deprotonated glyphosate were evaluated in the gas phase by infrared multiple photon dissociation spectroscopy to determine the exact nature of these coordination, protonation, and deprotonation states in the gas phase. In this context, Natural Bond Orbital analyses were carried out to unravel interactions that govern glyphosate (de)protonation states in the gas phase. The solvent effect on the protonation/deprotonation equilibria was also investigated by implicit (Solvation Model Based on Density and polarizable continuum models) and explicit solvation models (Monte Carlo and Molecular Dynamics simulations). These results show that glyphosate is protonated in the phosphonate group in the gas phase because of the strong hydrogen bond between the carboxylic oxygen (O7) and the protonated phosphonate group (O8-H19), while the most stable species in water is protonated at the amino group because of the preferential interaction of the NH₂ ⁺ group and the solvent water molecules. Similarly, deprotonated glyphosate [Glyp-H]^- was shown to be deprotonated at the phosphonate group in the gas phase but not in solution, also because of the preferential solvation of the NH₂ ⁺ group present in the other deprotomers. Therefore, these results show that the stabilization of the protonated amino group by the solvent molecules is the governing factor of the (de)protonation equilibrium of glyphosate in water.

5-Halogenation of Uridine Suppresses Protonation-Induced Tautomerization and Enhances Glycosidic Bond Stability of Protonated Uridine: Investigations via IRMPD Action Spectroscopy, ER-CID Experiments, and Theoretical Calculations

Uridine (Urd), a canonical nucleoside of RNA, is the most commonly modified nucleoside among those that occur naturally. Uridine has also been an important target for the development of modified nucleoside analogues for pharmaceutical applications. In this work, the effects of 5-halogenation of uracil on the structures and glycosidic bond stabilities of protonated uridine nucleoside analogues are examined using tandem mass spectrometry and computational methods. Infrared multiple photon dissociation (IRMPD) action spectroscopy experiments and theoretical calculations are performed to probe the structural influences of these modifications. Energy-resolved collision-induced dissociation experiments along with survival yield analyses are performed to probe glycosidic bond stability. The measured IRMPD spectra are compared to linear IR spectra predicted for the stable low-energy conformations of these species computed at the B3LYP/6-311+G(d,p) level of theory to determine the conformations experimentally populated. Spectral signatures in the IR fingerprint and hydrogen-stretching regions allow the 2,4-dihydroxy protonated tautomers (T) and O4- and O2-protonated conformers to be readily differentiated. Comparisons between the measured and predicted spectra indicate that parallel to findings for uridine, both T and O4-protonated conformers of the 5-halouridine nucleoside analogues are populated, whereas O2-protonated conformers are not. Variations in yields of the spectral signatures characteristic of the T and O4-protonated conformers indicate that the extent of protonation-induced tautomerization is suppressed as the size of the halogen substituent increases. Trends in the energy-dependence of the survival yield curves find that 5-halogenation strengthens the glycosidic bond and that the enhancement in stability increases with the size of the halogen substituent.

p-Aminobenzoic acid protonation dynamics in an evaporating droplet by ab initio molecular dynamics

Protonation equilibria are known to vary from the bulk to microdroplet conditions, which could induce many chemical and physical phenomena. Protonated p-aminobenzoic acid (PABA + H⁺) can be considered a model for probing the protonation dynamics in an evaporating droplet, as its protonation equilibrium is highly dependent on the formation conditions from solution via atmospheric pressure ionization sources. Experiments using diverse experimental techniques have shown that protic solvents allow formation of the O-protomer (PABA protonated in the carboxylic acid group) stable in the gas phase, while aprotic solvents yield the N-protomer (protonated in the amino group) that is the most stable protomer in solution. In this work, we explore the protonation equilibrium of PABA solvated by different numbers of water molecules (n = 0 to 32) using ab initio molecular dynamics. For n = 8-32, the protonation is either at the NH₂ group or in the solvent network. The solvent network interacts with the carboxylic acid group, but there is no complete proton transfer to form the O-protomer. For smaller clusters, however, solvent-mediated proton transfers to the carboxylic acid were observed, both via the Grotthuss mechanism and the vehicle or shuttle mechanism (for n = 1 and 2). Thermodynamic considerations allowed a description of the origins of the kinetic trapping effect, which explains the observation of the solution structure in the gas phase. This effect likely occurs in the final evaporation steps, which are outside the droplet size range covered by previous classical molecular dynamics simulations of charged droplets. These results may be considered relevant in determining the nature of the species observed in the ubiquitous ESI based mass spectrometry analysis, and in general for droplet chemistry, explaining how protonation equilibria are drastically changed from bulk to microdroplet conditions.

Structural and Energetic Effects of O2'-Ribose Methylation of Protonated Pyrimidine Nucleosides

The 2'-substituents distinguish DNA from RNA nucleosides. 2'-O-methylation occurs naturally in RNA and plays important roles in biological processes. Such 2'-modifications may alter the hydrogen-bonding interactions of the nucleoside and thus may affect the conformations of the nucleoside in an RNA chain. Structures of the protonated 2'-O-methylated pyrimidine nucleosides were examined by infrared multiple photon dissociation (IRMPD) action spectroscopy, assisted by electronic structure calculations. The glycosidic bond stabilities of the protonated 2'-O-methylated pyrimidine nucleosides, [Nuom+H]⁺, were also examined and compared to their DNA and RNA nucleoside analogues via energy-resolved collision-induced dissociation (ER-CID). The preferred sites of protonation of the 2'-O-methylated pyrimidine nucleosides parallel their canonical DNA and RNA nucleoside analogues, [dNuo+H]⁺ and [Nuo+H]⁺, yet their nucleobase orientation and sugar puckering differ. The glycosidic bond stabilities of the protonated pyrimidine nucleosides follow the order: [dNuo+H]⁺ < [Nuo+H]⁺ < [Nuom+H]⁺. The slightly altered structures help explain the stabilization induced by 2'-O-methylation of the pyrimidine nucleosides.

The effect of adenine protonation on RNA phosphodiester backbone bond cleavage elucidated by deaza-nucleobase modifications and mass spectrometry

The catalytic strategies of small self-cleaving ribozymes often involve interactions between nucleobases and the ribonucleic acid (RNA) backbone. Here we show that multiply protonated, gaseous RNA has an intrinsic preference for the formation of ionic hydrogen bonds between adenine protonated at N3 and the phosphodiester backbone moiety on its 5'-side that facilitates preferential phosphodiester backbone bond cleavage upon vibrational excitation by low-energy collisionally activated dissociation. Removal of the basic N3 site by deaza-modification of adenine was found to abrogate preferential phosphodiester backbone bond cleavage. No such effects were observed for N1 or N7 of adenine. Importantly, we found that the pH of the solution used for generation of the multiply protonated, gaseous RNA ions by electrospray ionization affects phosphodiester backbone bond cleavage next to adenine, which implies that the protonation patterns in solution are at least in part preserved during and after transfer into the gas phase. Our study suggests that interactions between protonated adenine and phosphodiester moieties of RNA may play a more important mechanistic role in biological processes than considered until now.

Structures and Relative Glycosidic Bond Stabilities of Protonated 2'-Fluoro-Substituted Purine Nucleosides

The 2'-substituent is the primary distinguishing feature between DNA and RNA nucleosides. Modifications to this critical position, both naturally occurring and synthetic, can produce biologically valuable nucleoside analogues. The unique properties of fluorine make it particularly interesting and medically useful as a synthetic nucleoside modification. In this work, the effects of 2'-fluoro modification on the protonated gas-phase purine nucleosides are examined using complementary tandem mass spectrometry and computational methods. Direct comparisons are made with previous studies on related nucleosides. Infrared multiple photon dissociation action spectroscopy performed in both the fingerprint and hydrogen-stretching regions allows for the determination of the experimentally populated conformations. The populated conformers of protonated 2'-fluoro-2'-deoxyadenosine, [Adofl+H]⁺, and 2'-fluoro-2'-deoxyguanosine, [Guofl+H]⁺, are highly parallel to their respective canonical DNA and RNA counterparts. Both N3 and N1 protonation sites are accessed by [Adofl+H]⁺, stabilizing syn and anti nucleobase orientations, respectively. N7 protonation and anti nucleobase orientation dominates in [Guofl+H]⁺. Spectroscopically observable intramolecular hydrogen-bonding interactions with fluorine allow more definitive sugar puckering determinations than possible for the canonical systems. [Adofl+H]⁺ adopts C2'-endo sugar puckering, whereas [Guofl+H]⁺ adopts both C2'-endo and C3'-endo sugar puckering. Energy-resolved collision-induced dissociation experiments with survival yield analyses provide relative glycosidic bond stabilities. The N-glycosidic bond stabilities of the protonated 2'-fluoro-substituted purine nucleosides are found to exceed those of their canonical analogues. Further, the N-glycosidic bond stability is found to increase with increasing electronegativity of the 2'-substituent, i.e., H < OH < F. The N-glycosidic bond stability is also greater for the adenine nucleoside analogues than the guanine nucleoside analogues.

Amino acid-linked platinum(II) compounds: non-canonical nucleoside preferences and influence on glycosidic bond stabilities

Abstract

Nucleobases serve as ideal targets where drugs bind and exert their anticancer activities. Cisplatin (cisPt) preferentially coordinates to 2′-deoxyguanosine (dGuo) residues within DNA. The dGuo adducts that are formed alter the DNA structure, contributing to inhibition of function and ultimately cancer cell death. Despite its success as an anticancer drug, cisPt has a number of drawbacks that reduce its efficacy, including repair of adducts and drug resistance. Some approaches to overcome this problem involve development of compounds that coordinate to other purine nucleobases, including those found in RNA. In this work, amino acid-linked platinum(II) (AAPt) compounds of alanine and ornithine (AlaPt and OrnPt, respectively) were studied. Their reactivity preferences for DNA and RNA purine nucleosides (i.e., 2′-deoxyadenosine (dAdo), adenosine (Ado), dGuo, and guanosine (Guo)) were determined. The chosen compounds form predominantly monofunctional adducts by reacting at the N1, N3, or N7 positions of purine nucleobases. In addition, features of AAPt compounds that impact the glycosidic bond stability of Ado residues were explored. The glycosidic bond cleavage is activated differentially for AlaPt-Ado and OrnPt-Ado isomers. Formation of unique adducts at non-canonical residues and subsequent destabilization of the glycosidic bonds are important features that could circumvent platinum-based drug resistance.

Graphic abstract

Electronic supplementary material

The online version of this article (10.1007/s00775-019-01693-y) contains supplementary material, which is available to authorized users.

Relative glycosidic bond stabilities of naturally occurring methylguanosines: 7-methylation is intrinsically activating

The frequency and diversity of posttranscriptional modifications add an additional layer of chemical complexity beyond canonical nucleic acid sequence. Methylations are particularly frequently occurring and often highly conserved throughout the kingdoms of life. However, the intricate functions of these modified nucleic acid constituents are often not fully understood. Systematic foundational research that reduces systems to their minimum constituents may aid in unraveling the complexities of nucleic acid biochemistry. Here, we examine the relative intrinsic N-glycosidic bond stabilities of guanosine and five naturally occurring methylguanosines (O2'-, 1-, 7-, N2,N2-di-, and N2,N2,O2'-trimethylguanosine) probed by energy-resolved collision-induced dissociation tandem mass spectrometry and complemented with quantum chemical calculations. Apparent glycosidic bond stability is generally found to increase with increasing methyl substitution (canonical < mono- < di- < trimethylated). Many biochemical transformations, including base excision repair mechanisms, involve protonation and/or noncovalent interactions to increase nucleobase leaving-group ability. The protonated gas-phase methylguanosines require less activation energy for glycosidic bond cleavage than their sodium cationized forms. However, methylation at the N7 position intrinsically weakens the glycosidic bond of 7-methylguanosine more significantly than subsequent cationization, and thus 7-methylguanosine is suggested to be under perpetually activated conditions. N7 methylation also alters the nucleoside geometric preferences relative to the other systems, including the nucleobase orientation in the neutral form, sugar puckering in the protonated form, and the preferred protonation and sodium cation binding sites. All of the methylated guanosines examined here are predicted to have proton affinities and gas-phase basicities that exceed that of canonical guanosine. Additionally, the proton affinity and gas-phase basicity trends exhibit a roughly inverse correlation with the apparent glycosidic bond stabilities.

Modified Quadrupole Ion Trap Mass Spectrometer for Infrared Ion Spectroscopy: Application to Protonated Thiated Uridines

Modifications to a Paul-type quadrupole ion trap mass spectrometer providing optical access to the trapped ion cloud as well as hardware and software for coupling to a table-top IR optical parametric oscillator laser (OPO) are detailed. Critical experimental parameters for infrared multiple photon dissociation (IRMPD) on this instrument are characterized. IRMPD action spectra, collected in the hydrogen-stretching region with this instrument, complemented by spectra in the IR fingerprint region acquired at the FELIX facility, are employed to characterize the structures of the protonated forms of 2-thiouridine, [s²Urd+H]⁺, and 4-thiouridine, [s⁴Urd+H]⁺. The measured spectra are compared with predicted linear IR spectra calculated at the B3LYP/6-311+G(d,p) level of theory to determine the conformers populated in the experiments. This comparison indicates that thiation at the 2- or 4-positions shifts the protonation preference between the 2,4-H tautomer and 4-protonation in opposite directions versus canonical uridine, which displays a roughly equal preference for the 2,4-H tautomer and O4 protonation. As found for canonical uridine, protonation leads to a mixture of conformers exhibiting C2'-endo and C3'-endo sugar puckering with an anti nucleobase orientation being populated for both 2- and 4-thiated uridine. Graphical Abstract ᅟ.

Structural and Energetic Effects of O2'-Ribose Methylation of Protonated Purine Nucleosides

The chemical difference between DNA and RNA nucleosides is their 2'-hydrogen versus 2'-hydroxyl substituents. Modification of the ribosyl moiety at the 2'-position and 2'-O-methylation in particular, is common among natural post-transcriptional modifications of RNA. 2'-Modification may alter the electronic properties and hydrogen-bonding characteristics of the nucleoside and thus may lead to enhanced stabilization or malfunction. The structures and relative glycosidic bond stabilities of the protonated forms of the 2'-O-methylated purine nucleosides, 2'-O-methyladenosine (Adom) and 2'-O-methylguanosine (Guom), were examined using two complementary tandem mass spectrometry approaches, infrared multiple photon dissociation action spectroscopy and energy-resolved collision-induced dissociation. Theoretical calculations were also performed to predict the structures and relative stabilities of stable low-energy conformations of the protonated forms of the 2'-O-methylated purine nucleosides and their infrared spectra in the gas phase. Low-energy conformations highly parallel to those found for the protonated forms of the canonical DNA and RNA purine nucleosides are also found for the protonated 2'-O-methylated purine nucleosides. Importantly, the preferred site of protonation, nucleobase orientation, and sugar puckering are preserved among the DNA, RNA, and 2'-O-methylated variants of the protonated purine nucleosides. The 2'-substituent does however influence hydrogen-bond stabilization as the 2'-O-methyl and 2'-hydroxyl substituents enable a hydrogen-bonding interaction between the 2'- and 3'-substituents, whereas a 2'-hydrogen atom does not. Further, 2'-O-methylation reduces the number of stable low-energy hydrogen-bonded conformations possible and importantly inverts the preferred polarity of this interaction versus that of the RNA analogues. Trends in the CID_50% values extracted from survival yield analyses of the 2'-O-methylated and canonical DNA and RNA forms of the protonated purine nucleosides are employed to elucidate their relative glycosidic bond stabilities. The glycosidic bond stability of Adom is found to exceed that of its DNA and RNA analogues. The glycosidic bond stability of Guom is also found to exceed that of its DNA analogue; however, this modification weakens this bond relative to its RNA counterpart. The glycosidic bond stability of the protonated purine nucleosides appears to be correlated with the hydrogen-bond stabilization of the sugar moiety.

Direct Experimental Characterization of Glycosyl Cations by Infrared Ion Spectroscopy

Glycosyl cations are crucial intermediates formed during enzymatic and chemical glycosylation. The intrinsic high reactivity and short lifetime of these reaction intermediates make them very challenging to characterize using spectroscopic techniques. Herein, we report the use of collision induced dissociation tandem mass spectrometry to generate glycosyl cations in the gas phase followed by infrared ion spectroscopy using the FELIX infrared free electron laser. The experimentally observed IR spectra were compared to DFT calculated spectra enabling the detailed structural elucidation of elusive glycosyl oxocarbenium and dioxolenium ions.

The intrinsic basicity of the phosphate backbone exceeds that of uracil and thymine residues: protonation of the phosphate moiety is preferred over the nucleobase for pdThd and pUrd

The gas-phase conformations of the protonated forms of thymidine-5'-monophosphate and uridine-5'-monophosphate, [pdThd+H]⁺ and [pUrd+H]⁺, are investigated by infrared multiple photon dissociation (IRMPD) action spectroscopy and electronic structure calculations. The IRMPD action spectra of [pdThd+H]⁺ and [pUrd+H]⁺ are measured over the IR fingerprint and hydrogen-stretching regions using the FELIX free electron laser and an OPO/OPA laser system. Low-energy conformations of [pdThd+H]⁺ and [pUrd+H]⁺ and their relative stabilities are computed at the MP2(full)/6-311+G(2d,2p)//B3LYP/6-311+G(d,p) and B3LYP/6-311+G(2d,2p)//B3LYP/6-311+G(d,p) levels of theory. Comparisons of the measured IRMPD action spectra and B3LYP/6-311+G(d,p) linear IR spectra computed for the low-energy conformers indicate that the dominant conformers of [pdThd+H]⁺ and [pUrd+H]⁺ populated in the experiments are protonated at the phosphate oxo oxygen atom, with a syn nucleobase orientation that is stabilized by strong P[double bond, length as m-dash]OH⁺O2 and P-OHO4' hydrogen-bonding interactions, and C2'-endo sugar puckering. Minor abundance of conformers protonated at the O2 carbonyl of the nucleobase residue may also contribute for [pdThd+H]⁺, but do not appear to be important for [pUrd+H]⁺. Comparisons to previous IRMPD spectroscopy investigations of the protonated forms of thymidine and uridine, [dThd+H]⁺ and [Urd+H]⁺, and the deprotonated forms of pdThd and pUrd, [pdThd-H]^- and [pUrd-H]^-, provide insight into the effects of the phosphate moiety and protonation on the conformational features of the nucleobase and sugar moieties. Most interestingly, the thymine and uracil nucleobases remain in their canonical forms for [pdThd+H]⁺ and [pUrd+H]⁺, unlike [dThd+H]⁺ and [Urd+H]⁺, where protonation occurs on the nucleobases and induces tautomerization of the thymine and uracil residues.

Influence of Linkage Stereochemistry and Protecting Groups on Glycosidic Bond Stability of Sodium Cationized Glycosyl Phosphates

Energy-resolved collision-induced dissociation (ER-CID) experiments of sodium cationized glycosyl phosphate complexes, [GP _x +Na]⁺, are performed to elucidate the effects of linkage stereochemistry (α versus β), the geometry of the leaving groups (1,2-cis versus 1,2-trans), and protecting groups (cyclic versus non-cyclic) on the stability of the glycosyl phosphate linkage via survival yield analyses. A four parameter logistic dynamic fitting model is used to determine CID_50% values, which correspond to the level of rf excitation required to produce 50% dissociation of the precursor ion complexes. Present results suggest that dissociation of 1,2-trans [GP _x +Na]⁺ occurs via a McLafferty-type rearrangement that is facilitated by a syn orientation of the leaving groups, whereas dissociation of 1,2-cis [GP_x+Na]⁺ is more energetic as it involves the formation of an oxocarbenium ion intermediate. Thus, the C1-C2 configuration plays a major role in determining the stability/reactivity of glycosyl phosphate stereoisomers. For 1,2-cis anomers, the cyclic protecting groups at the C4 and C6 positions stabilize the glycosidic bond, whereas for 1,2-trans anomers, the cyclic protecting groups at the C4 and C6 positions tend to activate the glycosidic bond. The C3 O-benzyl (3 BnO) substituent is key to determining whether the sugar or phosphate moiety retains the sodium cation upon CID. For 1,2-cis anomers, the 3 BnO substituent weakens the glycosidic bond, whereas for 1,2-trans anomers, the 3 BnO substituent stabilizes the glycosidic bond. The C2 O-benzyl substituent does not significantly impact the glycosidic bond stability regardless of its orientation. Graphical abstract ᅟ.

Influence of Transition Metal Cationization versus Sodium Cationization and Protonation on the Gas-Phase Tautomeric Conformations and Stability of Uracil: Application to [Ura+Cu]+ and [Ura+Ag]<sup/>

The gas-phase conformations of transition metal cation-uracil complexes, [Ura+Cu]⁺ and [Ura+Ag]⁺, were examined via infrared multiple photon dissociation (IRMPD) action spectroscopy and theoretical calculations. IRMPD action spectra were measured over the IR fingerprint and hydrogen-stretching regions. Structures and linear IR spectra of the stable tautomeric conformations of these complexes were initially determined at the B3LYP/6-31G(d) level. The four most stable structures computed were also examined at the B3LYP/def2-TZVPPD level to improve the accuracy of the predicted IR spectra. Two very favorable modes of binding are found for [Ura+Cu]⁺ and [Ura+Ag]⁺ that involve O2N3 bidentate binding to the 2-keto-4-hydroxy minor tautomer and O4 monodentate binding to the canonical 2,4-diketo tautomer of Ura. Comparisons between the measured IRMPD and calculated IR spectra enable elucidation of the conformers present in the experiments. These comparisons indicate that both favorable binding modes are represented in the experimental tautomeric conformations of [Ura+Cu]⁺ and [Ura+Ag]⁺. B3LYP suggests that Cu⁺ exhibits a slight preference for O4 binding, whereas Ag⁺ exhibits a slight preference for O2N3 binding. In contrast, MP2 suggests that both Cu⁺ and Ag⁺ exhibit a more significant preference for O2N3 binding. The relative band intensities suggest that O4 binding conformers comprise a larger portion of the population for [Ura+Ag]⁺ than [Ura+Cu]⁺. The dissociation behavior and relative stabilities of the [Ura+M]⁺ complexes, M⁺ = Cu⁺, Ag⁺, H⁺, and Na⁺) are examined via energy-resolved collision-induced dissociation experiments. The IRMPD spectra, dissociation behaviors, and binding preferences of Cu⁺ and Ag⁺ are compared with previous and present results for those of H⁺ and Na⁺. Graphical Abstract ᅟ.

IRMPD Action Spectroscopy, ER-CID Experiments, and Theoretical Studies of Sodium Cationized Thymidine and 5-Methyluridine: Kinetic Trapping During the ESI Desolvation Process Preserves the Solution Structure of [Thd+Na]<sup/>

Thymidine (dThd) is a fundamental building block of DNA nucleic acids, whereas 5-methyluridine (Thd) is a common modified nucleoside found in tRNA. In order to determine the conformations of the sodium cationized thymine nucleosides [dThd+Na]⁺ and [Thd+Na]⁺ produced by electrospray ionization, their infrared multiple photon dissociation (IRMPD) action spectra are measured. Complementary electronic structure calculations are performed to determine the stable low-energy conformations of these complexes. Geometry optimizations and frequency analyses are performed at the B3LYP/6-311+G(d,p) level of theory, whereas energies are calculated at the B3LYP/6-311+G(2d,2p) level of theory. As protonation preferentially stabilizes minor tautomers of dThd and Thd, tautomerization facilitated by Na⁺ binding is also considered. Comparisons of the measured IRMPD and computed IR spectra find that [dThd+Na]⁺ prefers tridentate (O2,O4',O5') coordination to the canonical 2,4-diketo form of dThd with thymine in a syn orientation. In contrast, [Thd+Na]⁺ prefers bidentate (O2,O2') coordination to the canonical 2,4-diketo tautomer of Thd with thymine in an anti orientation. Although 2,4-dihydroxy tautomers and O2 protonated thymine nucleosides coexist in the gas phase, no evidence for minor tautomers is observed for the sodium cationized species. Consistent with experimental observations, the computational results confirm that the sodium cationized thymine nucleosides exhibit a strong preference for the canonical form of the thymine nucleobase. Survival yield analyses based on energy-resolved collision-induced dissociation (ER-CID) experiments suggest that the relative stabilities of protonated and sodium cationized dThd and Thd follow the order [dThd+H]⁺ < [Thd+H]⁺ < [dThd+Na]⁺ < [Thd+Na]⁺. Graphical Abstract ᅟ.

Probing the geometry reorganization from solution to gas-phase in putrescine derivatives by IRMPD, 1H-NMR and theoretical calculations

Geometry reorganization of ESI formed ions are demonstrated and explicit calculations of the solution phase are shown to be relevant.