Drug-like Bioactive Structures and Conformational Coverage with the LigPrep/ConfGen Suite: Comparison to Programs MOE and Catalyst

Abstract

Computational conformational sampling underpins many aspects of small molecule modeling and design in pharmaceutical work. This work examined in detail the widely distributed LigPrep/ConfGen software suite and the conformational models it produces for drug-like compounds. We also compare LigPrep/ConfGen to MOE and Catalyst. Tests of the conformational sampling protocols included the reproduction of known bioactive structures of ligands, characterization of the size, coverage and diversity of the output conformational models, and relative computation times. The present tests will help the user to make informed choices among the predefined ConfGen protocols (Very fast, Fast, Intermediate, and Comprehensive), and the adjustable input parameters. The parameters governing the initial compound preparation (LigPrep) and the subsequent conformational sampling were explored. This analysis has led to a new protocol called "ConfGen Optimized", which improves upon the predefined protocols. ConfGen Optimized is computationally tractable and reproduced 80% of the bioactive structures within 1 A, versus 66% for the default ConfGen Fast protocol. We also addressed the issue of the reproduction of compact/folded bioactive structures by ConfGen. It involved the compilation of a new set of 50 folded diverse drug-like bioactive structures. This indicates that heuristics penalizing folded conformers hinder reproduction of some binding modes. Overall, ConfGen offers great flexibility of use and provides a valuable addition to the molecular modeling toolbox.

Energy-Based Pharmacophore Modeling, Virtual Screening, and Molecular Dynamics to Identify Potential Inhibitors for Glycogen Synthase Kinase 3 Beta

“…Water molecules which do not participate in interactions were removed, the bond order was assigned, and hydrogen atoms were added. Optimization and minimization of the protein were carried out by applying the OPLS-2005 force field [15].…”

“…All the ligands were optimized using OPLS_2005 force field, and conformations were generated using Epik in the range of pH 7.0±2.0 [15,16]. The structures were desalted and tautomerized.…”

Human Rab8b Protein as a Cancer Target - An In Silico Study

“…Ligand preparation: Schrödinger's LigPrep [41] is a tool used to generate the 3D structures from 2D representations, search for tautomers, steric isomers, and ionization states of different ligand molecules, followed by geometry optimization of ligands based upon the OPLS-2005 force field [42]. The process of ligand preparation involves conserving the specified chiralities to generate minimum 4 stereoisomers per ligand, using default conditions at pH 7.0 ± 2.0, to obtain several conformers.…”

Targeting the ubiquitin-conjugating enzyme E2D4 for cancer drug discovery–a structure-based approach

“…com). The small structural ligand molecules are subjected to the LigPrep module [55] (LigPrep, version 5.6, Schrödinger, LLC, New York, NY, 2010) of the Schrödinger suite for ligand preparation, which is set to generate five stereo isomeric structures per ligand using OPLS_2005 force field. For each ligand, multiple ionization states, tautomers, stereoisomers and five low energy conformers are generated at pH 7.0 ± 2.0.…”