A theoretical investigation on the conformation and the interaction of CHF2OCF2CHF2 (desflurane II) with one water molecule

The conformation and the interaction of CHF₂OCF₂CHF₂ (desflurane II) with one water molecule is investigated theoretically using the ab initio MP2/aug-cc-pvdz and DFT-based M062X/6-311++G(d,p) methods. The calculations include the optimized geometries, the harmonic frequencies of relevant vibrational modes along with a natural bond orbital (NBO) analysis including the NBO charges, the hybridization of the C atom and the intra- and intermolecular hyperconjugation energies. In the two most stable conformers, the CH bond of the F₂HCO- group occupies the gauche position. The hyperconjugation energies are about the same for both conformers and the conformational preference depends on the interaction between the non-bonded F and H atoms. The deprotonation enthalpies of the CH bonds are about the same for both conformers, the proton affinity of the less stable conformer being 3 kcal mol^?1 higher. Both conformers of desflurane II interact with water forming cyclic complexes characterized by CH…O and OH…F hydrogen bonds. The binding energies are moderate, ranging from ?2.4 to ?3.2 kcal mol^?1 at the MP2 level. The origin of the blue shifts of the ν(CH) vibrations is analyzed. In three of the complexes, the water molecule acts as an electron donor. Interestingly, in these cases a charge transfer is also directed to the non bonded OH group of the water molecule. This effect seems to be a property of polyfluorinated ethers.     

Intramolecular proton transfer in adenine imino tautomers

The electronic structural impact on intramolecular proton transfer in the cis- and trans-imino N7 and N9 tautomers of adenine (A) has been studied quantum mechanically, using density functional theory (B3LYP/TZVP, SAOP/TZ2P, LB94/TZ2P) and Green function (OVGF/TZVP) models. It is found that proton transfer does not significantly change isotropic properties but has profound impact on electron distributions of the species through anisotropic properties. The relative energies with respect to the canonical A tautomer (amino-9H), ΔE, for imino 7Hcis, imino 7Htrans, imino 9Hcis and imino 9Htrans are calculated as 16.15, 16.43, 18.46 and 13.80 kcal mol^? 1 (B3LYP/TZVP model) and some minor changes in perimeters of the purine ring is also observed. The Hirshfeld atomic charges indicate that whether a proton attached to N₍₇₎ or N₍₉₎ causes a significant local charge redistribution. However, these charges are insensitive to cis–trans proton transfer. Condensed Fukui function reveals N₍₁₀₎ and C₍₈₎ as the most electrophilic reactive site among N- and C-atom sites, respectively. We also found that proton transfer significantly alters in-plane σ orbitals, rather than out of plane π orbitals including the frontier orbital 6a″. Moreover, orbital based responses to various proton transfers are presented: the orbital 29a′ (HOMO-1) is a signature orbital differentiating all the four tautomers. Orbital 27a′ is a site (N₍₇₎ and N₍₉₎) sensitive orbital, whereas orbital 22a′ is only sensitive to proton orientation on the imino group = N–H.     

The physicochemical essence of the purine·pyrimidine transition mismatches with Watson-Crick geometry in DNA: A·C* versa A*·C. A QM and QTAIM atomistic understanding

It was established for the first time by DFT and MP2 quantum-mechanical (QM) methods either in vacuum, so in the continuum with a low dielectric constant (ε = 4), typical for hydrophobic interfaces of specific protein-nucleic acid interactions, that the repertoire for the tautomerisation of the biologically important adenine·cytosine* (A·C*) mismatched DNA base pair, formed by the amino tautomer of the A and the imino mutagenic tautomer of the C, into the A*·C base mispair (?G = 2.72 kcal?mol^?1 obtained at the MP2 level of QM theory in the continuum with ε = 4), formed by the imino mutagenic tautomer of the A and the amino tautomer of the C, proceeds via the asynchronous concerted double proton transfer along two antiparallel H-bonds through the transition state (TS_A·C*?A*·C). The limiting stage of the A·C*→A*·C tautomerisation is the final proton transfer along the intermolecular N6H···N4 H-bond. It was found that the A·C*/A*·C DNA base mispairs with Watson–Crick geometry are associated by the N6H?N4/N4H?N6, N3H?N1/N1H?N3 and C2H?O2 H-bonds, respectively, while the TS_A·C*?A*·C is joined by the N6–H–N4 covalent bridge and the N1H?N3 and C2H?O2 H-bonds. It was revealed that the A·C*?A*·C tautomerisation is assisted by the true C2H?O2 H-bond, that in contrast to the two others conventional H-bonds exists along the entire intrinsic reaction coordinate (IRC) range herewith becoming stronger at the transition from vacuum to the continuum with ε = 4. To better understand the behavior of the intermolecular H-bonds and base mispairs along the IRC of the A·C*?A*·C tautomerisation, the profiles of their electron-topological, energetical, geometrical, polar and charge characteristics are reported in this study. It was established based on the profiles of the H-bond energies that all three H-bonds are cooperative, mutually strengthening each other. The nine key points, providing a detailed physicochemical picture of the A·C*?A*·C tautomerisation, were revealed and thoroughly examined along the IRC. It was shown that the A*·C base mispair with the population ~1 % obtained at the MP2 level of QM theory in the continuum with ε = 4 is thermodynamically and dynamically stable structure. Its lifetime was calculated to be 5.76·10^?10 s at the MP2 level of QM theory in the continuum with ε = 4. This lifetime, from the one side, enables all six low-frequency intermolecular vibrations to develop, but, from the other side, it is by order less than the time (several ns) required for the replication machinery to forcibly dissociate a base pair into the monomers during DNA replication. This means that the A*·C base mispair “slips away from the hands” of the replication machinery into the A·C* mismatched base pair. Consequently, the authors came to the conclusion that exactly the A·C* base mispair is an active player of the point mutational events and is effectively dissociated by the replication machinery into the A and C* monomers in contrast to the A*·C base mispair, playing the mediated role of a provider of the A·C* base mispair in DNA that is synthesised.     

Kinetic and thermodynamic study of the substituent effect on the amino-Claisen rearrangement of para-substituted N-allyl-N-arylamine: a Hammett study via DFT

In order to find the susceptibility of the amino-Claisen rearrangement and the next proton shift reaction of N-allyl-N-arylamine to the substituent effects in the para position, the kinetic and thermodynamic parameters were calculated at the B3LYP level using the 6-31G^** basis set. The calculated activation energies for the rearrangements and proton shift reactions are close to 44.4 and 49.5 kcal mol^? 1, respectively. The transition states of the rearrangement with electron-donor substituents are more stable than those with electron-withdrawing substituent groups, but for the proton shift reaction, this situation is reversed (with the exception of fluorine atom for the rearrangement and fluorine and chlorine atoms for the proton shift reaction). Negative values for the activation entropy confirm the concerted mechanism for the amino-Claisen rearrangement and proton shift reaction. The Hammett ρ values of ? 2.4172 and ? 1.7791 are obtained for σ_p and σ^?  (enhanced sigma) in the amino-Claisen rearrangement, respectively. The correlation between log(k _X/k _H) and σ_p is weaker than that with σ^?  (enhanced sigma). A negative Hammett ρ value indicates that the electron-donating groups slightly increase the rate of amino-Claisen rearrangement. A positive Hammett ρ value (2.4921) for the proton shift reaction indicates that electron-withdrawing groups increase the rate of reaction.     

Prototropic tautomerism and basic molecular principles of hypoxanthine mutagenicity: an exhaustive quantum-chemical analysis

The molecular structures, relative stability order, and dipole moments of a complete family of 21 planar hypoxanthine (Hyp) prototropic molecular–zwitterionic tautomers including ylidic forms were computationally investigated at the MP2/6–311++G(2df,pd)//B3LYP/6–311++G(d,p) level of theory in vacuum and in three different surrounding environments: continuum with a low dielectric constant (??=?4) corresponding to a hydrophobic interface of protein–nucleic acid interactions, dimethylsulfoxide (DMSO), and water. The keto-N1HN7H tautomer was established to be the global minimum in vacuum and in continuum with ??=?4, while Hyp molecule exists as a mixture of the keto-N1HN9H and keto-N1HN7H tautomers in approximately equal amounts in DMSO and in water at T?=?298.15?K. We found out that neither intramolecular tautomerization by single proton transfer in the Hyp base, nor intermolecular tautomerization by double proton transfer in the most energetically favorable Hyp·Hyp homodimer (symmetry C _2h ), stabilized by two equivalent N1H…O6 H-bonds, induces the formation of the enol tautomer (marked with an asterisk) of Hyp with cis-oriented O6H hydroxyl group relative to neighboring N1C6 bond. We first discovered a new scenario of the keto–enol tautomerization of Hyp?·?Hyp homodimer (C _2h ) via zwitterionic near-orthogonal transition state (TS), stabilized by N1⁺H…N1^? and O6⁺H…N1^? H-bonds, to heterodimer Hyp^??·?Hyp (C _s ), stabilized by O6H…O6 and N1H…N1 H-bonds. We first showed that Hyp^??·?Thy mispair (C _s ), stabilized by O6H…O4, N3H…N1, and C2H…O2 H-bonds, mimicking Watson–Crick base pairing, converts to the wobble Hyp?·?Thy base pair (C _s ), stabilized by N3H…O6 and N1H…O2 H-bonds, via high- and low-energy TSs and intermediate Hyp?·?Thy^?, stabilized by O4H…O6, N1H…N3, and C2H…O2 H-bonds. The most energetically favorable TS is the zwitterionic pair Hyp⁺?·?Thy^? (C _s ), stabilized by O6⁺H…O4^?, O6⁺H…N3^?, N1⁺H…N3^?, and N1⁺H…O2^? H-bonds. The authors expressed and substantiated the hypothesis, that the keto tautomer of Hyp is a mutagenic compound, while enol tautomer Hyp^? does not possess mutagenic properties. The lifetime of the nonmutagenic tautomer Hyp^? exceeds by many orders the time needed to complete a round of DNA replication in the cell. For the first time purine–purine planar H-bonded mispairs containing Hyp in the anti-orientation with respect to the sugar moiety – Hyp?·?Ade _syn , Hyp?·?Gua^? _syn , and Hyp?·?Gua _syn , that closely resembles the geometry of the Watson–Crick base pairs, have been suggested as the source of transversions. An influence of the surrounding environment (??=?4) on the stability of studied complexes and corresponding TSs was estimated by means of the conductor-like polarizable continuum model. Electron-topological, structural, vibrational, and energetic characterictics of all conventional and nonconventional H-bonds in the investigated structures are presented. Presented data are key to understanding elementary molecular mechanisms of mutagenic action of Hyp as a product of the adenine deamination in DNA.     

Quantum chemical study of the structure and properties of citrinin

Detailed structures and electronic properties of three tautomeric forms of the toxin citrinin were investigated using several quantum calculation methods. Energetic preference of the predominant p- and o-quinone methide tautomeric forms is dependent on the method of calculation. A previously unstudied carboxylic acid enol tautomer was calculated to be surprisingly stable in vacuo, being within 2.5 kcal mol^? 1 at the B3LYP/6-311++G(2d,2p) level of theory. Despite differences in bond nature and connectivity of tautomers, the natural bond orbital analysis revealed that tautomeric forms share similar natural charges and natural electron configurations. Calculated bond lengths corresponded with experimentally observed values and assignments for the calculated infrared vibrational frequencies are reported.     

DPT tautomerization of the long A∙A* Watson-Crick base pair formed by the amino and imino tautomers of adenine: combined QM and QTAIM investigation

Combining quantum-mechanical (QM) calculations with quantum theory of atoms in molecules (QTAIM) and using the methodology of sweeps of the energetic, electron-topological, geometric and polar parameters, which describe the course of the tautomerization along the intrinsic reaction coordinate (IRC), we showed for the first time that the biologically important A?A* base pair (C_s symmetry) formed by the amino and imino tautomers of adenine (A) tautomerizes via asynchronous concerted double proton transfer (DPT) through a transition state (TS), which is the A⁺?A^? zwitterion with the separated charge, with C_s symmetry. The nine key points, which can be considered as electron-topological “fingerprints” of the asynchronous concerted A?A*?A*?A tautomerization process via the DPT, were detected and completely investigated along the IRC of the A?A*?A*?A tautomerization. Based on the sweeps of the H-bond energies, it was found that intermolecular antiparallel N6Н?N6 (7.01 kcal mol^?1) and N1H?N1 (6.88 kcal mol^?1) H-bonds are significantly cooperative and mutually reinforce each other. It was shown for the first time that the A?A*?A*?A tautomerization is assisted by the third C2H?HC2 dihydrogen bond (DHB), which, in contrast to the two others N6H?N6 and N1H?N1 H-bonds, exists within the IRC range from ?2.92 to 2.92 Å. The DHB cooperatively strengthens, reaching its maximum energy 0.42 kcal mol^?1 at IRC?=??0.52 Å and minimum energy 0.25 kcal mol^?1 at IRC?=??2.92 Å, and is accompanied by strengthening of the two other aforementioned classical H-bonds. We established that the C2H?HC2 DHB completely satisfies the electron-topological criteria for H-bonding, in particular Bader’s and all eight “two-molecule” Koch and Popelier’s criteria. The positive value of the Grunenberg’s compliance constant (5.203 Å/mdyn) at the TS_A?A*?A*?A proves that the C2H?HC2 DHB is a stabilizing interaction. NBO analysis predicts transfer of charge from σ(C2–H) bonding orbital to σ*(H–C2) anti-bonding orbital; at this point, the stabilization energy E⁽²⁾ is equal to 0.19 kcal mol^?1 at the TS_A?A*?A*?A.     

Binding of NDGA and morin with phospholipase A2: experimental and computational evidences

The effects of morin and nordihydroguaiaretic acid (NDGA), two plant secondary metabolites, on porcine pancreatic phospholipase A₂ (PLA₂) were investigated by isothermal titration calorimetry (ITC) and in silico docking analyses. The binding energies obtained for NDGA and morin from the ITC studies are ? 6.36 and ? 5.91 kcal mol^? 1, respectively. Similarly, the glide scores obtained for NDGA and morin towards PLA₂ were ? 7.32 and ? 7.23 kcal mol^? 1, respectively. Further the docked complexes were subjected to MD simulation in the presence of explicit water molecules to check the binding stability of the ligands in the active site of PLA₂. The bound ligands make hydrogen bonds with the active site residues of the enzyme and coordinate bonds with catalytically important Ca²⁺ ion. The binding of ligands at the active site of PLA₂ may also contribute to the reported anti-inflammatory properties of NDGA and morin.     

Density functional theory study on the reaction mechanism of synthesizing 1,3-dimethyl-2-imidazolidinone by urea method

We report a first-principles density functional theory investigation on tailoring the fundamental reaction mechanism of synthesizing 1,3-dimethyl-2-imidazolidinone (DMI) through the urea method with water serving as both solvent and catalyst. The nucleophilic cyclization reaction is implemented by two ammonia removal steps. One –NH group of dimethylethylenediamine (DMEDA) first attacks the carbon atom of urea, eliminating one –NH₃ group and forming an intermediate state CH₃NHC₂H₄N(CH₃)CONH₂ (IM^I). IM^I subsequently undergoes the cyclization process through a secondary ammonia removal via similar manner. Without water, the two ammonia removal steps are both slightly exothermic with high activation barriers (~50 kcal mol^-1). As water participated in the reaction, the kinetics of the two steps can be significantly improved, respectively. The role that water plays, beside as solvent, more importantly, is to serve as a proton exchange bridge. Due to the spatial configuration, the direct proton migration from the N atoms of ethylenediamine to urea is difficult to occur. The water bridge facilitates the proton migration by shortening the migration distance. As a consequence, the activation barriers are considerably lowered down to ~30 kcal mol^-1, indicating a strong catalytic effect from water. In contrast, the three possible side reactions of IM^I, even catalyzed by water, have higher activation barriers due to strong steric inhibitive effect and consequently become difficult to occur at the same condition. The current computational understanding on the prototypical reaction to DMI can be extended to guide developing more efficient routes to synthesize imidazolidinone derivatives through the urea method.     

Synthesis, 18F‐Radiolabelling and Biological Evaluations of C‐6 Alkylated Pyrimidine Nucleoside Analogues

Synthesis of pyrimidine derivatives with a side‐chain attached to the C‐6 of pyrimidine ring (6–14) is reported. Target compounds 8 and 12 were subjected to in vitro phosphorylation tests, determination of their binding affinities to herpes simplex virus (HSV‐1) thymidine kinase (TK) and catalytic turnover constants. Fluorinated pyrimidine derivative 12 (40 µM) exhibited better binding affinity for HSV‐1 TK than acyclovir (ACV, 170 µM) and ganciclovir (GCV, 48 µM). Catalytic turnover constant (k _cat) of 12 (0.08 s^? 1) was close to the k _cat values of ACV (0.10 s^? 1) and GCV (0.10 s^? 1). Furthermore, compounds 8 and 12 showed no cytotoxic effects in HSV‐1 TK‐transduced and non‐transduced cell lines. Besides, compounds 8 and 12 did not exhibit antiviral or cytostatic activities against several viruses and malignant tumor cell lines that were evaluated. The new fluorinated pyrimidine derivative 16 that is phosphorylated by HSV‐1 TK could be developed as non‐toxic PET‐tracer molecule. Thus, ¹⁸F labelling of the precursor 14 was performed by nucleophilic substitution using [¹⁸F] tetrabutylammonium fluoride as the fluorinating reagent.     

Structural and functional studies on a variant of cystatin purified from brain of Capra hircus

Cystatins, known for their ubiquitous presence in mammalian system are thiol protease inhibitors serving important physiological functions. Here, we present a variant of cystatin isolated from brain of Capra hircus (goat) which is glycosylated but lacks disulphide bonds. Caprine brain cystatin (CBC) was isolated using alkaline treatment, ammonium sulphate fractionation (40–60%) and gel filtration chromatography on Sephacryl S-100HR column with an overall yield of 26.29% and 322-fold purification. The inhibitor gave a molecular mass of ~44 kDa as determined by SDS-PAGE and gel filtration behaviour. The Stokes radius and diffusion coefficient of CBC were 27.14 Å and 8.18 × 10^?7 cm² s^?1, respectively. Kinetic data revealed that CBC inhibited thiol proteases reversibly and competitively, with the highest inhibition towards papain (K_i = 4.10 nM) followed by ficin and bromelain. CBC possessed 34.7% α-helical content as observed by CD spectroscopy. UV, fluorescence, CD and FTIR spectroscopy revealed significant conformational change upon CBC-papain complex formation. Isothermal titration calorimetry (ITC) was used to measure the thermodynamic parameters – ΔH, ΔS, ΔG along with N (binding stoichiometry) for CBC-papain complex formation. Binding stoichiometry (N = .97 ± .07 sites) for the CBC-papain complex indicates that cystatin is surrounded by nearly one papain molecule. Negative ΔH (?5.78 kcal mol^?1) and positive ΔS (11.01 cal mol^?1 deg^?1) values suggest that the interaction between CBC and papain is enthalpically as well as entropically favoured process. The overall negative ΔG (?9.19 kcal mol^?1) value implies a spontaneous CBC-papain interaction.     

Molecular aspects on adriamycin interaction with hmga1 regulatory region and its inhibitory effect on HMGA1 expression in human cervical cancer

High mobility group A1 (HMGA1), a non-histone chromosomal protein, is highly expressed in a wide range of human cancers including cervical, breast, and prostate cancers. Therefore, hmga1 gene is considered as an attractive potential target for anticancer drugs. We have chosen 27 bp DNA sequence from a regulatory region of hmga1 promoter and studied its interaction with adriamycin (ADM) and in vitro expression of HMGA1 in the presence of ADM in HeLa cell line. A variety of biophysical techniques were employed to understand the characteristics of [DNA–ADM] complex. Spectrophotometric titration data, DNA denaturation profiles, and quenching of fluorescence of ADM in the presence of DNA demonstrated a strong complexation between DNA and ADM with a high binding affinity (Ka) of 1.3 × 10⁶ M^?1 and a stoichiometry of 1:3 (drug:nucleotide). The energetics of binding obtained from isothermal titration calorimetry and differential scanning calorimetry suggest the binding to be exothermic and enthalpy (?H, ?6.7 ± 2.4 kcal M^?1) and entropy (TΔS, 18.5 ± 6.4 kcal M^?1) driven (20°C), which is typical of intercalative mode of binding. Further, results on decreased expression (by ~70%) of HMGA1 both at mRNA and protein levels in association with the observed cell death (by ~75%) in HeLa cell line, clearly confirm that ADM does target hmga1; however, the effect of ADM on genes other than hmga1 either directly or via hmga1-mediated pathways cannot be ruled out in the observed cytotoxicity. Therefore, hmga1 in general and particularly the regulatory region is a promising target for therapeutic strategy in combating cancer.     

Effect of dielectric constant of medium on protonation equilibria of ethylenediamine

Abstract

The effect of dielectric constant of medium on protonation equilibria has been studied by determining protonation constants of ethylenediamine pH metrically in various concentrations (0–60%v/v) of acetoni-trile– and ethylene glycol–water mixtures, at an ionic strength of 0.16mol L^?1 and at 303.0 K. MINIQUAD75 computer program has been used for the calculation of protonation constants. Linear and non-linear variations of step-wise protonation constants with reciprocal of dielectric constant of the solvent mixtures have been attributed to the dominance of the electrostatic and non-electrostatic forces, respectively. The trend is explained on the basis of solute–solute and solute–solvent interactions, solvation, proton transfer processes and dielectric constants of the media.     

DFT investigation on the decarboxylation mechanism of <Emphasis Type="Italic">ortho</Emphasis> hydroxy benzoic acids with acid catalysis

A density functional theory (DFT) study was performed to explore the mechanisms of the acid-catalyzed decarboxylation reaction of salicylic acids using the B3LYP method with 6-31++G(d,p) basis set in both gas phase and aqueous environment. The α-protonated cation of carboxylate acid was formed during the decarboxylation process in acidic conditions, and the presence of hydrogen ions promotes decarboxylation greatly by significantly decreasing the overall reaction energy barriers to 20.98 kcal mol^?1 in gas phase and 20.93 kcal mol^?1 in water, respectively. The hydrogen in the α-carbon came directly from the acid rather than from the carboxyl group in neutral state. Compared with the reaction in gas phase, water in aqueous state causes the reaction to occur more easily. Substituents of methyl group, chlorine and fluorine at the ortho-position to the carboxyl of salicylic acid could further lower the decarboxylation energy barriers and facilitate the reaction.     

Combination of Azidothymidine (AZT) and (E)-5-(2-Bromovinyl)-2′-deoxyuridine (BVDU) Inhibits the Replication of Herpes Simplex Virus Type 1 (HSV-1) and Type 2 (HSV-2) and Varicella Zoster Virus (VZV) Strains That Are Deficient in the Expression of the Viral Thymidine Kinase (tk)

Abstract

Combinations of high concentrations of AZT with BVDU, acyclovir (ACV) or ganciclovir (GCV) show antagonism against TK⁺ HSV-1, but not TK⁺ VZV strains, in cell cultures. When BVDU and AZT were used in combination against TK^? HSV-1, TK^? HSV-2 and TK^? VZV strains, a pronounced inhibition of viral replication was observed. This potentiating effect was not seen if AZT was combined with ACV or GCV.     

Binding of anti-cardiovascular drug to serum albumin: an insight in the light of spectroscopic and computational approaches

Isoprenaline hydrochloride is a potential cardiovascular drug helps in the smooth functioning of the heart muscles. So, we have performed the binding study of ISO with BSA. This study was investigated by UV absorption, fluorescence, synchronous fluorescence, circular dichroism, etc. The analysis of intrinsic fluorescence data showed the low binding affinity of ISO. The binding constant K_b was 2.8 × 103 M-1 and binding stoichiometry (n) was approximately one and the Gibb’s free energy change at 310 K was determined to be -8.69 kcal mol^?1. Negative Gibb’s free energy change shows the spontaneity of the BSA and ISO interaction. We have found ISO-induced alternation in the UV absorption, synchronous fluorescence and CD spectra in the absence and presence of the quencher indicates the complex formation. In synchronous fluorescence, red shift was obtained because of the complex formation of BSA and ISO. The distance (r) between the BSA (donor) and ISO (acceptor) was 2.89 nm, determined by FRET. DLS measurements interpreted complex formation due to the reduction in hydrodynamic radii of the protein in the presence of the drug. The binding site of ISO was found to be nearer to Trp 134 with the help of molecular docking and the ΔG° was found to be –10.2 kcal mol^?1. The esterase activity result suggests that ISO acts as competitive inhibitor. Thus, this study would help to determine the binding capacity of the drug to the protein which may indicate the efficiency of diffusion of ISO into the blood for the treatment of heart diseases.     

Quantum-chemical and empirical calculations on phospholipids. III. Hydration of the dimethylphosphate anion

The binding of water to the dimethylphosphate anion (DMP^?) was calculated using the PCILO method. We found binding energies of 25.95 kcal·mol^?1 in the O1-P-O3 plane bridging the anionic oxygen atoms and 19.3 kcal·mol^?1 for the one-site association of a water molecule to an anionic oxygen atom of DMP^?. In this range one water molecule added to DMP^? in the O1 … O3 bridged configuration has a significantly higher binding energy to DMP^? than water molecules added to other binding sites. The total binding energy of 5 water molecules to DMP^? is 92 kcal·mol^?1, a quantity which is about 10% less than the sum of the binding energies of the corresponding monohydrates.     

Local dynamics and structure of the solvated hydroxide ion in water

Details of the structural diffusion mechanism in proton transfer (PT) reactions involving a hydroxide ion in water and the structure of the solvated ions and transfer intermediates are in dispute. Here, we elucidate the mechanism of PT involving a hydroxide ion in water by molecular dynamics simulations using a dissociating water model based on ab initio calculations. We find that the hydroxide ion in bulk water is present as the four-coordinated OH^? (H₂O)₄ complex, which loses a water molecule before a PT occurs through the formation of proton sharing intermediates in general agreement with the previously disputed first principles studies of small systems.