Study of the hydrogen bond in different orientations of adenine-thymine base pairs: An <Emphasis Type="Italic">ab initio</Emphasis> study

In order to gain deeper insight into structure, charge distribution, and energies of A-T base pairs, we have performed quantum chemical ab initio and density functional calculations at the HF (Hartree-Fock) and B3LYP levels with 3-21G*, 6-31G*, 6-31G**, and 6-31++G** basis sets. The calculated donor-acceptor atom distances in the Watson-Crick A-T base pair are in good agreement with the experimental mean values obtained from an analysis of 21 high resolution DNA structures. In addition, for further correction of interaction energies between adenine and thymine, the basis set superposition error (BSSE) associated with the hydrogen bond energy has been computed via the counterpoise method using the individual bases as fragments. In the Watson-Crick A-T base pair there is a good agreement between theory and experimental results. The distances for (N2...H23-N19), (N8-H13...O24), and (C1...O18) are 2.84, 2.94, and 3.63 A, respectively, at B3LYP/6-31G** level, which is in good agreement with experimental results (2.82, 2.98, and 3.52 A). Interaction energy of the Watson-Crick A-T base pair is -13.90 and -10.24 kcal/mol at B3LYP/6-31G** and HF/6-31G** levels, respectively. The interaction energy of model (9) A-T base pair is larger than others, -18.28 and -17.26 kcal/mol, and for model (2) is the smallest value, -13.53 and -13.03 kcal/mol, at B3LYP/6-31G** and B3LYP/6-31++G** levels, respectively. The computed B3LYP/6-31G** bond enthalpies for Watson-Crick A-T pairs of -14.4 kcal/mol agree well with the experimental results of -12.1 kcal/mol deviating by as little as -2.3 kcal/mol. The BSSE of some cases is large (9.85 kcal/mol) and some is quite small (0.6 kcal/mol).     

Reaction pathway and free energy profiles for butyrylcholinesterase-catalyzed hydrolysis of acetylthiocholine

The catalytic mechanism for butyrylcholineserase (BChE)-catalyzed hydrolysis of acetylthiocholine (ATCh) has been studied by performing pseudobond first-principles quantum mechanical/molecular mechanical-free energy (QM/MM-FE) calculations on both acylation and deacylation of BChE. Additional quantum mechanical (QM) calculations have been carried out, along with the QM/MM-FE calculations, to understand the known substrate activation effect on the enzymatic hydrolysis of ATCh. It has been shown that the acylation of BChE with ATCh consists of two reaction steps including the nucleophilic attack on the carbonyl carbon of ATCh and the dissociation of thiocholine ester. The deacylation stage includes nucleophilic attack of a water molecule on the carboxyl carbon of substrate and dissociation between the carboxyl carbon of substrate and hydroxyl oxygen of Ser198 side chain. QM/MM-FE calculation results reveal that the acylation of BChE is rate-determining. It has also been demonstrated that an additional substrate molecule binding to the peripheral anionic site (PAS) of BChE is responsible for the substrate activation effect. In the presence of this additional substrate molecule at PAS, the calculated free energy barrier for the acylation stage (rate-determining step) is decreased by ~1.7 kcal/mol. All of our computational predictions are consistent with available experimental kinetic data. The overall free energy barriers calculated for BChE-catalyzed hydrolysis of ATCh at regular hydrolysis phase and substrate activation phase are ~13.6 and ~11.9 kcal/mol, respectively, which are in reasonable agreement with the corresponding experimentally derived activation free energies of 14.0 kcal/mol (for regular hydrolysis phase) and 13.5 kcal/mol (for substrate activation phase).     

Theoretical potential functions and vibrational analysis for halocarbonyl azides CXO-NNN (X=F,Cl and Br)

The structural stability of halocarbonyl azides CXO-NNN (X=F, Cl and Br) was investigated by DFT and MP2 calculations using the 6-311++G** basis set. From the calculations, the molecules were found to have an s-cis<--> s-trans conformational equilibrium with cis being the lower -energy form. Full energy optimizations were carried out for the transition states and the minima at the B3LYP/6 -311++G** and MP2/6 -311++G** levels, from which the rotational barriers were calculated to be of the order 8-10 kcal x mol(-1). The vibrational frequencies were computed at the DFT -B3LYP level and the vibrational assignments for the normal modes of the stable conformers were made on the basis of normal coordinate calculations.     

The Gibbs-Donnan near-equilibrium system of heart

The gradients of the major inorganic ions across the plasma membrane of heart were examined to determine the factors controlling the extent and direction of the changes induced during injury, certain diseases, and electrolyte disturbances. The ionic environment was altered by changing only the concentration of inorganic phosphate, [sigma Pi]o, from 0 to 1.2 to 5 mM in the Krebs-Henseleit buffer perfusing working rat hearts. Raising [sigma Pi]o from 1.2 to 5 mM resulted in a decrease in total Mg2+ content and calculated free cytosolic [Mg2+] from 0.44 to 0.04 mM, conversion of 4 mmol of MgATP2- to ATP4- and a decrease in measured intracellular [Cl-]i from 41 to 16 mM. At all levels of [sigma Pi]o, both the [Na+]i and [K+]i were invariant at about 3 mM and 130 mM, respectively, as was the energy of hydrolysis of the terminal phosphate bond of sigma ATP, delta GATP Hydr, of -13.2 kcal/mol. The relationship maintained between the ions on both sides of the plasma membrane by the 3Na+/2K(+)transporting ATPase (EC 3.6.1.37) and an open K+ channel was: (formula; see text) The energy of the gradients of the other inorganic ions across the plasma membrane, delta G[ion]o/i, exhibited three distinct quanta of energy derived from the prime quantum of delta GATP Hydr of -13.2 kcal/mol. The second quantum was about one-third of delta GATP Hydr or +/- 4.4 kcal/mol and comprised the delta G[Na+]o/i, delta G[Mg2+]o/i, and delta G[HPO42-]o/i. These results indicated near-equilibrium was achieved by the reactants of the 3Na+/2K(+)-ATPase, the K+ channel, the Na(+)-Pi co-transporter, and a postulated net Mg2+/H2PO4- exchanger. The third quantum was one-third of delta G[Na+]o/i or about +/- 1.5 kcal/mol and comprised delta G[H+]o/i, delta G[HCO3-]o/i, and delta G[Cl-]o/i. The delta G[K+]o/i was 0, indicating near-equilibrium between the chemical energy of [K+]o/i and the E across the plasma membrane of -83 mV. It is concluded that the gradients of the major inorganic ions across the plasma membrane and the potential across that membrane constitute a Gibbs-Donnan equilibrium system catalyzed by transport enzymes sharing common substrates. The chemical and electrical energies of those gradients are equal in magnitude and opposite in sign to the chemical energy of ATP hydrolysis.     

Theoretical investigation of (-)1-(benzofuran-2-yl)-2-propylaminopentane[(-)-BPAP] as a hydroxyl radical scavenger

The chemical reactions between (-)-BPAP and .OH were studied using molecular orbital theory, with several simplified models. (-)-BPAP was proved to be a good radical scavenger. It was found that every atom of the benzofuran ring, except carbon 3, was capable of easily trapping the radical, where the most active site was carbon 1 on the furan part. The activation energies for the trappings were ca. 10kcal/mol by the calculations at UHF/6-31G(d)//UHF/3-21G level. Since the single radical trapping products were still radicals, these could trap further radicals by way of cascade without any activation energy. Thus, the double radical trapping products were very stable, of which the lowest energy level was about 80kcal/mol lower than the starting reactants.     

Heat of PPi hydrolysis varies depending on the enzyme used. Yeast and corn vacuolar pyrophosphatase

With yeast-soluble inorganic pyrophosphatase, the heat released during PP(i) hydrolysis was -6.3 kcal/mol regardless of the KCl concentration in the medium. With the membrane-bound pyrophosphatase of corn vacuoles, the heat released varies between -23.5 and -7.5 kcal/mol depending on the KCl concentration in the medium and whether or not a H(+) gradient is formed across the vacuole membranes. The data support the proposal that enzymes are able to handle the energy derived from phosphate compound hydrolysis in such a way as to determine the parcel that is used for work and the fraction that is converted into heat.     

The thermodynamics of bovine and porcine insulin and proinsulin association determined by concentration difference spectroscopy

Difference spectroscopy was used to determine the equilibrium constants and thermodynamic parameters for the monomer-dimer association of bovine and porcine insulin and bovine proinsulin at pH 2.0 and 7.0. At pH 2 delta G degree 25, delta S degree, and delta H degree for dimerization of bovine insulin were found to be -6.6 kcal/mol, -18 cal/mol-deg, and -12 kcal/mol, respectively. Porcine insulin behaved similarly to bovine insulin in its dimerization properties in that delta G degree 25, delta S degree, and delta H degree were found to be -6.8 kcal/mol, -14 cal/mol-deg, and -11 kcal/mol, respectively. At pH 7 delta G degree 25, delta S degree, and delta H degree for dimerization of bovine insulin were found to be -7.2 kcal/mol, -16 cal/mol/deg, and -12 kcal/mol, respectively. At pH 7.0 delta G degree 25, delta S degree, and delta H degree for dimerization of porcine insulin were -6.7 kcal/mol, -11.6 cal/mol-deg, and -10 kcal/mol, respectively. The similarity in the thermodynamic parameters of both insulin species at the different pH's suggests that there are minimal structural changes at the monomer-monomer contact site over this pH range. The dimerization of both insulin species is under enthalpic control. This may suggest that the formation of the insulin dimer is not driven by hydrophobic bonding but, rather, is driven by the formation between subunits of four hydrogen bonds in an apolar environment. At pH 2 delta G degree 25, delta S degree, and delta H degree for dimerization of bovine proinsulin were found to be -5.3 kcal/mol, -26 cal/mol-deg, and -13 kcal/mol, respectively. At pH 7 delta G degree 25, delta S degree, and delta H degree for dimerization of proinsulin were -5.9 kcal/mol, -4.2 cal/mol-deg, and -7.2 kcal/mol, respectively. Although the presence of the C-peptide on proinsulin does not drastically affect the overall free energy change of dimer formation (as compared to insulin), the other thermodynamic parameters are rather drastically altered. This may be because of electrostatic interactions of groups on the C-peptide with groups on the B-chain which are near the subunit contact site in the insulin dimer.     

Theoretical insights into the protonation states of active site cysteine and citrullination mechanism of Porphyromonas gingivalis peptidylarginine deiminase

Porphyromonas gingivalis peptidylarginine deiminase (PPAD) catalyzes the citrullination of peptidylarginine, which plays a critical role in the rheumatoid arthritis (RA) and gene regulation. For a better understanding of citrullination mechanism of PPAD, it is required to establish the protonation states of active site cysteine, which is still a controversial issue for the members of guanidino‐group‐modifying enzyme superfamily. In this work, we first explored the transformation between the two states: State N (both C351 and H236 are neutral) and State I (both residues exist as a thiolate–imidazolium ion pair), and then investigated the citrullination reaction of peptidylarginine, using a combined QM/MM approach. State N is calculated to be more stable than State I by 8.46 kcal/mol, and State N can transform to State I via two steps of substrate‐assisted proton transfer. Citrullination of the peptidylarginine contains deamination and hydrolysis. Starting from State N, the deamination reaction corresponds to an energy barrier of 18.82 kcal/mol. The deprotonated C351 initiates the nucleophilic attack to the substrate, which is the key step for deamination reaction. The hydrolysis reaction contains two chemical steps. Both the deprotonated D238 and H236 can act as the bases to activate the hydrolytic water, which correspond to similar energy barriers (～17 kcal/mol). On the basis of our calculations, C351, D238, and H236 constitute a catalytic triad, and their protonation states are critical for both the deamination and hydrolysis processes. In view of the sequence similarity, these findings may be shared with human PAD1–PAD4 and other guanidino‐group‐modifying enzymes. Proteins 2017; 85:1518–1528. © 2017 Wiley Periodicals, Inc.     

Hybrid density functional theory with a specific reaction parameter: hydrogen abstraction reaction of difluoromethane by the hydroxyl radical

Accurate potential energy surfaces for the OH + CH2F2 --> H2O + CHF2 reaction are constructed using hybrid and hybrid meta density functional theory methods (mPW1PW91, B1B95, and mPW1B95) with specific reaction parameters in conjunction with the 6-31 + G(d,p) basis set. The accuracy of a surface is examined by comparing the calculated rate constants with the experimental ones. The rate constants are calculated over the temperature range 200-1,500 K using variational transition state theory with multidimensional tunneling contributions. The hybrid density functional theory methods with specific-reaction-parameter Hartree-Fock exchange contributions (39.2-41.0% for mPW1PW91, 41.0-42.2% for B1B95, and 44.9-46.3% for mPW1B95, respectively) provide accurate rate constants over an extended temperature range. The classical barrier height for the hydrogen abstraction reaction on these potential energy surfaces is determined to be 5.0-5.3 kcal mol(-1), and the best estimate value is 5.14 kcal mol(-1).     

Ca(2+) release and heat production by the endoplasmic reticulum Ca(2+)-ATPase of blood platelets. Effect of the platelet activating factor.

Different sarco/endoplasmic reticulum Ca(2+)-ATPases isoforms are found in blood platelets and in skeletal muscle. The amount of heat produced during ATP hydrolysis by vesicles derived from the endoplasmic reticulum of blood platelets was the same in the absence and presence of a transmembrane Ca(2+) gradient. Addition of platelets activating factor (PAF) to the medium promoted both a Ca(2+) efflux that was arrested by thapsigargin and an increase of the yield of heat produced during ATP hydrolysis. The calorimetric enthalpy of ATP hydrolysis (DeltaH(cal)) measured during Ca(2+) transport varied between -10 and -12 kcal/mol without PAF and between -20 and -24 kcal/mol with 4 microM PAF. Different from platelets, in skeletal muscle vesicles a thapsigargin-sensitive Ca(2+) efflux and a high heat production during ATP hydrolysis were measured without PAF and the DeltaH(cal) varied between -10 and -12 kcal/mol in the absence of Ca(2+) and between -22 up to -32 kcal/mol after formation of a transmembrane Ca(2+) gradient. PAF did not enhance the rate of thapsigargin-sensitive Ca(2+) efflux nor increase the yield of heat produced during ATP hydrolysis. These findings indicate that the platelets of Ca(2+)-ATPase isoforms are only able to convert osmotic energy into heat in the presence of PAF.     

Effect of carboxylate and Na ions on the tautomeric status of 9-methylguanine

Interactions of 9-methylguanine (m9Gua) with carboxylate ion of acetic acid (CH3COO-) and Na+ were studied by 1H NMR spectroscopy and ab initio quantum chemical calculations of the B3LYP/6-31++G(d,p) and B3LYP/6-311++G(d,p) levels of theory. Changes in the m9Gua 1H NMR spectrum in the presence of the equimolar amount of sodium acetate (NaAc), which in anhydrous DMSO dissociates to CH3COO- and Na+, were interpreted as a consequence of a complex formation of m9Gua in the amino-keto-N1H tautomeric form (m9GuaN1H) with carboxylate ions via two H-bonds involving amino and N1H-imino protons. Quantum chemical calculations of interactions of the m9GuaN1H ground-state tautomer and the m9GuaN3H high energy one with relative energy 20.01 kcal/mol show that the ground state tautomer forms the ground-state complex CH3COO-:mgGuaN1H, by 5.57 kcal/mol more stable than the CH3COO-:m9GuaN3H complex, and coordination of Na+ with the O6 and N7 atoms reduces this energy difference to 2.57 kcal/mol. Such a coordination of Na+ with tautomer m9GuaN3H therewith decreases its relative energy only to 13.31 kcal/mol. Non-additivity of the two ligands contributions to the 8-times reduction of the relative energy of the high energy tautomer in the CH3COO-:m9GuaN3H:Na+(O6,N7) triple complex was concluded, the role of CH3COO- being dominant. Besides, coordination with Na+ resulted in an iminoproton transfer from the base to CH3COO- in the triple complexes of both tautomers, according to calculations in vacuum. Biological significance of the results is noticed.     

The investigation of hydrocarbon cracking reaction energetics with composite energy methods

Hydrocarbo cracking reactions are one of the most commonly encountered reactions in the petroleum industry, and the energetics of the reactions are crucial in understanding the reaction mechanisms and predicting reaction rates. In this work, a modified composite energy method (CBS-RAD(MP2)) is created as a version of the CBS-RAD method which gives accurate energetics for hydrocarbon free radical reactions. It replaces the time consuming QCISD(fc)/6-31g* method in the geometry optimization and frequency calculation steps with MP2(full)/6-31g* level calculations. The accuracy of the new CBS-RAD(MP2) method is compared with the widely used G2, G3 and CBS-QB3 composite methods for predicting heats of reaction and activation barriers of 14 hydrocarbon cracking reactions. We find that the new CBS-RAD(MP2) method has the second least RMS error of 1.22 kcal/mol for heats of reaction calculations. For activation energy calculations, the new CBS-RAD(MP2) method has the least RMS error of 1.37 kcal/mol. Moreover, the CBS-RAD(MP2) method was found to require only 81% of the computational time required compared to the CBS-QB3 method, 32% of G3 and 15% of the G2 method, making it an attractive alternative for predicting hydrocarbon cracking reaction energetics.     

Evaluation of carbohydrate molecular mechanical force fields by quantum mechanical calculations

A quantitative evaluation of 20 second-generation carbohydrate force fields was carried out using ab initio and density functional methods. Geometry-optimized structures (B3LYP/6-31G(d)) and relative energies using augmented correlation consistent basis sets were calculated in gas phase for monosaccharide carbohydrate benchmark systems. Selected results are: (i). The interaction energy of the alpha-d-glucopyranose.H(2)O heterodimer is estimated to be 4.9 kcal/mol, using a composite method including terms at highly correlated (CCSD(T)) level. Most molecular mechanics force fields are in error in this respect; (ii). The (3)E envelope (south) pseudorotational conformer of methyl 5-deoxy-beta-d-xylofuranoside is 0.66 kcal/mol more stable than the (3)E envelope (north) conformer and the alpha-anomer of methyl d-glucopyranoside is 0.82 kcal/mol more stable than the beta-anomer; (iii). The relative energies of the (gg, gt and tg) rotamers of methyl alpha-d-glucopyranoside and methyl alpha-d-galactopyranoside are (0.13, 0.00, 0.15) and (0.64, 0.00, 0.77) kcal/mol, respectively. The results of the quantum mechanical calculations are compared with the results of calculations using the 20 second-generation carbohydrate force fields. No single force field is consistently better than the others for all the test cases. A statistical assessment of the performance of the force fields indicates that CHEAT(95), CFF, certain versions of Amber and of MM3 have the best overall performance, for these gas phase monosaccharide systems.     

Equilibrium constants under physiological conditions for the reactions of choline kinase and the hydrolysis of phosphorylcholine to choline and inorganic phosphate.

The observed equilibrium constants (Kobs) of the P-choline hydrolysis reaction have been determined under physiological conditions of temperature (38 degrees) and ionic strength (0.25 M) and physiological ranges of pH and free [Mg2+]. Using sigma and square brackets to indicate total concentrations: (see article.) The value of Kobs has been found to be relatively insensitive to variations in pH and free [Mg2+]. At pH 7.0 and taking the standard state of liquid water to have unit activity ([H2O] = 1), Kobs = 26.6 M at free [Mg2+] = 0 [epsilon G0obs = -2.03 kcal/mol(-8.48 kJ/mol)], 26.8 M at free [Mg2+] = 10(-3) M, and 28.4 M at free [Mg2+] = 10(-2) M. At pH 8.0, Kobs = 18.8 M at free [Mg2+] = 0, 19.2 M at free [Mg2+] = 10(-3), and 22.2 M at free [Mg2+] = 10(-2) M. These values apply only to situations where choline and Pi concentrations are both relatively low (such as the conditions found in most tissues). At higher concentrations of phosphate and choline, the value of Kobs becomes significantly increased since HPO42- complexes choline weakly (association constant = 3.3 M-1). The value of K at 38 degrees and I = 0.25 M is calculated to be 16.4 +/- 0.3 M [epsilonG0 = 1.73 kcal/mol (-7.23 kJ/mol)]. The K for the P-choline hydrolysis reaction has been combined with the K for the ATP hydrolysis reaction determined previously under physiological conditions to calculate a value of 4.95 X 10(-3 M [deltaG0 j.28 kcal/mol (13.7 kJ/mol] for the K of the choline kinase reaction (EC 2.7.1.32), an important step in phospholipid metabolism: (see article.) Likewise, values for Kobs for the choline kinase reaction at 38 degrees, pH 7.0, and I = 0.25 M have been calculated to be 5.76 X 10(4) [deltaG0OBS = -6.77 KCAL/MOL (-28.3 KJ/mol)] at [Mg2+] = 0; 1.24 X 10(4) [deltaG0obs = -5.82 kcal/mol (-24.4 kJ/mol)] at [Mg2+] = 10(-3) M and 8.05 X 10(3) [delta G0obs = -5.56 kcal/mol (-23.3 kJ/mol)] at [Mg2+ = 10(-2) M. Attempts to determine the Kobs of the choline kinase reaction directly were unsuccessful because of the high value of the constant. The results indicate that in contrast to the high deltaG0obs for the hydrolysis of the ester bond of acetylcholine, the deltaG0obs for the hydrolysis of the ester bond of P-choline is quite low, among the lowest known for phosphate ester bonds of biological interest.     

Theoretical Investigation of the Molecular Structure of the πκ DNA Base Pair

Abstract

The structure of the nonclassical πκ base pair (7–methyl-oxoformycin … 2,4-diaminopyrimidine) was studied at the ab initio Hartree-Fock (HF) and MP2 levels using the 6–31G* and 6–31G** basis sets. The πκ base pair is bound by three parallel hydrogen bonds with the donor-acceptor-donor recognition pattern. Recently, these bases were proposed as an extension of the genetic alphabet from four to six letters (Piccirilli et al. Nature 343, 33(1990)). By the HF/6- 31G* method with full geometry optimization we calculated the 12 degree propeller twist for the minimum energy structure of this complex. The linearity of hydrogen bonds is preserved in the twisted structure by virtue of the pyramidal arrangement of the κ-base amino groups. The rings of both the π and κ molecules remain nearly planar. This nonplanar structure of the πκ base pair is only 0.1 kcal/mol more stable than the planar (Cs) conformation. The HF/6- 31G* level gas-phase interaction energy of πκ (—13.5 kcal/mol) calculated by us turned out to be nearly the same as the interaction energy obtained previously for the adenine-thymine base pair (—13.4 kcal/mol) at the same computational level. The inclusion of p-polarization functions on hydrogens, electron correlation effects (MP2/6–31G** level), and the correction for the basis set superposition error (BSSE) increase this energy to -14.0 kcal/mol.     

Does Al4H 14 — cluster anion exist? High-level ab initio study

A comprehensive ab initio investigation using coupled cluster theory with the aug-cc-pVnZ, n = D,T basis sets is carried out to identify distinct structures of the Al(4)H(14)(-) cluster anion and to evaluate its fragmentation stability. Both thermodynamic and mechanistic aspects of the fragmentation reactions are studied. The observation of this so far the most hydrogenated aluminum tetramer was reported in the recent mass spectrometry study of Li et al. (2010) J Chem Phys 132:241103-241104. The four Al(4)H(14)(-) anion structures found are chain-like with the multiple-coordinate Al center and can be viewed approximately as comprising Al(2)H(7)(-) and Al(2)H(7) moieties. Locating computationally some of the Al(4)H(14)(-) minima on the correlated ab initio potential energy surfaces required the triple-zeta quality basis set to describe adequately the Al multi-coordinate bonding. For the two most stable Al(4)H(14)(-) isomers, the mechanism of their low-barrier interconversion is described. The dissociation of Al(4)H(14)(-) into the Al(2)H(7)(-) and Al(2)H(7) units is predicted to require 20-22 (10-13) kcal mol(-1) in terms of ΔH (ΔG) estimated at T = 298.15 K and p = 1 atm. However, Al(4)H(14)(-) is found to be a metastable species in the gas phase: the H(2) loss from the radical moiety of its most favorable isomer is exothermic by 18 kcal mol(-1) in terms of ΔH (298.15 K) and by 25 kcal mol(-1) in terms of ΔG(298.15 K), with the enthalpic/free energy barrier involved being less than 1 kcal mol(-1). By contrast with alane Al(4)H(14)(-), only a weakly bound complex between Ga(4)H(12)(-) and H(2) has been identified for the gallium analogue using the relativistic effective core potential.     

DFT study of alpha- and beta-D-allopyranose at the B3LYP/6-311++G ** level of theory

One hundred and two conformations of alpha- and beta-D-allopyranose, the C-3 substituted epimer of glucopyranose, were geometry optimized using the density functional, B3LYP, and the basis set, 6-311++G **. Full geometry optimization was performed on different ring geometries and on the hydroxymethyl rotamers (gg/gt/tg). Analytically derived Hessians were used to calculate zero point energy, enthalpy, and entropy. The lowest energy and free energy conformation found is the alpha-tg(g-)-4C1-c conformation, which is only slightly higher in electronic (approximately 0.2 kcal/mol) and free energy than the lowest energy alpha-D-glucopyranose. The in vacuo calculations showed a small (approximately 0.3 kcal/mol) energetic preference for the alpha- over the beta-anomer for allopyranose in the 4C1 conformation, whereas in the 1C4 conformation a considerable (approximately 1.6 kcal/mol) energetic preference for the beta- over the alpha-anomer for allopyranose was encountered. The results are compared to previous aldohexose calculations in vacuo. Boat and skew forms were found that remained stable upon gradient optimization although many starting boat conformations moved to other skew forms upon optimization. As found for glucose, mannose, and galactose the orientation and interaction of the hydroxyl groups make the most significant contributions to the conformation/energy relationship in vacuo. A comparison of different basis sets and density functionals is made in the Discussion section, confirming the appropriateness of the level of theory used here.     

The enthalpy change for the reduction of nicotinamide–adenine dinucleotide

The heat of the reaction NAD(+)+propan-2-ol=NADH+acetone+H(+) was determined to be 42.5+/-0.6kJ/mol (10.17+/-0.15kcal/mol) from equilibrium measurements at 9-42 degrees C catalysed by yeast alcohol dehydrogenase. With the aid of thermochemical data for acetone and propan-2-ol the values of DeltaH=-29.2kJ/mol (-6.99kcal/mol) and DeltaG(0)=22.1kJ/mol (5.28kcal/mol) are derived for the reduction of NAD (NAD(+)+H(2)=NADH+H(+)). These values are consistent with analogous but less accurate data for the ethanol-acetaldehyde reaction. Thermodynamic data for the reduction of NAD and NADP are summarized.     

Pyrophosphate of high and low energy. Contributions of pH, Ca2+, Mg2+, and water to free energy of hydrolysis

The equilibrium between inorganic pyrophosphate and inorganic orthophosphate was determined at pH values varying between 6.0 and 8.0, in the presence of different concentrations of MgCl2, mixtures of MgCl2 and CaCl2, and different organic solvents. The reactions were catalyzed by yeast inorganic pyrophosphatase. It was found that at 35 degrees C, depending on the conditions used, the observed equilibrium constant of pyrophosphate hydrolysis vary from a value higher than 4 X 10(3) M (delta Goobs more negative than -5.1 kcal/mol) to a value as low as 3 M (delta Goobs -0.7 kcal/mol). The experimental data were used to compute the equilibrium constants of the reactions involving different ionic species. The data presented are interpreted according to the concept that the Keq of hydrolysis of a high energy compound depends on the difference in solvation energy of reactants and products.     

Hydrolysis of nerve agents by model nucleophiles: a computational study

Density functional theory calculations were employed to study the reaction of five nerve agents with model nucleophiles, including EtX(-) and EtXH (X=O, S, Se) for serine, cysteine and selenocysteine, respectively. Calculations at the B3LYP/6-311++G(2d,p) level of theory predict an exothermic reaction between ethoxide and all of the nerve agents studied. As compared to EtO(-) as a nucleophile, these reactions become approximately 30kcal/mol more endothermic for EtS(-), and by approximately 40kcal/mol for EtSe(-). The equivalent reactions with the neutral nucleophiles (EtXH) were more endothermic. The effect of solvation on the reaction thermochemistry was determined using a polarizable continuum model simulating the dielectric constant of chloroform. While there was a large exothermic shift for reactions involving charged nucleophiles with solvation modeling, the corresponding shift was minimal for the reaction with neutral nucleophiles.