Hydrogen bonding mediated by key orbital interactions determines hydration enthalpy differences of phosphate water clusters

Structure and Ultrafast X-ray Diffraction of the Hydrated Metaphosphate

We study the pathway of metaphosphate hydration when a metaphosphate anion is dissolved in liquid water with an explicit water model. For this purpose, we propose a sequential Monte Carlo algorithm incorporated with the ab initio quantum mechanics/molecular mechanics (QM/MM) method, which can reduce the amount of ab initio QM/MM sampling while retaining the accuracy of the simulation. We demonstrate the numerical calculation of the standard enthalpy change for the successive addition reaction PO3-·2H2O + H2O ⇌ PO3-·3H2O in the liquid phase, which helps to clarify the hydration pathway of the metaphosphate. With the obtained hydrated structure of the metaphosphate anion, we perform ab initio calculations for its relaxation dynamics upon vibrational excitation and characterize the energy transfer process in solution with simulated ultrafast X-ray diffraction signals, which can be experimentally implemented with X-ray free-electron lasers.

Hyperconjugation-mediated solvent effects in phosphoanhydride bonds

Density functional theory and natural bond orbital analysis are used to explore the impact of solvent on hyperconjugation in methyl triphosphate, a model for "energy rich" phosphoanhydride bonds, such as found in ATP. As expected, dihedral rotation of a hydroxyl group vicinal to the phosphoanhydride bond reveals that the conformational dependence of the anomeric effect involves modulation of the orbital overlap between the donor and acceptor orbitals. However, a conformational independence was observed in the rotation of a solvent hydrogen bond. As one lone pair orbital rotates away from an optimal antiperiplanar orientation, the overall magnitude of the anomeric effect is compensated approximately by the other lone pair as it becomes more antiperiplanar. Furthermore, solvent modulation of the anomeric effect is not restricted to the antiperiplanar lone pair; hydrogen bonds involving gauche lone pairs also affect the anomeric interaction and the strength of the phosphoanhydride bond. Both gauche and anti solvent hydrogen bonds lengthen nonbridging O-P bonds, increasing the distance between donor and acceptor orbitals and decreasing orbital overlap, which leads to a reduction of the anomeric effect. Solvent effects are additive with greater reduction in the anomeric effect upon increasing water coordination. By controlling the coordination environment of substrates in an active site, kinases, phosphatases, and other enzymes important in metabolism and signaling may have the potential to modulate the stability of individual phosphoanhydride bonds through stereoelectronic effects.

Influence of the molecular environment on phosphorylated amino acid models: a density functional theory study

A protein environment can affect the structure and charge distribution of substrate molecules. Here, the structure and partial charges were studied for different phosphorylated amino acid models in varying environments using density functional theory. The three systems investigated, acetyl phosphate, methyl phosphate, and p-tolyl phosphate are representative models for aspartyl phosphate, serine or threonine phosphate, and tyrosine phosphate, respectively. Combined with the CPCM continuum model, explicit HF and H(2)O molecules were added in order to model environmental effects and interactions that may occur in a protein matrix. We show how the different interactions affect the scissile P-O(R) bond and that the elongation can be explained by an anomeric effect. An increasing scissile bond length will result in transfer of negative charge to the leaving group and in a widening of the angle between the terminal oxygens of the phosphate molecule, features that can expose the phosphate group to attacking nucleophiles. Lastly, calculations were performed on the active site of the Ca(2+)-ATPase E2P intermediate, which provide an example of how a protein environment facilitates the formation of a destabilized ground state.

Systematic study of hydration patterns of phosphoric(V) acid and its mono-, di-, and tripotassium salts in aqueous solution

Fourier transform infrared (FTIR) spectroscopy of the OD band of HDO molecules has been applied to perform a systematic study of various phosphate forms in the order of decreasing protonation: H3PO4, KH2PO4, K2HPO4, K3PO4. HDO isotopically diluted in H2O has been prepared by adding adequate amounts of D2O to aqueous solutions in ordinary water. The difference spectra procedure has been applied to remove the contribution of bulk water and thus to separate the spectra of solute-affected HDO. The position at maximum of the principal anion-affected HDO band for potassium phosphates moves in the order KH2PO4 (2478 cm(-1))>K2HPO4 (2363 cm(-1))>K3PO4 (2301 cm(-1)), that is, decreases with increasing solute basicity and charge. The number of moles of water affected by one mole of solute (N) equals 11.0, 13.8 and 16.2, respectively. Phosphoric acid affects statistically 13.9 water molecules and appears to be a "structure making" solute in water. The isotopic substitution with deuterium occurs also on the phosphate anions and phosphoric acid. The thus formed P-O-D groups interact with water molecules via strong hydrogen bonds and the relative strength of this interaction increases with increasing solute acidity. The plausible assignments of OD bands of HDO have been confirmed by calculating equilibrium structures of small aqueous clusters of the studied individual utilizing density functional theory. Further interpretation of the energetic and structural properties of hydrating water is enabled by calculating intermolecular interaction energy of water and probability distributions for interatomic oxygen-oxygen distance.