Structure of sodiated octa-glycine: IRMPD spectroscopy and molecular modeling

Perspective: intrinsic interactions of metal ions with biological molecules as studied by threshold collision-induced dissociation and infrared multiple photon dissociation

In this perspective, gas-phase studies of group 1 monocations and group 12 dications with amino acids and small peptides are highlighted. Although the focus is on two experimental techniques, threshold collision-induced dissociation and infrared multiple photon dissociation action spectroscopy, these methods as well as complementary approaches are summarized. The synergistic interplay with theory, made particularly powerful by the small sizes of the systems explored and the absence of solvent and support, is also elucidated. Importantly, these gas-phase methods permit quantitative insight into the structures and thermodynamics of metal cations interacting with biological molecules. Periodic trends in how these interactions vary as the metal cations get heavier are discussed as are quantitative trends with changes in the amino acid side chain and effects of hydration. Such trends allow these results to transcend the limitations associated with the biomimetic model systems.

The impact of the electric field of metal ions on the vibrations and internal hydrogen bond strength in alkali metal ion di- and triglycine complexes

Using infrared predissociation spectroscopy of cryogenic ions, we revisit the vibrational spectra of alkali metal ion (Li+, Na+, K+) di- and triglycine complexes. We assign their most stable conformation, which involves metal ion coordination to all C=O groups and an internal NH⋯NH2 hydrogen bond in the peptide backbone. An analysis of the spectral shifts of the OH and C=O stretching vibrations across the different metal ions and peptide chain lengths shows that these are largely caused by the electric field of the metal ion, which varies in strength as a function of the square of the distance. The metal ion–peptide interaction also remotely modulates the strength of internal hydrogen bonding in the peptide backbone via the weakening of the amide C=O bond, resulting in a decrease in internal hydrogen bond strength from Li+ > Na+ > K+.

Experimental and computational studies of monovalent metal cation-peptide interactions in [M·GlyGlyHis]+ (M = Li, Na, K, Rb, Cs, and Ag) complexes

[M·GlyGlyHis]⁺ (M = Li, Na, K, Rb, Cs, and Ag) complexes were generated using the electrospray ionization method and were subjected to collision-induced dissociation. Metal ion loss is the primary dissociation channel for [Cs·GlyGlyHis]⁺ whereas other complexes yield metal-bound peptide sequence ions and dehydrated ions as the main products. [Li·GlyGlyHis]⁺ and [Ag·GlyGlyHis]⁺ also generate product ions that are not observed for other complexes. Density functional theory calculations suggest that metal ion-peptide ligand interaction occurs through covalent interactions in [Li·GlyGlyHis]⁺ and [Ag·GlyGlyHis]⁺, and through electrostatic attraction in [Na·GlyGlyHis]⁺, [K·GlyGlyHis]⁺, [Rb·GlyGlyHis]⁺, and [Cs·GlyGlyHis]⁺. The calculations also suggest that fragmentation behavior of these complexes is affected by charge transfer to the ligand and ion-ligand interaction energy, and to a lesser extent by the ion size.

Gas-phase metal ion chelation investigated with IRMPD spectroscopy: A brief review of Robert Dunbar's contributions

With the passing of Prof. Robert C. Dunbar on 31 October 2017, the field of ion chemistry lost one of its modern heroes. Throughout his career in mass spectrometry, two of his main research interests involved the interaction of trapped ions with electromagnetic radiation and the chelation motifs of metal ions with organic ligands. The focus of his early career was on the fundamental processes that take place in molecules upon ultraviolet and infrared excitation. From 2003 to 2017, his scientific interests shifted to more structural questions, notably to resolving the structures and binding motifs of metal ion chelation complexes by application of infrared photodissociation spectroscopy. These experiments were carried out during numerous visits to the (Free Electron Laser for Infrared eXperiments) (FELIX) facility in the Netherlands and were complemented by extensive theoretical investigations by Rob. As a tribute to our friend, we present in this contribution a brief review of this work.

Replacing H+ by Na+ or K+ in phosphopeptide anions and cations prevents electron capture dissociation

By successively replacing H+ by Na+ or K+ in phosphopeptide anions and cations, we show that the efficiency of fragmentation into c and z˙ or c˙ and z fragments from N-Cα backbone bond cleavage by negative ion electron capture dissociation (niECD) and electron capture dissociation (ECD) substantially decreases with increasing number of alkali ions attached. In proton-deficient phosphopeptide ions with a net charge of 2-, we observed an exponential decrease in electron capture efficiency with increasing number of Na+ or K+ ions attached, suggesting that electrons are preferentially captured at protonated sites. In proton-abundant phosphopeptide ions with a net charge of 3+, the electron capture efficiency was not affected by replacing up to four H+ ions with Na+ or K+ ions, but the yield of c, z˙ and c˙, z fragments from N-Cα backbone bond cleavage generally decreased next to Na+ or K+ binding sites. We interpret the site-specific decrease in fragmentation efficiency as Na+ or K+ binding to backbone amide oxygen in competition with interactions of protonated sites that would otherwise lead to backbone cleavage into c, z˙ or c˙, z fragments. Our findings seriously challenge the hypothesis that the positive charge responsible for ECD into c, z˙ or c˙, z fragments can generally be a sodium or other metal ion instead of a proton.

Unfolding and Folding of the Three-Helix Bundle Protein KIX in the Absence of Solvent

Electron capture dissociation was used to probe the structure, unfolding, and folding of KIX ions in the gas phase. At energies for vibrational activation that were sufficiently high to cause loss of small molecules such as NH3 and H2O by breaking of covalent bonds in about 5% of the KIX (M + nH)(n+) ions with n = 7-9, only partial unfolding was observed, consistent with our previous hypothesis that salt bridges play an important role in stabilizing the native solution fold after transfer into the gas phase. Folding of the partially unfolded ions on a timescale of up to 10 s was observed only for (M + nH)(n+) ions with n = 9, but not n = 7 and n = 8, which we attribute to differences in the distribution of charges within the (M + nH)(n+) ions. Graphical Abstract ᅟ.

Complexes of Ni(ii) and Cu(ii) with small peptides: deciding whether to deprotonate

Infrared multiple photon dissociation (IRMPD) spectroscopy differentiates two binding modes (iminol versus charge solvated) for Ni(ii) bound to model peptides.

Chirality-dependent structuration of protonated or sodiated polyphenylalanines: IRMPD and ion mobility studies

Ion mobility experiments are combined with Infra-Red Multiple Photon Dissociation (IRMPD) spectroscopy and quantum chemical calculations for assessing the role of chirality in the structure of protonated and sodiated di- or tetra-peptides. Sodiated systems show a strong chirality dependence of the competition between Na(+)O and Na(+)π interactions. Chirality effects are more subtle in protonated systems and manifest themselves by differences in the secondary interactions such hydrogen bonds between neutral groups or those involving the aromatic rings.

Divalent Metal-Ion Complexes with Dipeptide Ligands Having Phe and His Side-Chain Anchors: Effects of Sequence, Metal Ion, and Anchor

Conformational preferences have been surveyed for divalent metal cation complexes with the dipeptide ligands AlaPhe, PheAla, GlyHis, and HisGly. Density functional theory results for a full set of complexes are presented, and previous experimental infrared spectra, supplemented by a number of newly recorded spectra obtained with infrared multiple photon dissociation spectroscopy, provide experimental verification of the preferred conformations in most cases. The overall structural features of these complexes are shown, and attention is given to comparisons involving peptide sequence, nature of the metal ion, and nature of the side-chain anchor. A regular progression is observed as a function of binding strength, whereby the weakly binding metal ions (Ba(2+) to Ca(2+)) transition from carboxylate zwitterion (ZW) binding to charge-solvated (CS) binding, while the stronger binding metal ions (Ca(2+) to Mg(2+) to Ni(2+)) transition from CS binding to metal-ion-backbone binding (Iminol) by direct metal-nitrogen bonds to the deprotonated amide nitrogens. Two new sequence-dependent reversals are found between ZW and CS binding modes, such that Ba(2+) and Ca(2+) prefer ZW binding in the GlyHis case but prefer CS binding in the HisGly case. The overall binding strength for a given metal ion is not strongly dependent on the sequence, but the histidine peptides are significantly more strongly bound (by 50-100 kJ mol(-1)) than the phenylalanine peptides.

The protonated and sodiated dimers of proline studied by IRMPD spectroscopy in the N-H and O-H stretching region and computational methods

IRMPD spectroscopy and computational chemistry techniques have been used to determine that the proton- and sodium-bound dimers of proline exist as a mixture of a number of different structures. Simulated annealing computations were found to be helpful in determining the unique structures of the protonated and sodiated dimers, augmenting chemical intuition. The experimental and computational results are consistent with the proton-bound dimer of N-protonated proline bound to zwitterionic proline. There was no spectroscopic evidence in the 3200-3800 cm(-1) region for a canonical structure which is predicted to have a weak N-H stretch at about 3440 cm(-1). A well resolved band at 1733 cm(-1) from a previous spectroscopic study (DOI: 10.1021/ja068715a ) was reassigned from a high energy canonical isomer to the C=O stretch of a lower energy zwitterionic structure. This band is a free carboxylate C=O stretch where protonated proline is hydrogen bonded to the other carboxylate oxygen which is also involved in an intramolecular hydrogen bond. Fifteen structures of the sodium bound proline dimer were computed to be within 10 kJ mol(-1) of Gibbs energy and eight structures were within 5 kJ mol(-1). None of these structures can be ruled out based on the experimental IRMPD spectrum. They all have an N-H stretching band predicted in a position that agrees with the experimental spectrum. However, only structures where one of the proline monomers is in the canonical form and having a free O-H bond can produce the band at ∼3600 cm(-1).

Unraveling non-covalent interactions within flexible biomolecules: from electron density topology to gas phase spectroscopy

The NCI (Non-Covalent Interactions) method, a recently-developed theoretical strategy to visualize weak non-covalent interactions from the topological analysis of the electron density and of its reduced gradient, is applied in the present paper to document intra- and inter-molecular interactions in flexible molecules and systems of biological interest in combination with IR spectroscopy. We first describe the conditions of application of the NCI method to the specific case of intramolecular interactions. Then we apply it to a series of stable conformations of isolated molecules as an interpretative technique to decipher the different physical interactions at play in these systems. Examples are chosen among neutral molecular systems exhibiting a large diversity of interactions, for which an extensive spectroscopic characterization under gas-phase isolation conditions has been obtained using state-of-the-art conformer-specific experimental techniques. The interactions presently documented range from weak intra-molecular H-bonds in simple amino-alcohols, to more complex patterns, with interactions of various strengths in model peptides, as well as in chiral bimolecular systems, where invaluable hints for the understanding of chiral recognition are revealed. We also provide a detailed technical appendix, which discusses the choices of cut-offs as well as the applicability of the NCI analysis to specific constrained systems, where local effects require attention. Finally, the NCI technique provides IR spectroscopists with an elegant visualization of the interactions that potentially impact their vibrational probes, namely the OH and NH stretching motions. This contribution illustrates the power and the conditions of use of the NCI technique, with the aim of providing an easy tool for all chemists, experimentalists and theoreticians, for the visualization and characterization of the interactions shaping complex molecular systems.

Probing the mechanisms and dynamics of gas phase hydrogen-deuterium exchange reactions of sodiated polyglycines

The rate constants for H-D exchange reactions of sodiated polyglycines (GnNa(+), n = 2-8) and polyalanines (AnNa(+), n = 2, 3 and 5) with ND3 have been measured in the cell of an FT-ICR mass spectrometer. All peptides except G2Na(+) are found to undergo three exchange reactions, all of which are consecutive with no sign of multiple exchanges within a single collision event. This information has been used to construct full mechanistic scenarios with the help of detailed quantum chemical calculations of the possible reaction paths for H-D exchange. The first exchange is always located at the C terminus however with different mechanisms depending upon whether the peptide termini can (larger peptides) or cannot (smaller peptides) interact directly without strong energy penalty. The most favourable mechanisms for the second and third exchanges of the N terminus protons, are found to be different from those for the first for all peptide sizes. The peptide distortions that are necessary in order for some of these reactions to occur are made possible by the energy reservoir provided by the favorable interaction of the peptide ion with ND3. Their occurrence and variety preclude any general relationship between H-D exchange kinetics and the most stable ion structures. There is however a break at G7Na(+) in the kinetics trend, with a first exchange rate which is much smaller than for all other peptide sizes. This break can be directly related to a different structural type in which the C terminus is neither free nor close to the N terminus.

Structure and further fragmentation of significant [a3 + Na - H]+ ions from sodium-cationized peptides

A good understanding of gas-phase fragmentation chemistry of peptides is important for accurate protein identification. Additional product ions obtained by sodiated peptides can provide useful sequence information supplementary to protonated peptides and improve protein identification. In this work, we first demonstrate that the sodiated a3 ions are abundant in the tandem mass spectra of sodium-cationized peptides although observations of a3 ions have rarely been reported in protonated peptides. Quantum chemical calculations combined with tandem mass spectrometry are used to investigate this phenomenon by using a model tetrapeptide GGAG. Our results reveal that the most stable [a3 + Na - H](+) ion is present as a bidentate linear structure in which the sodium cation coordinates to the two backbone carbonyl oxygen atoms. Due to structural inflexibility, further fragmentation of the [a3 + Na - H](+) ion needs to overcome several relatively high energetic barriers to form [b2 + Na - H](+) ion with a diketopiperazine structure. As a result, low abundance of [b2 + Na - H](+) ion is detected at relatively high collision energy. In addition, our computational data also indicate that the common oxazolone pathway to generate [b2 + Na - H](+) from the [a3 + Na - H](+) ion is unlikely. The present work provides a mechanistic insight into how a sodium ion affects the fragmentation behaviors of peptides.

Structure of the Pb²⁺-deprotonated dGMP complex in the gas phase: a combined MS-MS/IRMPD spectroscopy/ion mobility study

The structure of the Pb(2+)-deprotonated 2'-deoxyguanosine-5'-monophosphate (dGMP) complex, generated in the gas phase by electrospray ionization, was examined by combining tandem mass spectrometry, mid-infrared multiple-photon dissociation (IRMPD) spectroscopy and ion mobility. In the gas phase, the main binding site of Pb(2+) onto deprotonated dGMP is the deprotonated phosphate group, but the question is whether an additional stabilization of the metallic complex can occur via participation of the carbonyl group of guanine. Such macrochelates indeed correspond to the most stable structures according to theoretical calculations. A multiplexed experimental approach was used to characterize the gas-phase conformation of the metallic complex and hence determine the binding mode of Pb(2+) with [dGMP](-). MS/MS analysis, observation of characteristic bands by IRMPD spectroscopy, and measurement of the ion mobility collision cross section suggest that gaseous [Pb(dGMP)-H](+) complexes adopt a macrochelate folded structure, which consequently differs strongly from the zwitterionic forms postulated in solution from potentiometric studies.

Structural investigation of protonated azidothymidine and protonated dimer

Infrared multiple photon dissociation (IRMPD) spectroscopy experiments and quantum chemical calculations have been used to explore the possible structures of protonated azidothymidine and the corresponding protonated dimer. Many interesting differences between the protonated and neutral forms of azidothymidine were found, particularly associated with keto-enol tautomerization. Comparison of computational vibrational and the experimental IMRPD spectra show good agreement and give confidence that the dominant protonated species has been identified. The protonated dimer of azidothymidine exhibits three intramolecular hydrogen bonds. The IRMPD spectrum of the protonated dimer is consistent with the spectrum of the most stable computational structure. This work brings to light interesting keto-enol tautomerization and exocyclic hydrogen bonding involving azidothymidine and its protonated dimer. The fact that one dominant protonated species is observed in the gas phase, despite both the keto and enol structures being similar in energy, is proposed to be the direct result of the electrospray ionization process in which the dominant protonated dimer structure dissociates in the most energetically favorable way.

Consecutive fragmentation mechanisms of protonated ferulic acid probed by infrared multiple photon dissociation spectroscopy and electronic structure calculations

Protonated ferulic acid and its principle fragment ion have been characterized using infrared multiple photon dissociation spectroscopy and electronic structure calculations at the B3LYP/6-311 + G(d,p) level of theory. Due to its extensively conjugated structure, protonated ferulic acid is observed to yield three stable fragment ions in IRMPD experiments. It is proposed that two parallel fragmentation pathways of protonated ferulic acid are being observed. The first pathway involves proton transfer, resulting in the loss of water and subsequently carbon monoxide, producing fragment ions m/z 177 and 149, respectively. Optimization of m/z 177 yields a species containing an acylium group, which is supported by a diagnostic peak in the IRMPD spectrum at 2168 cm(-1). The second pathway involves an alternate proton transfer leading to loss of methanol and rearrangement to a five-membered ring.

Infrared multiple-photon dissociation spectroscopy of tripositive ions: lanthanum-tryptophan complexes

Collision-induced charge disproportionation limits the stability of triply charged metal ion complexes and has thus far prevented successful acquisition of their gas-phase IR spectra. This has curtailed our understanding of the structures of triply charged metal complexes in the gas phase and in biological environments. Herein we report the first gas-phase IR spectra of triply charged La(III) complexes with a derivative of tryptophan (N-acetyl tryptophan methyl ester), and an unusual dissociation product, a lanthanum amidate. These spectra are compared with those predicted using density functional theory. The best structures are those of the lowest energies that differ by details in the π-interaction between La(3+) and the indole rings. Other binding sites on the tryptophan derivative are the carbonyl oxygens. In the lanthanum amidate, La(3+) replaces an H(+) in the amide bond of the tryptophan derivative.

Determination of noncovalent docking by infrared spectroscopy of cold gas-phase complexes

Multidentate, noncovalent interactions between small molecules and biopolymer fragments are central to processes ranging from drug action to selective catalysis. We present a versatile and sensitive spectroscopic probe of functional groups engaged in hydrogen bonding in such contexts. This involves measurement of the frequency changes in specific covalent bonds upon complex formation, information drawn from otherwise transient complexes that have been extracted from solution and conformationally frozen near 10 kelvin in gas-phase clusters. Resonances closely associated with individual oscillators are easily identified through site-specific isotopic labeling, as demonstrated by application of the method to an archetypal system involving a synthetic tripeptide known to bind biaryl substrates through tailored hydrogen bonding to catalyze their asymmetric bromination. With such data, calculations readily converge on the plausible operative structures in otherwise computationally prohibitive, high-dimensionality landscapes.

Experimental and theoretical investigation of the proton-bound dimer of lysine

The structure of the proton-bound lysine dimer has been investigated by infrared multiple photon dissociation (IRMPD) spectroscopy and electronic structure calculations. The structures of different possible isomers of the proton-bound lysine dimer have been optimized at the B3LYP/6-31 + G(d) level of theory and IR spectra calculated using the same computational method. Based on relative Gibbs free energies (298 K) calculated at the MP2/aug-cc-pVTZ//B3LYP/6-31 + G(d) level of theory, LL-CS01, and followed closely (1.1 kJ mol(-1)) by LL-CS02 are the most stable non-zwitterionic isomers. At the MP2/aug-cc-pVTZ//6-31 + G(d) and MP2/aug-cc-pVTZ//6-31 + (d,p) levels of theory, isomer LL-CS02 is favored by 3.0 and 2.3 kJ mol(-1), respectively. The relative Gibbs free energies calculated by the aforementioned levels of theory for LL-CS01 and LL-CS02 are very close and strongly suggest that diagnostic vibrational signatures found in the IRMPD spectrum of the proton-bound dimer of lysine can be attributed to the existence of both isomers. LL-ZW01 is the most stable zwitterionic isomer, in which the zwitterionic structure of the neutral lysine is well stabilized by the protonated lysine moiety via a very strong intermolecular hydrogen bond. At the MP2/aug-cc-pVTZ//B3LYP/6-31 + G(d), MP2/aug-cc-pVTZ//6-31 + G(d) and MP2/aug-cc-pVTZ//6-31 + G(d,p) levels of theory, the most stable zwitterionic isomer (LL-ZW01) is less favored than LL-CS01 by 7.3, 4.1 and 2.3 kJ mol(-1), respectively. The experimental IRMPD spectrum also confirms that the proton-bound dimer of lysine largely exists as charge-solvated isomers. Investigation of zwitterionic and charge-solvated species of amino acids in the gas phase will aid in a further understanding of structure, property, and function of biological molecules.