Binding free energy, energy and entropy calculations using simple model systems

Free energy differences are calculated for a set of two model host molecules, binding acetone and methanol. Two active sites of different characteristics were constructed based on an artificially extended C60 fullerene molecule, possibly functionalised to include polar interactions in an otherwise apolar, spherical cavity. The model host systems minimise the necessary sampling of conformational space while still capturing key aspects of ligand binding. The estimates of the free energies are split up into energetic and entropic contributions, using three different approaches investigating the convergence behaviour. For these systems, a direct calculation of the total energy and entropy is more efficient than calculating the entropy from the temperature dependence of the free energy or from a direct thermodynamic integration formulation. Furthermore, the compensating surrounding–surrounding energies and entropies are split off by calculating reduced ligand-surrounding energies and entropies. These converge much more readily and lead to properties that are more straightforwardly interpreted in terms of molecular interactions and configurations. Even though not experimentally accessible, the reduced thermodynamic properties may prove highly relevant for computational drug design, as they may give direct insights into possibilities to further optimise ligand binding while optimisation in the surrounding–surrounding energy or entropy will exactly cancel and not lead to improved affinity.     

Calculation of solvation free energies of charged solutes using mixed cluster/continuum models

We derive a consistent approach for predicting the solvation free energies of charged solutes in the presence of implicit and explicit solvents. We find that some published methodologies make systematic errors in the computed free energies because of the incorrect accounting of the standard state corrections for water molecules or water clusters present in the thermodynamic cycle. This problem can be avoided by using the same standard state for each species involved in the reaction under consideration. We analyze two different thermodynamic cycles for calculating the solvation free energies of ionic solutes: (1) the cluster cycle with an n water cluster as a reagent and (2) the monomer cycle with n distinct water molecules as reagents. The use of the cluster cycle gives solvation free energies that are in excellent agreement with the experimental values obtained from studies of ion-water clusters. The mean absolute errors are 0.8 kcal/mol for H(+) and 2.0 kcal/mol for Cu(2+). Conversely, calculations using the monomer cycle lead to mean absolute errors that are >10 kcal/mol for H(+) and >30 kcal/mol for Cu(2+). The presence of hydrogen-bonded clusters of similar size on the left- and right-hand sides of the reaction cycle results in the cancellation of the systematic errors in the calculated free energies. Using the cluster cycle with 1 solvation shell leads to errors of 5 kcal/mol for H(+) (6 waters) and 27 kcal/mol for Cu(2+) (6 waters), whereas using 2 solvation shells leads to accuracies of 2 kcal/mol for Cu(2+) (18 waters) and 1 kcal/mol for H(+) (10 waters).     

Predicting aqueous free energies of solvation as functions of temperature

This work introduces a model, solvation model 6 with temperature dependence (SM6T), to predict the temperature dependence of aqueous free energies of solvation for compounds containing H, C, and O in the range 273-373 K. In particular, we extend solvation model 6 (SM6), which was previously developed (Kelly, C. P.; Cramer, C. J.; Truhlar, D. G. J. Chem. Theory Comput. 2005, 1, 1133) for predicting aqueous free energies of solvation at 298 K, to predict the variation of the free energy of solvation relative to 298 K. Also, we describe the database of experimental aqueous free energies of solvation for compounds containing H, C, and O that was used to parametrize and test the new model. SM6T partitions the temperature dependence of the free energy of solvation into two components: the temperature dependence of the bulk electrostatic contribution to the free energy of solvation, which is computed using the generalized Born equation, and the temperature dependence of first-solvation-shell effects which is modeled using a parametrized solvent-exposed surface-area-dependent term. We found that SM6T predicts the temperature dependence of aqueous free energies of solvation with a mean unsigned error of 0.08 kcal/mol over our entire database, whereas using the experimental value at 298 K produces a mean unsigned error of 0.53 kcal/mol.     

Calculation of the solvation free energy of neutral and ionic molecules in diverse solvents

The solvation free energy density (SFED) model was modified to extend its applicability and predictability. The parametrization process was performed with a large, diverse set of solvation free energies that included highly polar and ionic molecules. The mean absolute error for 1200 solvation free energies of the 379 neutral molecules in 9 organic solvents and water was 0.40 kcal/mol, and for 90 hydration free energies of ions was 1.7 kcal/mol. Overall, the calculated solvation free energies of a wide range of solute functional groups in diverse solvents were consistent with experimental data.     

Aqueous solvation free energies of ions and ion-water clusters based on an accurate value for the absolute aqueous solvation free energy of the proton

Thermochemical cycles that involve pKa, gas-phase acidities, aqueous solvation free energies of neutral species, and gas-phase clustering free energies have been used with the cluster pair approximation to determine the absolute aqueous solvation free energy of the proton. The best value obtained in this work is in good agreement with the value reported by Tissandier et al. (Tissandier, M. D.; Cowen, K. A.; Feng, W. Y.; Gundlach, E.; Cohen, M. J.; Earhart, A. D.; Coe, J. V. J. Phys. Chem. A 1998, 102, 7787), who applied the cluster pair approximation to a less diverse and smaller data set of ions. We agree with previous workers who advocated the value of -265.9 kcal/mol for the absolute aqueous solvation free energy of the proton. Considering the uncertainties associated with the experimental gas-phase free energies of ions that are required to use the cluster pair approximation as well as analyses of various subsets of data, we estimate an uncertainty for the absolute aqueous solvation free energy of the proton of no less than 2 kcal/mol. Using a value of -265.9 kcal/mol for the absolute aqueous solvation free energy of the proton, we expand and update our previous compilation of absolute aqueous solvation free energies; this new data set contains conventional and absolute aqueous solvation free energies for 121 unclustered ions (not including the proton) and 147 conventional and absolute aqueous solvation free energies for 51 clustered ions containing from 1 to 6 water molecules. When tested against the same set of ions that was recently used to develop the SM6 continuum solvation model, SM6 retains its previously determined high accuracy; indeed, in most cases the mean unsigned error improves when it is tested against the more accurate reference data.     

Calculation of solvation and binding free energy differences between VX-478 and its analogs by free energy perturbation and AMSOL methods

Summary VX-478 belongs to a novel class of HIV-1 protease inhibitors that are based on N,N-disubstituted benzene sulfonamides. Force field parameters for the N,N-dialkyl benzene sulfonamide moiety have been assembled from the literature and from our own ab initio calculations. These parameters were employed to calculate solvation and binding free energy differences between VX-478 and two analogs. The free energy perturbation method has been used to determine these differences using two approaches. In the first approach, intergroup interaction terms only were included in the calculation of free energies (as in most reports of free energy calculations using AMBER). In the second approach, both the inter- and intragroup interaction terms were included. The results obtained with the two approaches are in excellent agreement with each other and are also in close agreement with the experimental results. The solvation free energies of N,N-dimethyl benzene sulfonamide derivatives (truncated models of the inhibitors), calculated using continuum solvation (AMSOL) methods, are found to be in qualitative agreement with the experimental and free energy perturbation results. The binding and solvation free energy results are discussed in the context of structure-based drug design to show how physicochemical properties (for example aqueous solubilities and bioavailabilities) of these HIV-1 protease inhibitors were improved, while maintaining their inhibitory potency.     

Predicting the energy of the water exchange reaction and free energy of solvation for the uranyl ion in aqueous solution

The structures and vibrational frequencies of UO2(H2O)4(2+) and UO2(H2O)5(2+) have been calculated using density functional theory and are in reasonable agreement with experiment. The energies of various reactions were calculated at the density functional theory (DFT) and MP2 levels; the latter provides the best results. Self-consistent reaction field calculations in the PCM and SCIPCM approximations predicted the free energy of the water exchange reaction, UO2(H2O)4(2+) + H2O <--> UO2(H2O)5(2+). The calculated free energies of reaction are very sensitive to the choice of radii (O and H) and isodensity values in the PCM and SCIPCM models, respectively. Results consistent with the experimental HEXS value of -1.19 +/- 0.42 kcal/mol (within 1-3 kcal/mol) are obtained with small cavities. The structures and vibrational frequencies of the clusters with second solvation shell waters: UO2(H2O)4(H2O)8(2+), UO2(H2O)4(H2O)10(2+), UO2(H2O)4(H2O)11(2+), UO2(H2O)5(H2O)7(2+), and UO2(H2O)5(H2O)10(2+), were calculated and are in better agreement with experiment as compared to reactions involving only UO2(H2O)4(2+) and UO2(H2O)5(2+). The MP2 reaction energies for water exchange gave gas-phase results that agreed with experiment in the range -5.5 to +3.3 kcal/mol. The results were improved by inclusion of a standard PCM model with differences of -1.2 to +2.7 kcal/mol. Rearrangement reactions based on an intramolecular isomerization leading to a redistribution of water in the two shells provide good values in comparison to experiment with values of Delta G(exchange) from -2.2 to -0.5 kcal/mol so the inclusion of a second hydration sphere accounts for most solvation effects. Calculation of the free energy of solvation of the uranyl cation yielded an upper bound to the solvation energy of -410 +/- 5 kcal/mol, consistent with the best experimental value of -421 +/- 15 kcal/mol.     

Thermodynamic study of alkali sulphates: Calculation of enthalpy,entropy and Gibbs free energy

Joined with measures of specific heat at low temperature , a calorimetric study of alkali sulphates from ordinary temperature to 1500 K has allowed the calculation of enthalpy, entropy and Gibbs free energy of these salts. The value of entropy at melting point is compared with the value assessed by an acoustical method.     

Calculation of solvation free energy using RISM theory for peptide in salt solution

We developed a robust, highly efficient algorithm for solving the full reference interaction site model (RISM) equations for salt solutions near a solute molecule with many atomic sites. It was obtained as an extension of our previously reported algorithm for pure water near the solute molecule. The algorithm is a judicious hybrid of the Newton–Raphson and Picard methods. The most striking advantage is that the Jacobian matrix is just part of the input data and need not be recalculated at all. To illustrate the algorithm, we solved the full RISM equations for a dipeptide (NH₂(SINGLE BOND)CHCH₃(SINGLE BOND)CONH(SINGLE BOND)CHCH₃(SINGLE BOND)COOH) in a 1 M NaCl solution. The extended simple point charge (SPC/E) model was employed for water molecules. Two different conformations of the dipeptide were considered. It was assumed for each conformation that the dipeptide was present either as an un-ionized form or as a zwitterion. The structure of the salt solution near the dipeptide and salt effects on the solvation free energy were also discussed. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1724–1735, 1998     

Calculation of solvation free energy from quantum mechanical charge density and continuum dielectric theory

We have combined ultrasoft pseudopotential density functional theory utilizing plane wave basis with a Poisson-Boltzmann/solvent-accessible surface area (PB/SA) model to calculate the solvation free energy of small neutral organic compounds in water. The solute charge density obtained from density functional theory was directly used in solving the Poisson-Boltzmann equation to obtain the reaction field. The polarized electronic wave function of the solute in the solvent was solved by including the reaction field in the density functional Hamiltonian. The quantum mechanical and Poisson-Boltzmann equations were solved self-consistently until the charge density and reaction field converged. Using the solute charge density directly instead of a point-charge representation permitted asymmetric distortion and spreading out of the electron cloud. Because the electron density could leave the van der Waals surface to penetrate into the high-dielectric solvent, the reaction field generated by this density was generally smaller than that obtained by using the point-charge representation. In applying this model to calculate the solvation free energy of 31 small neutral organic molecules spanning a range of 25 kcal/mol, we obtained a root-mean-square error of only 1.3 kcal/mol if we allowed one adjustable parameter to shift the calculated solvation free energy.     

Calculation of the free energy of solvation for neutral analogs of amino acid side chains

The ability of the GROMOS96 force field to reproduce partition constants between water and two less polar solvents (cyclohexane and chloroform) for analogs of 18 of the 20 naturally occurring amino acids has been investigated. The estimations of the solvation free energies in water, in cyclohexane solution, and chloroform solution are based on thermodynamic integration free energy calculations using molecular dynamics simulations. The calculations show that while the force field reproduces the experimental solvation free energies of nonpolar analogs with reasonable accuracy the solvation free energies of polar analogs in water are systematically overestimated (too positive). The dependence of the calculated free energies on the atomic partial charges was also studied.     

Calculation of relative solvation free energy differences by thermodynamic perturbation method: Dependence of free energy results on simulation length

Molecular dynamics (MD) simulations in conjunction with the thermodynamic cycle perturbation approach has been used to calculate relative solvation free energies for acetone to acetaldehyde, acetone to pyruvic acid, acetone to 1,1,1-trifluoroacetone, acetone to 1,1,1-trichloroacetone, acetone to 2,3-butanedione, acetone to cyclopropanone, and formaldehyde hydrate to formaldehyde. To evaluate the dependence of relative solvation free energy convergence on MD simulation length and starting configuration two studies were performed. In the first study, each simulation started from the same well-equilibrated configuration and the length was varied from 153 to 1530 ps. In the second study, the relative solvation free energy differences were calculated starting from three different configurations and using 510 ps of MD simulation for each mutation. These results clearly indicate that, even for molecules with limited conformational flexibility, a simulation length of 510 ps or greater is required to obtain satisfactory convergence and, for the mutations of large structural changes between reactant and product, such as cyclopropanone to acetone, require much longer simulation lengths to achieve satisfactory convergence. These results also show that performing one long simulation is better than averaging results from three shortest simulations of the same length using different starting conformations. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1018–1027, 1999     

Box size effects are negligible for solvation free energies of neutral solutes

Hydration free energy calculations in explicit solvent have become an integral part of binding free energy calculations and a valuable test of force fields. Most of these simulations follow the conventional norm of keeping edge length of the periodic solvent box larger than twice the Lennard-Jones (LJ) cutoff distance, with the rationale that this should be sufficient to keep the interactions between copies of the solute to a minimum. However, for charged solutes, hydration free energies can exhibit substantial box size-dependence even at typical box sizes. Here, we examine whether similar size-dependence affects hydration of neutral molecules. Thus, we focused on two strongly polar molecules with large dipole moments, where any size-dependence should be most pronounced, and determined how their hydration free energies vary as a function of simulation box size. In addition to testing a variety of simulation box sizes, we also tested two LJ cut-off distances, 0.65 and 1.0 nm. We show from these simulations that the calculated hydration free energy is independent of the box-size as well as the LJ cut-off distance, suggesting that typical hydration free energy calculations of neutral compounds indeed need not be particularly concerned with finite-size effects as long as standard good practices are followed.     

Redox entropy of plastocyanin: developing a microscopic view of mesoscopic polar solvation

We report applications of analytical formalisms and molecular dynamics (MD) simulations to the calculation of redox entropy of plastocyanin metalloprotein in aqueous solution. The goal of our analysis is to establish critical components of the theory required to describe polar solvation at the mesoscopic scale. The analytical techniques include a microscopic formalism based on structure factors of the solvent dipolar orientations and density and continuum dielectric theories. The microscopic theory employs the atomistic structure of the protein with force-field atomic charges and solvent structure factors obtained from separate MD simulations of the homogeneous solvent. The MD simulations provide linear response solvation free energies and reorganization energies of electron transfer in the temperature range of 280-310 K. We found that continuum models universally underestimate solvation entropies, and a more favorable agreement is reported between the microscopic calculations and MD simulations. The analysis of simulations also suggests that difficulties of extending standard formalisms to protein solvation are related to the inhomogeneous structure of the solvation shell at the protein-water interface combining islands of highly structured water around ionized residues along with partial dewetting of hydrophobic patches. Quantitative theories of electrostatic protein hydration need to incorporate realistic density profile of water at the protein-water interface.     

Calculation of acidic dissociation constants in water: solvation free energy terms. Their accuracy and impact

Three polarizable continuum models, DPCM, CPCM, and IEFPCM, have been applied to calculate free energy differences for nine neutral compounds and their anions. On the basis of solvation free energies, the pK_a values were obtained for the compounds in question by using three thermodynamic cycles: one, involving the combined experimental and calculated data, as well as two other cycles solely with calculated data. This paper deals with the influence of factors such as the SCRF model applied, choice of a particular thermodynamic cycle, atomic radii used to build a cavity in the solvent (water), optimization of geometry in water, inclusion of electron correlation, and the dimension of the basis set on the solvation free energies and on the calculated pK_a values. Electronic supplementary material The online version of this article (doi: ) contains supplementary material, which is available to authorized users.     

Prediction of SAMPL2 aqueous solvation free energies and tautomeric ratios using the SM8, SM8AD,and SMD solvation models

We applied the solvation models SM8, SM8AD, and SMD in combination with the Minnesota M06-2X density functional to predict vacuum-water transfer free energies (Task 1) and tautomeric ratios in aqueous solution (Task 2) for the SAMPL2 test set. The bulk-electrostatic contribution to the free energy of solvation is treated as follows: SM8 employs the generalized Born model with the Coulomb field approximation, SM8AD employs the generalized Born approximation with asymmetric descreening, and SMD solves the nonhomogeneous Poisson equation. The non-bulk-electrostatic contribution arising from short-range interactions between the solute and solvent molecules in the first solvation shell is treated as a sum of terms that are products of geometry-dependent atomic surface tensions and solvent-accessible surface areas of the individual atoms of the solute. On average, three models tested in the present work perform similarly. In particular, we achieved mean unsigned errors of 1.3 (SM8), 2.0 (SM8AD), and 2.6 kcal/mol (SMD) for the aqueous free energies of 30 out of 31 compounds with known reference data involved in Task 1 and mean unsigned errors of 2.7 (SM8), 1.8 (SM8AD), and 2.4 kcal/mol (SMD) in the free energy differences (tautomeric ratios) for 21 tautomeric pairs in aqueous solution involved in Task 2.     

Calculation of the entropy and free energy by the hypothetical scanning Monte Carlo method: application to peptides

A new approach, the hypothetical scanning Monte Carlo (HSMC), for calculating the absolute entropy, S, and free energy, F, has been introduced recently and applied first to fluids (argon and water) and later to peptides. In this paper the method is further developed for peptide chains in vacuum. S is calculated from a given MC sample by reconstructing each sample conformation i step-by-step, i.e., calculating transition probabilities (TPs) for the dihedral and bond angles and fixing the related atoms at their positions. At step k of the process the chain's coordinates that have already been determined are kept fixed (the "frozen past") and TP(k) is obtained from a MC simulation of the "future" part of the chain whose TPs as yet have not been determined; when the process is completed the contribution of conformation i to the entropy is, S(i) approximately -ln Pi(k) TP(k). In a recent paper we studied polyglycine chains, modeled by the AMBER force field with constant bond lengths and bond angles (the rigid model). Decaglycine [(Gly)(10)] was studied in the helical, extended, and hairpin microstates, while (Gly)(16) was treated only in the first two microstates. In this paper the samples are increased and restudied, (Gly)(16) is also investigated in the hairpin microstate, and for (Gly)(10) approximations are tested where only part of the future is considered for calculating the TPs. We calculate upper and lower bounds for F and demonstrate that like for fluids, F can be obtained from multiple reconstructions of a single conformation. We also test a more realistic model of (Gly)(10) where the bond angles are allowed to move (the flexible model). Very accurate results for S and F are obtained which are compared to results obtained by the quasiharmonic approximation and the local states method. Thus, differences in entropy and free energy between the three microstates are obtained within errors of 0.1-0.3 kcal/mol. The HSMC method can be applied to a macromolecule with any degree of flexibility, ranging from local fluctuations to a random coil. The present results demonstrate that the difference in stability, DeltaF(mn)=F(m)-F(n) between significantly different microstates m and n, can be obtained from two simulations only without the need to resort to thermodynamic integration. Our long-term goal is to extend this method to any peptide and apply it to a peptide immersed in a box with explicit water.