Experimental validation of Gaussian-3 lithium cation affinities of amides: implications for the gas-phase lithium cation basicity scale.

Using a refined Gaussian-3 (G3) protocol, the highest level of ab initio calculations reported so far, we have established the Li+ cation binding enthalpy (affinity) at 0 K (in kJ mol-1) for formamide (195.7), N-methylformamide (209.2), N,N'-dimethylformamide (220.0), acetamide (211.7), N-methylacetamide (222.5), and N,N'-dimethylacetamide (230.1), with an estimated maximum uncertainty of +/-8 kJ mol-1. With these six theoretical lithium cation binding affinities as reference values, the absolute Li+ affinities of imidazole and dimethoxyethane were determined by the extended kinetic method, and by adopting the statistical data treatment protocol recently proposed by Armentrout. The Li+ affinities obtained for these two ligands are in good agreement (within 6 kJ mol-1) with recent values determined by the threshold collision-induced dissociation method, and consistent with the Li+ basicity values first reported by Taft and co-workers in 1990. Our study confirms that the previously suggested, and recently implemented, downward revision of Taft's original basicity scale by 10.9 kJ mol-1 is justified for ligands with revised basicities less than 151 kJ mol-1. However, for selected ligands with Li+ basicities greater than 151 kJ mol-1, including some of the six amides studied in this work, the reported discrepancy between theoretical and experimental estimates in the revised Li+ basicity scale of Burk et al. is likely to arise from experimental uncertainties.     

Experimental validation of theoretical potassium and sodium cation affinities of amides by mass spectrometric kinetic method measurements

In this study the theoretical Gaussian-2 K(+)/Na(+) binding affinities (enthalpies) at 0 K (in kJ mol(-1)) for six amides in the order: formamide (109.2/138.5) < N-methylformamide (117.7/148.6) < acetamide (118.7/149.5) < N,N-dimethylformamide (123.9/156.4) < N-methylacetamide (125.6/157.7) < N,N-dimethylacetamide (129.2/162.6), reported previously (Siu et al., J. Chem. Phys. 2001; 114: 7045-7051), were validated experimentally by mass spectrometric kinetic method measurements. By monitoring the collision-induced dissociation (CID) of K(+)/Na(+)-bound heterodimers of the amides, the relative affinities were shown to be accurate to within +/-2 kJ mol(-1). With these six theoretical K(+)/Na(+) binding affinities as reference values, the absolute K(+)/Na(+) affinities of imidazole, 1-methylimidazole, pyridazine and 1,2-dimethoxyethane were determined by the extended kinetic method, and found to be consistent (to within +/-9 kJ mol(-1)) with literature experimental values obtained by threshold-CID, equilibrium high-pressure mass spectrometry, and Fourier transform ion cyclotron resonance/ligand-exchange equilibrium methods. A self-consistent resolution is proposed for the inconsistencies in the relative order of K(+)/Na(+) affinities of amides reported in the literature. These two sets of validated K(+) and Na(+) affinity values are useful as reference values in kinetic method measurements of K(+)/Na(+) affinity of model biological ligands, such as the K(+) affinities of aliphatic amino acids.     

The sodium ion affinities of simple Di-, Tri-, and tetrapeptides

The sodium ion affinities (binding energies) of nineteen peptides containing 2-4 residues have been determined by experimental and computational approaches. Na(+)-bound heterodimers with amino acid and peptide ligands (Pep(1), Pep(2)) were produced by electrospray ionization. The dissociations of these Pep(1)-Na(+)-Pep(2) ions to Pep(1)-Na(+) and Pep(2)-Na(+) were examined by collisionally activated dissociation to construct a ladder of relative affinities via the kinetic method. The accuracy of this ladder was subsequently ascertained by experiments using several excitation energies for four peptide pairs. The relative scale was converted to absolute affinities by anchoring the relative values to the known Na(+) affinity of GlyGly. The Na(+) affinities of AlaAla, HisGly, GlyHis, GlyGlyGly, AlaAlaAla, GlyGlyGlyGly, and AlaAlaAlaAla were also calculated at the MP2(full)/6-311 + G(2d,2p) level of ab initio theory using geometries that were optimized at the MP2(full)/6-31G(d) level for AlaAla or HF/6-31G(d) level for the other peptides; the resulting values agree well with experimental Na(+) affinities. Increasing the peptide size is found to dramatically augment the Na(+) binding energy. The calculations show that in nearly all cases, all available carbonyl oxygens are sodium binding sites in the most stable structures. Whenever side chains are available, as in HisGly and GlyHis, specific additional binding sites are provided to the cation. Oligoglycines and oligoalanines have similar binding modes for the di- and tripeptides, but differ significantly for the tetrapeptides: while the lowest energy structure of GlyGlyGlyGly-Na(+) has the peptide folded around the ion with all four carbonyl oxygens in close contact with Na(+), that of AlaAlaAlaAla-Na(+) involves a pseudo-cyclic peptide in which the C and N termini interact via hydrogen bonding, while Na(+) sits on top of the oxygens of three nearly parallel C=O bonds.     

Experimental and theoretical studies of sodium cation interactions with D-arabinose, xylose, glucose, and galactose

The binding of Na (+) to arabinose (Ara), xylose (Xyl), glucose (Glc), and galactose (Gal) is examined in detail by studying the collision-induced dissociation (CID) of the four sodiated monosaccharide complexes with Xe using a guided ion beam tandem mass spectrometer (GIBMS). Analysis of the energy-dependent CID cross-sections provides 0 K sodium cation affinities for experimental complexes after accounting for unimolecular decay rates, internal energy of reactant ions, and multiple ion-neutral collisions. Quantum chemical calculations for a number of geometric conformations of each Na (+)(L) complex with a comprehensive analysis of the alpha and beta anomeric forms are determined at the B3LYP/6-311+G(d,p) level with single-point energies calculated at MP2(full), B3LYP, and B3P86 levels using a 6-311+G(2d,2p) basis set. This coordinated examination of both experimental work and quantum chemical calculations allows for determination of the bond energy for both the alpha and beta forms of each monosaccharide studied here. An understanding of the energetic contributions of individual structural characteristics as well as the energetic trends in binding among the monosaccharides is developed. Structural characteristics that affect the energetics of binding involve multidentate sodium cation coordination, ring sterics, and hydrogen bonding schemes. The overall trend in sodium binding affinities for the eight ligands follows beta-Ara < alpha-Ara < beta-Xyl < beta-Glc < alpha-Glc < alpha;-Xyl < alpha-Gal < beta-Gal.     

Experimental and theoretical studies of sodium cation interactions with the acidic amino acids and their amide derivatives

The binding of Na+ to aspartic acid (Asp), glutamic acid (Glu), asparagine (Asn), and glutamine (Gln) is examined in detail by studying the collision-induced dissociation (CID) of the four sodiated amino acid complexes with Xe using a guided ion beam tandem mass spectrometer (GIBMS). Analysis of the energy-dependent CID cross sections provides 0 K sodium cation affinities for the complexes after accounting for unimolecular decay rates, internal energy of the reactant ions, and multiple ion-molecule collisions. Quantum chemical calculations for a number of geometric conformations of each Na+(L) complex are determined at the B3LYP/6-311+G(d,p) level with single-point energies calculated at MP2(full), B3LYP, and B3P86 levels using a 6-311+G(2d,2p) basis set. This coordinated examination of both experimental work and quantum chemical calculations allows the energetic contributions of individual functionalities as well as steric influences of relative chain lengths to be thoroughly explored. Na+ binding affinities for the amide complexes are systematically stronger than those for the acid complexes by 14 +/- 1 kJ/mol, which is attributed to an inductive effect of the OH group in the carboxylic acid side chain. Additionally, the Na+ binding affinity for the longer-chain amino acids (Glx) is enhanced by 4 +/- 1 kJ/mol compared to the shorter-chain Asx because steric effects are reduced.     

Competition between pi and non-pi cation-binding sites in aromatic amino acids: a theoretical study of alkali metal cation (Li+, Na+, K+)-phenylalanine complexes

To understand the cation-pi interaction in aromatic amino acids and peptides, the binding of M(+) (where M(+) = Li(+), Na(+), and K(+)) to phenylalanine (Phe) is studied at the best level of density functional theory reported so far. The different modes of M(+) binding show the same order of binding affinity (Li(+)>Na(+)>K(+)), in the approximate ratio of 2.2:1.5:1.0. The most stable binding mode is one in which the M(+) is stabilized by a tridentate interaction between the cation and the carbonyl oxygen (O[double bond]C), amino nitrogen (--NH(2)), and aromatic pi ring; the absolute Li(+), Na(+), and K(+) affinities are estimated theoretically to be 275, 201, and 141 kJ mol(-1), respectively. Factors affecting the relative stabilities of various M(+)-Phe binding modes and conformers have been identified, with ion-dipole interaction playing an important role. We found that the trend of pi and non-pi cation bonding distances (Na(+)-pi>Na(+)-N>Na(+)-O and K(+)-pi>K(+)-N>K(+)-O) in our theoretical Na(+)/K(+)-Phe structures are in agreement with the reported X-ray crystal structures of model synthetic receptors (sodium and potassium bound lariat ether complexes), even though the average alkali metal cation-pi distance found in the crystal structures is longer. This difference between the solid and the gas-phase structures can be reconciled by taking the higher coordination number of the cations in the lariat ether complexes into account.     

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.     

Gas-phase conformations of cationized poly(styrene) oligomers

The gas-phase conformations of poly(styrene) oligomers cationized by Li+, Na+, Cu+, and Ag+ (M+ PSn) were examined using ion mobility experiments and molecular mechanics/ dynamics calculations. M+PSn ions were formed by MALDI and their ion-He collision cross-sections were measured by ion mobility methods. The experimental collision cross-sections of each M+PS n-mer were similar for all four metal cations and increased linearly with n. Molecular modeling of selected M+ PS oligomers cationized by Li+ and Na+ yielded quasi-linear structures with the metal cation sandwiched between two phenyl groups. The relative energies of the structures were approximately 2-3 kcal/mol more stable when the metal cation was sandwiched near the middle of the oligomer chain than when it was near the ends of the oligomer. The cross-sections of these theoretical structures agree well with the experimental values with deviations typically around 1-2%. The calculations also show that the metal cation tends to align the phenyl groups on the same side of the -CH2-CH- backbone. Calculations on neutral poly(styrene), on the other hand, showed structures in which the phenyl groups were more randomly positioned about the oligomer backbone. The conformations and metal-oligomer binding energies of M+PS are also used to help explain CID product distributions and fragmentation mechanisms of cationized PS oligomers. etry     

Calix[4]quinones derived from double calix[4]arenes: synthesis, complexation, and electrochemical properties toward alkali metal ions

Bis(calix[4]diquinones) 1 and 2 and double calix[4]diquinone 3 have been synthesized from their corresponding double calix[4]arenes 4, 5, and 6, respectively. Compounds 4-6 have been prepared from one-pot and stepwise syntheses under high pressure. Complexation studies of ligands 1-3 with alkali metal ions such as Li+, Na+, K+, and Cs+ were carried out by 1H NMR titrations. Receptors 1 can selectively form 1:1 complexes with Na+. Ligand 2 prefers to form 1:1 complexes with K+ and Cs+. Receptor 3 retained the cone conformation of the calix[4]arene unit upon binding K+ but changed the conformation when complexing Li+ and Na+. Electrochemical studies using cyclic voltammetry and square wave voltammetry showed significant changing of voltammograms of 2 and 3 in the presence of alkali metal ions. Receptor 3 showed the electrochemically switched binding property toward Na+ and K+.     

Cation-pi interactions with a model for the side chain of tryptophan: structures and absolute binding energies of alkali metal cation-indole complexes

Threshold collision-induced dissociation techniques are employed to determine bond dissociation energies (BDEs) of mono- and bis-complexes of alkali metal cations, Li+, Na+, K+, Rb+, and Cs+, with indole, C8H7N. The primary and lowest energy dissociation pathway in all cases is endothermic loss of an intact indole ligand. Sequential loss of a second indole ligand is observed at elevated energies for the bis-complexes. Density functional theory calculations at the B3LYP/6-31G level of theory are used to determine the structures, vibrational frequencies, and rotational constants of these complexes. Theoretical BDEs are determined from single point energy calculations at the MP2(full)/6-311+G(2d,2p) level using the B3LYP/6-31G* geometries. The agreement between theory and experiment is very good for all complexes except Li+ (C8H7N), where theory underestimates the strength of the binding. The trends in the BDEs of these alkali metal cation-indole complexes are compared with the analogous benzene and naphthalene complexes to examine the influence of the extended pi network and heteroatom on the strength of the cation-pi interaction. The Na+ and K+ binding affinities of benzene, phenol, and indole are also compared to those of the aromatic amino acids, phenylalanine, tyrosine, and tryptophan to elucidate the factors that contribute to the binding in complexes to the aromatic amino acids. The nature of the binding and trends in the BDEs of cation-pi complexes between alkali metal cations and benzene, phenol, and indole are examined to help understand nature's preference for engaging tryptophan over phenylalanine and tyrosine in cation-pi interactions in biological systems.     

Experimental and theoretical investigation of alkali metal cation interactions with hydroxyl side-chain amino acids

Absolute bond dissociation energies of serine (Ser) and threonine (Thr) to alkali metal cations are determined experimentally by threshold collision-induced dissociation of M+AA complexes, where M+=Li+, Na+, and K+ and AA=Ser and Thr, with xenon in a guided ion beam tandem mass spectrometer. Experimental results show that the binding energies of both amino acids to the alkali metal cations are very similar to one another and follow the order of Li+>Na+>K+. Quantum chemical calculations at three different levels, B3LYP, B3P86, and MP2(full), using the 6-311+G(2d,2p) basis set with geometries and zero-point energies calculated at the B3LYP/6-311+G(d,p) level show good agreement with the experimental bond energies. Theoretical calculations show that all M+AA complexes have charge-solvated structures (nonzwitterionic) with [CO, N, O] tridentate coordination.     

Alkali metal ion binding to amino acids versus their methyl esters: affinity trends and structural changes in the gas phase

The relative alkali metal ion (M(+)) affinities (binding energies) between seventeen different amino acids (AA) and the corresponding methyl esters (AAOMe) were determined in the gas phase by the kinetic method based on the dissociation of AA-M(+)-AAOMe heterodimers (M=Li, Na, K, Cs). With the exception of proline, the Li(+), Na(+), and K(+) affinities of the other aliphatic amino acids increase in the order AAAAOMe is already observed for K(+). Proline binds more strongly than its methyl ester to all M(+) except Li(+). Ab initio calculations on the M(+) complexes of alanine, beta-aminoisobutyric acid, proline, glycine methyl ester, alanine methyl ester, and proline methyl ester show that their energetically most favorable complexes result from charge solvation, except for proline which forms salt bridges. The most stable mode of charge solvation depends on the ligand (AA or AAOMe) and, for AA, it gradually changes with metal ion size. Esters chelate all M(+) ions through the amine and carbonyl groups. Amino acids coordinate Li(+) and Na(+) ions through the amine and carbonyl groups as well, but K(+) and Cs(+) ions are coordinated by the O atoms of the carboxyl group. Upon consideration of these differences in favored binding geometries, the theoretically derived relative M(+) affinities between aliphatic AA and AAOMe are in good overall agreement with the above given experimental trends. The majority of side chain functionalized amino acids studied show experimentally the affinity order AAAAOMe. The latter ranking is attributed to salt bridge formation.     

Ionization energies of hypervalent Li2F,Li2Cl and Na2Cl molecules obtained by surface ionization electron impact neutralization mass spectrometry

Ionization energies of hypervalent Li(2)F, Li(2)Cl and Na(2)Cl molecules detected by surface ionization electron impact neutralization mass spectrometry are reported. The ionization energies were 3.78 +/- 0.2 eV for Li(2)F, 4.93 +/- 0.2 eV for Li(2)Cl, and 4.21 +/- 0.2 eV for Na(2)Cl. The ionization energies (IE) agree with theoretical ionization energies calculated by ab initio methods, supporting the theoretical prediction that Li(2)F has a hyperlithiated configuration in which the odd electron delocalizes over the two lithiums and with photoionization measurement. The first ionization energy of Na(2)Cl was experimentally confirmed earlier and for Li(2)Cl as well.8 We have developed and used this new approach for the problem--in the present work ions were first formed by surface ionization, followed by electron attachment (neutralization).     

pH-triggered assembly of organometallic receptors for lithium ions

The reaction of half-sandwich complexes of ruthenium, rhodium, and iridium with amino-substituted 3-hydroxy-2-pyridone ligands in aqueous solution gives monomeric O,O'-chelate complexes. Upon addition of base, the complexes assemble to form trimeric metallamacrocycles, as evidenced by NMR spectroscopy and single-crystal X-ray analyses. The macrocycles are able to act as highly selective receptors for lithium ions. The binding constants depend on the nature of the half-sandwich complex, the ligand, and the pH. With a commercially available (cymene)Ru complex, a receptor with a Li+ binding constant of K(a) = 5.8 (+/-1.0) x 10(4) M(-1) and a Li+-Na+ selectivity of 10 000:1 can be obtained. The fact that the assembly process of the receptor is pH-dependent can be used to detect the presence of lithium ions by a pH measurement. Furthermore, it is possible to transduce the binding of Li+ into a change of color by means of a chemical reaction with FeCl(3). This allows the detection of Li+ in the pharmacologically relevant concentration range of 0.5-1.5 mM by the "naked eye".     

A PGSE diffusion and electrophoretic NMR study of Cs+ and Na+ dynamics in aqueous crown ether systems

Multinuclear pulsed gradient spin-echo (PGSE) NMR diffusion and linewidth measurements were used to probe binding and transport in aqueous Na+-15-crown-5, Na+-18-crown-6, Cs+-15-crown-5 and Cs+-18-crown-6 systems. Since direct PGSE observation of many alkali cations is precluded by either low inherent sensitivity or rapid relaxation (or both), the feasibility of proton-detected electrophoretic NMR (ENMR) measurements to complement PGSE data was investigated. ENMR measurements were performed on aqueous Cs+-, Li+-, Na+-, K+-, and Rb+- 18-crown-6 systems. The data analysis is based on a two-site binding model and its corresponding association constants. Cs+ was found to bind considerably more tightly to 18-crown-6 (K=8 M-1) than to 15-crown-5 (K approximately 2 M-1), whereas Na+ had almost equal affinity (K approximately 4.5 M-1) for 15-crown-5 and 18-crown-6. The difficulties encountered in analysing the NMR parameters, methodological limitations and the implied need for more complicated binding models are discussed.     

Experimental and theoretical studies of potassium cation interactions with the acidic amino acids and their amide derivatives

The binding of K(+) to aspartic acid (Asp), glutamic acid (Glu), asparagine (Asn), and glutamine (Gln) is examined in detail by studying the collision-induced dissociation (CID) of the four potassium cation-bound amino acid complexes with Xe using a guided ion beam tandem mass spectrometer (GIBMS). Formed by electrospray ionization, these complexes have energy-dependent CID cross sections that are analyzed to provide 0 K bond energies after accounting for unimolecular decay rates, internal energy of reactant ions, and multiple ion-molecule collisions. Quantum chemical calculations for a number of geometric conformations of each K(+)(L) complex are determined at the B3LYP/6-311+G(d,p) level with single-point energies calculated at B3LYP, B3P86, and MP2(full) levels using a 6-311+G(2d,2p) basis set. Theoretical bond dissociation energies are in good agreement with the experimental values. This coordinated examination of both experimental work and quantum chemical calculations allows for a comprehensive understanding of the molecular interactions of K(+) with the Asx and Glx amino acids. K(+) binding affinities for the amide complexes are systematically stronger than those for the acid complexes by 9+/-1 kJ/mol, which is attributed to an inductive effect of the OH group in the carboxylic acid side chain. Additionally, the K(+) binding affinity for the longer-chain amino acids (Glx) is enhanced by 5+/-1 kJ/mol compared to the shorter-chain Asx because steric effects are reduced. Further, a detailed comparison between experimental and theoretical results reveals interesting differences in the binding of K(+) and Na(+) to these amino acids.     

Experimental and theoretical studies of sodium cation complexes of the deamidation and dehydration products of asparagine, glutamine, aspartic acid, and glutamic acid

The deamidation and dehydration products of Na+(L), where L = asparagine (Asn), glutamine (Gln), aspartic acid (Asp), and glutamic acid (Glu), are examined in detail utilizing collision-induced dissociation (CID) with Xe in a guided ion beam tandem mass spectrometer (GIBMS). Results establish that the Na+(L) complexes decompose upon formation in our dc discharge/flow tube ion source to form a bis-ligand complex, Na+(L-HX)(HX), composed of a sodium cation, the (L-HX) decomposition product, and HX, where HX = NH3 for the amides and H2O for the acids. Analysis of the energy-dependent CID cross sections for the Na+(L-HX)(HX) complexes provides unambiguous identification of the (L-HX) fragmentation products as 3-amino succinic anhydride (a-SA) for Asx and oxo-proline (O-Pro) for Glx. Furthermore, these experiments establish the 0 K sodium cation affinities for these five-membered ring decomposition products and the H2O and NH3 binding affinities of the Na+(a-SA) and Na+(O-Pro) complexes after accounting for unimolecular decay rates, the internal energy of reactant ions, and multiple ion-molecule collisions. Quantum chemical calculations are determined for a number of geometric conformations of all reaction species as well as a number of candidate species for (L-HX) at the B3LYP/6-311+G(d,p) level with single-point energies calculated at MP2(full), B3LYP, and B3P86 levels using a 6-311+G(2d,2p) basis set. This coordinated examination of both the experimental work and quantum chemical calculations allows for a complete characterization of the products of deamidation and dehydration of Asx and Glx, as well as the details of Na+, H2O, and NH3 binding to the decomposition species.     

A recipe for designing molecules with ever-increasing electron affinities

Halogens possess among the highest electron affinities of elements in the periodic table. Superhalogen molecules with electron affinities higher than those of halogen atoms have been known to form when a metal atom is surrounded with halogen atoms. Recently, it was discovered that a new class of molecules called hyperhalogens with electron affinities higher than those of superhalogens can form when the latter serve as the building block. By use of density functional theory and B3LYP hybrid exchange-correlation functional we show that molecules with electron affinities even higher can be formed by using hyperhalogens as building blocks. We demonstrate this by using Na and Li as metal atoms and F, BF(4), and Na(BF(4))(2) as halogen, superhalogen, and hyperhalogen building blocks. The predicted electron affinities of Na[Na(BF(4))(2)](2) and Li[Li(BF(4))(2)](2) are 9.18 and 9.01 eV, which are, respectively, 0.85 and 0.5 eV higher than those of their hyperhalogen [Na(BF(4))(2) and Li(BF(4))(2)] counterparts.     

Influence of halogenation on the properties of uracil and its noncovalent interactions with alkali metal ions. Threshold collision-induced dissociation and theoretical studies

The influence of halogenation on the properties of uracil and its noncovalent interactions with alkali metal ions is investigated both experimentally and theoretically. Bond dissociation energies of alkali metal ion-halouracil complexes, M+(XU), are determined using threshold collision-induced dissociation techniques in a guided ion beam mass spectrometer, where M+ = Li+, Na+, and K+ and XU = 5-fluorouracil, 5-chlorouracil, 6-chlorouracil, 5-bromouracil, and 5-iodouracil. The structures and theoretical bond dissociation energies of these complexes are determined from ab initio calculations. Theoretical calculations are also performed to examine the influence of halogenation on the acidities, proton affinities, and Watson-Crick base pairing energies. Halogenation of uracil is found to produce a decrease in the proton affinity, an increase in the alkali metal ion binding affinities, an increase in the acidity, and stabilization of the A::U base pair. In addition, alkali metal ion binding is expected to lead to an increase in the stability of nucleic acids by reducing the charge on the nucleic acid in a zwitterion effect as well as through additional noncovalent interactions between the alkali metal ion and the nucleobases.