uu.seUppsala University Publications
Change search
ReferencesLink to record
Permanent link

Direct link
Generalized Born and Explicit Solvent Models for Free Energy Calculations in Organic Solvents: Cyclodextrin Dimerization
Univ Sci & Technol Beijing, Sch Chem & Biol Engn, Dept Biol Sci & Engn, Beijing 100083, Peoples R China.;Beijing Univ Chem Technol, Dept Biochem Engn, Beijing Key Lab Bioproc, Beijing 100029, Peoples R China..
Beijing Univ Chem Technol, Dept Biochem Engn, Beijing Key Lab Bioproc, Beijing 100029, Peoples R China..
Uppsala University, Disciplinary Domain of Science and Technology, Biology, Department of Cell and Molecular Biology, Computational and Systems Biology.
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
Abstract [en]

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.
National Category
Physical Chemistry
URN: urn:nbn:se:uu:diva-269249DOI: 10.1021/acs.jctc.5b00620ISI: 000364614000006PubMedID: 26574308OAI: oai:DiVA.org:uu-269249DiVA: diva2:884977
Swedish Research Council, 2013-5947
Available from: 2015-12-17 Created: 2015-12-15 Last updated: 2015-12-17Bibliographically approved

Open Access in DiVA

No full text

Other links

Publisher's full textPubMed

Search in DiVA

By author/editor
Van der Spoel, David
By organisation
Computational and Systems Biology
In the same journal
Journal of Chemical Theory and Computation
Physical Chemistry

Search outside of DiVA

GoogleGoogle Scholar
The number of downloads is the sum of all downloads of full texts. It may include eg previous versions that are now no longer available

Altmetric score

Total: 57 hits
ReferencesLink to record
Permanent link

Direct link