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 451 records · Page 25Linked to original sources
PubMed2006

Enzyme isoselective inhibitors: a tool for binding-trend analysis.

A transition-state analogue inhibitor that covalently reversibly binds to an enzyme formally consists of two parts: the chemical site, CS and the recognition site, RS. We have experimentally and theoretically demonstrated that the trend of binding affinity in a series of isoselective inhibitors (with identical RS and different CS fragments) depends mainly on their CS fragments. Isoselective inhibitors have the same affinity trend toward different enzymes of the same family with a common catalytic mechanism. Thus, very good correlation between experimentally determined and theoretically calculated Ki values was demonstrated. A practical outcome is the application of the described method as a tool for an expert analysis in virtual screening of inhibitor libraries and in the design of new enzyme inhibitors.

Binding Sites↗
PubMed2002

Pattern recognition and massively distributed computing.

A feature of Peter Kollman's research was his exploitation of the latest computational techniques to devise novel applications of the free energy perturbation method. He would certainly have seized upon the opportunities offered by massively distributed computing. Here we describe the use of over a million personal computers to perform virtual screening of 3.5 billion druglike molecules against protein targets by pharmacophore pattern matching, together with other applications of pattern recognition such as docking ligands without any a priori knowledge about the binding site location.

Algorithms↗
PubMed2006

HIV-1 integrase inhibitors: 2003-2004 update.

The integration of viral cDNA into the host genome is an essential step in the HIV-1-life cycle and is mediated by the virally encoded enzyme, integrase (IN). Inhibition of this process provides an attractive strategy for antiviral drug design. The discovery of beta-diketo acid inhibitors played a major role in validating IN as a legitimate antiretroviral drug target. Over a decade of research, a plethora of IN inhibitors have been discovered and some showed antiviral activity consistent with their effect on IN. To date, at least two compounds have been tested in human but none are close to the FDA approval. In this review, we provide a comprehensive report of all small-molecule IN inhibitors discovered during the years 2003 and 2004. Compilation of such data will prove beneficial in developing QSAR, virtual screening, pharmacophore hypothesis generation, and validation.

Anti-Retroviral Agents↗
PubMed2002

ConsDock: A new program for the consensus analysis of protein-ligand interactions.

Protein-based virtual screening of chemical libraries is a powerful technique for identifying new molecules that may interact with a macromolecular target of interest. Because of docking and scoring limitations, it is more difficult to apply as a lead optimization method because it requires that the docking/scoring tool is able to propose as few solutions as possible and all of them with a very good accuracy for both the protein-bound orientation and the conformation of the ligand. In the present study, we present a consensus docking approach (ConsDock) that takes advantage of three widely used docking tools (Dock, FlexX, and Gold). The consensus analysis of all possible poses generated by several docking tools is performed sequentially in four steps: (i) hierarchical clustering of all poses generated by a docking tool into families represented by a leader; (ii) definition of all consensus pairs from leaders generated by different docking programs; (iii) clustering of consensus pairs into classes, represented by a mean structure; and (iv) ranking the different means starting from the most populated class of consensus pairs. When applied to a test set of 100 protein-ligand complexes from the Protein Data Bank, ConsDock significantly outperforms single docking with respect to the docking accuracy of the top-ranked pose. In 60% of the cases investigated here, ConsDock was able to rank as top solution a pose within 2 A RMSD of the X-ray structure. It can be applied as a postprocessing filter to either single- or multiple-docking programs to prioritize three-dimensional guided lead optimization from the most likely docking solution.

Animals↗
PubMed2003

Pharmacophore-based molecular docking to account for ligand flexibility.

Rapid computational mining of large 3D molecular databases is central to generating new drug leads. Accurate virtual screening of large 3D molecular databases requires consideration of the conformational flexibility of the ligand molecules. Ligand flexibility can be included without prohibitively increasing the search time by docking ensembles of precomputed conformers from a conformationally expanded database. A pharmacophore-based docking method whereby conformers of the same or different molecules are overlaid by their largest 3D pharmacophore and simultaneously docked by partial matches to that pharmacophore is presented. The method is implemented in DOCK 4.0.

Algorithms↗
PubMed2003

Improved protein-ligand docking using GOLD.

The Chemscore function was implemented as a scoring function for the protein-ligand docking program GOLD, and its performance compared to the original Goldscore function and two consensus docking protocols, "Goldscore-CS" and "Chemscore-GS," in terms of docking accuracy, prediction of binding affinities, and speed. In the "Goldscore-CS" protocol, dockings produced with the Goldscore function are scored and ranked with the Chemscore function; in the "Chemscore-GS" protocol, dockings produced with the Chemscore function are scored and ranked with the Goldscore function. Comparisons were made for a "clean" set of 224 protein-ligand complexes, and for two subsets of this set, one for which the ligands are "drug-like," the other for which they are "fragment-like." For "drug-like" and "fragment-like" ligands, the docking accuracies obtained with Chemscore and Goldscore functions are similar. For larger ligands, Goldscore gives superior results. Docking with the Chemscore function is up to three times faster than docking with the Goldscore function. Both combined docking protocols give significant improvements in docking accuracy over the use of the Goldscore or Chemscore function alone. "Goldscore-CS" gives success rates of up to 81% (top-ranked GOLD solution within 2.0 A of the experimental binding mode) for the "clean list," but at the cost of long search times. For most virtual screening applications, "Chemscore-GS" seems optimal; search settings that give docking speeds of around 0.25-1.3 min/compound have success rates of about 78% for "drug-like" compounds and 85% for "fragment-like" compounds. In terms of producing binding energy estimates, the Goldscore function appears to perform better than the Chemscore function and the two consensus protocols, particularly for faster search settings. Even at docking speeds of around 1-2 min/compound, the Goldscore function predicts binding energies with a standard deviation of approximately 10.5 kJ/mol.

Algorithms↗
PubMed2004

Rapid protein-ligand docking using soft modes from molecular dynamics simulations to account for protein deformability: binding of FK506 to FKBP.

Most current docking methods to identify possible ligands and putative binding sites on a receptor molecule assume a rigid receptor structure to allow virtual screening of large ligand databases. However, binding of a ligand can lead to changes in the receptor protein conformation that are sterically necessary to accommodate a bound ligand. An approach is presented that allows relaxation of the protein conformation in precalculated soft flexible degrees of freedom during ligand-receptor docking. For the immunosuppressant FK506-binding protein FKBP, the soft flexible modes are extracted as principal components of motion from a molecular dynamics simulation. A simple penalty function for deformations in the soft flexible mode is used to limit receptor protein deformations during docking that avoids a costly recalculation of the receptor energy by summing over all receptor atom pairs at each step. Rigid docking of the FK506 ligand binding to an unbound FKBP conformation failed to identify a geometry close to experiment as favorable binding site. In contrast, inclusion of the flexible soft modes during systematic docking runs selected a binding geometry close to experiment as lowest energy conformation. This has been achieved at a modest increase of computational cost compared to rigid docking. The approach could provide a computationally efficient way to approximately account for receptor flexibility during docking of large numbers of putative ligands and putative docking geometries.

Binding Sites↗
PubMed2004

Probing flexibility and "induced-fit" phenomena in aldose reductase by comparative crystal structure analysis and molecular dynamics simulations.

Aldose reductase is a promising target for the treatment of diabetic complications, and as such, has become the focus of various drug design projects. As revealed by a survey of available crystal structures, the protein shows pronounced induced-fit effects upon ligand binding. Although helping to explain the enzyme's substrate promiscuity, phenomena of this kind are still responsible for significant complications in structure-based design efforts directed to aldose reductase. Accordingly, a deeper understanding of the principles governing conformational alterations in this enzyme would be of utmost practical importance. As a first step in addressing this issue, molecular dynamics (MD) simulations have been carried out. The ultrahigh resolution crystal structure of aldose reductase complexed with inhibitor IDD594 served as ideal starting point for a set of different simulations of nanosecond time scale: the native complexed state with bound inhibitor, the uncomplexed state (after removal of the inhibitor) at standard temperature, and the uncomplexed state at elevated temperature. The reference simulation of the complex exhibits extraordinary stability of the overall fold, whereas two distinct conformational substates are found for the binding-site region. In contrast, already at standard temperature pronounced changes are observed in the binding region during the simulation of the uncomplexed state. Leu300, for example, closes the access to the pocket opened by IDD594. On the other hand, conformations around the catalytic site are highly conserved, with the His110-Tyr48-NADP+ orientation being stabilized by a water molecule. Detailed analysis of the trajectories allows to reveal a set of distinct conformational substates that may prove useful as alternative structural templates in virtual screening for new aldose reductase inhibitors.

Acetates↗
PubMed2004

Toward the active conformations of rhodopsin and the beta2-adrenergic receptor.

Using sets of experimental distance restraints, which characterize active or inactive receptor conformations, and the X-ray crystal structure of the inactive form of bovine rhodopsin as a starting point, we have constructed models of both the active and inactive forms of rhodopsin and the beta2-adrenergic G-protein coupled receptors (GPCRs). The distance restraints were obtained from published data for site-directed crosslinking, engineered zinc binding, site-directed spin-labeling, IR spectroscopy, and cysteine accessibility studies conducted on class A GPCRs. Molecular dynamics simulations in the presence of either "active" or "inactive" restraints were used to generate two distinguishable receptor models. The process for generating the inactive and active models was validated by the hit rates, yields, and enrichment factors determined for the selection of antagonists in the inactive model and for the selection of agonists in the active model from a set of nonadrenergic GPCR drug-like ligands in a virtual screen using ligand docking software. The simulation results provide new insights into the relationships observed between selected biochemical data, the crystal structure of rhodopsin, and the structural rearrangements that occur during activation.

Animals↗
PubMed2005

General and targeted statistical potentials for protein-ligand interactions.

We present a novel atom-atom potential derived from a database of protein-ligand complexes. First, we clarify the similarities and differences between two statistical potentials described in the literature, PMF and Drugscore. We highlight shortcomings caused by an important factor unaccounted for in their reference states, and describe a new potential, which we name the Astex Statistical Potential (ASP). ASP's reference state considers the difference in exposure of protein atom types towards ligand binding sites. We show that this new potential predicts binding affinities with an accuracy similar to that of Goldscore and Chemscore. We investigate the influence of the choice of reference state by constructing two additional statistical potentials that differ from ASP only in this respect. The reference states in these two potentials are defined along the lines of Drugscore and PMF. In docking experiments, the potential using the new reference state proposed for ASP gives better success rates than when these literature reference states were used; a success rate similar to the established scoring functions Goldscore and Chemscore is achieved with ASP. This is the case both for a large, general validation set of protein-ligand structures and for small test sets of actives against four pharmaceutically relevant targets. Virtual screening experiments for these targets show less discrimination between the different reference states in terms of enrichment. In addition, we describe how statistical potentials can be used in the construction of targeted scoring functions. Examples are given for cdk2, using four different targeted scoring functions, biased towards increasingly large target-specific databases. Using these targeted scoring functions, docking success rates as well as enrichments are significantly better than for the general ASP scoring function. Results improve with the number of structures used in the construction of the target scoring functions, thus illustrating that these targeted ASP potentials can be continuously improved as new structural data become available.

Binding Sites↗
PubMed2006

Effective handling of induced-fit motion in flexible docking.

For structure-based drug design, where various ligand structures need to be docked to a target protein structure, a docking method that can handle conformational flexibility of not only the ligand, but also the protein, is indispensable. We have developed a simple and effective approach for dealing with the local induced-fit motion of the target protein, and implemented it in our docking tool, ADAM. Our approach efficiently combines the following two strategies: a vdW-offset grid in which the protein cavity is enlarged uniformly, and structure optimization allowing the motion of ligand and protein atoms. To examine the effectiveness of our approach, we performed docking validation studies, including redocking in 18 test cases and foreign-docking, in which various ligands from foreign crystal structures of complexes are docked into a target protein structure, in 22 cases (on five target proteins). With the original ADAM, the correct docking modes (RMSD < 2.0 A) were not present among the top 20 models in one case of redocking and four cases of foreign-docking. When the handling of induced-fit motion was implemented, the correct solutions were acquired in all 40 test cases. In foreign-docking on thymidine kinase, the correct docking modes were obtained as the top-ranked solutions for all 10 test ligands by our combinatorial approach, and this appears to be the best result ever reported with any docking tool. The results of docking validation have thus confirmed the effectiveness of our approach, which can provide reliable docking models even in the case of foreign-docking, where conformational change of the target protein cannot be ignored. We expect that this approach will contribute substantially to actual drug design, including virtual screening.

Ligands↗
PubMed2002

Efficient macrocyclization for cyclicpeptide using solid-phase reaction.

Cyclicpeptides are important targets in peptide synthesis because of their interesting biological properties. Constraining highly flexible linear peptides by cyclization is one of the mostly widely used approaches to define the bioactive conformation of peptides. Cyclic peptides often have increased receptor affinity and metabolic stability over their linear counterparts. We carried out virtual screening experiment via docking in order to understand the interaction between HLE-Human Leukocyte Elastase and ligand peptide and to identify the sequence that can be a target in various ligand peptides. We made cyclic peptides as a target base on Met-Ile-Phe sequence having affinity for ligand and receptor active site docking. There are three ways to cyclize certain sequences of amino acids such as Met-Ile-Phe-Gly-Ile. First is head-to-tail cyclization method, linking between N-terminal and C-terminal. Second method utilizes amino acid side chain such as thiol functional group in Cys, making a thioether bond. The last one includes an application of resin-substituted amino acids in solid phase reaction. Among the three methods, solid phase reaction showed the greatest yield. Macrocyclization of Fmoc-Met-Ile-Phe-Gly-Ile-OBn after cleavage of Fmoc protection in solution phase was carried out to give macrocyclic compound 5 in about 7% yield. In the contrast with solution phase reaction, solid phase reaction for macrocyclization of Met-Ile-Phe-Gly-Ile-Asp-Tentagel in normal concentrated condition gave macrocyclic compound 7 in more than 35% yield.

Cyclization↗
PubMed1999

Expanding our understanding of immunoglobulin, T-cell antigen receptor, and novel immune-type receptor genes: a subset of the immunoglobulin gene superfamily.

The immunoglobulin superfamily (IgSF) is an extensively diversified multigene family whose members share a common structural feature, the Ig fold. Members of the Ig/T-cell antigen receptor (TCR) subset of the IgSF mediate antigen-specific recognition in adaptive immune responses. Antigen-binding receptors belonging to this subset are present in all species of jawed vertebrates. To explore whether there are additional structurally related but otherwise distinct members of this subset, we have developed a technique termed the short-primer polymerase chain reaction (PCR) that targets structurally conserved short motifs in the Ig fold. Large-scale sequencing efforts and recent advances in information biotechnology, including "electronic PCR," provide additional computational means to implement similarly directed searches within databases. The use of these approaches has led to the discoveries of Ig/TCR homologues in a variety of phylogenetically diverse organisms, a diversified family of novel immune-type receptor genes, as well as a novel human IgSF member. The potential of random sequencing efforts and virtual screening of databases is described in the context of two novel genes in bony fish. The various methodologies that are discussed and the examples shown provide means for further investigating, and/or elucidating novel, IgSF receptors as well as components of pathways that are involved in immune responses in both traditional and nontraditional model systems.

Amino Acid Sequence↗
PubMed2004

Predicting anti-HIV activity of 2,3-diaryl-1,3-thiazolidin-4-ones: computational approach using reformed eccentric connectivity index.

The eccentric connectivity index, which has recently been employed successfully for the development of numerous mathematical models for the prediction of biological activities of diverse nature, has been reformed to overcome its limitations caused by degeneracy and insensitivity towards heteroatoms. The reformed eccentric connectivity index, termed the eccentric connectivity topochemical index, overcomes the limitations of the eccentric connectivity index by exhibiting very low degeneracy and displaying sensitivity to both the presence and relative position of heteroatoms without compromizing the discriminating power of the eccentric connectivity index. The relationship of the eccentric connectivity topochemical index, eccentric connectivity index and Wiener's index with regard to the anti-HIV activity of 2, 3-diaryl-1, 3-thiazolidin-4-one derivatives was subsequently investigated. The values of the eccentric connectivity topochemical index, the eccentric connectivity index and Wiener's index of each of 31 analogues comprizing the data set were computed using in-house computer program. Resultant data was analyzed and suitable models developed after identification of active ranges. Subsequently, each derivative was assigned a biological activity using these models, which was then compared with the reported anti-HIV activity. The accuracy of prediction using these models was found to vary from 81 to 90%. The proposed index offers a vast potential for virtual screening of combinatorial libraries, structure property/activity studies and drug design. [figure]. Basic structure of 2,3-diaryl-1, 3-thiazoidin-4-ones.

Anti-HIV Agents↗
PubMed2006

A fast surface-matching procedure for protein-ligand docking.

A very simple, fast, and efficient scheme is proposed for performing preliminary protein-ligand docking as the first step of intensive high-throughput virtual screening. The procedure acts as a surface-complementarity filter that first calculates the 2D-contour maps of both the protein cavity and of the ligands using a spherical harmonics description of the associated molecular surfaces. Next, the obtained 2D-fingerprint images are compared to detect their complementarity. This scheme was tested on three typical cases of protein cavities, namely, a well-closed pocket, a small open pocket, and a large open one. For that purpose, for each case, a sample of 101 ligand conformers was generated (the X-ray one and 100 different conformers generated using simulated annealing), and these conformational samples were ranked according to the complementarity with the protein cavity surface. Compared to traditional docking procedures such as FRED (considered as typical of a very fast rigid body docking algorithms) and GOLD (considered as typical of the more accurate flexible docking algorithms), our procedure was much faster and more successful in detecting the right X-ray conformation. We did, however, identify a certain weakness in the case of the very large pocket where results were not as expected. In general, our method could be used for incorporating indirectly flexibility in protein-ligand docking calculations as such a scheme can easily handle several conformational states of both the protein and the ligand.

Binding Sites↗
PubMed2004

Unsupervised guided docking of covalently bound ligands.

An approach for docking covalently bound ligands in protein enzymes or receptors was implemented in MacDOCK, a similarity-driven docking program based on DOCK 4.0. This approach was tested with a small number of covalent ligand-protein structures, using both native and non-native protein structures. In all cases, MacDOCK was able to generate orientations consistent with the known covalent binding mode of these complexes, with a performance similar to that of other docking programs. This method was also applied to search for known covalent thrombin inhibitors in a medium-sized molecular database (ca. 11,000 compounds). Detection of functional groups suitable for covalent docking was carried out automatically. A significant enrichment in known active molecules in the first 5% of the database was obtained, showing that MacDOCK can be used efficiently for the virtual screening of covalently bound ligands.

Computer Simulation↗
PubMed2005

Characterization of beta3-adrenergic receptor: determination of pharmacophore and 3D QSAR model for beta3 adrenergic receptor agonism.

The beta3-adrenoreceptor (beta3-AR) has been shown to mediate various pharmacological and physiological effects such as lipolysis, thermogenesis, and intestinal smooth muscle relaxation. It also plays an important role in glucose homeostasis and energy balance. Molecular modeling studies were undertaken to develop predictive pharmacophoric hypothesis and 3D-QSAR model, which may explain variations in beta3-AR agonistic activity in terms of chemical features and physicochemical properties. The two softwares, CATALYST for pharmacophoric alignment and APEX-3D for 3D-QSAR modeling were used to establish the structure activity relationships for beta3-AR agonistic activity. Among the several statistically significant models, the selection of the best pharmacophore and 3D-QSAR model was based on its ability to estimate the activity of external test sets of similar and different structural types along with the reasonable consistency of the model with the limited information of the active site of beta3-AR. The final 3D-QSAR model was derived using the pharmacophoric alignments from the hypothesis which consisted of four chemical features: basic or positive ionizable feature on the nitrogen of the aryloxypropylamino group, two ring aromatic features corresponding to the phenyl ring of the phenoxide and the benzenesulphonamido groups and a hydrogen-bond donor (HBD) in the vicinity of the nitrogen atom of the benzenesulphonamido group with the most active molecule mapping in an energetically favorable extended conformation. This hypothesis was in agreement with the site directed mutagenesis studies on human beta3-AR and correlated well the observed and estimated activity both in, training and both the external test sets. It also mapped reasonably well to six beta3-AR agonists of different structural classes under clinical development and thus this hypothesis may have a universal applicability in providing a powerful template for virtual screening and also for designing new chemical entities (NCEs) as beta3-AR agonists.

Adrenergic beta-3 Receptor Agonists↗
PubMed2006

Permuting input for more effective sampling of 3D conformer space.

SMILES strings and other classic 2D structural formats offer a convenient way to represent molecules as a simplistic connection table, with the inherent advantages of ease of handling and storage. In the context of virtual screening, chemical databases to be screened are often initially represented by canonicalised SMILES strings that can be filtered and pre-processed in a number of ways, resulting in molecules that occupy similar regions of chemical space to active compounds of a therapeutic target. A wide variety of software exists to convert molecules into SMILES format, namely, Mol2smi (Daylight Inc.), MOE (Chemical Computing Group) and Babel (Openeye Scientific Software). Depending on the algorithm employed, the atoms of a SMILES string defining a molecule can be ordered differently. Upon conversion to 3D coordinates they result in the production of ostensibly the same molecule. In this work we show how different permutations of a SMILES string can affect conformer generation, affecting reliability and repeatability of the results. Furthermore, we propose a novel procedure for the generation of conformers, taking advantage of the permutation of the input strings--both SMILES and other 2D formats, leading to more effective sampling of conformation space in output, and also implementing fingerprint and principal component analyses step to post process and visualise the results.

Algorithms↗