{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,9,8]],"date-time":"2026-09-08T17:31:15Z","timestamp":1788888675053,"version":"build-2803163510"},"reference-count":51,"publisher":"MDPI AG","issue":"6","license":[{"start":{"date-parts":[[2021,3,17]],"date-time":"2021-03-17T00:00:00Z","timestamp":1615939200000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/501100001659","name":"Deutsche Forschungsgemeinschaft","doi-asserted-by":"publisher","award":["56091768"],"award-info":[{"award-number":["56091768"]}],"id":[{"id":"10.13039\/501100001659","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100008530","name":"European Regional Development Fund","doi-asserted-by":"publisher","award":["ZS\/2016\/04\/78155"],"award-info":[{"award-number":["ZS\/2016\/04\/78155"]}],"id":[{"id":"10.13039\/501100008530","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100004189","name":"Max-Planck-Gesellschaft","doi-asserted-by":"publisher","award":["MSD"],"award-info":[{"award-number":["MSD"]}],"id":[{"id":"10.13039\/501100004189","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["IJMS"],"abstract":"<jats:p>Small molecule receptor-binding is dominated by weak, non-covalent interactions such as van-der-Waals hydrogen bonding or electrostatics. Calculating these non-covalent ligand-receptor interactions is a challenge to computational means in terms of accuracy and efficacy since the ligand may bind in a number of thermally accessible conformations. The conformational rotamer ensemble sampling tool (CREST) uses an iterative scheme to efficiently sample the conformational space and calculates energies using the semi-empirical \u2018Geometry, Frequency, Noncovalent, eXtended Tight Binding\u2019 (GFN2-xTB) method. This combined approach is applied to blind predictions of the modes and free energies of binding for a set of 10 drug molecule ligands to the cucurbit[n]urils CB[8] receptor from the recent \u2018Statistical Assessment of the Modeling of Proteins and Ligands\u2019 (SAMPL) challenge including morphine, hydromorphine, cocaine, fentanyl, and ketamine. For each system, the conformational space was sufficiently sampled for the free ligand and the ligand-receptor complexes using the quantum chemical Hamiltonian. A multitude of structures makes up the final conformer-rotamer ensemble, for which then free energies of binding are calculated. For those large and complex molecules, the results are in good agreement with experimental values with a mean error of 3 kcal\/mol. The GFN2-xTB energies of binding are validated by advanced density functional theory calculations and found to be in good agreement. The efficacy of the automated QM sampling workflow allows the extension towards other complex molecular interaction scenarios.<\/jats:p>","DOI":"10.3390\/ijms22063078","type":"journal-article","created":{"date-parts":[[2021,3,17]],"date-time":"2021-03-17T11:48:22Z","timestamp":1615981702000},"page":"3078","update-policy":"https:\/\/doi.org\/10.3390\/mdpi_crossmark_policy","source":"Crossref","is-referenced-by-count":30,"title":["Accurate Receptor-Ligand Binding Free Energies from Fast QM Conformational Chemical Space Sampling"],"prefix":"10.3390","volume":"22","author":[{"given":"Esra","family":"Boz","sequence":"first","affiliation":[{"name":"Max Planck Institute for Dynamics of Complex Technical Systems, Molecular Simulations and Design Group, 39106 Magdeburg, Germany"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-7793-0052","authenticated-orcid":false,"given":"Matthias","family":"Stein","sequence":"additional","affiliation":[{"name":"Max Planck Institute for Dynamics of Complex Technical Systems, Molecular Simulations and Design Group, 39106 Magdeburg, Germany"}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"1968","published-online":{"date-parts":[[2021,3,17]]},"reference":[{"key":"ref_1","doi-asserted-by":"crossref","first-page":"531","DOI":"10.1146\/annurev-biophys-070816-033654","article-title":"Predicting binding free energies: Frontiers and benchmarks","volume":"46","author":"Mobley","year":"2017","journal-title":"Annu. Rev. Biophys."},{"key":"ref_2","doi-asserted-by":"crossref","first-page":"884","DOI":"10.1038\/s41573-020-0086-4","article-title":"Targeting receptor complexes: A new dimension in drug discovery","volume":"19","author":"Rosenbaum","year":"2020","journal-title":"Nat. Rev. Drug Discov."},{"key":"ref_3","doi-asserted-by":"crossref","unstructured":"Yang, S.Q., Ye, Q., Ding, J.J., Yin, M.Z., Lu, A.P., Chen, X., Hou, T.J., and Cao, D.S. (2020). Current advances in ligand-based target prediction. WIREs Comput. Mol. Sci., e1504.","DOI":"10.1002\/wcms.1504"},{"key":"ref_4","first-page":"902","article-title":"Recent advances in the design of choline kinase \u03b1 inhibitors and the molecular basis of their inhibition","volume":"41","year":"2020","journal-title":"Med. Res. Rev."},{"key":"ref_5","doi-asserted-by":"crossref","first-page":"3331","DOI":"10.1021\/ct5000296","article-title":"Binding free energy calculations for lead optimization: Assessment of their accuracy in an industrial drug design context","volume":"10","author":"Homeyer","year":"2014","journal-title":"J. Chem. Theory Comput."},{"key":"ref_6","doi-asserted-by":"crossref","first-page":"S1","DOI":"10.1002\/chir.20009","article-title":"Theoretical reassessment of Whelk-O1 as an enantioselective receptor for 1-(4-halogeno-phenyl)-1-ethylamine derivatives","volume":"16","author":"Hayes","year":"2004","journal-title":"Chirality"},{"key":"ref_7","doi-asserted-by":"crossref","first-page":"3572","DOI":"10.1021\/jp0373797","article-title":"Accurate calculations of ligand binding free energies: Chiral separation with enantioselective receptors","volume":"108","author":"Hayes","year":"2004","journal-title":"J. Phys. Chem. A"},{"key":"ref_8","doi-asserted-by":"crossref","first-page":"5674","DOI":"10.1002\/chem.200500444","article-title":"Novel enantioselective receptors for N-protected glutamate and aspartate","volume":"11","author":"Ragusa","year":"2005","journal-title":"Chem. A Eur. J."},{"key":"ref_9","doi-asserted-by":"crossref","first-page":"21","DOI":"10.1146\/annurev.biophys.36.040306.132550","article-title":"Calculation of protein-ligand binding affinities","volume":"36","author":"Gilson","year":"2007","journal-title":"Annu. Rev. Biophys. Biomol. Struct."},{"key":"ref_10","doi-asserted-by":"crossref","first-page":"449","DOI":"10.1016\/S1367-5931(97)80038-5","article-title":"Computational approaches to molecular recognition","volume":"1","author":"Lamb","year":"1997","journal-title":"Curr. Opin. Chem. Biol."},{"key":"ref_11","doi-asserted-by":"crossref","first-page":"1813","DOI":"10.1126\/science.1096361","article-title":"The many roles of computation in drug discovery","volume":"303","author":"Jorgensen","year":"2004","journal-title":"Science"},{"key":"ref_12","doi-asserted-by":"crossref","first-page":"889","DOI":"10.1021\/ar000033j","article-title":"Calculating structures and free energies of complex molecules: Combining molecular mechanics and continuum models","volume":"33","author":"Kollman","year":"2000","journal-title":"Acc. Chem. Res."},{"key":"ref_13","doi-asserted-by":"crossref","first-page":"601","DOI":"10.1007\/s10822-020-00290-5","article-title":"The SAMPL6 sampling challenge: Assessing the reliability and efficiency of binding free energy calculations","volume":"34","author":"Rizzi","year":"2020","journal-title":"J. Comput. Aided Mol. Des."},{"key":"ref_14","doi-asserted-by":"crossref","first-page":"937","DOI":"10.1007\/s10822-018-0170-6","article-title":"Overview of the SAMPL6 host-guest binding affinity prediction challenge","volume":"32","author":"Rizzi","year":"2018","journal-title":"J. Comput. Aided Mol. Des."},{"key":"ref_15","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1007\/s10822-016-9974-4","article-title":"Overview of the SAMPL5 host-guest challenge: Are we doing better?","volume":"31","author":"Yin","year":"2017","journal-title":"J. Comput. Aided Mol. Des."},{"key":"ref_16","doi-asserted-by":"crossref","first-page":"475","DOI":"10.1007\/s10822-012-9554-1","article-title":"Blind prediction of host\u2013guest binding affinities: A new SAMPL3 challenge","volume":"26","author":"Muddana","year":"2012","journal-title":"J. Comput. Aided Mol. Des."},{"key":"ref_17","doi-asserted-by":"crossref","first-page":"540","DOI":"10.1021\/ja993376p","article-title":"New cucurbituril homologues:\u2009 Syntheses, isolation, characterization, and X-ray crystal structures of cucurbit[n]uril (n = 5, 7, and 8)","volume":"122","author":"Kim","year":"2000","journal-title":"J. Am. Chem. Soc."},{"key":"ref_18","doi-asserted-by":"crossref","first-page":"12320","DOI":"10.1021\/acs.chemrev.5b00341","article-title":"Cucurbituril-based molecular recognition","volume":"115","author":"Barrow","year":"2015","journal-title":"Chem. Rev."},{"key":"ref_19","doi-asserted-by":"crossref","first-page":"2842","DOI":"10.1021\/jp2110067","article-title":"Cucurbit[8]uril mediated donor-acceptor ternary complexes: A model system for studying charge-transfer interactions","volume":"116","author":"Biedermann","year":"2012","journal-title":"J. Phys. Chem. B"},{"key":"ref_20","doi-asserted-by":"crossref","first-page":"15959","DOI":"10.1021\/ja055013x","article-title":"The cucurbit[n]uril family: Prime components for self-sorting systems","volume":"127","author":"Liu","year":"2005","journal-title":"J. Am. Chem. Soc."},{"key":"ref_21","doi-asserted-by":"crossref","first-page":"150","DOI":"10.1080\/10610278.2018.1516885","article-title":"Cucurbit 8 uril center dot guest complexes: Blinded dataset for the SAMPL6 challenge","volume":"31","author":"Murkli","year":"2019","journal-title":"Supramol. Chem."},{"key":"ref_22","doi-asserted-by":"crossref","first-page":"511","DOI":"10.1007\/s10822-020-00278-1","article-title":"LogP prediction performance with the SMD solvation model and the M06 density functional family for SAMPL6 blind prediction challenge molecules","volume":"34","author":"Guan","year":"2020","journal-title":"J. Comput. Aided Mol. Des."},{"key":"ref_23","doi-asserted-by":"crossref","first-page":"589","DOI":"10.1007\/s10822-020-00294-1","article-title":"SAMPL6 host-guest binding affinities and binding poses from spherical-coordinates-biased simulations","volume":"34","author":"Sun","year":"2020","journal-title":"J. Comput. Aided Mol. Des."},{"key":"ref_24","doi-asserted-by":"crossref","first-page":"1047","DOI":"10.1007\/s10822-018-0154-6","article-title":"Blinded predictions of standard binding free energies: Lessons learned from the SAMPL6 challenge","volume":"32","author":"Papadourakis","year":"2018","journal-title":"J. Comput. Aided Mol. Des."},{"key":"ref_25","doi-asserted-by":"crossref","first-page":"1097","DOI":"10.1007\/s10822-018-0159-1","article-title":"SAMPL6 host-guest challenge: Binding free energies via a multistep approach","volume":"32","author":"Eken","year":"2018","journal-title":"J. Comput. Aided Mol. Des."},{"key":"ref_26","doi-asserted-by":"crossref","first-page":"1013","DOI":"10.1007\/s10822-018-0153-7","article-title":"Detailed potential of mean force studies on host-guest systems from the SAMPL6 challenge","volume":"32","author":"Song","year":"2018","journal-title":"J. Comput. Aided Mol. Des."},{"key":"ref_27","doi-asserted-by":"crossref","first-page":"527","DOI":"10.1007\/s10822-011-9524-z","article-title":"Binding affinities in the SAMPL3 trypsin and host-guest blind tests estimated with the MM\/PBSA and LIE methods","volume":"26","author":"Mikulskis","year":"2012","journal-title":"J. Comput. Aided Mol. Des."},{"key":"ref_28","doi-asserted-by":"crossref","first-page":"9955","DOI":"10.1002\/chem.201200497","article-title":"Supramolecular binding thermodynamics by dispersion-corrected density functional theory","volume":"18","author":"Grimme","year":"2012","journal-title":"Chem. A Eur. J."},{"key":"ref_29","doi-asserted-by":"crossref","first-page":"3431","DOI":"10.1021\/jp411616b","article-title":"Blind prediction of binding affinities for charged supramolecular host\u2013guest systems: Achievements and shortcomings of DFT-D3","volume":"118","author":"Sure","year":"2014","journal-title":"J. Phys. Chem. B"},{"key":"ref_30","doi-asserted-by":"crossref","first-page":"7169","DOI":"10.1039\/C9CP06869D","article-title":"Automated exploration of the low-energy chemical space with fast quantum chemical methods","volume":"22","author":"Pracht","year":"2020","journal-title":"Phys. Chem. Chem. Phys."},{"key":"ref_31","doi-asserted-by":"crossref","first-page":"3397","DOI":"10.1039\/b704725h","article-title":"Double-hybrid density functionals with long-range dispersion corrections: Higher accuracy and extended applicability","volume":"9","author":"Schwabe","year":"2007","journal-title":"Phys. Chem. Chem. Phys."},{"key":"ref_32","doi-asserted-by":"crossref","first-page":"3098","DOI":"10.1002\/chem.202004380","article-title":"In vitro and in vivo sequestration of phencyclidine by Me4Cucurbit [8] uril","volume":"27","author":"Murkli","year":"2021","journal-title":"Chem. A Eur. J."},{"key":"ref_33","unstructured":"Isaacs, L. (2021, March 13). CB8-DOA-SAMPL-Answer-Sheet-20201014.Pdf. Available online: https:\/\/github.com\/samplchallenges\/SAMPL8\/blob\/master\/host_guest\/Analysis\/ExperimentalMeasurements\/CB8-DOA-SAMPL-Answer-Sheet-20201014.pdf."},{"key":"ref_34","doi-asserted-by":"crossref","first-page":"1652","DOI":"10.1021\/acs.jctc.8b01176","article-title":"GFN2-xTB\u2014an accurate and broadly parametrized self-consistent tight-binding quantum chemical method with multipole electrostatics and density-dependent dispersion contributions","volume":"15","author":"Bannwarth","year":"2019","journal-title":"J. Chem. Theory Comput."},{"key":"ref_35","unstructured":"Mobley, D.L., and Amezcua, M. (2021, February 21). The SAMPL8 CB8 \u201cDrugs of Abuse\u201d Challenge. Available online: https:\/\/github.com\/samplchallenges\/SAMPL8\/blob\/master\/host_guest\/CB8\/README.md."},{"key":"ref_36","unstructured":"Grimme, S., and Pracht, P. (2020). CREST Conformer-Rotamer Ensemble Sampling Tool Based on the GFN-, Mulliken Center for Theoretical Chemistry, University of Bonn."},{"key":"ref_37","doi-asserted-by":"crossref","first-page":"e1493","DOI":"10.1002\/wcms.1493","article-title":"Extended tight-binding quantum chemistry methods","volume":"11","author":"Bannwarth","year":"2020","journal-title":"WIREs Comput. Mol. Sci."},{"key":"ref_38","unstructured":"Grimme, S. (2020). Semi-empirical Extended Tight-Binding Program Package xtb.v- 6.3.2, Mulliken Center for Theoretical Chemistry, University of Bonn."},{"key":"ref_39","doi-asserted-by":"crossref","first-page":"3636","DOI":"10.1021\/acs.jpcb.0c00549","article-title":"Quantum chemical calculation of molecular and periodic peptide and protein structures","volume":"124","author":"Schmitz","year":"2020","journal-title":"J. Phys. Chem. B"},{"key":"ref_40","doi-asserted-by":"crossref","first-page":"2847","DOI":"10.1021\/acs.jctc.9b00143","article-title":"Exploration of chemical compound, conformer, and reaction space with meta-dynamics simulations based on tight-binding quantum chemical calculations","volume":"15","author":"Grimme","year":"2019","journal-title":"J. Chem. Theory Comput."},{"key":"ref_41","doi-asserted-by":"crossref","first-page":"1463","DOI":"10.1002\/jcc.20078","article-title":"Accurate description of van der Waals complexes by density functional theory including empirical corrections","volume":"25","author":"Grimme","year":"2004","journal-title":"J. Comput. Chem."},{"key":"ref_42","doi-asserted-by":"crossref","first-page":"5029","DOI":"10.1063\/1.478401","article-title":"Assessment of the Perdew\u2013Burke\u2013Ernzerhof exchange-correlation functional","volume":"110","author":"Ernzerhof","year":"1999","journal-title":"J. Chem. Phys."},{"key":"ref_43","doi-asserted-by":"crossref","first-page":"034108","DOI":"10.1063\/1.2148954","article-title":"Semiempirical hybrid density functional with perturbative second-order correlation","volume":"124","author":"Grimme","year":"2006","journal-title":"J. Chem. Phys."},{"key":"ref_44","doi-asserted-by":"crossref","first-page":"291","DOI":"10.1021\/ct100466k","article-title":"Efficient and accurate double-hybrid-meta-gga density functionals\u2014evaluation with the extended GMTKN30 database for general main group thermochemistry, kinetics, and noncovalent interactions","volume":"7","author":"Goerigk","year":"2011","journal-title":"J. Chem. Theory Comput."},{"key":"ref_45","doi-asserted-by":"crossref","first-page":"73","DOI":"10.1002\/wcms.81","article-title":"The Orca program system","volume":"2","author":"Neese","year":"2012","journal-title":"WIREs Comput. Mol. Sci."},{"key":"ref_46","doi-asserted-by":"crossref","first-page":"e1327","DOI":"10.1002\/wcms.1327","article-title":"Software update: The Orca program system, version 4.0","volume":"8","author":"Neese","year":"2018","journal-title":"WIREs Comput. Mol. Sci."},{"key":"ref_47","doi-asserted-by":"crossref","first-page":"1740","DOI":"10.1002\/jcc.10318","article-title":"An improvement of the resolution of the identity approximation for the formation of the Coulomb matrix","volume":"24","author":"Neese","year":"2003","journal-title":"J. Comput. Chem."},{"key":"ref_48","doi-asserted-by":"crossref","first-page":"3297","DOI":"10.1039\/b508541a","article-title":"Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy","volume":"7","author":"Weigend","year":"2005","journal-title":"Phys. Chem. Chem. Phys."},{"key":"ref_49","doi-asserted-by":"crossref","first-page":"283","DOI":"10.1016\/0009-2614(95)00621-A","article-title":"Auxiliary basis sets to approximate coulomb potentials","volume":"240","author":"Eichkorn","year":"1995","journal-title":"Chem. Phys. Lett."},{"key":"ref_50","doi-asserted-by":"crossref","first-page":"1456","DOI":"10.1002\/jcc.21759","article-title":"Effect of the damping function in dispersion corrected density functional theory","volume":"32","author":"Grimme","year":"2011","journal-title":"J. Comput. Chem."},{"key":"ref_51","doi-asserted-by":"crossref","first-page":"6378","DOI":"10.1021\/jp810292n","article-title":"Universal solvation model based on solute electron density and on a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions","volume":"113","author":"Marenich","year":"2009","journal-title":"J. Phys. Chem. B"}],"container-title":["International Journal of Molecular Sciences"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/www.mdpi.com\/1422-0067\/22\/6\/3078\/pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,10,11]],"date-time":"2025-10-11T05:37:10Z","timestamp":1760161030000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.mdpi.com\/1422-0067\/22\/6\/3078"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2021,3,17]]},"references-count":51,"journal-issue":{"issue":"6","published-online":{"date-parts":[[2021,3]]}},"alternative-id":["ijms22063078"],"URL":"https:\/\/doi.org\/10.3390\/ijms22063078","relation":{},"ISSN":["1422-0067"],"issn-type":[{"value":"1422-0067","type":"electronic"}],"subject":[],"published":{"date-parts":[[2021,3,17]]}}}