Solvation free energies can now be calculated precisely from molecular simulations, providing a valuable test of the energy functions underlying these simulations. Here, we briefly review “alchemical” approaches for calculating the solvation free energies of small, neutral organic molecules from molecular simulations, and illustrate by applying them to calculate aqueous solvation free energies (hydration free energies). These approaches use a non-physical pathway to compute free energy differences from a simulation or set of simulations and appear to be a particularly robust and general-purpose approach for this task. We also present an update (version 0.5) to our FreeSolv database of experimental and calculated hydration free energies of neutral compounds and provide input files in formats for several simulation packages. This revision to FreeSolv provides calculated values generated with a single protocol and software version, rather than the heterogeneous protocols used in the prior version of the database. We also further update the database to provide calculated enthalpies and entropies of hydration and some experimental enthalpies and entropies, as well as electrostatic and nonpolar components of solvation free energies.

The hydration structure of Li+ has been studied by means of hybrid quantum-mechanical molecular mechanical molecular dynamics simulations at Hartree–Fock and density-functional level of theory. The size of the quantum-mechanical region and the form of the potential function are shown to be of crucial importance for reliable results. Radial distribution functions, coordination number distributions, and various angular distributions have been used to discuss details of the hydration structure, together with bond lengths and bond angles of the water molecules in the first hydration shell. The lithium ion is found to be mainly fourfold coordinated with some smaller amounts of fivefold coordination. The lithium–water cluster exhibits a nearly perfect tetrahedral but still very flexible structure, in which coordinated water molecules are considerably tilted away from planarity. Water molecules in the first hydration shell are shown to be considerably polarized compared to gas-phase structures.

Alchemical free energy calculations are an increasingly important modern simulation technique to calculate free energy changes on binding or solvation. Contemporary molecular simulation software such as AMBER, CHARMM, GROMACS, and SOMD include support for the method. Implementation details vary among those codes, but users expect reliability and reproducibility, i.e., for a given molecular model and set of force field parameters, comparable free energy differences should be obtained within statistical bounds regardless of the code used. Relative alchemical free energy (RAFE) simulation is increasingly used to support molecule discovery projects, yet the reproducibility of the methodology has been less well tested than its absolute counterpart. Here we present RAFE calculations of hydration free energies for a set of small organic molecules and demonstrate that free energies can be reproduced to within about 0.2 kcal/mol with the aforementioned codes. Absolute alchemical free energy simulations have been carried out as a reference. Achieving this level of reproducibility requires considerable attention to detail and package-specific simulation protocols, and no universally applicable protocol emerges. The benchmarks and protocols reported here should be useful for the community to validate new and future versions of software for free energy calculations.

FESetup is a new pipeline tool which can be used flexibly within larger workflows. The tool aims to support fast and easy setup of alchemical free energy simulations for molecular simulation packages such as AMBER, GROMACS, Sire, or NAMD. Post-processing methods like MM-PBSA and LIE can be set up as well. Ligands are automatically parametrized with AM1-BCC, and atom mappings for a single topology description are computed with a maximum common substructure search (MCSS) algorithm. An abstract molecular dynamics (MD) engine can be used for equilibration prior to free energy setup or standalone. Currently, all modern AMBER force fields are supported. Ease of use, robustness of the code, and automation where it is feasible are the main development goals. The project follows an open development model, and we welcome contributions.

A classical molecular dynamics simulation including three-body corrections was compared with combined ab initio quantum mechanics/molecular mechanics molecular dynamics simulations (QM/MM–MD), which were carried out at Hartree–Fock (HF) and density functional theory (DFT) level for Ca2+ in water. In the QM approach the region of primary interest—the first hydration sphere of the calcium ion—was treated by Born–Oppenheimer quantum mechanics, while the rest of the system was described by classical pair potentials. Coordination numbers of 7.1, 7.6, and 8.1 were found in the classical, the HF, and the DFT simulation, respectively, using the same double-ζ basis set in both QM methods. The CPU time for one DFT step was about 50% above the time for a HF step, but due to a smaller number of steps needed for equilibration in the DFT case, there was no significant difference in the overall simulation time.