An accurate description of solvation effects is of high importance in modeling biomolecular systems. Our main interest is to find an accurate yet efficient solvation model for semiempirical quantum-mechanical methods applicable to large protein-ligand complexes in the context of computer-aided drug design. We present a survey of readily available methods and a new reparametrization of the COSMO solvent model for PM6 and PM7 calculations in MOPAC. We have tested the reparametrized method on validation data sets of small drug-like molecules for which experimental solvation free energies are available as well as on a set of large model systems of the active site of carbonic anhydrase II interacting with a series of ligands for which experimental affinity values are known. In both cases, there is a significant improvement in accuracy after the reparametrization and the addition of a nonpolar term to the COSMO solvent model.

Department of Physical and Macromolecular Chemistry Faculty of Science Charles University Hlavova 8 128 40 Prague 2 Czech Republic

Institute of Organic Chemistry and Biochemistry Czech Academy of Sciences 166 10 Prague Czech Republic

References provided by Crossref.org

End-Point Affinity Estimation of Galectin Ligands by Classical and Semiempirical Quantum Mechanical Potentials

Nature-Inspired Gallinamides Are Potent Antischistosomal Agents: Inhibition of the Cathepsin B1 Protease Target and Binding Mode Analysis

SQM2.20: Semiempirical quantum-mechanical scoring function yields DFT-quality protein-ligand binding affinity predictions in minutes

Highly potent inhibitors of cathepsin K with a differently positioned cyanohydrazide warhead: structural analysis of binding mode to mature and zymogen-like enzymes

Structural and Thermodynamic Analysis of the Resistance Development to Pimodivir (VX-787), the Clinical Inhibitor of Cap Binding to PB2 Subunit of Influenza A Polymerase