Generalized Born and Explicit Solvent Models for Free Energy Calculations in Organic Solvents: Cyclodextrin Dimerization
2015 (English)In: Journal of Chemical Theory and Computation, ISSN 1549-9618, E-ISSN 1549-9626, Vol. 11, no 11, 5103-5113 p.Article in journal (Refereed) PublishedText
Evaluation of solvation (binding) free energies with implicit solvent models in different dielectric environments for biological simulations as well as high throughput ligand screening remain challenging endeavors. In order to address how well implicit solvent models approximate explicit ones we examined four generalized Born models (GB(still), GB(HCT), GB(OBC)I, and GB(OBC)II) for determining the dimerization free energy (Delta G(0)) of beta-cyclodextrin monomers in 17 implicit solvents with dielectric constants (D) ranging from 5 to 80 and compared the results to previous free energy calculations with explicit solvents (Zhang et al. J. Phys. Chem. B 2012, 116, 12684-12693). The comparison indicates that neglecting the environmental dependence of Born radii appears acceptable for such calculations involving cyclodextrin and that the GB(still) and GB(OBC)I models yield a reasonable estimation of Delta G(0), although the details of binding are quite different from explicit solvents. Large discrepancies between implicit and explicit solvent models occur in high-dielectric media with strong hydrogen bond (HB) interruption properties. Delta G(0) with the GB models is shown to correlate strongly to 2(D-1)/(2D+1) (R-2 similar to 0.90) in line with the Onsager reaction field (Onsager J. Am. Chem. Soc. 1936, 58, 1486-1493) but to be very sensitive to D (D < 10) as well. Both high-dielectric environments where hydrogen bonds are of interest and low-dielectric media such as protein binding pockets and membrane interiors therefore need to be considered with caution in GB-based calculations. Finally, a literature analysis of Gibbs energy of solvation of small molecules in organic liquids shows that the Onsager relation does not hold for real molecules since the correlation between Delta G(0) and 2(D-1)/(2D+1) is low for most solutes. Interestingly, explicit solvent calculations of the solvation free energy (Zhang et al. J. Chem. Inf. Model. 2015, SS, 1192-1201) reproduce the weak experimental correlations with 2(D-1)/(2D+1) very well.
Place, publisher, year, edition, pages
2015. Vol. 11, no 11, 5103-5113 p.
IdentifiersURN: urn:nbn:se:uu:diva-269249DOI: 10.1021/acs.jctc.5b00620ISI: 000364614000006PubMedID: 26574308OAI: oai:DiVA.org:uu-269249DiVA: diva2:884977
FunderSwedish Research Council, 2013-5947