No abstract

No abstract

We present self-consistent reaction field (SCRF) calculations, utilizing correlated ab initio quantum mechanics, of aqueous solvation free energies for a large data base of molecular solutes. We identify a subset of chemical functional groups for which there are systematic deviations in the comparison of theory and experiment; furthermore, for one case which has been extensively investigated, methylated amines, similar deviations appear in explicit solvent free energy perturbation calculations employing several commonly used molecular mechanics potential functions. By carrying out high-level ab initio quantum chemical calculations of hydrogen-bonding energies of the solutes to a water molecule, we arrive at a coherent explanation of the disagreements between theory and experiment, namely, that hydrogen-bonding energies are in some cases poorly correlated with classical electrostatic interaction energies. We show that the deviation in hydrogen-bonding energies of a solute from a reference molecule (for which there is good agreement between the SCRF calculations and experiment) is an excellent predictor of the errors made for that solute in the SCRF calculations. A new SCRF model is developed in which short-range empirical corrections, based upon solvent accessibility, are made for these chemical functional groups; this reduces the mean error of the calculated solvation free energies for the entire data base by a factor of ∼2, to 0.37 kcal/mol. These results have significant implications for the accuracy of explicit solvent potential functions as well as dielectric continuum models. Finally, we also identify cases where the observed discrepancies in solvation free energies cannot be explained by pair hydrogen-bonding results and suggest problems here that may be specific to dielectric continuum theory.

No abstract

The energy of dimerization of two N-methylacetamide (NMA) molecules in vacuum is calculated using density functional theory. Natural orbital analysis suggests that the dimerization energy of -6.6 kcal/mol is predominantly due to the (NsH‚‚‚OdC) donor-acceptor interaction. The gas phase to water hydration free energies and the free energies of transfer from the aqueous phase to liquid alkane of hydrogen bonded, (NsH‚‚‚OdC), and nonbonded, (NsH,OdC), groups are calculated using a continuum solvent model. On the basis of these calculations, we estimate the free energy of forming an amide hydrogen bond in the context of the NMA dimer in water and in liquid alkane as ∼-1 and ∼-5 kcal/mol, respectively. The relevance of these calculations to processes such as protein folding and membrane insertion of proteins is discussed.