Multiply charged anions in the gas phase

Molecular anions

The experimental and theoretical study of molecular anions has undergone explosive growth over the past 40 years. Advances in techniques used to generate anions in appreciable numbers as well as new ion-storage, ion-optics, and laser spectroscopic tools have been key on the experimental front. Theoretical developments on the electronic structure and molecular dynamics fronts now allow one to achieve higher accuracy and to study electronically metastable states, thus bringing theory in close collaboration with experiment in this field. In this article, many of the experimental and theoretical challenges specific to studying molecular anions are discussed. Results from many research groups on several classes of molecular anions are overviewed, and both literature citations and active (in online html and pdf versions) links to numerous contributing scientists' Web sites are provided. Specific focus is made on the following families of anions: dipole-bound, zwitterion-bound, double-Rydberg, multiply charged, metastable, cluster-based, and biological anions. In discussing each kind of anion, emphasis is placed on the structural, energetic, spectroscopic, and chemical-reactivity characteristics that make these anions novel, interesting, and important.

Imaging intramolecular coulomb repulsions in multiply charged anions

The properties of multiply charged anions are dominated by intramolecular Coulomb repulsion (ICR). Using photoelectron imaging, we show the effect of ICR on photoelectron angular distributions for a series of dianions, -O2C(CH2)_{n}CO_{2};{-} (D_{n};{2-}). The observed photoemission band of D_{n};{2-} was due to a perpendicular transition from the charged end group. However, photoemission intensities were observed to peak along the laser polarization for smaller n due to the strong ICR that forces electrons to be emitted along the molecular axis. This emission pattern weakens with increasing n and at D112- the angular distribution reverses back to peak at the perpendicular direction due to the reduced ICR.

Stabilization of amino acid zwitterions with varieties of anionic species: the intrinsic mechanism

With high-level ab initio theoretical methods, varieties of novel anions (two monoanions and dianions with two binding sites) were explored to stabilize the glycine zwitterions. Unlike the malonic and oxalic dianions (J. Am. Chem. Soc. 2005, 127, 13098.), the presently found anions are self-stable and widely available, ensuring the direct and convenient applications to stabilize the zwitterions. Some of the complexes formed with the anions and glycine zwitterions have very large vertical dissociation energies (>500.0 kJ mol (-1)), implying that the contained anions can be used to stabilize much more unstable zwitterions than the glycine zwitterions. Further studies revealed that the stabilization effects are closely related with the proton-capturing capacities of the anions. In order to stabilize the glycine zwitterions, the proton affinity (PA) of the anionic species should fall within the range of 1567.0-1983.6 kJ mol (-1) and, meanwhile, the proton-affinity differences of the two binding sites (DeltaPA) should be less than 104.2 kJ mol (-1). The present results can be used to direct the efficient designs of the stabilizers to other amino acid zwitterions as well as other types of zwitterions. In addition, the density functional theory was used and compared with the default MP2 theory, with the details given in the discussions.

Dianions of 7,7,8,8-tetracyano-p-quinodimethane and perfluorinated tetracyanoquinodimethane: Information on excited states from lifetime measurements in an electrostatic storage ring and optical absorption spectroscopy

We have developed an experimental technique that allows us to study the physics of short lived molecular dianions in the gas phase. It is based on the formation of monoanions via electrospray ionization, acceleration of these ions to keV energies, and subsequent electron capture in a sodium vapor cell. The dianions are stored in an electrostatic ion storage ring in which they circulate with revolution times on the order of 100 micros. This enables lifetime studies in a time regime covering five orders of magnitude, 10(-5)-1 s. We have produced dianions of 7,7,8,8-tetracyano-p-quinodimethane and 2,3,5,6-tetrafluoro-7,7,8,8-tetracyano-p-quinodimethane (TCNQ-F(4)) and measured their lifetimes with respect to electron autodetachment. Our data indicate that most of the dianions were initially formed in electronically excited states in the electron transfer process. Two levels of excitation were identified by spectroscopy on the dianion of TCNQ-F(4), and the absorption spectrum was compared with spectra obtained from spectroelectrochemistry of TCNQ-F(4) in acetonitrile solution.

Dimensional scaling treatment of stability of atomic anions induced by superintense, high-frequency laser fields

We show that dimensional scaling, combined with the high-frequency Floquet theory, provides useful means to evaluate the stability of gas phase atomic anions in a superintense laser field. At the large-dimension limit (D-->infinity), in a suitably scaled space, electrons become localized along the polarization direction of the laser field. We find that calculations at large D are much simpler than D=3, yet yield similar results for the field strengths needed to bind an "extra" one or two electrons to H and He atoms. For both linearly and circularly polarized laser fields, the amplitude of quiver motion of the electrons correlates with the detachment energy. Despite large differences in scale, this correlation is qualitatively like that found between internuclear distances and dissociation energies of chemical bonds.

Stabilization of excess charge in isolated adenosine 5'-triphosphate and adenosine 5'-diphosphate multiply and singly charged anions

Multiply charged anions (MCAs) represent highly energetic species in the gas phase but can be stabilized through formation of molecular clusters with solvent molecules or counterions. We explore the intramolecular stabilization of excess negative charge in gas-phase MCAs by probing the intrinsic stability of the [adenosine 5'-triphosphate-2H](2-) ([ATP-2H](2-)), [adenosine 5'-diphosphate-2H](2-) ([ADP-2H](2-)), and H(3)P(3)O(10)(2-) dianions and their protonated monoanionic analogues. The relative activation barriers for decay of the dianions via electron detachment or ionic fragmentation are investigated using resonance excitation of ions isolated within a quadrupole trap. All of the dianions decayed via ionic fragmentation demonstrating that the repulsive Coulomb barriers (RCB) for ionic fragmentation lie below the RCBs for electron detachment. Both the electrospray ionization mass spectra (ESI-MS) and total fragmentation energies for [ATP-2H](2-), [ADP-2H](2-), and H(3)P(3)O(10)(2-) indicate that the multiply charged H(3)P(3)O(10)(2-) phosphate moiety is stabilized by the presence of the adenosine group and the stability of the dianions increases in the order H(3)P(3)O(10)(2-) < [ADP-2H](2-) < [ATP-2H](2-). Fully optimized, B3LYP/6-31+G* minimum energy structures illustrate that the excess charges in all of the phosphate anions are stabilized by intramolecular hydrogen bonding either within the phosphate chain or between the phosphate and the adenosine. We develop a model to illustrate that the relative magnitudes of the RCBs and hence the stability of these ions is dominated by the extent of intramolecular hydrogen bonding.

Computational study on the negative electron affinities of NO2−∙(H2O)n clusters (n=0–30)

Here we report negative electron affinities of NO(2)(-).(H2O)n clusters (n=0-30) obtained from density functional theory calculations and a simple correction to Koopmans' theorem. The method relies on the calculation of the detachment energy of the monoanion and its highest occupied molecular orbital and lowest unoccupied molecular orbital energies, and explicit calculations on the dianion itself are avoided. A good agreement with resonances in the cross section for neutral production in electron scattering experiments is found for n=0, 1, and 2. We find several isomeric structures of NO(2)(-).(H2O)2 of similar energy that elucidate the interplay between water-water and ion-water interactions. The topology is predicted to influence the electron affinity by 0.5 and 0.4 eV for NO(2)(-).(H2O) and NO(2)(-).(H2O)2, respectively. The electron affinity of larger clusters is shown to follow a (n+delta)-1/3 dependence, where delta=3 represents the number of water molecules that in volume, could replace NO(2) (-).

PtF62− dianion and its detachment spectrum: A fully relativistic study

In this work we calculate the photoelectron spectrum of the PtF(6)2- dianion by application of the third-order Dirac-Hartree-Fock one-particle propagator technique. Relativistic effects and electron correlation are hereby treated on a consistent theoretical basis which is mandatory for systems containing heavy elements. A PtF6(2-) gas phase photoelectron spectrum is not yet available and our calculations therefore have predictive character. As it is characteristic for dianionic systems a strong dependence on basis set size and molecular geometry is observed. In contrast to the already calculated PtCl(6)2- photoelectron spectrum no valence orbital inversion due to strong interplay of spin-orbit coupling and electron correlation is observed. Furthermore an unusually strong spin-orbit splitting was found for the sigma-type subvalence 1t1u molecular spinor despite its very small platinum p population. The double ionization threshold is strongly lowered by relativistic effects now enabling an interatomic Coulombic decay process after ionization from the sigma-bonding orbitals. The results stress the importance of spin-orbit coupling for the understanding of the spectral structure which cannot be reproduced by a scalar-relativistic treatment only.

Base-dependent electron photodetachment from negatively charged DNA strands upon 260-nm laser irradiation

DNA multiply charged anions stored in a quadrupole ion trap undergo one-photon electron ejection (oxidation) when subjected to laser irradiation at 260 nm (4.77 eV). Electron photodetachment is likely a fast process, given that photodetachment is able to compete with internal conversion or radiative relaxation to the ground state. The DNA [6-mer]3- ions studied here show a marked sequence dependence of electron photodetachment yield. Remarkably, the photodetachment yield (dG6 > dA6 > dC6 > dT6) is inversely correlated with the base ionization potentials (G < A < C < T). Sequences with guanine runs show increased photodetachment yield as the number of guanine increases, in line with the fact that positive holes are the most stable in guanine runs. This correlation between photodetachment yield and the stability of the base radical may be explained by tunneling of the electron through the repulsive Coulomb barrier. Theoretical calculations on dinucleotide monophosphates show that the HOMO and HOMO-1 orbitals are localized on the bases. The wavelength dependence of electron detachment yield was studied for dG63-. Maximum electron photodetachment is observed in the wavelength range corresponding to base absorption (260-270 nm). This demonstrates the feasibility of gas-phase UV spectroscopy on large DNA anions. The calculations and the wavelength dependence suggest that the electron photodetachment is initiated at the bases and not at the phosphates. This also indicates that, although direct photodetachment could also occur, autodetachment from excited states, presumably corresponding to base excitation, is the dominant process at 260 nm. Excited-state dynamics of large DNA strands still remains largely unexplored, and photo-oxidation studies on trapped DNA multiply charged anions can help in bridging the gap between gas-phase studies on isolated bases or base pairs and solution-phase studies on full DNA strands.

Characterizing the intrinsic stability of gas-phase clusters of transition metal complex dianions with alkali metal counterions: Counterion perturbation of multiply charged anions

The authors report the gas-phase generation and characterization of a series of cation-dianion clusters, e.g., M(+).PtCl(6) (2-), M(+).PtCl(4) (2-), M(+).Pt(CN)(6) (2-), and M(+).Pd(CN)(4) (2-), where M(+)=Na(+),K(+),Rb(+), as model systems for investigating gas-phase contact ionpairs. Low-energy collisional excitation of these systems isolated within a quadrupole ion trap reveals that the fragmentation products are determined by the dianion and are independent of the counterion. This indicates that cation-dianion clusters represent gaseous ion-pair complexes, in line with recent findings for K(+).Pt(CN)(n) (2-), n=4,6 [Burke et al., J. Chem. Phys. 125, 021105 (2006)]. The relative fragmentation energies of several cation-dianion systems are obtained as a function of the counterion to explore the nature of ion-pair binding. For most of the systems studied, e.g., M(+).PtCl(6) (2-), the fragmentation energy increases as the cation size decreases, in line with a simple electrostatic description of the cation-dianion binding. However, the M(+).Pt(CN)(4) (2-) clusters displayed the reverse trend with the fragmentation energy increasing as the cation size increases. Density functional theory calculations of the cation-dianion fragmentation potential energy surfaces reveal the existence of a novel double-minima surface, separated by a repulsive Coulomb barrierlike feature at short range. The experimentally observed trends in the fragmentation energies can be fully understood with reference to the computed surfaces, hence providing strong support for the existence of the double-minima surface.

Small gas-phase dianions of Zn3O42−, Zn4O52−, CuZn2O42−, Si2GeO62−, Ti2O52−and Ti3O72−

We have searched for new species of small oxygen-containing gas-phase dianions produced in a secondary ion mass spectrometer by Cs+ ion bombardment of solid samples with simultaneous exposure of their surfaces to O2 gas. The targets were a pure zinc metal foil, a copper-contaminated zinc-based coin, two silicon-germanium samples (Si(1-x)Ge(x)(with x= 6.5% or 27%)) and a piece of titanium metal. The novel dianions Zn3O(4)(2-), Zn4O(5)(2-), CuZn2O(4)(2-), Si2GeO(6)(2-), Ti2O(5)(2-) and Ti3O(7)(2-) have been observed at half-integer m/z values in the negative ion mass spectra. The heptamer dianions Zn3O(4)(2-) and Ti2O(5)(2-) have been unambiguously identified by their isotopic abundances. Their flight times through the mass spectrometer are approximately 20 micros and approximately 17 micros, respectively. The geometrical structures of the two heptamer dianions Ti2O(5)(2-), and Zn3O(4)(2-) are investigated using ab initio methods, and the identified isomers are compared to those of the novel Ge2O(5)(2-) and the known Si2O(5)(2-) and Be3O(4)(2-) dianions.

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.

Dissolving metal reduction of acetylenes: a computational study

The two-electron, two-proton reduction of alkynes to trans-alkenes has been studied computationally using the polarizable continuum regime to model liquid ammonia, the solvent in which such reductions are generally carried out. Two computational approaches have been used. In one, the energies of species (alkyne, radical anion, vinyl radical, vinyl anion, dianion, and alkene) that are implicated as possible participants in the reduction are obtained using high level ab initio single-point computations under the polarizable continuum model (PCM) conditions. In the other approach, the same species are surrounded by ten explicit ammonia molecules before undergoing the same single-point PCM analysis. It has been shown that the two methods provide nearly identical results in terms of relative energies. Other findings include the probable bent nature of the radical anion species in ammonia, the likelihood that the trans stereochemistry of the reduction is determined at the vinyl anion stage, and the elimination of a dianion as a possible species that determines the stereochemical result. Various observations relating the solvent effects of ammonia are made relative to known gas-phase properties of the species studied.

Hydrogen loss from nucleobase nitrogens upon electron attachment to isolated DNA and RNA nucleotide anions

Electron transfer to isolated nucleotide monoanions in collisions with Na vapor induces hydrogen loss from nitrogen of the transient nucleobase anion. The cross section for this process is linearly correlated with the number of N-H hydrogens and is highest for guanine. The process is much faster than microseconds since only dehydrogenated dianions survived for mass spectrometric detection. The lifetime of the adenosine 5(')-monophospate dianions was measured to be 0.2 ms in an electrostatic ion storage ring but also a longer-lived component with a lifetime of at least 10 ms was identified. Implications of dissociation along the N-H coordinate for a nucleotide in DNA are briefly discussed in terms of Watson-Crick base pairs.

Infrared spectroscopy of hydrated sulfate dianions

We report the first infrared spectra of multiply-charged anions in the gas phase. The spectra of SO(4) (2-)(H(2)O)(n), with n=3-24, show four main bands assigned to two vibrations of the dianionic core, the water bending mode, and solvent libration. The triply degenerate SO(4) (2-) antisymmetric stretch vibration probes the local solvent symmetry, while the solvent librational band is sensitive to the hydrogen bonding network. The spectra and accompanying electronic structure calculations indicate a highly symmetric structure for the n=6 cluster and closure of the first solvation shell at n=12.

Photoelectron Spectroscopy of Free Multiply Charged Keggin Anions α-[PM12O40]3- (M = Mo, W) in the Gas Phase

Two polyoxometalate Keggin-type anions, alpha-PM12O40(3-) (M = Mo, W), were transferred to the gas phase by electrospray; their electronic structure and stability were probed by photoelectron spectroscopy. These triply charged anions were found to be highly stable in the gas phase with large adiabatic electron detachment energies of 1.7 and 2.1 eV for M = Mo and W, respectively. The magnitude of the repulsive Coulomb barrier was measured as approximately 3.4 eV for both anions, providing an experimental estimate for the intramolecular Coulomb repulsion present in these highly charged anions. Density functional theory calculations were carried out and compared with the experimental data, providing insight into the electronic structure and valence molecular orbitals of the two Keggin anions. The calculations indicated that the highest occupied molecular orbital and other frontier orbitals for PM12O40(3-) are localized primarily on the mu2-oxo bridging ligands of the polyoxometalate framework, consistent with the reactivity on the mu2-oxo sites observed in solution. It was shown that the HOMO of PW12O40(3-) is stabilized relative to that of PMo12O40(3-) by approximately 0.35 eV. The experimental adiabatic electron detachment energies of PM12O40(3-) (i.e., the electron affinities of PM12O40(2-)) are combined with recent calculations on the proton affinity of PM12O40(3-) to yield O-H bond dissociation energies in PM12O39(OH)2- as approximately 5.1 eV.

Counter-ion perturbation of the fragmentation pathways of multiply charged anions: Evidence for exit channel complexes on the fragmentation potential energy surfaces

We report the first low-energy collisional excitation measurements and density functional theory calculations to characterize the ground state potential energy surfaces of contact ion-pair complexes that contain multiply charged anions (MCAs). Excitation of K+.Pt(CN)(4) (2-) and K+.Pt(CN)(6) (2-) result in fragmentation products associated with decay of the isolated constituent dianions, revealing that the ground state ion-pair surfaces are dominated by the intrinsic characteristics of the MCA. This observation is important since it indicates that counter-ion complexation only weakly perturbs the electronic structure of an MCA. For K+.Pt(CN)(4) (2-), where the Pt(CN)(4) (2-) dianion decays with production of two ionic fragments, we observe evidence for the existence of a novel exit-channel complex corresponding to a polar KCN salt unit bound to the Pt(CN)(3) (-) anion. The results described provide a basis for understanding the potential energy surfaces and fragmentation characteristics of other ion-pair complexes that involve MCAs.

Photoelectron Spectroscopy of the Bis(dithiolene) Anions [M(mnt)2]n- (M = Fe − Zn; n = 1, 2): Changes in Electronic Structure with Variation of Metal Center and with Oxidation

A detailed understanding of the electronic structures of transition metal bis(dithiolene) centers is important in the context of their interesting redox, magnetic, and optical properties. The electronic structures of the series [M(mnt)2]n- (M = Fe - Zn; mnt = 1,2-S2C2(CN)2; n = 1, 2) were examined by a combination of photodetachment photoelectron spectroscopy and density functional theory calculations, providing insights into changes in electronic structure with variation of the metal center and with oxidation. Significant changes were observed for the dianions [M(mnt)2]2- due to stabilization of the metal 3d levels from Fe to Zn and the transition from square-planar to tetrahedral coordination about the metal center (Fe-Ni, D(2h) --> Cu D2 --> Zn, D(2d). Changes with oxidation from [M(mnt)2]2- to [M(mnt)2]1- were largely dependent on the nature of the redox-active orbital in the couple [M(mnt)2](2-/1-). In particular, the first detachment feature for [Fe(mnt)2]2- originated from a metal-based orbital (Fe(II) --> Fe(III)) while that for [Fe(mnt)2]1- originated from a ligand-based orbital, a consequence of stabilization of Fe 3d levels in the latter. In contrast, the first detachment feature for both of [Ni(mnt)2]2- and [Ni(mnt)2]1- originated from the same ligand-based orbital in both cases, a result of occupied Ni 3d levels being stabilized relative those of Fe 3d and occurring below the highest energy occupied ligand-based orbital for both of [Ni(mnt)2]2- and [Ni(mnt)2]1- . The combined data illustrate the subtle interplay between metal- and ligand-based redox chemistry in these species and demonstrate changes in their electronic structures with variation of metal center, oxidation, and coordination geometry.

New stable multiply charged negative atomic ions in linearly polarized superintense laser fields

Singly charged negative atomic ions exist in the gas phase and are of fundamental importance in atomic and molecular physics. However, theoretical calculations and experimental results clearly exclude the existence of any stable doubly-negatively-charged atomic ion in the gas phase, only one electron can be added to a free atom in the gas phase. In this report, using the high-frequency Floquet theory, we predict that in a linear superintense laser field one can stabilize multiply charged negative atomic ions in the gas phase. We present self-consistent field calculations for the linear superintense laser fields needed to bind extra one and two electrons to form He-, He2-, and Li2-, with detachment energies dependent on the laser intensity and maximal values of 1.2, 0.12, and 0.13 eV, respectively. The fields and frequencies needed for binding extra electrons are within experimental reach. This method of stabilization is general and can be used to predict stability of larger multiply charged negative atomic ions.

Differences in electronic properties of fluorinated and trifluoromethylated fullerenes revealed by their propensity for dianion formation

The cross sections for electron transfer from sodium to C(60)F(n) (-) and C(60)(CF(3))(n) (-) anions in 50-keV collisions as a function of the number of functional groups n are reported. There are clear differences between derivatives of fluorine and trifluoromethyl due to the different electron withdrawing properties of F and CF(3). The role of inductive effects and pi electron delocalization on the electron affinity is discussed, assuming a correlation between the cross section and the electron affinity of the anion.

A Density Functional Study on Beryllium-Doped Carbon Dianion Clusters CnBe2- (n = 4−14)

Making use of the software of molecular graphics, we designed numerous models of C(n)()Be(2-) (n = 4-14). We carried out geometry optimization and calculation on vibration frequency by means of the B3LYP density functional method. After comparison of structure stability, we found that the ground-state isomers of C(n)()Be(2-) (n = 4-14) are linear with the beryllium atom located inside the C(n)() chain. When a side carbon chain is with an even number of carbon atoms, it is polyacetylene-like, whereas when a side chain is with an odd number of carbon atoms, it is cumulene-like. The C(n)Be(2-) (n = 4-14) clusters with an even number of carbon atoms are more stable than that with an odd number of carbon atoms, matching the peak pattern observed in accelerator mass spectrometry (AMS) and Coulomb Explosion Imaging (CEI) investigations of C(n)()Be(2-) (n = 4-14). The trend of such odd/even alternation is explained based on concepts of bonding characteristics, electronic configuration, electron detachment, and incremental binding energy.

Probing the Intrinsic Electronic Structure of the Bis(dithiolene) Anions [M(mnt)2]2- and [M(mnt)2]1- (M = Ni, Pd, Pt; mnt = 1,2-S2C2(CN)2) in the Gas Phase by Photoelectron Spectroscopy

A detailed understanding of the electronic structure of transition metal bis(dithiolene) complexes is important because of their interesting redox, magnetic, optical, and conducting properties and their relevance to enzymes containing molybdenum and tungsten bis(dithiolene) centers. The electronic structures of the bis(dithiolene) anions [M(mnt)(2)](n-) (M = Ni, Pd, Pt; mnt = 1,2-S(2)C(2)(CN)(2); n = 0-2) were examined by a combination of photodetachment photoelectron spectroscopy (PES) and density functional theory calculations. The combined experimental and theoretical data provide insight into the molecular orbital energy levels of [M(mnt)(2)](2-) and the ground and excited states of [M(mnt)(2)](1-) and [M(mnt)(2)]. Detachment features from ligand-based orbitals of [M(mnt)(2)](2-) occur at similar energies for each species, independent of the metal center, while those arising from metal-based orbitals occur at higher energies for the heavier congeners. Electronic excitation energies inferred for [M(mnt)(2)](1-) from the PES experiments agree well with those obtained in optical absorption experiments in solution, with the PES experiments providing additional insight into the changes in energy of these transitions as a function of metal. The singly charged anions [M(mnt)(2)](1-) were also prepared and studied independently. Electron detachment from the ground states of these doublet anions accessed the lowest singlet and triplet states of neutral [M(mnt)(2)], thereby providing a direct experimental measure of their singlet-triplet splitting.

Evaluation of Two-Center, Two-Electron Integrals

We present a new analytic treatment of two-electron integrals over two-center integrals including correlation (interelectronic distance) explicitly in the wave function. All the integrals needed for the evaluation of the matrix elements of any diatomic two-electron molecule are obtained as analytic recursion expressions. As an application of this method in molecular physics, we calculate the value of the ground-state energy and equilibrium internuclear distance of the hydrogen molecule in the Born-Oppenheimer approximation.

Mechanism of Triphosphate Hydrolysis in Aqueous Solution: QM/MM Simulations in Water Clusters

The mechanism of the hydrolysis reaction of the unprotonated methyl triphosphate (MTP) ester in water clusters has been modeled. The effective fragment potential based quantum mechanical-molecular mechanical (QM/MM) approach has been applied in the simulations. It is shown that the minimum energy reaction path is consistent with an assumption of a two-step dissociative-type process similar to the case of the guanosine triphosphate (GTP) hydrolysis in the Ras-GAP protein complex (Grigorenko, B. L.; Nemukhin, A. V.; Topol, I. A.; Cachau, R. E.; Burt, S. K. Proteins: Struct., Funct., Bioinf. 2005, 60, 495). At the first stage, a unified action of environmental molecular groups and the catalytic water molecule leads to a substantial spatial separation of the gamma-phosphate group from the rest of the molecule. At the second stage, inorganic phosphate H2PO4- is formed from water and the metaphosphate anion PO3- through the chain of proton transfers along hydrogen bonds. The estimated activation barriers for MTP in aqueous solution at both stages (20 and 14 kcal/mol) are substantially higher than the corresponding barriers for the GTP hydrolysis in the protein.