Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy.

Unlike other methods for docking ligands to the rigid 3D structure of a known protein receptor, Glide approximates a complete systematic search of the conformational, orientational, and positional space of the docked ligand. In this search, an initial rough positioning and scoring phase that dramatically narrows the search space is followed by torsionally flexible energy optimization on an OPLS-AA nonbonded potential grid for a few hundred surviving candidate poses. The very best candidates are further refined via a Monte Carlo sampling of pose conformation; in some cases, this is crucial to obtaining an accurate docked pose. Selection of the best docked pose uses a model energy function that combines empirical and force-field-based terms. Docking accuracy is assessed by redocking ligands from 282 cocrystallized PDB complexes starting from conformationally optimized ligand geometries that bear no memory of the correctly docked pose. Errors in geometry for the top-ranked pose are less than 1 A in nearly half of the cases and are greater than 2 A in only about one-third of them. Comparisons to published data on rms deviations show that Glide is nearly twice as accurate as GOLD and more than twice as accurate as FlexX for ligands having up to 20 rotatable bonds. Glide is also found to be more accurate than the recently described Surflex method.

Novel inhibitors of dengue virus methyltransferase: discovery by in vitro-driven virtual screening on a desktop computer grid.

Dengue fever is a viral disease that affects 50-100 million people annually and is one of the most important emerging infectious diseases in many areas of the world. Currently, neither specific drugs nor vaccines are available. Here, we report on the discovery of new inhibitors of the viral NS5 RNA methyltransferase, a promising flavivirus drug target. We have used a multistage molecular docking approach to screen a library of more than 5 million commercially available compounds against the two binding sites of this enzyme. In 263 compounds chosen for experimental verification, we found 10 inhibitors with IC(50) values of <100 microM, of which four exhibited IC(50) values of <10 microM in in vitro assays. The initial hit list also contained 25 nonspecific aggregators. We discuss why this likely occurred for this particular target. We also describe our attempts to use aggregation prediction to further guide the study, following this finding.

De novo structural prediction of transition metal complexes: application to technetium.

De novo structural prediction of transition metal complexes is investigated. Technetium complexes are chosen given their importance in medical imaging and nuclear waste remediation and for the chemical diversity they display. A new conformational searching algorithm (LIGB) for transition metals is described that allows one to search for different conformational and geometric isomers within a single simulation. In the preponderance of cases, both conformational searching techniques (LIGB and high-temperature molecular dynamics/simulated annealing) provide comparable results, while LIGB is superior for macrocyclic complexes. A genetic algorithm-optimized PM3(tm) parametrization for Tc is compared with the standard implementation and found to yield a significant improvement in predictive ability for the most prevalent Tc structural motifs. The utility of a coupled molecular mechanics-semiempirical quantum mechanics protocol is demonstrated for very rapid, efficient, and effective de novo prediction of transition metal complex geometries.