Monte Carlo simulations of the solution structure of simple alcohols in water-acetonitrile mixtures

Monte Carlo simulations have been performed to explore the solution structure of ethyl, isopropyl, isobutyl, and tertiary butyl alcohols in pure water, pure acetonitrile, and different mixtures of the two solvents. The explicit solvent studies in NpT ensembles at T = 298 K illustrate that the solute "discriminates" the solvent's components and that the composition of the first solvation shell differs from that of the bulk solution. Since the polarizable continuum dielectric method (PCM) does not presently model the solvation of molecules with both polar and apolar sites in mixed protic solvents, we suggest a direction for further program development wherein a continuum dielectric method would accept more than one solvent and the solute sites would be solvated by user-defined solvent components. The prevailing solvation model will be determined upon the lowest free energy calculated for a particular solvation pattern of the solute having a specific conformational/tautomeric state. Characterization of equilibrium hydrogen-bond formation becomes a complicated problem that depends on the chemical properties of the solute and its conformation, as well as upon the varying nature of the first solvation shell. For example, while the number of hydrogen bonds to secondary and tertiary alcohol solutes are nearly constant in pure water and in water-acetonitrile mixtures with at least 50% water content, the number of hydrogen bonds to primary alcohols gradually decreases for most of their conformations when acetonitrile content is increased. Nonetheless, the calculations indicate that O-H...O(water) hydrogen bonds are still possible in a small fraction of the arrangements for the solution models with water content of 30% or less. The isopentene solute does not form any observable hydrogen bonds, despite having an electron-rich, double-bond site.     

Correlation effects on energy curves for proton transfer. The cation [H5O2]+

Potential curves for proton transfer in [H₅O₂]⁺ and for the dissociation of one OH bond in [H₃O]⁺ were calculated by both ab initio and semi-empirical LCAO MO SCF CI methods. The energy barrier of the symmetric double minimum potential in [H₅O₂]⁺ is very sensitive to electron correlation. At an OO distance of 2.74 Å it decreases from the HF value of 9.5 kcal/mole to about 7.0 kcal/mole. The results of the semi-empirical calculations agree well with the ab initio data as long as only relative effects are regarded. The partitioning of correlation energy into contributions of individual electron pairs is very similar for proton transfer in [H₅O₂]⁺ and for the dissociation of one OH bond in [H₃O]⁺. In this example the proton transfer appears as a superposition of two “contracted ionic dissociation” processes. An interpretation of the behaviour of correlation during these processes is presented.     

Ab initio calculations on the interaction of oxygen containing ligands with alkali cations. The system H2CO/Li+

Ab initio calculations on formaldehyde/Li⁺ complexes are presented. The most stable arrangement is characterized by an energy of interaction of 43.2 kcal/mole, C_2v symmetry and an oxygen—lithium distance of R_OLi = 1.77 Å. A detailed analysis of the electron density function gives proof of the electrostatic nature of the complexes H₂O/Li⁺ and H₂Co/Li⁺ and shows extensive mutual polarization. The failure of the semi-empirical method to predict the changes in electron density at the Li⁺ cation correctly is explained.     

Molecular orbital theoretical prediction of the isomeric products formed from reactions of arene oxides and related metabolites of polycyclic aromatic hydrocarbons

MO theoretical calculations based on the perturbational method of Dewar provide good correlation between predicted and observed structures of products formed during: (1) isomerization of arene oxides to phenols; (2) hydration and nucleophilic addition to arene oxides; and (3) dehydration of arene dihydrodiols. The method is equally applicable to the arene oxides, dihydrodiols, etc. derived from carcinogenic and noncarcinogenic polycyclic hydrocarbons. Extension to the related enzymatic reactions occurring during metabolism of carcinogenic hydrocarbons and to the reactions of the biologically active arene diolepoxides and aryloxirenes suggests the potential utility of this approach in predicting (a) metabolite structure and (b) the structural requirements for carcinogenic and mutagenic activity.     

Free energy perturbation calculations on models of active sites: Applications to adenosine deaminase inhibitors

Free energy perturbation calculations were performed to determine the free energy of binding associated with the presence of perhaps an unusual hydroxyl group in the transition state analog of nebularine, an inhibitor of the enzyme adenosine deaminase. The presence of a single hydroxyl group in this inhibitor has been found to contribute ?9.8 kcal/mol to the free energy of binding, with a 10⁸-fold increase in the binding affinity by the enzyme. In this work, we calculate the difference in solvation free energy for the 1,6-dihydropurine complex versus that of the 6-hydroxyl-1,6-dihydropurine complex to determine if this marked increase in binding affinity is attributed to an unusually hydrophobic hydroxyl group. The calculated ΔG associated for the solvation free energy is ?11.8 kcal/mol. This large change in the solvation free energy suggests that this hydroxyl is instead unusually hydrophilic and that the difference in free energy of interaction for the two inhibitors to the enzyme must be at least ca. 20 kcal/mol. Although the crystal structure for adenosine deaminase is currently not known, we attempt to mimic the nature of the active site by constructing models which simulate the enzyme-inhibitor complex. We present a first attempt at determining the change in free energy of binding for a system in which structural data for the enzyme is incomplete. To do this, we construct what we believe is a minimal model of the binding between adenosine deaminase and an inhibitor. The active site is simulated as a single charged carboxyl group which can form a hydrogen bond with the hydroxyl group of the analog. Two different carboxyl anion models are used. In the first model, the association is modeled between an acetic acid anion and the modified inhibitor. The second model consists of a hydrophobic amino acid pocket with an interior Glu residue in the active site. From these models we calculate the change in free energy of association and the overall change in free energy of binding. We calculate the free energies of interaction both in the absence and presence of water. We conclude from this that the presence of a single suitably placed-CO^?₂ group probably cannot explain the binding effect of the-OH group and that additional interactions will be found in the adenosine deaminase active site.     

The structure of the H3O+ (hydronium) ion

Near Hartree-Fock level ab initio molecular orbital calculations on H₃O⁺ and a minimum energy structure with θ(HOH) = 112.5° and r(OH) = 0.963 Å and an inversion barrier of 1.9 kcal/mole. By comparing these results to calculations on NH₃ and H₂O, where precise experimental geometries are known, we estimate the “true” geometry of isolated H₃O⁺ to have a structure with θ(HOH) = 110-112°, r(OH) = 0.97–0.98 Å and an inversion barrier of 2–3 kcal/mole. Our prediction for the proton affinity of water is ≈ 170 kcal/mole, which is somewhat smaller than the currently accepted value.     

On the oscillatory modes of quarks in baryons

The color bond structure of a quark-antiquark system is extended, in the long-range approximation, self-consistently to the baryonic three-quark bond structure for SU(3)_c and generally to the N-quark bond structure for SU(N)_c. The universal (N-independent) mass square eigenvalues for massless quarks are

$\text{M}{}^2\text{=(H}{}_{\mathrm{N}}\text{)}{}^2\text{?2m}{}_{\mathrm{\rho }}{}^2\text{∑}\text{α=1}\text{3N?3}\text{ν}{}_{\mathrm{\alpha }}\text{+}\text{constant}\text{, ν}{}_{\mathrm{\alpha }}\text{=0,1,2,…}$

.     

On uniform distribution of double sequences

In this paper a new notion of uniform distribution of double sequences is introduced. In the most elementary case a double sequence (x_mn)_{m, nN} is called uniformly distributed in this sense if all finite matrices of equal size appear with the same frequency as submatrices of the infinite matrix (x_mn). Some special double sequences with such a distribution property are studied and metric theorems are proved.     

“Deuterated” Raney nickel: deuteration (reduction) of alkenes,carbonyl compounds and aromatic rings. Proton-deuterium exchange of “activated” aliphatic and aromatic ring hydrogens.

In the absence of deuterium gas, the reaction of deuterated Raney nickel with benzene derivatives, carbonyl compounds, alkenes and activated methylene groups gave either fully deuterated (reduced) products or led to hydrogen-deuterium exchange.