The role of proton affinity, acidity, and electrostatics on the stability of neutral versus ion-pair forms of molecular dimers Eric F. Strittmatter, Evan R. Williams* Department of Chemistry, University of California, Berkeley, CA 94720-1460, USA Received 21 September 2000, accepted 17 May 2001 Abstract Ion-pair formation via proton transfer from an acidic hydrogen of one functional group to a basic functional group plays an important role in the structure and reactivity of biomolecules in the gas phase. The relative stabilities of the ion-pair and the neutral-pair forms of five dimers composed of a basic molecule and a trifluoroacetic acid molecule were compared using density-functional calculations. The proton affinity of the basic molecules investigated ranged from 246 to 254 kcal/mol. The gas phase acidity of trifluoroacetic acid is 323.8 kcal/mol. The results of the B3LYP (6-311++G**) calculations indicate that the structures of the dimers change from a neutral pair to an ion pair as the proton affinity of the bases increases. This result is consistent with previous blackbody infrared radiative dissociation experiments on protonated trimolecular complexes (or trimers) consisting of two basic molecules and trifluoroacetic acid, which indicates that the predominant structure of the trimer changes from a charge-solvation structure to a salt bridge structure with the increasing gas phase basicity of the base. The electrostatic character of the interaction between the basic molecule and the trifluoroacetic acid molecule was determined using the natural energy-decomposition analysis (NEDA) program. In the ion pair, a majority (69%–77%) of the attractive energy of the dimer is comprised of the electrostatic component. Two models are derived that include the acidity of the acidic molecule, the proton affinity of the basic molecule, and an electrostatic binding term for both the ion pair and the neutral pair. Several nonelectrostatic interaction terms can be replaced by a single correction or constant term so that both models, one using NEDA electrostatic terms and the other using integration of point-charge interactions, provide reasonably accurate results. This indicates that electrostatic models similar to the ones used here may be useful in studying salt bridge formation in larger molecules. (Int J Mass Spectrom 212 (2001) 287–300) © 2001 Elsevier Science B.V. 1. Introduction For over a century, it has been recognized that amino acids exhibit amphoteric properties (thus can exist as acid or base) in aqueous solution [1,2]. The structure of amino acids was deduced largely from the solution behavior of these compounds under the influence of electric fields. In 1894, Georg Bredig suggested that betaine exists as an “inneres salz” (inner salt) containing both a positive and a negative charge in the same molecule. A few years later, Ku ¨ster [3] recognized that the pH-dependant color change of methylorange arises from the transition between the anionic form and inner salt form (H + (CH 3 ) 2 N-C 6 H 4 -N 2 -C 6 H 4 -SO 3 + ). The inner salt com- pound is ionized in aqueous solution, yet it does not conduct current under an electric field. Ku ¨ster termed this nonconducting ion a “zwitterion” [3]. All amino acids were later found to exist in their zwitterionic state under physiological conditions [2]. * Corresponding author. E-mail: williams@cchem.berkeley.edu Dedicated to R. Graham Cooks on the occasion of his sixtieth birthday. 1387-3806/01/$20.00 © 2001 Elsevier Science B.V. All rights reserved PII S1387-3806(01)00475-4 International Journal of Mass Spectrometry 212 (2001) 287–300 www.elsevier.com/locate/ijms