Use of Solvent Mapping for Characterizing the Binding Site and for Predicting the Inhibition of the Human Ether-á-Go-Go-Related K+ Channel.

Molecular dynamics was used to optimize the droperidol-hERG complex obtained from docking. To accommodate the inhibitor, residues T623, S624, V625, G648, Y652, and F656 did not move significantly during the simulation, while F627 moved significantly. Binding sites in cryo-EM structures and in structures obtained from molecular dynamics simulations were characterized using solvent mapping and Atlas ligands, which were negative images of the binding site, were generated. Atlas ligands were found to be useful for identifying human ether-á-go-go-related potassium channel (hERG) inhibitors by aligning compounds to them or by guiding the docking of compounds in the binding site. A molecular dynamics optimized structure of hERG led to improved predictions using either compound alignment to the Atlas ligand or docking. The structure was also found to be suitable to define a strategy for lowering inhibition based on the proposed binding mode of compounds in the channel.

The influence of calculated physicochemical properties of compounds on their ADMET profiles.

We analyzed the influence of calculated physicochemical properties of more than 20,000 compounds on their P-gp and BCRP mediated efflux, microsomal stability, hERG inhibition, and plasma protein binding. Our goal was to provide guidance for designing compounds with desired pharmacokinetic profiles. Our analysis showed that compounds with ClogP less than 3 and molecular weight less than 400 will have high microsomal stability and low plasma protein binding. Compounds with logD less than 2.2 and/or basic pKa larger than 5.3 are likely to be BCRP substrates and compounds with basic pKa less than 5.2 and/or acidic pKa less than 13.4 are less likely to inhibit hERG. Based on these results, compounds with MW < 400, ClogP < 3, basic pKa < 5.2 and acidic pKa < 13.4 are likely to have good bioavailability and low hERG inhibition.

The statistics of virtual screening and lead optimization.

Analytic formulae are used to estimate the error for two virtual screening metrics, enrichment factor and area under the ROC curve. These analytic error estimates are then compared to bootstrapping error estimates, and shown to have excellent agreement with respect to area under the ROC curve and good agreement with respect to enrichment factor. The major advantage of the analytic formulae is that they are trivial to calculate and depend only on the number of actives and inactives and the measured value of the metric, information commonly reported in papers. In contrast to this, the bootstrapping method requires the individual compound scores. Methods for converting the error, which is calculated as a variance, into more familiar error bars are also discussed.

FRED and HYBRID docking performance on standardized datasets.

The docking performance of the FRED and HYBRID programs are evaluated on two standardized datasets from the Docking and Scoring Symposium of the ACS Spring 2011 national meeting. The evaluation includes cognate docking and virtual screening performance. FRED docks 70 % of the structures to within 2 Å in the cognate docking test. In the virtual screening test, FRED is found to have a mean AUC of 0.75. The HYBRID program uses a modified version of FRED's algorithm that uses both ligand- and structure-based information to dock molecules, which increases its mean AUC to 0.78. HYBRID can also implicitly account for protein flexibility by making use of multiple crystal structures. Using multiple crystal structures improves HYBRID's performance (mean AUC 0.80) with a negligible increase in docking time (~15 %).

FRED pose prediction and virtual screening accuracy.

Results of a previous docking study are reanalyzed and extended to include results from the docking program FRED and a detailed statistical analysis of both structure reproduction and virtual screening results. FRED is run both in a traditional docking mode and in a hybrid mode that makes use of the structure of a bound ligand in addition to the protein structure to screen molecules. This analysis shows that most docking programs are effective overall but highly inconsistent, tending to do well on one system and poorly on the next. Comparing methods, the difference in mean performance on DUD is found to be statistically significant (95% confidence) 61% of the time when using a global enrichment metric (AUC). Early enrichment metrics are found to have relatively poor statistical power, with 0.5% early enrichment only able to distinguish methods to 95% confidence 14% of the time.

Gaussian docking functions.

A shape-based Gaussian docking function is constructed which uses Gaussian functions to represent the shapes of individual atoms. A set of 20 trypsin ligand-protein complexes are drawn from the Protein Data Bank (PDB), the ligands are separated from the proteins, and then are docked back into the active sites using numerical optimization of this function. It is found that by employing this docking function, quasi-Newton optimization is capable of moving ligands great distances [on average 7 A root mean square distance (RMSD)] to locate the correctly docked structure. It is also found that a ligand drawn from one PDB file can be docked into a trypsin structure drawn from any of the trypsin PDB files. This implies that this scoring function is not limited to more accurate x-ray structures, as is the case for many of the conventional docking methods, but could be extended to homology models.