MD simulation of the Na+-phenylalanine complex in water: competition between cation-pi interaction and aqueous solvation

The competition between cation-pi interaction and aqueous solvation for the Na+ ion has been investigated by molecular dynamics simulations, using the phenylalanine amino acid as the test pi system. Starting from one of the best standard force fields, we have developed new parameters that significantly improve the agreement with experimental and high quality quantum mechanical results for the complexes of Na+ with phenylalanine, benzene, and water. The modified force field performs very well in forecasting energy and geometry of cation coordination for the complexes. Next, analysis of MD trajectories and steered MD simulations indicate that the Na+-phenylalanine complex survives for a significant time in aqueous solution and that the free energy barrier opposing dissociation of the complex is sizable. Finally, we analyze the role of different intermolecular interactions in determining the preference for cation-pi bonding with respect to aqueous solvation. We thus confirm that the Na+-phenylalanine stabilization energy may overcome the interactions with water.     

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.     

How does ammonium dynamically interact with benzene in aqueous media? A first principle study using the Car-Parrinello molecular dynamics method

The Car-Parrinello molecular dynamics (CPMD) method was used to study the dynamic characteristics of the cation-pi interaction between ammonium and benzene in gaseous and aqueous media. The results obtained from the CPMD calculation on the cation-pi complex in the gaseous state were very similar to those calculated from the Gaussian98 program with DFT and MP2 algorithms, demonstrating that CPMD is a valid approach for studying this system. Unlike the interaction in the gaseous state, our 12-ps CPMD simulation showed that the geometry of the complex in aqueous solution changes frequently in terms of the interaction angles and distances. Furthermore, the simulation revealed that the ammonium is constantly oscillating above the benzene plane in an aqueous environment and interacts with benzene mostly through three of its hydrogen atoms. In contrast, the interaction of the cation with the aromatic molecule in the gaseous state involves two hydrogen atoms. In addition, the free energy profile in aqueous solution was studied using constrained CPMD simulations, resulting in a calculated binding free energy of -5.75 kcal/mol at an optimum interaction distance of approximately 3.25 A, indicating that the cation-pi interaction between ammonium and benzene is stable even in aqueous solution. Thus, this CPMD study suggested that the cation-pi interaction between an ammonium (group) and an aromatic structure could take place even on surfaces of protein or nucleic acids in solution.     

Cation-pi interactions and the gas-phase thermochemistry of the Na(+)/phenylalanine complex

The complex of Na(+) with phenylalanine (Phe) is a prototype for the participation of cation-pi interactions in metal-ion binding to biological molecules. A recent comparison of this complex with the Na(+)/alanine (Na(+)/Ala) counterpart suggested only a small contribution of the phenyl ring interaction to binding, casting doubt on the extent of the cation-pi effect. The present work reexamines this thermochemistry using ligand-exchange equilibrium measurements in the Fourier transform ion cyclotron resonance (FT-ICR) ion trapping mass spectrometer. An increment of 7 +/- 2 kcal mol(-1) was found in the Ala/Phe comparison of binding enthalpies, confirming the importance of cation-pi binding enhancement in the Phe case. Absolute Na(+) binding enthalpies of 38 +/- 2 and 45 +/- 2 kcal mol(-1) were assigned for Ala and Phe, respectively, using pyridine as the thermochemical reference ligand. All of these results were supported by quantum calculations using both density functional and Hartree-Fock/MP2 methods, improved in several respects over previous calculations. Alanine methyl ester (AlaMe) was also observed, and found to have an Na(+) ion affinity larger by 2.3 kcal mol(-1) than Ala. New, lower energy conformations of neutral Phe were discovered in the computations.     

Fullerenols revisited as stable radical anions

The first exhaustive purification and characterization of the much-studied "fullerenols", prepared by reaction of C(60) in toluene with an oxygenated, aqueous NaOH solution using tetrabutylammonium hydroxide as a phase transfer catalyst, has been performed. The resulting fullerenol is not simply polyhydroxylated C(60) but rather is a structurally and electronically complex C(60) radical anion with a molecular formula of Na(+)(n)[C(60)O(x)(OH)(y)](n)(-) (where n = 2-3, x = 7-9, and y = 12-15) for three different, but identical, preparations. Surprisingly, Na(+)-fullerenol is paramagnetic, exhibiting mu(B) values in aqueous solution of 1.9-2.1 B.M. at 0.5 T and 300 K and R(1) proton relaxivities of 0.55-0.77 mM(-1)s(-1) at 20 MHz and 40 degrees C, values both slightly higher than those expected for a pure S = 1/2 spin system. ESR studies (ESE-FS and 2D nutation) of frozen aqueous solutions at 1.5 and 5.0 K establish that Na(+)-fullerenol is mainly S = 1/2 with a minor, but significant, component of S = 1. Thus, this is the first report to characterize these widely studied, water-soluble fullerenols as stable radical anions. The stability of the S = 1/2 Na(+)-fullerenol radical is likely due to a highly derivatized C(60) surface that protects a cyclopentadienyl radical center on the fullerene.     

Glycosidic linkage conformation of methyl-alpha-mannopyranoside

We study the preferred conformation of the glycosidic linkage of methyl-alpha-mannopyranoside in the gas phase and in aqueous solution. Results obtained utilizing Car-Parrinello molecular dynamics (CPMD) simulations are compared to those obtained from classical molecular dynamics (MD) simulations. We describe classical simulations performed with various water potential functions to study the impact of the chosen water potential on the predicted conformational preference of the glycosidic linkage of the carbohydrate in aqueous solution. In agreement with our recent studies, we find that results obtained with CPMD simulations differ from those obtained from classical simulations. In particular, this study shows that the trans (t) orientation of the glycosidic linkage of methyl-alpha-mannopyranoside is preferred over its gauche anticlockwise (g-) orientation in aqueous solution. CPMD simulations indicate that this preference is due to intermolecular hydrogen bonding with surrounding water molecules, whereas no such information could be demonstrated by classical MD simulations. This study emphasizes the importance of ab initio MD simulations for studying the structural properties of carbohydrates in aqueous solution.     

Preferred conformation of the glycosidic linkage of methyl-beta-mannose

The conformational preference of the glycosidic linkage of methyl-beta-mannose was studied in the gas phase and in aqueous solution by ab initio calculations, and by molecular dynamics (MD) and Car-Parrinello molecular dynamics (CPMD) simulations. MD simulations were performed with various water potential functions to study the impact of the chosen water potential on the predicted conformational preference of the glycosidic linkage of the carbohydrate in solution. This study shows that the trans (t) orientation of the glycosidic linkage of methyl-beta-mannose is preferred over its gauche clockwise (g+) orientation in solution. CPMD simulations clearly indicate that this preference is due to intermolecular hydrogen bonding with surrounding water molecules, whereas no such information could be demonstrated by MD simulations. This study demonstrates the importance of ab initio molecular dynamics simulations in studying the structural properties of carbohydrate-water interactions.     

Infrared multiple photon dissociation spectroscopy of cationized histidine: effects of metal cation size on gas-phase conformation

The gas phase structures of cationized histidine (His), including complexes with Li(+), Na(+), K(+), Rb(+), and Cs(+), are examined by infrared multiple photon dissociation (IRMPD) action spectroscopy utilizing light generated by a free electron laser, in conjunction with quantum chemical calculations. To identify the structures present in the experimental studies, measured IRMPD spectra are compared to spectra calculated at B3LYP/6-311+G(d,p) (Li(+), Na(+), and K(+) complexes) and B3LYP/HW*/6-311+G(d,p) (Rb(+) and Cs(+) complexes) levels of theory, where HW* indicates that the Hay-Wadt effective core potential with additional polarization functions was used on the metals. Single point energy calculations were carried out at the B3LYP, B3P86, and MP2(full) levels using the 6-311+G(2d,2p) basis set. On the basis of these experiments and calculations, the only conformation that reproduces the IRMPD action spectra for the complexes of the smaller alkali metal cations, Li(+)(His) and Na(+)(His), is a charge-solvated, tridentate structure where the metal cation binds to the backbone carbonyl oxygen, backbone amino nitrogen, and nitrogen atom of the imidazole side chain, [CO,N(α),N(1)], in agreement with the predicted ground states of these complexes. Spectra of the larger alkali metal cation complexes, K(+)(His), Rb(+)(His), and Cs(+)(His), have very similar spectral features that are considerably more complex than the IRMPD spectra of Li(+)(His) and Na(+)(His). For these complexes, the bidentate [CO,N(1)] conformer in which the metal cation binds to the backbone carbonyl oxygen and nitrogen atom of the imidazole side chain is a dominant contributor, although features associated with the tridentate [CO,N(α),N(1)] conformer remain, and those for the [COOH] conformer are also clearly present. Theoretical results for Rb(+)(His) and Cs(+)(His) indicate that both [CO,N(1)] and [COOH] conformers are low-energy structures, with different levels of theory predicting different ground conformers.     

Molecular structures of cation...pi(arene) interactions for alkali metals with pi- and sigma-modalities

The monovalent cations of Na(+), K(+), Rb(+), and Cs(+) derived from the highly electropositive alkali metals represent prototypical charged spheres that are mainly subject to relatively simple electrostatic and solvation (hydration) forces. We now find that the largest of these Rb(+) and Cs(+) are involved in rather strong cation...pi(arene) interactions when they are suitably disposed with the ambifunctional hexasubstituted benzene C(6)E(6). The ether tentacles (E = methoxymethyl) allow these cations to effect eta(1)-bonding to the benzene center in a manner strongly reminiscent of the classical sigma-arene complexes with positively charged electrophiles where Z(+) = CH(3)(+), Br(+), Cl(+), Et(3)Si(+), etc. The somewhat smaller potassium cation is involved in a similar M(+)...pi(arene) interaction that leads to eta(2)-bonding with the aromatic center in the pi-mode previously defined in the well-known series of silver(I)/arene complexes. We can find no evidence for significant Na(+)... pi(arene) interaction under essentially the same conditions. As such, the sigma-structure of the Rb(+) and Cs(+) complexes and pi-structure of the K(+) complex are completely integrated into the continuum of sigma-pi bondings of various types of electrophilic (cationic) acceptors with arene donors that were initially identified by Mulliken as charge-transfer.     

Enhanced Li+ binding energies in alkylbenzene derivatives: the scorpion effect

The gas-phase lithium cation basicities (LCBs; Gibbs free energy of binding) of ethyl-, n-butyl-, and n-heptylbenzene have been measured by Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometry. The structures of the corresponding complexes and their relative stabilities were investigated through the use of B3LYP/6-311G(+)(3df,2p)//B3LYP/6-31G(d) density functional theory calculations. For n-butylbenzene and n-heptylbenzene, the most stable adducts correspond to pi complexes in which the alkyl chain coils toward the aromatic ring to favor its interaction with the metal cation. The extra stabilization provided by the flexible alkyl chain polarized by the charge on Li(+) is named the "scorpion effect". Conversely, these coiled conformations are among the least stable in the neutral system; they are not all stationary points on the potential-energy surface. The formation of complexes with a coiled alkyl chain leads to a significant enhancement of the Li(+) bonding energies (LBEs), which are approximately 20-30 kJ mol(-1) higher than those calculated for alkylbenzene pi complexes in which an uncoiled chain remains distant from the cation and thus minimizes the scorpion effect. This enhancement is less significant when LCBs are concerned, because the scorpion effect is entropically disfavored. There is very good agreement between the experimental Li(+) gas-phase basicities and the calculated values, provided that the statistical distribution of the conformers present in the gas phase is taken into account in this calculation.     

Monovalent ion adsorption at the muscovite (001)-solution interface: relationships among ion coverage and speciation, interfacial water structure, and substrate relaxation

The interfacial structure between the muscovite (001) surface and aqueous solutions containing monovalent cations (3 × 10(-3) m Li(+), Na(+), H(3)O(+), K(+), Rb(+), or Cs(+), or 3 × 10(-2) m Li(+) or Na(+)) was measured using in situ specular X-ray reflectivity. The element-specific distribution of Rb(+) was also obtained with resonant anomalous X-ray reflectivity. The results demonstrate complex interdependencies among adsorbed cation coverage and speciation, interfacial hydration structure, and muscovite surface relaxation. Electron-density profiles of the solution near the surface varied systematically and distinctly with each adsorbed cation. Observations include a broad profile for H(3)O(+), a more structured profile for Li(+) and Na(+), and increasing electron density near the surface because of the inner-sphere adsorption of K(+), Rb(+), and Cs(+) at 1.91 ± 0.12, 1.97 ± 0.01, and 2.26 ± 0.01 ?, respectively. Estimated inner-sphere coverages increased from ~0.6 to 0.78 ± 0.01 to ~0.9 per unit cell area with decreasing cation hydration strength for K(+), Rb(+), and Cs(+), respectively. Between 7 and 12% of the Rb(+) coverage occurred as an outer-sphere species. Systematic trends in the vertical displacement of the muscovite lattice were observed within ~40 ? of the surface. These include a <0.1 ? shift of the interlayer K(+) toward the interface that decays into the crystal and an expansion of the tetrahedral-octahedral-tetrahedral layers except for the top layer in contact with solution. The distortion of the top tetrahedral sheet depends on the adsorbed cation, ranging from an expansion (by ~0.05 ? vertically) in 3 × 10(-3)m H(3)O(+) to a contraction (by ~0.1 ?) in 3 × 10(-3) m Cs(+). The tetrahedral tilting angle in the top sheet increases by 1 to 4° in 3 × 10(-3) m Li(+) or Na(+), which is similar to that in deionized water where the adsorbed cation coverages are insufficient for full charge compensation.     

Alternative low-energy mechanisms for isotopic exchange in gas-phase D2O-H+(H2O)n reactions.

Molecular-dynamics (MD) trajectories and high-level ab initio methods have been used to study the low-energy mechanism for D(2)O-H(+)(H(2)O)(n) reactions. At low collisional energies, MD simulations show that the collisional complexes are long-lived and undergo fast monomolecular isomerization, converting between different isomers within 50-500 ps. Such processes, primarily involving water-molecule shifts along a water chain, require the surmounting of very-low-energy barriers and present sizable non- Rice-Ramsperger-Kassel-Marcus (RRKM) effects, which are interpreted as a lack of randomization of the internal kinetic energy. Interestingly, the rate of water shifts was found to increase upon increasing the size of the cluster. Based on these findings, we propose to incorporate the following steps into the mechanism for low-energy isotopic scrambling these D(2)O-H(+)(H(2)O)(n) reactions: a) formation of the collisional complex [H(+)(H(2)O)(n)D(2)O]* in a vibro-rotational excited state; b) incorporation of the heavy-water molecule in the cluster core as HD(2)O(+) by means of isomerization involving molecular shifts; c) displacement of solvation molecules from the first shell of HD(2)O(+) inducing de-deuteration (shift of a D(+) to a neighbor water molecule); d) reorganization of the clusters and/or expulsion of one of the isotopic variants of water (H(2)O, HDO or D(2)O) from the periphery of the complex.     

Sodium cation affinities of MALDI matrices determined by guided ion beam tandem mass spectrometry: application to benzoic acid derivatives

Threshold collision-induced dissociation of Na(+)(xBA) complexes with Xe is studied using guided ion beam mass spectrometry. The xBA ligands studied include benzoic acid and all of the mono- and dihydroxy-substituted benzoic acids: 2-, 3-, and 4-hydroxybenzoic acid and 2,3-, 2,4-, 2,5-, 2,6-, 3,4-, and 3,5-dihydroxybenzoic acid. In all cases, the primary product corresponds to endothermic loss of the intact xBA ligand. The cross section thresholds are interpreted to yield 0 and 298 K bond dissociation energies (BDEs) for Na(+)-xBA after accounting for the effects of multiple ion-neutral collisions, internal and kinetic energy distributions of the reactants, and dissociation lifetimes. Density functional theory calculations at the B3LYP/6-31G* level of theory are used to determine the structures of these complexes and provide the molecular constants necessary for the thermodynamic analysis of the experimental data. Theoretical BDEs are determined at the B3LYP/6-311+G(2d,2p) and MP2(full)/6-311+G(2d,2p) levels using the B3LYP/6-31G* optimized geometries. The trends in the measured BDEs suggest two very different binding modes for the Na(+)(xBA) complexes, while theory finds four. In general, the most stable binding conformation involves the formation of a six-membered chelation ring via interaction with the carbonyl and 2-hydroxyl oxygen atoms. The ground state geometries of the Na(+)(xBA) complexes in which the ligand does not possess a 2-hydroxyl group generally involve binding of Na(+) to either the carbonyl oxygen atom or to both oxygen atoms of the carboxylic acid group. These binding modes tend to be competitive because the enhancement in binding associated with the chelation interactions in the latter is mediated by steric repulsion between the hydroxyl and ortho hydrogen atoms. When possible, hydrogen bonding interactions with the ring hydroxyl group(s) enhance the stability of these complexes. The agreement between the theoretical and experimental BDEs is quite good for B3LYP and somewhat less satisfactory for MP2(full).     

A computational study on the interaction of the nitric oxide ions NO+ and NO- with the side groups of the aromatic amino acids

The interaction of the nitric oxide ions NO+ and NO- with benzene (C6H6) and the aromatic R-groups of the amino acids phenylalanine (Phe), tyrosine (Tyr), histidine (His), and tryptophan (Trp) have been examined using the DFT method B3LYP and the conventional electron correlation method MP2. In particular, the structures and complexation energies of the resulting half-sandwich Ar...NO+/- and sandwich [Ar...NO...Ar]+/- complexes have been considered. For the Ar...NO+ complexes, the presence of an electron rich heteroatom within or attached to the ring is found to not preclude the cation...pi bound complex from being the most stable. Furthermore, unlike the anionic complexes, the pi...cation...pi ([Ar...NO...Ar]+) complexes do not correspond to a "doubling" of the parent half-sandwich.     

Chiral recognition between a substituted cyclooctatetraene dianion and a half crown ether substituted with ibuprofen

(S)-Verbenol was substituted onto cyclooctatetraene (COT) via an ether linkage. In tetrahydrofuran (THF), Cs(+) or Na(+) counterions are tightly ion associated with the verbenoxy-COT dianion. A cosolvent, consisting of an ibuprofen unit connected to a half crown ether, was added to the verbenoxy-COT(2)(-),M(+)(2) solutions. The intimate interaction between the chiral cosolvent (ibuprofoxymethoxyethoxyethane) and the ion-associated counterion (either Na(+) or Cs(+)) forces a chiral recognition between the verbenoxy moiety and the ibuprofoxy moiety. When a molar excess of the cosolvent is present in the dianion THF solution, separation of the cosolvent associated with the verbenoxy-COT(2)(-),M(+)(2) complex from the uncomplexed cosolvent allows partial resolution of the enantiomers of ibuprofoxymethoxyethoxyethane.     

Hydration of potassiated amino acids in the gas phase

The thermochemistry of stepwise hydration of several potassiated amino acids was studied by measuring the gas-phase equilibria, AAK(+)(H(2)O)(n-1) + H(2)O = AAK(+)(H(2)O)(n) (AA = Gly, AL, Val, Met, Pro, and Phe), using a high-pressure mass spectrometer. The AAK(+) ions were obtained by electrospray and the equilibrium constants K(n-1,n) were measured in a pulsed reaction chamber at 10 mbar bath gas, N(2), containing a known partial pressure of water vapor. Determination of the equilibrium constants at different temperatures was used to obtain the DeltaH(n)(o), DeltaS(n)(o), and DeltaG(n)(o) values. The results indicate that the water binding energy in AAK(+)(H(2)O) decreases as the K(+) affinity to AA increases. This trend in binding energies is explained in terms of changes in the side-chain substituent, which delocalize the positive charge from K(+) to AA in AAK(+) complexes, varying the AAK(+)-H(2)O electrostatic interaction.     

Infrared predissociation spectroscopy of M+ (C6H6)(1-4)(H2O)(1-2)Ar(0-1) cluster ions, M = Li, Na

Infrared predissociation (IRPD) spectra of Li(+)(C(6)H(6))(1-4)(H(2)O)(1-2)Ar(0-1) and Na(+)(C(6)H(6))(2-4)(H(2)O)(1-2)Ar(1) are presented along with ab initio calculations. The results indicate that the global minimum energy structure for Li(+)(C(6)H(6))(2)(H(2)O)(2) has each water forming a π-hydrogen bond with the same benzene molecule. This bonding motif is preserved in Li(+)(C(6)H(6))(3-4)(H(2)O)(2)Ar(0-1) with the additional benzene ligands binding to the available free OH groups. Argon tagging allows high-energy Li(+)(C(6)H(6))(2-4)(H(2)O)(2)Ar isomers containing water-water hydrogen bonds to be trapped and detected. The monohydrated, Li(+) containing clusters contain benzene-water interactions with varying strength as indicated by shifts in OH stretching frequencies. The IRPD spectra of M(+)(C(6)H(6))(1-4)(H(2)O)(1-2)Ar are very different for lithium-bearing versus sodium-bearing cluster ions emphasizing the important role of ion size in determining the most favorable balance of competing noncovalent interactions.     

Controlling cesium cation recognition via cation metathesis within an ion pair receptor

Ion pair receptor 3 bearing an anion binding site and multiple cation binding sites has been synthesized and shown to function in a novel binding-release cycle that does not necessarily require displacement to effect release. The receptor forms stable complexes with the test cesium salts, CsCl and CsNO(3), in solution (10% methanol-d(4) in chloroform-d) as inferred from (1)H NMR spectroscopic analyses. The addition of KClO(4) to these cesium salt complexes leads to a novel type of cation metathesis in which the "exchanged" cations occupy different binding sites. Specifically, K(+) becomes bound at the expense of the Cs(+) cation initially present in the complex. Under liquid-liquid conditions, receptor 3 is able to extract CsNO(3) and CsCl from an aqueous D(2)O layer into nitrobenzene-d(5) as inferred from (1)H NMR spectroscopic analyses and radiotracer measurements. The Cs(+) cation of the CsNO(3) extracted into the nitrobenzene phase by receptor 3 may be released into the aqueous phase by contacting the loaded nitrobenzene phase with an aqueous KClO(4) solution. Additional exposure of the nitrobenzene layer to chloroform and water gives 3 in its uncomplexed, ion-free form. This allows receptor 3 to be recovered for subsequent use. Support for the underlying complexation chemistry came from single-crystal X-ray diffraction analyses and gas-phase energy-minimization studies.     

Cation [M = H+, Li+, Na+, K+, Ca2+, Mg2+, NH4+, and NMe4+] interactions with the aromatic motifs of naturally occurring amino acids: a theoretical study

Ab initio (HF, MP2, and CCSD(T)) and DFT (B3LYP) calculations were done in modeling the cation (H(+), Li(+), Na(+), K(+), Ca(2+), Mg(2+), NH(4)(+), and NMe(4)(+)) interaction with aromatic side chain motifs of four amino acids (viz., phenylalanine, tyrosine, tryptophan and histidine). As the metal ion approaches the pi-framework of the model systems, they form strongly bound cation-pi complexes, where the metal ion is symmetrically disposed with respect to all ring atoms. In contrast, proton prefers to bind covalently to one of the ring carbons. The NH(4)(+) and NMe(4)(+) ions have shown N-H...pi interaction and C-H...pi interaction with the aromatic motifs. The interaction energies of N-H...pi and C-H...pi complexes are higher than hydrogen bonding interactions; thus, the orientation of aromatic side chains in protein is effected in the presence of ammonium ions. However, the regioselectivity of metal ion complexation is controlled by the affinity of the site of attack. In the imidazole unit of histidine the ring nitrogen has much higher metal ion (as well as proton) affinity as compared to the pi-face, facilitating the in-plane complexation of the metal ions. The interaction energies increase in the order of 1-M < 2-M < 3-M < 4-M < 5-M for all the metal ion considered. Similarly, the complexation energies with the model systems decrease in the following order: Mg(2+) > Ca(2+) > Li(+) > Na(+) > K(+) congruent with NH(4)(+) > NMe(4)(+). The variation of the bond lengths and the extent of charge transfer upon complexation correlate well with the computed interaction energies.     

Photoionization of C(4) molecular beam: ab initio calculations

Large computations are performed on the C(4) (+) cation in order to characterize its stable isomers and its lowest electronic excited states using configuration interaction methods and large basis sets. Several stable isomers are found including a linear C(4) (+)(l-C(4) (+)), a rhombic C(4) (+)(r-C(4) (+)) (or cyclic), and a branched (d-C(4) (+)) structure. Our calculations show a high density of electronic states for all of these isomers favoring their interactions. By combining the present ab initio data and those on neutral C(4), the l-C(4)(X)+hnu-->l-C(4) (+)(X(+))+e(-), d-C(4)(X)+hnu-->d-C(4) (+)(X(+))+e(-), and r-C(4)(X)+hnu-->r-C(4) (+)(X(+))+e(-) vertical photoionization transition energies are computed at 10.87, 10.92, and 10.77 eV, respectively. Photoionizing a C(4) molecular beam results on an onset at 10.4-10.5 eV and then to a linear increase of the signal due to the opening of several ionization channels involving most of the C(4) and C(4) (+) isomers and electronic states.