Toward an Understanding of the Hydrolysis of Mg−PPi. An ab Initio Study of the Isomerization Reactions of Neutral and Anionic Mg−Pyrophosphate Complexes

Structural characterization of nucleotide 5'-triphosphates by infrared ion spectroscopy and theoretical studies

The molecular family of nucleotide triphosphates (NTPs), with adenosine 5'-triphosphate (ATP) as its best-known member, is of high biochemical importance as their phosphodiester bonds form Nature's main means to store and transport energy. Here, gas-phase IR spectroscopic studies and supporting theoretical studies have been performed on adenosine 5'-triphosphate, cytosine 5'-triphosphate and guanosine 5'-triphosphate to elucidate the intrinsic structural properties of NTPs, focusing on the influence of the nucleobase and the extent of deprotonation. Mass spectrometric studies involving collision induced dissociation showed similar fragmentation channels for the three studied NTPs within a selected charge state. The doubly charged anions exhibit fragmentation similar to the energy-releasing hydrolysis reaction in nature, while the singly charged anions show different dominant fragmentation channels, suggesting that the charge state plays a significant role in the favorability of the hydrolysis reaction. A combination of infrared ion spectroscopy and quantum-chemical computations indicates that the singly charged anions of all NTPs are preferentially deprotonated at their β-phosphates, while the doubly-charged anions are dominantly αβ-deprotonated. The assigned three-dimensional structure differs for ATP and CTP on the one hand and GTP on the other, in the sense that ATP and CTP show no interaction between nucleobase and phosphate tail, while in GTP they are hydrogen bonded. This can be rationalized by considering the structure and geometry of the NTPs where the final three dimensional structure depends on a subtle balance between hydrogen bond strength, flexibility and steric hindrance.

Spectroscopic and theoretical investigations of adenosine 5'-diphosphate and adenosine 5'-triphosphate dianions in the gas phase

Doubly deprotonated adenosine 5'-diphosphate ([ADP-2H](2-)) and adenosine 5'-triphosphate ([ATP-2H](2-)) dianions were investigated using infrared multiple photon dissociation (IR-MPD) and photoelectron spectroscopy. Vibrational spectra acquired in the X-H stretch region (X = C, N, O) and augmented by isotope-labelling were compared to density functional theory (DFT) calculations at the B3LYP/TZVPP level. This suggests that in [ATP-2H](2-) the two phosphate groups adjacent to the ribose ring are preferentially deprotonated. Photoelectron spectra recorded at 4.66 and 6.42 eV photon energies revealed adiabatic detachment energies of 1.35 eV for [ADP-2H](2-) and 3.35 eV for [ATP-2H](2-). Repulsive Coulomb barriers were estimated at ~2.2 eV for [ADP-2H](2-) and ~1.9 eV for [ATP-2H](2-). Time-dependent DFT calculations have been used to simulate the photoelectron spectra. Photodetachment occurs primarily from lone pair orbitals on oxygen atoms within the phosphate chain.

Quantum and classical dynamics simulations of ATP hydrolysis in solution

ATP hydrolysis is a key reaction in living cells that drives many cellular processes. The reaction, which involves gamma phosphate cleavage from ATP, converting it to ADP, has been suggested to occur via an associative or dissociative mechanism dependent upon the surrounding environment. Prior quantum chemical studies suffered from short simulation timescales failing to capture free energy contributions due to relaxation of the surrounding aqueous environment. We have developed a highly parallelized QM/MM implementation in the NAMD and OpenAtom simulation packages, using the dual grid, dual length scale method for combined plane-wave and Eular exponential spline-based QM/MM simulations. This approach, using message-driven parallel quantum and classical dynamics, permits sufficient timescale simulations for quantum chemical events such as ATP hydrolysis, and is found to accurately and reliably include the free energy contributions of solvent relaxation to hydrolysis. In this paper we describe the application of the dual grid, dual length plane-wave-based QM/MM method to study both the associative and dissociative mechanisms of ATP hydrolysis, accounting for the free energy contribution from solvent relaxation, as well as for the key role of Mg(2+) in the reaction.

Where is electronic energy stored in adenosine triphosphate?

The gas-phase electronic energy of the hydrolysis of methyl triphosphate, a model of adenosine 5'-triphosphate (ATP), is partitioned into local (atomic and group) contributions. A modified definition of Lipmann's "group transfer potential" is proposed on the basis of the partitioning of the total electronic energy into atomic contributions within the framework of the quantum theory of atoms in molecules (QTAIM). The group transfer potential is defined here as the sum of the atomic energies forming the group in ATP minus the sum of the energies of the same atoms in inorganic phosphate. It is found that the transfer potential of the terminal phosphate group in ATP is significantly reduced, from +241.7 to +73.1 kcal/mol, as a result of complexation with magnesium. This is accompanied by a concomitant change in the energy of reaction from -168.6 to -24.9 kcal/mol. Regions within ATP where the electronic energy changes the most upon hydrolysis are identified. The study is conducted at the DFT/B3LYP/6-31+G(d,p) level of theory.

Atomic Partitioning of the Dissociation Energy of the P−O(H) Bond in Hydrogen Phosphate Anion (HPO42-): Disentangling the Effect of Mg2+

This paper has three goals: (1) to provide a first step in understanding the atomic basis of the role of magnesium in facilitating the dissociation of the P-O bond in phosphorylated biochemical fuel molecules (such as ATP or GTP), (2) to compare second-order Møller-Plesset perturbation theory (MP2) results with those obtained at the more economical density functional theory (DFT) level for a future study of larger more realistic models of ATP/GTP, and (3) to examine the calculation of atomic total energies from atomic kinetic energies within a Kohn-Sham implemention of DFT, as compared to ab initio methods. A newly described method based on the quantum theory of atoms in molecules (QTAIM), which is termed the "atomic partitioning of the bond dissociation energy" (APBDE), is applied to a simple model of phosphorylated biological molecules (HPO42-). The APBDE approach is applied in the presence and in the absence of magnesium. It is found that the P-O(H) bond in the magnesium complex is shorter, exhibits a higher stretching frequency, and has a higher electron density at the bond critical point than in the magnesium-free hydrogen phosphate anion. Though these data would seem to suggest a stronger P-O(H) bond in the magnesium complex compared to the magnesium-free case, the homolytic breaking of the P-O(H) bond in the complex is found to be easier, i.e., has a lower BDE. This effect is the result of the balance of several atomic contributions to the BDE induced by the magnesium cation, which stabilizes the dissociation product more than it stabilizes the intact model molecule.

New Nonsymmetric P(OH)3 Species. Comparison with the C3 Isomer and Themochemistry at the DFT, MP2, and CCSD(T) Levels of Theory

Two new less-symmetric P(OH)3 isomers that are more stable than the C3 structure are found at the density functional theory (B3PW91, B3LYP), MP2, and CCSD(T) levels with the large aug-cc-pvdz/pvtz basis sets. The C1 and C3 structures are qualitatively different from those found for the As(OH)3 molecule. An additional lower lying P(OH)3 structure with Cs symmetry has been obtained. With the largest basis set the Cs isomer is predicted to be the most stable. However, the inclusion of zero-point-energy corrections induces an inversion between the Cs and C1 isomers, with the latter becoming the lowest energy structure at the highest correlated level. Increasing inclusion of electronic correlation effects reduces the energy difference between the C1 and Cs structures while the C1-C3 energy difference and C1-Cs interconversion barrier become larger. In all cases, energy differences and barrier heights are around 1 kcal/mol.

Gas-Phase Chemistry of Diphosphate Anions as a Tool To Investigate the Intrinsic Requirements of Phosphate Ester Enzymatic Reactions: The [M1M2HP2O7]− Ions

Experimental studies on gaseous inorganic phosphate ions are practically nonexistent, yet they can prove helpful for a better understanding of the mechanisms of phosphate ester enzymatic processes. The present contribution extends our previous investigations on the gas-phase ion chemistry of diphosphate species to the [M(1)M(2)HP(2)O(7)](-) ions where M(1) and M(2) are the same or different and correspond to the Li, Na, K, Cs, and Rb cations. The diphosphate ions are formed by electrospray ionization of 10(-4) M solutions of Na(5)P(3)O(10) in CH(3)CN/H(2)O (1/1) and MOH bases or M salts as a source of M(+) cations. The joint application of mass spectrometric techniques and quantum-mechanical calculations makes it possible to characterize the gaseous [M(1)M(2)HP(2)O(7)](-) ions as a mixed ionic population formed by two isomeric species: linear diphosphate anion coordinated to two M(+) cations (group I) and [PO(3)M(1)M(2)HPO(4)](-) clusters (group II). The relative gas-phase stabilities and activation barriers for the isomerization I-->II, which depend on the nature of the M(+) cations, highlight the electronic susceptibility of P-O-P bond breaking in the active site of enzymes. The previously unexplored gas-phase reactivity of [M(1)M(2)HP(2)O(7)](-) ions towards alcohols of different acidity was investigated by Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR/MS). The reaction proceeds by addition of the alcohol molecule followed by elimination of a water molecule.

Phosphate transfer from inositol pyrophosphates InsP5PP and InsP4(PP)2: A semi-empirical investigation

A novel phosphate transfer process involving the non-enzymatic transfer of a phosphate group from inositol pyrophosphates to serine residues in proteins has been recently reported. Semi-empirical calculations at the PM3/SM5.2 level were undertaken to explore the effect of inositol pyrophosphate structure and overall charge on the thermodynamics of this phosphate transfer.

Molecular dynamics DFT:B3LYP study of guanosinetriphosphate conversion into guanosinemonophosphate upon Mg2+ chelation of alpha and beta phosphate oxygens of the triphosphate tail

A molecular dynamics DFT:B3LYP (6-31G(**) basis set) study is used to elucidate the mechanism of guanosinetriphosphate (GTP) conversion into guanosinemonophosphate (GMP) upon the action of Mg(2+) (magnesium cofactor). The computations are carried out at 310 K in a volume of 178 water molecules, which surround the Mg(2+)-GTP complex and imitate the effect of solution. Over 5 ps, Mg(2+)-GTP appears to be fully decomposed, yielding five final products: two hydrated molecules of inorganic phosphate Pi, a hydrated Mg(2+), atomic oxygen (which in the course of a couple of subsequent reactions gains two hydrogens and converts into a water molecule) and a highly active *GMP radical. The radical production is linked to presence of Mg(2+), which initiates a radical mechanism of GTP cleavage. At the initial stage, Mg(2+) undergoes reduction to Mg(+), accompanied by the formation of an ion-radical pair with GTP, (+)Mg*-*GTP(3-). Without Mg(2+), an inert form of GMP (the ionic mechanism of GTP hydrolytic cleavage) rather than GMP is produced. *GMP production, which is similar to that of *AMP (adenosinemonophosphate), *CMP (cytidinemonophosphate), TMP (thymidinemonophosphate) and *UMP (uridinemonophosphate), plays a crucial role in DNA and RNA single chain synthesis.

Effect of Alkali Metal Coordination on Gas-Phase Chemistry of the Diphosphate Ion: The MH2P2O7− Ions

Systematic experimental and theoretical studies on anionic phosphate species in the gas phase are almost nonexistent, even though they could provide a benchmark for enhanced comprehension of their liquid-phase chemical behavior. Gaseous MH(2)P(2)O(7) (-) ions (M=Li, Na, K, Rb, Cs), obtained from electrospray ionization of solutions containing H(4)P(2)O(7) and MOH or M salts as a source of M(+) ions were structurally assayed by collisionally activated dissociation (CAD) mass spectrometry and theoretical calculations at the B3LYP/6-31+G* level of theory. The joint application of mass spectrometric techniques and theoretical methods allowed the MH(2)P(2)O(7) (-) ions to be identified as having a structure in which the linear diphosphate anion is coordinated to the M(+) ion (I) and provides information on gas-phase isomerization processes in the [PO(3)...MH(2)PO(4)](-) clusters II and the [P(2)O(6)...M...H(2)O](-) clusters IV. Studies of gas-phase reactivity by Fourier transform ion cyclotron resonance (FTICR) and triple quadrupole (TQ) mass spectrometry revealed that the MH(2)P(2)O(7) (-) ions react with selected nucleophiles by clustering, proton transfer and addition-elimination mechanisms. The influence of the coordination of alkali metal ions on the chemical behavior of pyrophosphate is discussed.

Gaseous H5P2O8? Ions: A Theoretical and Experimental Study on the Hydrolysis and Synthesis of Diphosphate Ion

The structure and reactivity of gaseous H5P2O8- ions obtained from the chemical ionization (CI) of an H4P2O7/H2O mixture and from electrospray ionization (ESI) of CH3CN/H2O/H4P2O7 solutions were investigated by Fourier transform ion cyclotron (FTICR) and triple quadrupole mass spectrometry. Theoretical calculations performed at the B3LYP/6-31+G* level of theory and collisionally activated dissociation (CAD) mass spectrometric results allowed the ionic population obtained in the CI conditions to be structurally characterized as a mixture of gaseous [H3P2O7...H2O]-, [H3PO4...H2PO4]-, and [PO3...H3PO4...H2O]- clusters. The energy profile emerging from theoretical calculations affords insight into the mechanism of diphosphate ion hydrolysis and synthesis.

The structure and stability of biological metaphosphate, phosphate, and phosphorane compounds in the gas phase and in solution

Density functional calculations of a series of metaphosphates, acyclic and cyclic phosphates and phosphoranes relevant to RNA catalysis are presented. Solvent effects calculated with three well-established solvation models are analyzed and compared. The structure and stability of the compounds are characterized in terms of thermodynamic quantities for isomerization and ligand substitution reactions, gas-phase proton affinities, and microscopic solution pK(a)() values. The large dataset of compounds allows the estimation of bond energies to determine the relative strengths of axial and equatorial P-O phosphorane single bonds and P-O single and double bonds in metaphosphates and phosphates. The relative apicophilicty of hydroxyl and methoxy ligands in phosphoranes are characterized. The results presented here provide quantitative insight into RNA catalysis and serve as a first step toward the construction of a high-level quantum database for development of new semiempirical Hamiltonian models for biological reactions

Influence of key residues on the reaction mechanism of the cAMP-dependent protein kinase

The reaction mechanism of the catalytic phosphoryl transfer of cAMP-dependent protein kinase (cAPK) was investigated by semi-empirical AM1 molecular orbital computations of an active site model system derived from the crystal structure of the catalytic subunit of the enzyme. The activation barrier is calculated as 20.7 kcal mol(-1) and the reaction itself to be exothermic by 12.2 kcal mol(-1). The active site residue Asp166, which was often proposed to act as a catalytic base, does not accept a proton in any of the reaction steps. Instead, the hydroxyl hydrogen of serine is shifted to the simultaneously transferred phosphate group of ATP. Although the calculated transition state geometry indicates an associative phosphoryl transfer, no concentration of negative charge is found. To study the influence of protein mutations on the reaction mechanism, we compared two-dimensional energy hypersurfaces of the protein kinase wild-type model and a corresponding mutant in which Asp166 was replaced by alanine. Surprisingly, they show similar energy profiles despite the experimentally known decrease of catalytic activity for corresponding mutants. Furthermore, a model structure was examined, where the charged NH3 group of Lys168 was replaced by a neutral methyl group. The energetic hypersurface of this hypothetical mutant shows two possible pathways for phosphoryl transfer, which both require significantly higher activation energies than the other systems investigated, while the energetic stabilization of the reaction product is similar in all systems. As the position of the amino acid side chains and the substrate peptide is virtually unchanged in all model systems, our results suggest that the exchange of Asp166 by other amino acid is less important to the phosphoryl transfer itself, but crucial to maintain the configuration of the active site in vivo. The positively charged side chain of Lys168, however, is necessary to stabilize the intermediate reaction states, particularly the side chain of the substrate peptide.