Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “virtual screening”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18Linked to original sources

On evaluating molecular-docking methods for pose prediction and enrichment factors.

Four of the most well-known, commercially available docking programs, FlexX, GOLD, GLIDE, and ICM, have been examined for their ligand-docking and virtual-screening capabilities. The relative performance of the programs in reproducing the native ligand conformation from starting SMILES strings for 164 high-resolution protein-ligand complexes is presented and compared. Applying only the native scoring functions, the latest versions of these four docking programs were also used to conduct virtual screening for 12 protein targets of therapeutic interest, involving both publicly available structures and AstraZeneca in-house structures. The capability of the four programs to correctly rank-order target-specific active compounds over alternative binders and nonbinders (decoys plus randomly selected compounds) and thereby enrich a small subset of a screening library is compared. Enrichments from the virtual-screening experiments are contrasted with those obtained with alternative 3D shape-matching and 2D similarity database-search methods.

Journal Article↗

A new method to estimate ligand-receptor energetics.

In the discovery of new drugs, lead identification and optimization have assumed critical importance given the number of drug targets generated from genetic, genomics, and proteomic technologies. High-throughput experimental screening assays have been complemented recently by "virtual screening" approaches to identify and filter potential ligands when the characteristics of a target receptor structure of interest are known. Virtual screening mandates a reliable procedure for automatic ranking of structurally distinct ligands in compound library databases. Computing a rank score requires the accurate prediction of binding affinities between these ligands and the target. Many current scoring strategies require information about the target three-dimensional structure. In this study, a new method to estimate the free binding energy between a ligand and receptor is proposed. We extend a central idea previously reported (Bock, J. R., and Gough, D. A. (2001) Predicting protein-protein interactions from primary structure. Bioinformatics 17, 455-460; Bock, J. R., and Gough, D. A. (2002) Whole-proteome interaction mining. Bioinformatics, in press) that uses simple descriptors to represent biomolecules as input examples to train a support vector machine (Smola, A. J., and Schölkopf, B. (1998) A Tutorial on Support Vector Regression, Neuro-COLT Technical Report NC-TR-98-030, Royal Holloway College, University of London, UK) and the application of the trained system to previously unseen pairs, estimating their propensity for interaction. Here we seek to learn the function that maps features of a receptor-ligand pair onto their equilibrium free binding energy. These features do not comprise any direct information about the three-dimensional structures of ligand or target. In cross-validation experiments, it is demonstrated that objective measurements of prediction error rate and rank-ordering statistics are competitive with those of several other investigations, most of which depend on three-dimensional structural data. The size of the sample (n = 2,671) indicates that this approach is robust and may have widespread applicability beyond restricted families of receptor types. It is concluded that newly sequenced proteins, or those for which three-dimensional crystal structures are not easily obtained, can be rapidly analyzed for their binding potential against a library of ligands using this methodology.

Drug Design↗

Methods for compound selection focused on hits and application in drug discovery.

In the context of virtual screening calculations, a multiple fingerprint-based metric is applied to generate focused compound libraries by database searching. Different fingerprints are used to facilitate a similarity step for database mining, followed by a diversity step to assemble the final library. The method is applied, for example, to build libraries of limited size for hit-to-lead development efforts. In studies designed to inhibit a therapeutically relevant protein-protein interaction, small molecular hits were initially obtained by combined fingerprint- and structure-based virtual screening and used for the design of focused libraries. We review the applied virtual screening approach and report the statistics and results of screening as well as focused library design. While the structures of lead compounds cannot be disclosed, the analysis is thought to provide an example of the interplay of different methods applied in practical lead identification.

Binding Sites↗

[Receptor-ligand docking simulation for membrane proteins].

Protein structure-based molecular design using the computational techniques of protein structure prediction, ligand docking, and virtual screening is an integral part of drug discovery for limiting the application of the structure-based approach to target proteins such as G-protein-coupled receptors (GPCRs). GPCRs play an important role in living organisms and are of major interest to the pharmaceutical industry. However, structural data on ligand-binding forms for GPCRs from experiments to elucidate structural templates for docking simulations are lacking due to the difficulties associated with crystallization and crystallography. Therefore structural prediction of GPCRs in the ligand-bound state using computational methods has been introduced, but the prediction of ligand conformation onto target GPCRs is still constructed manually by human experts. We developed a molecular modeling technique for the prediction of ligand-receptor binding using comparative ligand-binding analysis (CoLBA) that not only considers interaction energy but also the similarity of interaction profiles among ligands. The advantage of CoLBA is that it can facilitate intuitive and flexible screening based on docking results when protein structures with low resolution (or theoretical models) are targeted. We applied CoLBA to ligand-binding prediction in several GPCRs. The predicted ligand-binding models were evaluated in site-directed mutagenesis experiments in collaborative research, and the enrichment rate of activated ligands was compared with random compounds in virtual screening simulations. We propose that CoLBA can be applied in large-scale modeling of ligand-receptor complexes and virtual screening for GPCRs.

Amino Acid Sequence↗

Receptor flexibility in the in silico screening of reagents in the S1' pocket of human collagenase.

A major difficulty in structure-based molecular design is the prediction of the structure of the protein-ligand complex because of the enormous number of degrees of freedom. Commonly, the target protein is kept rigid in a single low-energy conformation. However, this does not reflect the dynamic nature of protein structures. In this work, we investigate the influence of receptor flexibility in virtual screening of reagents on a common scaffold in the S1' pocket of human collagenase (matrix metalloproteinase-1). We compare screening using a single-crystal structure and multiple NMR structures, both apo and holo forms. We also investigate two computational methods of addressing receptor flexibility that can be used when NMR data are not available. The results from virtual screening using the experimental structures are compared to those obtained using the two computational methods. From the results, we draw conclusions about the impact of target flexibility on the identification of active and diverse reagents in a virtual screening protocol.

Apoenzymes↗

BRUTUS: optimization of a grid-based similarity function for rigid-body molecular superposition. II. Description and characterization.

Finding novel lead molecules is one of the primary goals in early phases of drug discovery projects. However, structurally dissimilar compounds may exhibit similar biological activity, and finding new and structurally diverse lead compounds is difficult for computer algorithms. Molecular energy fields are appropriate for finding structurally novel molecules, but they are demanding to calculate and this limits their usefulness in virtual screening of large chemical databases. In our approach, energy fields are computed only once per superposition and a simple interpolation scheme is devised to allow coarse energy field lattices having fewer grid points to be used without any significant loss of accuracy. The resulting processing speed of about 0.25 s per conformation on a 2.4 GHz Intel Pentium processor allows the method to be used for virtual screening on commonly available desktop machines. Moreover, the results indicate that grid-based superposition methods could be efficiently used for the virtual screening of compound libraries.

Algorithms↗

Efficient molecular docking of NMR structures: application to HIV-1 protease.

Docking ligands into an ensemble of NMR conformers is essential to structure-based drug discovery if only NMR structures are available for the target. However, sequentially docking ligands into each NMR conformer through standard single-receptor-structure docking, referred to as sequential docking, is computationally expensive for large-scale database screening because of the large number of NMR conformers involved. Recently, we developed an efficient ensemble docking algorithm to consider protein structural variations in ligand binding. The algorithm simultaneously docks ligands into an ensemble of protein structures and achieves comparable performance to sequential docking without significant increase in computational time over single-structure docking. Here, we applied this algorithm to docking with NMR structures. The HIV-1 protease was used for validation in terms of docking accuracy and virtual screening. Ensemble docking of the NMR structures identified 91% of the known inhibitors under the criterion of RMSD < 2.0 A for the best-scored conformation, higher than the average success rate of single docking of individual crystal structures (66%). In the virtual screening test, on average, ensemble docking of the NMR structures obtained higher enrichments than single-structure docking of the crystal structures. In contrast, docking of either the NMR minimized average structure or a single NMR conformer performed less satisfactorily on both binding mode prediction and virtual screening, indicating that a single NMR structure may not be suitable for docking calculations. The success of ensemble docking of the NMR structures suggests an efficient alternative method for standard single docking of crystal structures and for considering protein flexibility.

Algorithms↗

Discovery of diverse thyroid hormone receptor antagonists by high-throughput docking.

Treatment of hyperthyroidism, a common clinical condition that can have serious manifestations in the elderly, has remained essentially unchanged for >30 years. Directly antagonizing the effect of the thyroid hormone at the receptor level may be a significant improvement for the treatment of hyperthyroid patients. We built a computer model of the thyroid hormone receptor (TR) ligand-binding domain in its predicted antagonist-bound conformation and used a virtual screening algorithm to select 100 TR antagonist candidates out of a library of >250,000 compounds. We were able to obtain 75 of the compounds selected in silico and studied their ability to act as antagonists by using cultured cells that express TR. Fourteen of these compounds were found to antagonize the effect of T3 on TR with IC50s ranging from 1.5 to 30 microM. A small virtual library of compounds, derived from the highest affinity antagonist (1-850) that could be rapidly synthesized, was generated. A second round of virtual screening identified new compounds with predicted increased antagonist activity. These second generation compounds were synthesized, and their ability to act as TR antagonists was confirmed by transfection and receptor binding experiments. The extreme structural diversity of the antagonist compounds shows how receptor-based virtual screening can identify diverse chemistries that comply with the structural rules of TR antagonism.

Drug Design↗

"In silico" design of new uranyl extractants based on phosphoryl-containing podands: QSPR studies, generation and screening of virtual combinatorial library, and experimental tests.

This paper is devoted to computer-aided design of new extractants of the uranyl cation involving three main steps: (i) a QSPR study, (ii) generation and screening of a virtual combinatorial library, and (iii) synthesis of several predicted compounds and their experimental extraction studies. First, we performed a QSPR modeling of the distribution coefficient (logD) of uranyl extracted by phosphoryl-containing podands from water to 1,2-dichloroethane. Two different approaches were used: one based on classical structural and physicochemical descriptors (implemented in the CODESSA PRO program) and another one based on fragment descriptors (implemented in the TRAIL program). Three statistically significant models obtained with TRAIL involve as descriptors either sequences of atoms and bonds or atoms with their close environment (augmented atoms). The best models of CODESSA PRO include its own molecular descriptors as well as fragment descriptors obtained with TRAIL. At the second step, a virtual combinatorial library of 2024 podands has been generated with the CombiLib program, followed by the assessment of logD values using developed QSPR models. At the third step, eight of these hypothetical compounds were synthesized and tested experimentally. Comparison with experiment shows that developed QSPR models successfully predict logD values for 7 of 8 compounds from that "blind test" set.

Journal Article↗

Representing receptor flexibility in ligand docking through relevant normal modes.

Inspired by the current representation of the ligand-receptor binding process, a normal-mode-based methodology is presented to incorporate receptor flexibility in ligand docking and virtual screening. However, the systematic representation of the deformation space grows geometrically with the number of modes, and furthermore, midscale loop rearrangements like those found in protein kinase binding pockets cannot be accounted for with the first lowest-frequency modes. We thus introduced a measure of relevance of normal modes on a given region of interest and showed that only very few modes in the low-frequency range are necessary and sufficient to describe loop flexibility in cAMP-dependent protein kinase. We used this approach to generate an ensemble of representative receptor backbone conformations by perturbing the structure along a combination of relevant modes. Each ensemble conformation is complexed with known non-native binders to optimize the position of the binding-pocket side chains through a full flexible docking procedure. The multiple receptor conformations thus obtained are used in a small-scale virtual screening using receptor ensemble docking. We evaluated this algorithm on holo and apo structures of cAMP-dependent protein kinase that exhibit backbone rearrangements on two independent loop regions close to the binding pocket. Docking accuracy is improved, since the ligands considered in the virtual screening docked within 1.5 A to at least one of the structures. The discrimination between binders and nonbinders is also enhanced, as shown by the improvement of the enrichment factor. This constitutes a new step toward the systematic integration of flexible ligand-flexible receptor docking tools in structure-based drug discovery.

Algorithms↗

From genome to drug lead: identification of a small-molecule inhibitor of the SARS virus.

Virtual screening, a fast, computational approach to identify drug leads [Perola, E.; Xu, K.; Kollmeyer, T. M.; Kaufmann, S. H.; Prendergast, F. G. J. Med. Chem.2000, 43, 401; Miller, M. A. Nat. Rev. Drug Disc.2002, 1 220], is limited by a known challenge in crystallographically determining flexible regions of proteins. This approach has not been able to identify active inhibitors of the severe acute respiratory syndrome-associated coronavirus (SARS-CoV) using solely the crystal structures of a SARS-CoV cysteine proteinase with a flexible loop in the active site [Yang, H. T.; Yang, M. J.; Ding, Y.; Liu, Y. W.; Lou, Z. Y. Proc. Natl. Acad. Sci. U.S.A.2003, 100, 13190; Jenwitheesuk, E.; Samudrala, R. Bioorg. Med. Chem. Lett.2003, 13, 3989; Rajnarayanan, R. V.; Dakshanamurthy, S.; Pattabiraman, N. Biochem. Biophys. Res. Commun.2004, 321, 370; Du, Q.; Wang, S.; Wei, D.; Sirois, S.; Chou, K. Anal. Biochem.2005, 337, 262; Du, Q.; Wang, S.; Zhu, Y.; Wei, D.; Guo, H. Peptides2004, 25, 1857; Lee, V.; Wittayanarakul, K.; Remsungenen, T.; Parasuk, V.; Sompornpisut, P. Science (Asia)2003, 29, 181; Toney, J.; Navas-Martin, S.; Weiss, S.; Koeller, A. J. Med. Chem.2004, 47, 1079; Zhang, X. W.; Yap, Y. L. Bioorg. Med. Chem.2004, 12, 2517]. This article demonstrates a genome-to-drug-lead approach that uses terascale computing to model flexible regions of proteins, thus permitting the utilization of genetic information to identify drug leads expeditiously. A small-molecule inhibitor of SARS-CoV, exhibiting an effective concentration (EC50) of 23 microM in cell-based assays, was identified through virtual screening against a computer-predicted model of the cysteine proteinase. Screening against two crystal structures of the same proteinase failed to identify the 23-microM inhibitor. This study suggests that terascale computing can complement crystallography, broaden the scope of virtual screening, and accelerate the development of therapeutics to treat emerging infectious diseases such as SARS and Bird Flu.

Aminobenzoates↗

A kinase-focused compound collection: compilation and screening strategy.

Lead identification by high-throughput screening of large compound libraries has been supplemented with virtual screening and focused compound libraries. To complement existing approaches for lead identification at Biogen Idec, a kinase-focused compound collection was designed, developed and validated. Two strategies were adopted to populate the compound collection: a ligand shape-based virtual screening and a receptor-based approach (structural interaction fingerprint). Compounds selected with the two approaches were cherry-picked from an existing high-throughput screening compound library, ordered from suppliers and supplemented with specific medicinal compounds from internal programs. Promising hits and leads have been generated from the kinase-focused compound collection against multiple kinase targets. The principle of the collection design and screening strategy was validated and the use of the kinase-focused compound collection for lead identification has been added to existing strategies.

Combinatorial Chemistry Techniques↗

Novel scoring methods in virtual ligand screening.

Several different approaches have been proposed in the last decade to assess the binding affinity of a virtual small molecule ligand to a target protein, particularly with respect to screening large compound databases. Here we review the methods that have been proposed, and discuss techniques for optimizing scoring functions that have been applied in industrial settings.

Computer Simulation↗

Drug screening for influenza neuraminidase inhibitors.

Neuraminidase (NA) is one of the most important targets to screen the drugs of anti-influenza virus A and B. After virtual screening approaches were applied to a compound database which possesses more than 10000 compound structures, 160 compounds were selected for bioactivity assay, then a High Throughput Screening (HTS) model established for influenza virus NA inhibitors was applied to detect these compounds. Finally, three compounds among them displayed higher inhibitory activities, the range of their IC50 was from 0.1 micromol/L to 3 micromol/L. Their structural scaffolds are novel and different from those of NA inhibitors approved for influenza treatment, and will be useful for the design and research of new NA inhibitors. The result indicated that the combination of virtual screening with HTS was very significant to drug screening and drug discovery.

Antiviral Agents↗

Methodological developments and strategies for a fast flexible superposition of drug-size molecules.

An alternative to experimental high through-put screening is the virtual screening of compound libraries on the computer. In absence of a detailed structure of the receptor protein, candidate molecules are compared with a known reference by mutually superimposing their skeletons and scoring their similarity. Since molecular shape highly depends on the adopted conformation, an efficient conformational screening is performed using a knowledge-based approach. A comprehensive torsion library has been compiled from crystal data stored in the Cambridge Structural Database. For molecular comparison a strategy is followed considering shape associated physicochemical properties in space such as steric occupancy, electrostatics, lipophilicity and potential hydrogen-bonding. Molecular shape is approximated by a set of Gaussian functions not necessarily located at the atomic positions. The superposition is performed in two steps: first by a global alignment search operating on multiple rigid conformations and then by conformationally relaxing the best scored hits of the global search. A normalized similarity scoring is used to allow for a comparison of molecules with rather different shape and size. The approach has been implemented on a cluster of parallel processors. As a case study, the search for ligands binding to the dopamine receptor is given.

Algorithms↗

[Virtual colonoscopy: a new hope for colorectal cancer screening].

Colon cancer is the second most common cause of cancer death and the first one in non-smoking persons. None of the current colorectal caner screening options is ideal. As for mammography in breast cancer screening, virtual colonoscopy, although still in clinical evaluation, promises to become a primary screening method as it is performant, quick, non-invasive and cheaper than conventional colonoscopy.

Colonic Polyps↗

Determination and mapping of activity-specific descriptor value ranges for the identification of active compounds.

MAD (Mapping to Activity class-specific Descriptor value ranges) is a novel molecular similarity method that relies on the identification of activity-specific descriptors. Applying a categorical descriptor scoring function, value ranges of molecular descriptors in screening databases are compared with those in classes of active compounds and descriptors displaying significant deviations are selected. In order to identify new actives, database molecules are mapped to class-specific value ranges and ranked using a similarity function. As a mapping algorithm, MAD is distinct from many other molecular similarity and virtual screening methods. In systematic virtual screening trials, for small selection sets of only 30 database compounds, average hit and recovery rates over six activity classes ranged from about 10% to 25% and about 25% to 75%, respectively. Moreover, when mining a database of bioactive molecules many similar compounds were selected (with hit rates between about 15% and 79%). Our findings suggest that it is possible to generate compound class-directed descriptor reference spaces for molecular similarity analysis.

Algorithms↗

Structure-based discovery of a new affinity ligand to pancreatic alpha-amylase.

A ligand useful for affinity capture of porcine pancreatic alpha-amylase was found by virtual screening of the commercially available compound data base MDL Available Chemicals Directory. Hits from the virtual screening were investigated for binding by nuclear magnetic resonance (NMR) and surface plasmon resonance. Selected compounds were tested for inhibition of the enzyme using a NMR-based assay. One of the binders found was covalently coupled to a chromatographic resin and a column, packed with this resin, could retain alpha-amylase, which subsequently was eluted by introduction of the known inhibitor acarbose to the elution buffer.

Animals↗