Insights from first principles molecular dynamics studies toward infrared multiple-photon and single-photon action spectroscopy: Case study of the proton-bound dimethyl ether dimer

Resource Optimization for Quantum Dynamics with Tensor Networks: Quantum and Classical Algorithms

The exponential scaling of the quantum degrees of freedom with the size of the system is one of the biggest challenges in computational chemistry and particularly in quantum dynamics. We present a tensor network approach for the time-evolution of the nuclear degrees of freedom of multiconfigurational chemical systems at a reduced storage and computational complexity. We also present quantum algorithms for the resultant dynamics. To preserve the compression advantage achieved via tensor network decompositions, we present an adaptive algorithm for the regularization of nonphysical bond dimensions, preventing the potentially exponential growth of these with time. While applicable to any quantum dynamical problem, our method is particularly valuable for dynamical simulations of nuclear chemical systems. Our algorithm is demonstrated using ab initio potentials obtained for a symmetric hydrogen-bonded system, namely, the protonated 2,2'-bipyridine, and compared to exact diagonalization numerical results.

A spectroscopic test suggests that fragment ion structure annotations in MS/MS libraries are frequently incorrect

Modern untargeted mass spectrometry (MS) analyses quickly detect and resolve thousands of molecular compounds. Although features are readily annotated with a molecular formula in high-resolution small-molecule MS applications, the large majority of them remains unidentified in terms of their full molecular structure. Collision-induced dissociation tandem mass spectrometry (CID-MS2) provides a diagnostic molecular fingerprint to resolve the molecular structure through a library search. However, for de novo identifications, one must often rely on in silico generated MS2 spectra as reference. The ability of different in silico algorithms to correctly predict MS2 spectra and thus to retrieve correct molecular structures is a topic of lively debate, for instance in the CASMI contest. Underlying the predicted MS2 spectra are the in silico generated product ion structures, which are normally not used in de novo identification, but which can serve to critically assess the fragmentation algorithms. Here we evaluate in silico generated MSn product ion structures by comparison with structures established experimentally by infrared ion spectroscopy (IRIS). For a set of three dozen product ion structures from five precursor molecules, we find that virtually all fragment ion structure annotations in three major in silico MS2 libraries (HMDB, METLIN, mzCloud) are incorrect and caution the reader against their use for structure annotation of MS/MS ions.

Correlated proton dynamics in hydrogen bonding networks: the benchmark case of 3-hydroxyglutaric acid

Proton and hydrogen-bonded networks sustain a broad range of structural and charge transfer processes in supramolecular materials. The modelling of proton dynamics is however challenging and demands insights from prototypical benchmark systems. The intramolecular H-bonding networks induced by either protonation or deprotonation of 3-hydroxyglutaric acid provide intriguing case studies of correlated proton dynamics. The vibrational signatures associated with the fluxional proton bonding and its coupling with the hydroxyglutaric backbone are investigated here with infrared action ion spectroscopy experiments and Born-Oppenheimer molecular dynamics (BOMD) computations. Despite the formally similar symmetry of protonated and deprotonated hydroxyglutaric acid, the relative proton affinities of the oxygen centers of the carboxylic and carboxylate groups with respect to that of the central hydroxyl group lead to distinct proton dynamics. In the protonated acid, a tautomeric arrangement of the type HOCO·[HOH]+·OCOH is preferred with the proton binding tighter to the central hydroxyl moiety and the electronic density being shared between the two nearly symmetric H-bonds with the carboxylic end groups. In the deprotonated acid, the asymmetric [OCO]-·HO·HOCO configuration is more stable, with a stronger H-bonding on the bare carboxylate end. Both systems display active backbone dynamics and concerted Grothuss-like proton motions, leading to diffuse band structures in their vibrational spectra. These features are accurately reproduced by the BOMD computations.

A Dynamic Proton Bond: MH+·H2O ⇌ M·H3O+ Interconversion in Loosely Coordinated Environments

The interaction of organic molecules with oxonium cations within their solvation shell may lead to the emergence of dynamic supramolecular structures with recurrently changing host-guest chemical identity. We illustrate this phenomenon in benchmark proton-bonded complexes of water with polyether macrocyles. Despite the smaller proton affinity of water versus the ether group, water in fact retains the proton in the form of H₃O⁺, with increasing stability as the coordination number increases. Hindrance in many-fold coordination induces dynamic reversible (ether)·H₃O⁺ ⇌ (etherH⁺)·H₂O interconversion. We perform infrared action ion spectroscopy over a broad spectral range to expose the vibrational signatures of the loose proton bonding in these systems. Remarkably, characteristic bands for the two limiting proton bonding configurations are observed in the experimental vibrational spectra, superimposed onto diffuse bands associated with proton delocalization. These features cannot be described by static equilibrium structures but are accurately modeled within the framework of ab initio molecular dynamics.

Graph Theoretic Molecular Fragmentation for Multidimensional Potential Energy Surfaces Yield an Adaptive and General Transfer Machine Learning Protocol

Over a series of publications we have introduced a graph-theoretic description for molecular fragmentation. Here, a system is divided into a set of nodes, or vertices, that are then connected through edges, faces, and higher-order simplexes to represent a collection of spatially overlapping and locally interacting subsystems. Each such subsystem is treated at two levels of electronic structure theory, and the result is used to construct many-body expansions that are then embedded within an ONIOM-scheme. These expansions converge rapidly with many-body order (or graphical rank) of subsystems and have been previously used for ab initio molecular dynamics (AIMD) calculations and for computing multidimensional potential energy surfaces. Specifically, in all these cases we have shown that CCSD and MP2 level AIMD trajectories and potential surfaces may be obtained at density functional theory cost. The approach has been demonstrated for gas-phase studies, for condensed phase electronic structure, and also for basis set extrapolation-based AIMD. Recently, this approach has also been used to derive new quantum-computing algorithms that enormously reduce the quantum circuit depth in a circuit-based computation of correlated electronic structure. In this publication, we introduce (a) a family of neural networks that act in parallel to represent, efficiently, the post-Hartree-Fock electronic structure energy contributions for all simplexes (fragments), and (b) a new k-means-based tessellation strategy to glean training data for high-dimensional molecular spaces and minimize the extent of training needed to construct this family of neural networks. The approach is particularly useful when coupled cluster accuracy is desired and when fragment sizes grow in order to capture nonlocal interactions accurately. The unique multidimensional k-means tessellation/clustering algorithm used to determine our training data for all fragments is shown to be extremely efficient and reduces the needed training to only 10% of data for all fragments to obtain accurate neural networks for each fragment. These fully connected dense neural networks are then used to extrapolate the potential energy surface for all molecular fragments, and these are then combined as per our graph-theoretic procedure to transfer the learning process to a full system energy for the entire AIMD trajectory at less than one-tenth the cost as compared to a regular fragmentation-based AIMD calculation.

Mapping Quantum Chemical Dynamics Problems to Spin-Lattice Simulators

The accurate computational determination of chemical, materials, biological, and atmospheric properties has a critical impact on a wide range of health and environmental problems, but is deeply limited by the computational scaling of quantum mechanical methods. The complexity of quantum chemical studies arises from the steep algebraic scaling of electron correlation methods and the exponential scaling in studying nuclear dynamics and molecular flexibility. To date, efforts to apply quantum hardware to such quantum chemistry problems have focused primarily on electron correlation. Here, we provide a framework that allows for the solution of quantum chemical nuclear dynamics by mapping these to quantum spin-lattice simulators. Using the example case of a short-strong hydrogen-bonded system, we construct the Hamiltonian for the nuclear degrees of freedom on a single Born-Oppenheimer surface and show how it can be transformed to a generalized Ising model Hamiltonian. We then demonstrate a method to determine the local fields and spin-spin couplings needed to identically match the molecular and spin-lattice Hamiltonians. We describe a protocol to determine the on-site and intersite coupling parameters of this Ising Hamiltonian from the Born-Oppenheimer potential and nuclear kinetic energy operator. Our approach represents a paradigm shift in the methods used to study quantum nuclear dynamics, opening the possibility to solve both electronic structure and nuclear dynamics problems using quantum computing systems.

Proton in the ring: spectroscopy and dynamics of proton bonding in macrocycle cavities

The proton bond is a paradigmatic quantum molecular interaction and a major driving force of supramolecular chemistry. The ring cavities of crown ethers provide an intriguing environment, promoting competitive proton sharing with multiple coordination anchors. This study shows that protons confined in crown ether cavities form dynamic bonds that migrate to varying pairs of coordinating atoms when allowed by the flexibility of the macrocycle backbone. Prototypic native crown ethers (12-crown-4, 15-crown-5 and 18-crown-6) and aza-crown ethers (cyclen, 1-aza-18-crown-6 and hexacyclen) are investigated. For each system, Infrared action spectroscopy experiments and ab initio Molecular Dynamics computations are employed to elucidate the structural effects associated with proton diffusion and its entanglement with the conformational and vibrational dynamics of the protonated host.

Weighted-Graph-Theoretic Methods for Many-Body Corrections within ONIOM: Smooth AIMD and the Role of High-Order Many-Body Terms

We present a weighted-graph-theoretic approach to adaptively compute contributions from many-body approximations for smooth and accurate post-Hartree-Fock (pHF) ab initio molecular dynamics (AIMD) of highly fluxional chemical systems. This approach is ONIOM-like, where the full system is treated at a computationally feasible quality of treatment (density functional theory (DFT) for the size of systems considered in this publication), which is then improved through a perturbative correction that captures local many-body interactions up to a certain order within a higher level of theory (post-Hartree-Fock in this publication) described through graph-theoretic techniques. Due to the fluxional and dynamical nature of the systems studied here, these graphical representations evolve during dynamics. As a result, energetic "hops" appear as the graphical representation deforms with the evolution of the chemical and physical properties of the system. In this paper, we introduce dynamically weighted, linear combinations of graphs, where the transition between graphical representations is smoothly achieved by considering a range of neighboring graphical representations at a given instant during dynamics. We compare these trajectories with those obtained from a set of trajectories where the range of local many-body interactions considered is increased, sometimes to the maximum available limit, which yields conservative trajectories as the order of interactions is increased. The weighted-graph approach presents improved dynamics trajectories while only using lower-order many-body interaction terms. The methods are compared by computing dynamical properties through time-correlation functions and structural distribution functions. In all cases, the weighted-graph approach provides accurate results at a lower cost.

A theoretical study on the infrared signatures of proton-bound rare gas dimers (Rg-H+-Rg), Rg = {Ne, Ar, Kr, and Xe}

The infrared spectrum of proton-bound rare gas dimers has been extensively studied via matrix isolation spectroscopy. However, little attention has been paid on their spectrum in the gas phase. Most of the Rg₂H⁺ has not been detected outside the matrix environment. Recently, Ar_nH⁺ (n = 3-7) has been first detected in the gas-phase [D. C. McDonald et al., J. Chem. Phys. 145, 231101 (2016)]. In that work, anharmonic theory can reproduce the observed vibrational structure. In this paper, we extend the existing theory to examine the vibrational signatures of Rg₂H⁺, Rg = {Ne, Ar, Kr, and Xe}. The successive binding of Rg to H⁺ was investigated through the calculation of stepwise formation energies. It was found that this binding is anti-cooperative. High-level full-dimensional potential energy surfaces at the CCSD(T)/aug-cc-pVQZ//MP2/aug-cc-pVQZ were constructed and used in the anharmonic calculation via discrete variable representation. We found that the potential coupling between the symmetric and asymmetric Rg-H⁺ stretch (ν₁ and ν₃ respectively) causes a series of bright n₁ν₁ + ν₃ progressions. From Ne₂H⁺ to Xe₂H⁺, an enhancement of intensities for these bands was observed.

Intra-cavity proton bonding and anharmonicity in the anionophore cyclen

Proton bonding drives the supramolecular chemistry of a broad range of materials with polar moieties. Proton delocalization and electronic charge redistribution have a profound impact on the structure of proton-bound molecular frameworks, and pose fundamental challenges to quantum chemical modelling. This study provides insights into the structural and spectral signatures of the intramolecular proton bond formed in a benchmark polyazamacrocycle anionophore (cyclen, 1,4,7,10-tetraazacyclododecane). Infrared action spectroscopy is employed to characterize the macrocycle, isolated in protonated form. In its most stable configuration, protonated cyclen adopts an open arrangement of Cs symmetry with a particularly strong NHδ+N bond across the cavity. The quantum chemical analysis of the infrared spectrum reveals intrinsic difficulties for the accurate description of the vibrational modes of the system. The reconciliation of the computational predictions with experiment demands a careful anharmonic treatment of the proton motion, which exposes the limitations of current methods. Best results are obtained with the incorporation of anharmonicity only to the fundamental modes directly related to motions of the proton. However, the full anharmonic treatment of the system fails to describe correctly the vibrations related to the macrocycle backbone. The results should serve as motivation for new developments in the modelling of proton bonded systems.

Proton relays in anomalous carbocations dictate spectroscopy, stability, and mechanisms: case studies on C2H5+ and C3H3+

We present a detailed exposition of “Grothuss-like” proton shuttles in C₂H₅⁺ and C₃H₃⁺ that result in anomalous structural and spectral behavior.

A Grotthuss-like proton shuttle in the anomalous C2H3+ carbocation: energetic and vibrational properties for isotopologues

The proton shuttle in C₂H₃⁺ is reminiscent of the Grotthuss proton transfer mechanism.

Active Site Dynamical Effects in the Hydrogen Transfer Rate-limiting Step in the Catalysis of Linoleic Acid by Soybean Lipoxygenase-1 (SLO-1): Primary and Secondary Isotope Contributions

Using ab initio molecular dynamics (AIMD) simulations that facilitate the treatment of rare events, we probe the active site participation in the rate-determining hydrogen transfer step in the catalytic oxidation of linoleic acid by soybean lipoxygenase-1 (SLO-1). The role of two different active site components is probed. (a) On the hydrogen atom acceptor side of the active site, the hydrogen bonding propensity between the acceptor side hydroxyl group, which is bound to the iron cofactor, and the backbone carboxyl group of isoleucine (residue number 839) is studied toward its role in promoting the hydrogen transfer event. Primary and secondary (H/D) isotope effects are also probed and a definite correlation with subtle secondary H/D isotope effects is found. With increasing average nuclear kinetic energy, the increase in transfer probability is enhanced due to the presence of the hydrogen bond between the backbone carbonyl of I839 and the acceptor oxygen. Further increase in average nuclear kinetic energy reduces the strength of this secondary hydrogen bond which leads to a deterioration in hydrogen transfer rates and finally embrances an Arrhenius-like behavior. (b) On the hydrogen atom donor side, the coupling between vibrational modes predominantly localized on the donor-side linoleic acid group and the reactive mode is probed. There appears to be a qualitative difference in the coupling between modes that belong to linoleic acid and the hydrogen transfer mode, for hydrogen and deuterium transfer. For example, the donor side secondary hydrogen atom is much more labile (by nearly a factor of 5) during deuterium transfer as compared to the case for hydrogen transfer. This appears to indicate a greater coupling between the modes belonging to the linoleic acid scaffold and the deuterium transfer mode and also provides a new rationalization for the abnormal (nonclassical) secondary isotope effect results obtained by Knapp, Rickert, and Klinman in J. Am. Chem. Soc. , 2002 , 124 , 3865 . To substantiate our findings noted in point a above, we have suggested an I839 → A839 or I839 → V839 mutation. This will modify the bulkiness of hydrogen the bonding residue, allowing greater flexibility in the secondary hydrogen bond formation highlighted above and adversely affecting the reaction rate.

Proton-bound dimers of 1-methylcytosine and its derivatives: vibrational and NMR spectroscopy

Vibrational spectroscopy and NMR demonstrate that the proton-bound dimer of 1-methylcytosine, 1, has an unsymmetrical structure at room temperature. In the gas phase, investigation of isolated homodimer 1 reveals five fundamental NH vibrations by IR Multiple Photon Dissociation (IRMPD) action spectroscopy. The NH···N stretching vibration between the two ring nitrogens exhibits a frequency of 1570 cm(-1), as confirmed by examination of the proton-bound homodimers of 5-fluoro-1-methycytosine, 2, and of 1,5-dimethylcytosine, 3, which display absorptions in the same region that disappear upon deuterium substitution. (13)C, and (15)N NMR of the solid iodide salt of 1 confirm the nonequivalence of the two rings in the anhydrous proton-bound homodimer at room temperature. IRMPD spectra of the three possible heterodimers also show NH···N stretches in the same domain, and at least one of the heterodimers, the proton-bound dimer of 1,5-dimethylcytosine with 1-methylcytosine, exhibits two bands suggestive of the presence of two tautomers close in energy.

Extension of the AMBER Force Field for Nitroxide Radicals and Combined QM/MM/PCM Approach to the Accurate Determination of EPR Parameters of DMPO-H in Solution

A computational strategy that combines both time-dependent and time-independent approaches is exploited to accurately model molecular dynamics and solvent effects on the isotropic hyperfine coupling constants of the DMPO-H nitroxide. Our recent general force field for nitroxides derived from AMBER ff99SB is further extended to systems involving hydrogen atoms in β-positions with respect to NO-moiety. The resulting force-field has been employed in a series of classical molecular dynamics simulations, comparing the computed EPR parameters from selected molecular configurations to the corresponding experimental data in different solvents. The effect of vibrational averaging on the spectroscopic parameters is also taken into account, by second-order vibrational perturbation theory involving semidiagonal third energy derivatives together first and second property derivatives.

Gauging the flexibility of the active site in soybean lipoxygenase-1 (SLO-1) through an atom-centered density matrix propagation (ADMP) treatment that facilitates the sampling of rare events

We present a computational methodology to sample rare events in large biological enzymes that may involve electronically polarizing, reactive processes. The approach includes simultaneous dynamical treatment of electronic and nuclear degrees of freedom, where contributions from the electronic portion are computed using hybrid density functional theory and the computational costs are reduced through a hybrid quantum mechanics/molecular mechanics (QM/MM) treatment. Thus, the paper involves a QM/MM dynamical treatment of rare events. The method is applied to probe the effect of the active site elements on the critical hydrogen transfer step in the soybean lipoxygenase-1 (SLO-1) catalyzed oxidation of linoleic acid. It is found that the dynamical fluctuations and associated flexibility of the active site are critical toward maintaining the electrostatics in the regime where the reactive process can occur smoothly. Physical constraints enforced to limit the active site flexibility are akin to mutations and, in the cases studied, have a detrimental effect on the electrostatic fluctuations, thus adversely affecting the hydrogen transfer process.

"Pump-probe" atom-centered density matrix propagation studies to gauge anharmonicity and energy repartitioning in atmospheric reactive adducts: case study of the OH + isoprene and OH + butadiene reaction intermediates

Time-resolved "pump-probe" ab initio molecular dynamics studies are constructed to probe the stability of reaction intermediates, the mechanism of energy transfer, and energy repartitioning, for moieties involved during the interaction of volatile organic compunds with hydroxyl radical. These systems are of prime importance in the atmosphere. Specifically, the stability of reaction intermediates of hydroxyl radical adducts to isoprene and butadiene molecules is used as a case study to develop novel computational techniques involving "pump-probe" ab initio molecular dynamics. Starting with the various possible hydroxyl radical adducts to isoprene and butadiene, select vibrational modes of each of the adducts are populated with excess energy to mimic the initial conditions of an experiment. The flow of energy into the remaining modes is then probed by subjecting the excited adducts to ab initio molecular dynamics simulations. It is found that the stability of the adducts arises directly due to the anhormonically driven coupling of the modes to facilitate repartitioning of the excess vibrational energy. This kind of vibrational repartitioning has a critical influence on the energy density.

Quantum wavepacket ab initio molecular dynamics: generalizations using an extended Lagrangian treatment of diabatic states coupled through multireference electronic structure

We present a generalization to our previously developed quantum wavepacket ab initio molecular dynamics (QWAIMD) method by using multiple diabatic electronic reduced single particle density matrices, propagated within an extended Lagrangian paradigm. The Slater determinantal wavefunctions associated with the density matrices utilized may be orthogonal or nonorthogonal with respect to each other. This generalization directly results from an analysis of the variance in electronic structure with quantum nuclear degrees of freedom. The diabatic electronic states are treated here as classical parametric variables and propagated simultaneously along with the quantum wavepacket and classical nuclei. Each electronic density matrix is constrained to be N-representable. Consequently two sets of new methods are derived: extended Lagrangian-QWAIMD (xLag-QWAIMD) and diabatic extended Lagrangian-QWAIMD (DxLag-QWAIMD). In both cases, the instantaneous potential energy surface for the quantum nuclear degrees of freedom is constructed from the diabatic states using an on-the-fly nonorthogonal multireference formalism. By introducing generalized grid-based electronic basis functions, we eliminate the basis set dependence on the quantum nucleus. Subsequent reuse of the two-electron integrals during the on-the-fly potential energy surface computation stage yields a substantial reduction in computational costs. Specifically, both xLag-QWAIMD and DxLag-QWAIMD turn out to be about two orders of magnitude faster than our previously developed time-dependent deterministic sampling implementation of QWAIMD. Energy conservation properties, accuracy of the associated potential surfaces, and vibrational properties are analyzed for a family of hydrogen bonded systems.

Infrared consequence spectroscopy of gaseous protonated and metal ion cationized complexes

In this article, the new and exciting techniques of infrared consequence spectroscopy (sometimes called action spectroscopy) of gaseous ions are reviewed. These techniques include vibrational predissociation spectroscopy and infrared multiple photon dissociation spectroscopy and they typically complement one another in the systems studied and the information gained. In recent years infrared consequence spectroscopy has provided long-awaited direct evidence into the structures of gaseous ions from organometallic species to strong ionic hydrogen bonded structures to large biomolecules. Much is being learned with respect to the structures of ions without their stabilizing solvent which can be used to better understand the effect of solvent on their structures. This review mainly covers the topics with which the author has been directly involved in research: structures of proton-bound dimers, protonated amino acids and DNA bases, amino acid and DNA bases bound to metal ions and, more recently, solvated ionic complexes. It is hoped that this review reveals the impact that infrared consequence spectroscopy has had on the field of gaseous ion chemistry.

Vibrational spectroscopy of bare and solvated ionic complexes of biological relevance

The low density of ions in mass spectrometers generally precludes direct infrared (IR) absorption measurements. The IR spectrum of an ion can nonetheless be obtained by inducing photodissociation of the ion using a high-intensity tunable laser. The emergence of free electron lasers (FELs) and recent breakthroughs in bench-top lasers based on nonlinear optics have now made it possible to routinely record IR spectra of gas-phase ions. As the energy of one IR photon is insufficient to cause dissociation of molecules and strongly bound complexes, two main experimental strategies have been developed to effect photodissociation. In infrared multiple-photon dissociation (IR-MPD) many photons are absorbed resonantly and their energy is stored in the bath of vibrational modes, leading to dissociation. In the "messenger" technique a weakly bound van der Waals atom is detached upon absorption of a single photon. Fundamental, historical, and practical aspects of these methods will be presented. Both of these approaches make use of very different methods of ion preparation and manipulation. While in IR-MPD ions are irradiated in trapping mass spectrometers, the "messenger" technique is generally carried out in molecular beam instruments. The main focus of this review is the application of IR spectroscopy to biologically relevant molecular systems (amino acids, peptides, proteins). Particular issues that will be addressed here include gas-phase zwitterions, the (chemical) structures of peptides and their collision-induced dissociation (CID) products, IR spectra of gas-phase proteins, and the chelation of metal-ligand complexes. Another growing area of research is IR spectroscopy on solvated clusters, which offer a bridge between the gas-phase and solution environments. The development of state-of-the-art computational approaches has gone hand-in-hand with advances in experimental techniques. The main advantage of gas-phase cluster research, as opposed to condensed-phase experiments, is that the systems of interest can be understood in detail and structural effects can be studied in isolation. It will be shown that IR spectroscopy of mass-selected (bio)molecular systems is now well-placed to address specific questions on the individual effect of charge carriers (protons and metal ions), as well as solvent molecules on the overall structure.

Proton affinity and zwitterion stability: new results from infrared spectroscopy and theory of cationized lysine and analogues in the gas phase

The gas-phase structures of alkali metal cationized lysine (Lys), alpha-N-methyllysine (NMeLys), and epsilon-N,N-dimethyllysine (Lys(Me)(2)) are investigated using infrared multiple photon dissociation (IRMPD) spectroscopy utilizing light generated by a free electron laser and ab initio calculations. The proton affinities of the compounds span a range of approximately 20 kJ/mol. For NMeLys x M(+), experiment and theory indicate that NMeLys is nonzwitterionic for M = Li and zwitterionic for M = Na and K. For Lys(Me)(2) x M(+), experiment and theory indicate that Lys(Me)(2) is zwitterionic for M = Li, Na, and K. This is the first spectroscopic observation of the zwitterionic form of an amino acid complexed with Li(+). The results are compared with IRMPD spectra reported previously for Lys and -N-methyllysine (Lys(Me)) complexed with Li, Na, and K, and new calculations performed at higher levels of theory for those ions. The combined experimental and theoretical results indicate that protonation in the zwitterionic forms of the these amino acids is favored at the more basic methylated amine site, but that any relationship between the proton affinity of the amino acid and the relative zwitterion stability of the alkali metal cationized amino acid is only indirect. These results provide additional evidence that proton affinities are not a reliable indicator of zwitterion stability for cationized amino acids because side chains can have very different effects on the stability of different conformers in the neutral, protonated, and metal cationized forms.