Competitive Gas-Phase Solvation of Alkali Metal Ions by Water and Methanol

Preferential Ion Microsolvation in Mixed-Modifier Environments Observed Using Differential Mobility Spectrometry

The preferential solvation behavior for eight different derivatives of protonated quinoline was measured in a tandem differential mobility spectrometer mass spectrometer (DMS-MS). Ion-solvent cluster formation was induced in the DMS by the addition of chemical modifiers (i.e., solvent vapors) to the N₂ buffer gas. To determine the effect of more than one modifier in the DMS environment, we performed DMS experiments with varying mixtures of water, acetonitrile, and isopropyl alcohol solvent vapors. The results show that doping the buffer gas with a binary mixture of modifiers leads to the ions binding preferentially to one modifier over another. We used density functional theory to calculate the ion-solvent binding energies, and in all cases, calculations show that the quinolinium ions bind most strongly with acetonitrile, then isopropyl alcohol, and most weakly with water. Computational results support the hypothesis that the quinolinium ions bind exclusively to whichever solvent they have the strongest interaction with, regardless of the presence of other modifier gases.

LiCl solvation in N-methyl-acetamide (NMA) as a model for understanding Li(+) binding to an amide plane

The thermodynamics of ion solvation in non-aqueous solvents remains of great significance for understanding cellular transport and ion homeostasis for the design of novel ion-selective materials and applications in molecular pharmacology. Molecular simulations play pivotal roles in connecting experimental measurements to the microscopic structures of liquids. One of the most useful and versatile mimetic systems for understanding biological ion transport is N-methyl-acetamide (NMA). A plethora of theoretical studies for ion solvation in NMA have appeared recently, but further progress is limited by two factors. One is an apparent lack of experimental data on solubility and thermodynamics of solvation for a broad panel of 1 : 1 salts over an appropriate temperature and concentration range. The second concern is more substantial and has to do with the limitations hardwired in the additive (fixed charge) approximations used for most of the existing force-fields. In this submission, we report on the experimental evaluation of LiCl solvation in NMA over a broad range of concentrations and temperatures and compare the results with those of MD simulations with several additive and one polarizable force-field (Drude). By comparing our simulations and experimental results to density functional theory computations, we discuss the limiting factors in existing potential functions. To evaluate the possible implications of explicit and implicit polarizability treatments on ion permeation across biological channels, we performed potential of mean force (PMF) computations for Li(+) transport through a model narrow ion channel with additive and polarizable force-fields.

Cationic Noncovalent Interactions: Energetics and Periodic Trends

In this review, noncovalent interactions of ions with neutral molecules are discussed. After defining the scope of the article, which excludes anionic and most protonated systems, methods associated with measuring thermodynamic information for such systems are briefly recounted. An extensive set of tables detailing available thermodynamic information for the noncovalent interactions of metal cations with a host of ligands is provided. Ligands include small molecules (H2, NH3, CO, CS, H2O, CH3CN, and others), organic ligands (O- and N-donors, crown ethers and related molecules, MALDI matrix molecules), π-ligands (alkenes, alkynes, benzene, and substituted benzenes), miscellaneous inorganic ligands, and biological systems (amino acids, peptides, sugars, nucleobases, nucleosides, and nucleotides). Hydration of metalated biological systems is also included along with selected proton-based systems: 18-crown-6 polyether with protonated peptides and base-pairing energies of nucleobases. In all cases, the literature thermochemistry is evaluated and, in many cases, reanchored or adjusted to 0 K bond dissociation energies. Trends in these values are discussed and related to a variety of simple molecular concepts.

Role of methyl-induced polarization in ion binding

The chemical property of methyl groups that renders them indispensable to biomolecules is their hydrophobicity. Quantum mechanical studies undertaken here to understand the effect of point substitutions on potassium (K-) channels illustrate quantitatively how methyl-induced polarization also contributes to biomolecular function. K- channels regulate transmembrane salt concentration gradients by transporting K(+) ions selectively. One of the K(+) binding sites in the channel's selectivity filter, the S4 site, also binds Ba(2+) ions, which blocks K(+) transport. This inhibitory property of Ba(2+) ions has been vital in understanding K-channel mechanism. In most K-channels, the S4 site is composed of four threonine amino acids. The K channels that carry serine instead of threonine are significantly less susceptible to Ba(2+) block and have reduced stabilities. We find that these differences can be explained by the lower polarizability of serine compared with threonine, because serine carries one less branched methyl group than threonine. A T→S substitution in the S4 site reduces its polarizability, which, in turn, reduces ion binding by several kilocalories per mole. Although the loss in binding affinity is high for Ba(2+), the loss in K(+) binding affinity is also significant thermodynamically, which reduces channel stability. These results highlight, in general, how biomolecular function can rely on the polarization induced by methyl groups, especially those that are proximal to charged moieties, including ions, titratable amino acids, sulfates, phosphates, and nucleotides.

Determination of apparent decomposition threshold energies of lithium adducts of acylglycerols using tandem mass spectrometry and a novel derived effective reaction path length approach

Apparent decomposition threshold energies for the fragmentation pathways of lithiated acylglycerols were experimentally determined by collisional activation in a quadrupole-hexapole-quadrupole (QhQ) mass spectrometer. A previously developed 'derived effective reaction path length' approach for predicting bond dissociation energies (BDEs) of simple dissociations of electrostatic complexes such as alkali metal adducts (Li+), or halide adducts (Cl(-)) of acylglycerols, was extended to predict covalent bond apparent decomposition threshold energies of lithium adducts of a mono-acylglycerol, a 1,2-diacylglycerol, and a 1,3-diacylglycerol. The ability of the model to treat relatively large ionic systems (e.g. more than 100 atoms) represents a huge advantage of this approach. The model's calculated apparent decomposition threshold energies (Ea) are used in conjunction with the method of energy-resolved mass spectrometry, employing breakdown graphs, to give a more complete quantitative description of the fragmentation processes. Calculated Ea values allowed ranking of the 1,2-diacylglycerol as more reactive than the 1,3-diacylglycerol; the mono-acylglycerol was ranked the least reactive. The method was applied to the low molecular weight product ions generally associated with the hydrocarbon series CnH2n+1+, where two separate pathways are deduced as contributing to the production of the abundant m/z 81 fragment ion. The favored ranking of the neutral losses of fatty acyl substituents for the 1,2-diacylglycerol was determined as: loss of lithium fatty acetate > loss of fatty acid > loss of fatty acyl chain as ketene. For the 1,3-diacylglycerol, the descending order of ease of neutral loss was: loss of fatty acyl ketene > loss of lithium fatty acetate > loss of fatty acid. The results of this study demonstrate that the newly developed method is general in nature, and it can be used for the measurement of covalent bond decomposition threshold energies, as well as for the previously documented electrostatic (noncovalent) bond energies.

Infrared spectroscopy of Cr+(H2O) and Cr2+(H2O): the role of charge in cation hydration

Singly and doubly charged chromium-water ion-molecule complexes are produced by laser vaporization in a pulsed-nozzle cluster source. These species are detected and mass-selected in a specially designed time-of-flight mass spectrometer. Vibrational spectroscopy is measured for these complexes in the O-H stretching region using infrared photodissociation spectroscopy and the method of rare gas atom predissociation. Infrared excitation is not able to break the ion-water bonds in these systems, but it leads to elimination of argon, providing an efficient mechanism for detecting the spectrum. The O-H stretches for both singly and doubly charged complexes are shifted to frequencies lower than those for the free water molecule, and the intensity of the symmetric stretch band is strongly enhanced relative to the asymmetric stretch. Partially resolved rotational structure for both complexes shows that the H-O-H bond angle is greater than it is in the free water molecule. These polarization-induced effects are enhanced in the doubly charged ion relative to its singly charged analog.

A gas-phase study of the preferential solvation of Mn(2+) in mixed water/methanol clusters

The kinetic shift that exists between two competing unimolecular fragmentation processes has been used to establish whether or not gas-phase Mn(2+) exhibits preferential solvation when forming mixed clusters with water and methanol. Supported by molecular orbital calculations, these first results for a metal dication demonstrate that Mn(2+) prefers to be solvated by methanol in the primary solvation shell.

Study of doubly charged alkaline earth metal and 3-azidopropionitrile complexes by electrospray ionization mass spectrometry

The present work describes a study of the complexation of calcium and magnesium by 3-azidopropionitrile by means of electrospray ionization mass spectrometry (ESI-MS). Complexes were obtained from solutions of calcium and magnesium salts of the type CaX2 and MgX2 (where X = Cl or NO3) in water and methanol/water. The complexes detected were mainly double positively charged, with various stoichiometries not depending on the solvent, since water and 3-azidopropionitrile were always the main ligands. Solvation with methanol was not observed unlike in a previous study of complexation of nickel and cobalt by 3-azidopropionitrile. The complex ions [M(II)Az4(H2O)](2+), [M(II)Az5](2+) (where M = Ca and Mg) are the most abundant for both metals, and both counter ions. Tandem mass spectrometric (MS/MS) analysis showed that, under collision-induced dissociation (CID) conditions, the most important processes occurring were loss of neutral ligands and the replacement of 3-azidopropionitrile by water. A complex species containing reduced alkaline earth metal was due to radical loss, resulting from homolytic cleavage in the azide ligand. Some terminal ions, in the fragmentation sequences, point to the nitrile group as the coordination site in the 3-azidopropionitrile. Density functional theory (DFT) calculations confirmed this coordination site and proved that 3-azidopropionitrile behaves as a monodentate ligand in the systems under study. Moreover, the theoretical study proved that the presence of water ligand introduces stability through a hydrogen bond established between the water molecule and one nitrogen atom of the azido group. In addition, the strong dipole moment of 3-azidopropionitrile (4.76 D), which is mainly related to presence of the nitrile group, favors the stabilization of the metal-ligand complexes through charge-dipole interactions and the coordination of the metal to the nitrile group.

Properties of microsolvated ions: From the microenvironment of chromophore and alkali metal ions in proteins to negative ions in water clusters

Here we discuss the fascinating chemistry and physics of microsolvated ions that bridge the transition from bare ions in gas phase to ions in solution. Such ions occur in many situations in biochemistry and are crucial for several functions; metal ions, for example, must remove their water shell to pass through ion pumps in membranes. Furthermore, only a few water molecules are buried in the hydrophobic pockets of proteins where they are bound to charged amino acid residues or ionic chromophores. Another aspect is the reactivity of microsolvated ions and the importance in atmospheric, organic and inorganic chemistry. We close by a discussion of the stability of molecular dianions, and how hydration affects the electronic binding energy. There is a vast literature on microsolvated ions, and in this review we are far from being comprehensive, rather we mainly bring examples of our own work.

Solvation dynamics in Ni+ (H2O)n clusters probed with infrared spectroscopy

Infrared photodissociation spectroscopy is reported for mass-selected Ni+ (H2O)n complexes in the O-H stretching region up to cluster sizes of n = 25. These clusters fragment by the loss of one or more intact water molecules, and their excitation spectra show distinct bands in the region of the symmetric and asymmetric stretches of water. The first evidence for hydrogen bonding, indicated by a broad band strongly red-shifted from the free OH region, appears at the cluster size of n = 4. At larger cluster sizes, additional red-shifted structure evolves over a broader wavelength range in the hydrogen-bonding region. In the free OH region, the symmetric stretch gradually diminishes in intensity, while the asymmetric stretch develops into a closely spaced doublet near 3700 cm(-1). The data indicate that essentially all of the water molecules are in a hydrogen-bonded network by the size of n = 10. However, there is no evidence for the formation of clathrate structures seen recently via IR spectroscopy of protonated water clusters.

Extending, and Repositioning, a Thermochemical Ladder: High-Level Quantum Chemical Calculations on the Sodium Cation Affinity Scale

High-level ab initio quantum chemical calculations, at the CP-dG2thaw level of theory, are reported for coordination of Na+ to a wide assortment of small organic and inorganic ligands. The ligands range in size from H to C6H6, and include 22 of the ligands for which precise relative sodium ion binding free energies have been determined by recent Fourier transform ion cyclotron resonance and guided ion beam studies. Agreement with the relative experimental values is excellent (+/-1.1 kJ mol(-1)), and agreement with the absolute scale (obtained when these relative values are pegged to the CH3NH2 "anchor" value measured in a high-pressure mass spectrometric study) is only marginally poorer, with CP-dG2thaw values exceeding the absolute experimental DeltaG(298) values by an average of 2.1 kJ mol(-1). The excellent agreement between experiment and the CP-dG2thaw technique also suggests that the additional 97 ligands surveyed here (which, in many cases, are not readily susceptible to laboratory investigation) can also be reliably fitted to the existing experimental scale. However, while CP-dG2thaw and the experimental ladder are in close accord, a small set of higher level ab initio calculations on sodium ion/ligand complexes (including several values obtained here using the W1 protocol) suggests that the CP-dG2thaw values are themselves too low by approximately 2.5 kJ mol(-1), thereby implying that the accepted laboratory values are typically 4.6 kJ mol(-1) too low. The present work also highlights the importance of Na+/ligand binding energy determinations (whether by experimental or theoretical approaches) on a case-by-case basis: trends in increasing binding energy along homologous series of compounds are not reliably predictable, nor are binding site preferences or chelating tendencies in polyfunctional compounds.

Determination of Bond Dissociation Energies Using Electrospray Tandem Mass Spectrometry and a Derived Effective Reaction Path Length Approach

A new approach for calculating bond dissociation energies (BDEs) from ES-MS/MS measurements has been developed. The new method features a "derived effective reaction path length" that has been applied to measure BDEs of alkali metal (Li+) adducts and halide (Cl-) adducts of monoacylglycerol, 1,2-diacylglycerol, and 1,3-diacylglycerol lipids. Also studied were lithium-bound dimers of monoacylglycerols, 1,2-diacylglycerols, and 1,3-diacylglycerols. BDEs for the adducts and dimers of the lipids were derived from collision-induced dissociation experiments using a triple quadrupole mass spectrometer with electrospray as the ionization source. Mass spectral data were used to empirically derive a single-exponential growth equation that relates product cross section to collision energy. From these single-exponential equations, a general second-order polynomial was derived using a multivariate growth curve model that enables prediction of BDEs of unknown complexes. Mass spectral results were compared to computer-generated bond dissociation energies using Becke-style three-parameter density functional theory (B3LYP, employing the Lee-Yang-Parr correlation functional), with excellent agreement between experimental and theoretical energy values. The newly developed method is general in nature and can be used for the measurement of metal or halide ionic adduct bond dissociation energies and for the measurement of bond energies of noncovalent interactions such as dimer dissociation energies. The validity of the method has been rigorously established using a triple quadrupole, but it may also be applied to other mass spectrometers that allow user control of the collision cell potential.

Electron Binding Energies of Aqueous Alkali and Halide Ions: EUV Photoelectron Spectroscopy of Liquid Solutions and Combined Ab Initio and Molecular Dynamics Calculations

Photoelectron spectroscopy combined with the liquid microjet technique enables the direct probing of the electronic structure of aqueous solutions. We report measured and calculated lowest vertical electron binding energies of aqueous alkali cations and halide anions. In some cases, ejection from deeper electronic levels of the solute could be observed. Electron binding energies of a given aqueous ion are found to be independent of the counterion and the salt concentration. The experimental results are complemented by ab initio calculations, at the MP2 and CCSD(T) level, of the ionization energies of these prototype ions in the aqueous phase. The solvent effect was accounted for in the electronic structure calculations in two ways. An explicit inclusion of discrete water molecules using a set of snapshots from an equilibrium classical molecular dynamics simulations and a fractional charge representation of solvent molecules give good results for halide ions. The electron binding energies of alkali cations computed with this approach tend to be overestimated. On the other hand, the polarizable continuum model, which strictly provides adiabatic binding energies, performs well for the alkali cations but fails for the halides. Photon energies in the experiment were in the EUV region (typically 100 eV) for which the technique is probing the top layers of the liquid sample. Hence, the reported energies of aqueous ions are closely connected with both structures and chemical reactivity at the liquid interface, for example, in atmospheric aerosol particles, as well as fundamental bulk solvation properties.

Competing Non-covalent Interactions in Alkali Metal Ion−Acetonitrile−Water Clusters

Competitive ion-dipole, ion-water, and water-water interactions were investigated at the molecular level in M+ (CH3CN)n(H2O)m cluster ions for M = Na and K. Different [n,m] combinations for two different n + m cluster sizes were characterized with infrared predissociation spectroscopy in the O-H stretch region and MP2 calculations. In all cases, no differences were observed between the two alkali metal ions. The results showed that at the n + m = 4 cluster size, the solvent molecules interact only with the ion, and that the interaction between the ion and the large dipole moment of CH3CN decreases the ion-water electrostatic interactions. At the n + m = 5 cluster size, at least two different hydrogen-bonded structures were identified. In these structures, the ion-dipole interaction weakens the ability of the ion to polarize the hydrogen bonds and thus decreases the strength of the water-water interactions in the immediate vicinity of the alkali metal ion.

Formation of solvated ions in the atmospheric interface of an electrospray ionization triple-quadrupole mass spectrometer

A simple method capable of generating and investigating various solvent clusters and solvated ions was developed. The technique opens a door to studying these complexes on commercially available instruments. Formation of the desired solvated ion in the gas phase was achieved by introducing the appropriate volatile solvent vapour into the curtain gas stream. Capabilities of the technique are illustrated by generating alkali, alkaline earth and transition metal cations solvated by various volatile compounds such as water, methanol and acetonitrile. Depending on the ligands and on the experimental conditions, clusters of 2-100 molecules may be observed. Isotope labelling suggests that these are formed by a re-solvation process in the curtain gas region.

Coordination chemistry of silver cations

While in pure solvents Ag(+) is known to be tetrahedrally coordinated, in the presence of ligands such as ammonia it forms linear complexes, usually explained by the ion's tendency toward sd-hybridization. To explore this disparity, we have investigated the reaction of ammoniated silver cations Ag(+)(NH(3))(n)(), n = 11-23, with H(2)O as well as the complementary process, the reaction of Ag(+)(H(2)O)(n)(), n = 25-45, with NH(3) by means of FT-ICR mass spectrometry. In both cases, ligand exchange reactions take place, leading to clusters with a limited number of NH(3) ligands. The former reaction proceeds very rapidly until only three NH(3) ligands are left, followed by a much slower loss of an additional ligand to form Ag(+)(NH(3))(2)(H(2)O)(m)() clusters. In the complementary process, the reaction of Ag(+)(H(2)O)(n)() with NH(3) five ammonia ligands are very rapidly taken up by the clusters, with a much less efficient uptake of a sixth one. The accompanying DFT calculations reveal a delicate balance between competing effects where not only the preference of Ag(+) for sd-hybridization, but also its ability to polarize the ligands and thus affect the strength of their hydrogen bonding, as well as the ability of the solvent to form extended hydrogen-bonded networks are important.