Glide's ability to identify active compounds in a database screen is characterized by applying Glide to a diverse set of nine protein receptors. In many cases, two, or even three, protein sites are employed to probe the sensitivity of the results to the site geometry. To make the database screens as realistic as possible, the screens use sets of “druglike” decoy ligands that have been selected to be representative of what we believe is likely to be found in the compound collection of a pharmaceutical or biotechnology company. Results are presented for releases 1.8, 2.0, and 2.5 of Glide. The comparisons show that average measures for both “early” and “global” enrichment for Glide 2.5 are 3 times higher than for Glide 1.8 and more than 2 times higher than for Glide 2.0 because of better results for the least well-handled screens. This improvement in enrichment stems largely from the better balance of the more widely parametrized GlideScore 2.5 function and the inclusion of terms that penalize ligand−protei...

/pdf/glide-a-new-approach-for-rapid-accurate-docking-and-scoring-5ewth20tdn.pdf

Glide: a new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening.

A novel scoring function to estimate protein-ligand binding affinities has been developed and implemented as the Glide 4.0 XP scoring function and docking protocol. In addition to unique water desolvation energy terms, protein-ligand structural motifs leading to enhanced binding affinity are included: (1) hydrophobic enclosure where groups of lipophilic ligand atoms are enclosed on opposite faces by lipophilic protein atoms, (2) neutral-neutral single or correlated hydrogen bonds in a hydrophobically enclosed environment, and (3) five categories of charged-charged hydrogen bonds. The XP scoring function and docking protocol have been developed to reproduce experimental binding affinities for a set of 198 complexes (RMSDs of 2.26 and 1.73 kcal/mol over all and well-docked ligands, respectively) and to yield quality enrichments for a set of fifteen screens of pharmaceutical importance. Enrichment results demonstrate the importance of the novel XP molecular recognition and water scoring in separating active and inactive ligands and avoiding false positives.

/pdf/extra-precision-glide-docking-and-scoring-incorporating-a-2o3e7pl9pg.pdf

Extra Precision Glide: Docking and Scoring Incorporating a Model of Hydrophobic Enclosure for Protein-Ligand Complexes

Structure-based virtual screening plays an important role in drug discovery and complements other screening approaches. In general, protein crystal structures are prepared prior to docking in order to add hydrogen atoms, optimize hydrogen bonds, remove atomic clashes, and perform other operations that are not part of the x-ray crystal structure refinement process. In addition, ligands must be prepared to create 3-dimensional geometries, assign proper bond orders, and generate accessible tautomer and ionization states prior to virtual screening. While the prerequisite for proper system preparation is generally accepted in the field, an extensive study of the preparation steps and their effect on virtual screening enrichments has not been performed. In this work, we systematically explore each of the steps involved in preparing a system for virtual screening. We first explore a large number of parameters using the Glide validation set of 36 crystal structures and 1,000 decoys. We then apply a subset of protocols to the DUD database. We show that database enrichment is improved with proper preparation and that neglecting certain steps of the preparation process produces a systematic degradation in enrichments, which can be large for some targets. We provide examples illustrating the structural changes introduced by the preparation that impact database enrichment. While the work presented here was performed with the Protein Preparation Wizard and Glide, the insights and guidance are expected to be generalizable to structure-based virtual screening with other docking methods.

Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments

Computational approaches that 'dock' small molecules into the structures of macromolecular targets and 'score' their potential complementarity to binding sites are widely used in hit identification and lead optimization Indeed, there are now a number of drugs whose development was heavily influenced by or based on structure-based design and screening strategies, such as HIV protease inhibitors Nevertheless, there remain significant challenges in the application of these approaches, in particular in relation to current scoring schemes Here, we review key concepts and specific features of small-molecule-protein docking methods, highlight selected applications and discuss recent advances that aim to address the acknowledged limitations of established approaches

/pdf/docking-and-scoring-in-virtual-screening-for-drug-discovery-4nrii464b0.pdf

Docking and scoring in virtual screening for drug discovery: methods and applications.

A new coronavirus, known as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is the aetiological agent responsible for the 2019–2020 viral pneumonia outbreak of coronavirus disease 2019 (COVID-19)1–4. Currently, there are no targeted therapeutic agents for the treatment of this disease, and effective treatment options remain very limited. Here we describe the results of a programme that aimed to rapidly discover lead compounds for clinical use, by combining structure-assisted drug design, virtual drug screening and high-throughput screening. This programme focused on identifying drug leads that target main protease (Mpro) of SARS-CoV-2: Mpro is a key enzyme of coronaviruses and has a pivotal role in mediating viral replication and transcription, making it an attractive drug target for SARS-CoV-25,6. We identified a mechanism-based inhibitor (N3) by computer-aided drug design, and then determined the crystal structure of Mpro of SARS-CoV-2 in complex with this compound. Through a combination of structure-based virtual and high-throughput screening, we assayed more than 10,000 compounds—including approved drugs, drug candidates in clinical trials and other pharmacologically active compounds—as inhibitors of Mpro. Six of these compounds inhibited Mpro, showing half-maximal inhibitory concentration values that ranged from 0.67 to 21.4 μM. One of these compounds (ebselen) also exhibited promising antiviral activity in cell-based assays. Our results demonstrate the efficacy of our screening strategy, which can lead to the rapid discovery of drug leads with clinical potential in response to new infectious diseases for which no specific drugs or vaccines are available. A programme of structure-assisted drug design and high-throughput screening identifies six compounds that inhibit the main protease of SARS-CoV-2, demonstrating the ability of this strategy to isolate drug leads with clinical potential.

/pdf/structure-of-mpro-from-sars-cov-2-and-discovery-of-its-2hyiph66s6.pdf

Structure of Mpro from SARS-CoV-2 and discovery of its inhibitors

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 ha

/pdf/glide-a-new-approach-for-rapid-accurate-docking-and-scoring-lnxrmd5zbm.pdf

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

A coarse grain model for phospholipids was systematically parametrized to mimic structural properties obtained from an atomistic simulation of a dimyristoylphosphatidylcholine bilayer. The model semiquantitatively reproduces the cross-sectional structure of a preassembled phospholipid bilayer obtained from an atomistic simulation; a property that was not directly fit. The model is sufficiently fast to permit the simulation of the self-assembly of the bilayer starting from a random configuration.

A Coarse Grain Model for Phospholipid Simulations

Major advances have been made at several levels of computer simulation of surfactant solutions. Atomistic level studies of preassembled surfactant structures have become fairly routine. The development of structure in surfactant solutions has now been studied using atomistic, coarse grain and mesoscopic models. Coarse grain and mesoscopic simulations have been used to determine phase diagrams. The challenges involved in treating complex surfactant solutions will continue to drive this field forward.

Computer simulation of surfactant solutions

A computationally efficient coarse grain model designed to closely mimic specific phospholipids is used to study a number of phospholipid systems to demonstrate its strengths and weaknesses. A study of a membrane containing an anesthetic, halothane, illustrates the shortcomings of this model in treating systems without extensive parametrization. In contrast, the power of the model is demonstrated by its ability to successfully simulate the self-assembly of two phospholipid phases from random initial configurations:  a lamellar phase and a reverse hexagonal phase in a ternary system containing water, a hydrocarbon, and a phospholipid. The aqueous columns in the reverse hexagonal phase tend to adopt polygonal cross sections and the local structure of phospholipids is still bilayer-like. Molecular dynamics was found to be much more efficient at simulating self-assembly in the current systems than Monte Carlo.

Simulations of Phospholipids Using a Coarse Grain Model

Glide is a ligand docking program for predicting protein-ligand binding modes and ranking ligands via high-throughput virtual screening. Glide utilizes two different scoring functions, SP and XP GlideScore, to rank-order compounds. Three modes of sampling ligand conformational and positional degrees of freedom are available to determine the optimal ligand orientation relative to a rigid protein receptor geometry. This unit presents protocols for flexible ligand docking with Glide, optionally including ligand constraints or ligand molecular similarities.

Mee Shelley

Papers

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

A Coarse Grain Model for Phospholipid Simulations

Computer simulation of surfactant solutions

Simulations of Phospholipids Using a Coarse Grain Model

Flexible ligand docking with Glide.