SEARCH · PubMed Health

Results for “Molecular Docking Simulation”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

At least 1,459 records · Page 81Linked to original sources

PubMed2005

Human topoisomerase I poisoning: docking protoberberines into a structure-based binding site model.

Using the X-ray crystal structure of the human topoisomerase I (top1) - DNA cleavable complex and the Sybyl software package, we have developed a general model for the ternary cleavable complex formed with four protoberberine alkaloids differing in the substitution on the terminal phenyl rings and covering a broad range of the top1-poisoning activities. This model has the drug intercalated with its planar chromophore between the -1 and +1 base pairs flanking the cleavage site, with the nonplanar portion pointing into the minor groove. The ternary complexes were geometry-optimized and relative interaction energies, computed by using the Tripos force field, were found to rank in correct order the biological potency of the compounds; in addition, the model is also consistent with the top1-poisoning inactivity of berberine, a major prototype of the protoberberine alkaloids. The model might serve as a rational basis for elaboration of the most active compound as a lead structure, in order to develop more potent top1 poisons as next generation anti-cancer drugs.

Berberine Alkaloids↗

PubMed2005

Computational study of the dynamics of mannose disaccharides free in solution and bound to the potent anti-HIV virucidal protein cyanovirin.

In this paper, we present a computational study of the dynamics of the potent anti-HIV virucidal protein cyanovirin in complex with mannose disaccharides. Recently, it has been experimentally demonstrated that cyanovirin binds mannose oligomers on the surface of glycoprotein gp120. gp120, a protein on the surface of the HIV virus, is key in the process of viral docking and transfer of genetic material into human cells. Cyanovirin prevents the transfer of viral RNA into human cells. In this study, we found that, among all residues that show nuclear Overhauser effects in the solution NMR experiments, residues Glu41 and Arg76 appear to interact with the sugar at the high-affinity binding site through stronger Coulombic interactions. In particular, Arg76 participates in a dynamical mechanism that caps and locks the sugar once it is bound to the protein. We also studied the distribution of glycosidic torsional angles of mannose disaccharides in solution and compared it with those when bound at the high- and low-affinity sites of the protein. Throughout our 20 ns simulations, we find that the sugar bound to the high-affinity site preserves the most favorable conformation in solution while the sugar bound at the low-affinity site does not. The sugar at the low-affinity site can adopt both conformations, but we find it most predominantly on the one that is least probable for the free sugar in solution. We also carried out a detailed study of the interactions between the disaccharides and different amino acids as well as between the disaccharide and the solvent at both binding locations.

Anti-HIV Agents↗

PubMed2003

Docking of protein-protein complexes on the basis of highly ambiguous intermolecular distance restraints derived from 1H/15N chemical shift mapping and backbone 15N-1H residual dipolar couplings using conjoined rigid body/torsion angle dynamics.

A simple and reliable method for docking protein-protein complexes using (1)H(N)/(15)N chemical shift mapping and backbone (15)N-(1)H residual dipolar couplings is presented and illustrated with three complexes (EIN-HPr, IIA(Glc)-HPr, and IIA(Mtl)-HPr) of known structure. The (1)H(N)/(15)N chemical shift mapping data are transformed into a set of highly ambiguous, intermolecular distance restraints (comprising between 400 and 3000 individual distances) with translational and some degree of orientational information content, while the dipolar couplings provide information on relative protein-protein orientation. The optimization protocol employs conjoined rigid body/torsion angle dynamics in simulated annealing calculations. The target function also comprises three nonbonded interactions terms: a van der Waals repulsion term to prevent atomic overlap, a radius of gyration term (E(rgyr)) to avoid expansion at the protein-protein interface, and a torsion angle database potential of mean force to bias interfacial side chain conformations toward physically allowed rotamers. For the EIN-HPr and IIA(Glc)-HPr complexes, all structures satisfying the experimental restraints (i.e., both the ambiguous intermolecular distance restraints and the dipolar couplings) converge to a single cluster with mean backbone coordinate accuracies of 0.7-1.5 A. For the IIA(Mtl)-HPr complex, twofold degeneracy remains, and the structures cluster into two distinct solutions differing by a 180 degrees rotation about the z axis of the alignment tensor. The correct and incorrect solutions which have mean backbone coordinate accuracies of approximately 0.5 and approximately 10.5 A, respectively, can readily be distinguished using a variety of criteria: (a) examination of the overall (1)H(N)/(15)N chemical shift perturbation map (because the incorrect cluster predicts the presence of residues at the interface that experience only minimal chemical shift perturbations; this information is readily incorporated into the calculations in the form of ambiguous intermolecular repulsion restraints); (b) back-calculation of dipolar couplings on the basis of molecular shape; or (c) the E(rgyr) distribution which, because of its global nature, directly reflects the interfacial packing quality. This methodology should be particularly useful for high throughput, NMR-based, structural proteomics.

Bacterial Proteins↗

PubMed2006

Predicting the three-dimensional structure of human P-glycoprotein in absence of ATP by computational techniques embodying crosslinking data: insight into the mechanism of ligand migration and binding sites.

P-glycoprotein is a membrane protein involved in the phenomenon of multidrug resistance. Its activity and transport function have been largely characterized by various biochemical studies and a low-resolution image has been obtained by electron microscopy. Obtaining a high-resolution structure is, however, still remote due to the inherent difficulties in the experimental determination of membrane protein structures. We present here a three-dimensional (3D) atomic model of P-glycoprotein in absence of ATP. This model was obtained using a combination of computational techniques including comparative modeling and rigid body dynamics simulations that embody all available cysteine disulfide crosslinking data characterizing the whole protein in absence of ATP. The model features rather well most of the experimental interresidue distances derived both in the transmembrane domains and in the nucleotide binding domains. The model is also in good agreement with electron microscopy data, particularly in terms of size and topology. It features a large cavity detected in the protein core into which seven ligands were successfully docked. Their predicted affinity correlates well with experimental values. Locations of docked ligands compare favorably with those suggested by cysteine-scanning data. The finding of different positions both for a single ligand and for different ligands corroborates the experimental evidence indicating the existence of multiple drug binding sites. The interactions identified between P-glycoprotein and the docked ligands reveal that different types of interactions such as H-bonds, pi-pi and cation-pi interactions occur in agreement with a recently proposed pharmacophore model of P-glycoprotein ligands. Furthermore, the model also displays a lateral opening located in the transmembrane domains connecting the lipid bilayer to the central cavity. This feature supports rather well the commonly admitted mechanism of substrate uptake from the lipid bilayer. We propose that this 3D model may be an important tool to understand the structure-function relationship of P-glycoprotein.

ATP Binding Cassette Transporter, Subfamily B, Mem↗

PubMed2005

Phospho-pivot modeling predicts specific interactions of protein phosphatase-1 with a phospho-inhibitor protein CPI-17.

Phospho-amino acids in proteins are directly associated with phospho-receptor proteins, including protein phosphatases. Here we produced and tested a scheme for docking together interacting phospho-proteins whose monomeric 3D structures were known. The phosphate of calyculin A, an inhibitor for protein phosphatase-1 and 2A (PP1 and PP2A), or phospho-CPI-17, a PP1-specific inhibitor protein, was docked at the active site of PP1. First, a library of 186,624 virtual complexes was generated in silico, by pivoting the phospho-ligand at the phosphorus atom by step every 5 degrees on three rotational axes. These models were then graded for probability according to atomic proximity between two molecules. The predicted structure of PP1 x calyculin A complex fitted to the crystal structure with r.m.s.d. of 0.23 A, providing a validate test of the modeling method. Modeling of PP1 x phospho-CPI-17 complex yielded one converged structure. The segment of CPI-17 around phospho-Thr38 is predicted to fit in the active site of PP1. Positive charges at Arg33/36 of CPI-17 are in close proximity to Glu274 of PP1, where the sequence is unique among Ser/Thr phosphatases. Single mutations of these residues in PP1 reduced the affinity against phospho-CPI-17. Thus, the interface of the PP1 x CPI-17 complex predicted by the phospho-pivot modeling accounts for the specificity of CPI-17 against PP1.

Animals↗

PubMed2003

Complete reversal of coenzyme specificity by concerted mutation of three consecutive residues in alcohol dehydrogenase.

Gastric tissues from amphibian Rana perezi express the only vertebrate alcohol dehydrogenase (ADH8) that is specific for NADP(H) instead of NAD(H). In the crystallographic ADH8-NADP+ complex, a binding pocket for the extra phosphate group of coenzyme is formed by ADH8-specific residues Gly223-Thr224-His225, and the highly conserved Leu200 and Lys228. To investigate the minimal structural determinants for coenzyme specificity, several ADH8 mutants involving residues 223 to 225 were engineered and kinetically characterized. Computer-assisted modeling of the docked coenzymes was also performed with the mutant enzymes and compared with the wild-type crystallographic binary complex. The G223D mutant, having a negative charge in the phosphate-binding site, still preferred NADP(H) over NAD(H), as did the T224I and H225N mutants. Catalytic efficiency with NADP(H) dropped dramatically in the double mutants, G223D/T224I and T224I/H225N, and in the triple mutant, G223D/T224I/H225N (kcat/KmNADPH = 760 mm-1 min-1), as compared with the wild-type enzyme (kcat/KmNADPH = 133330 mm-1 min-1). This was associated with a lower binding affinity for NADP+ and a change in the rate-limiting step. Conversely, in the triple mutant, catalytic efficiency with NAD(H) increased, reaching values (kcat/KmNADH = 155000 mm-1 min-1) similar to those of the wild-type enzyme with NADP(H). The complete reversal of ADH8 coenzyme specificity was therefore attained by the substitution of only three consecutive residues in the phosphate-binding site, an unprecedented achievement within the ADH family.

Alcohol Dehydrogenase↗

PubMed2004

Evaluation and application of multiple scoring functions for a virtual screening experiment.

In order to identify novel chemical classes of factor Xa inhibitors, five scoring functions (FlexX, DOCK, GOLD, ChemScore and PMF) were engaged to evaluate the multiple docking poses generated by FlexX. The compound collection was composed of confirmed potent factor Xa inhibitors and a subset of the LeadQuest screening compound library. Except for PMF the other four scoring functions succeeded in reproducing the crystal complex (PDB code: 1FAX). During virtual screening the highest hit rate (80%) was demonstrated by FlexX at an energy cutoff of -40 kJ/mol, which is about 40-fold over random screening (2.06%). Limited results suggest that presenting more poses of a single molecule to the scoring functions could deteriorate their enrichment factors. A series of promising scaffolds with favorable binding scores was retrieved from LeadQuest. Consensus scoring by pair-wise intersection failed to enrich the hit rate yielded by single scorings (i.e. FlexX). We note that reported successes of consensus scoring in hit rate enrichment could be artificial because their comparisons were based on a selected subset of single scoring and a markedly reduced subset of double or triple scoring. The findings presented in this report are based upon a single biological system and support further studies.

Animals↗

PubMed2002

Concept and prototype of protein-ligand docking simulator with force feedback technology.

A novel concept for a protein-ligand docking simulator using Virtual Reality (VR) technologies, in particular the tactile sense technology, was designed and a prototype was developed. Most conventional docking simulators are based on numerical differential calculations of the total energy between a protein and a ligand. However, the basic concept of our method differs from that of conventional simulators. Our design utilizes the force between a ligand and a protein instead of the total energy. The most characteristic function of the system is its ability to enable the user to 'touch' and sense the electrostatic potential field of a protein molecule. The user can scan the surface of a protein using a globular probe, which is given an electrostatic charge, and is controlled by a force feedback device. The electrostatic force between the protein and the probe is calculated in real time and immediately fed back into the force feedback device. The user can easily search interactively for positions where the probe is strongly attracted to the force field. Such positions can be regarded as candidate sites where functional groups of ligands corresponding to the probe can bind to the target protein. Certain limitations remain; for example, only twenty protein atoms can be used to generate the electrostatic field. Furthermore, the system can only use globular probes, preventing drug molecules or small chemical groups from being simulated. These limitations are the result of our insufficient computer resources. However, our prototype system has the potential to become a novel application method as well as being applicable to conventional VR technologies, especially to force feedback technologies.

Computational Biology↗

PubMed2003

Fold recognition, homology modeling, docking simulations, kinetics analysis and mutagenesis of ATP/CTP:tRNA nucleotidyltransferase from Methanococcus jannaschii.

ATP/CTP:tRNA nucleotidyltransferases (NTases) and poly(A) polymerases (PAPs) belong to the same superfamily and their catalytic domains are remotely related. Based on the results of fold-recognition analysis and comparison of secondary structure patterns, we predicted that these two NTase families share three domains, corresponding to "palm," "fingers," and "fingernails" in the PAP crystal structure. A homology model of tRNA NTase from Methanococcus jannaschii was constructed. Energy minimization calculations of enzyme-nucleotide complexes and computer-aided docking of nucleotides onto the enzyme's surface were carried out to explore possible ATP and CTP binding sites. Theoretical models were used to guide experimental analysis. Recombinant His-tagged enzyme was expressed in Escherichia coli, and kinetic properties were characterized. The apparent K(M) for CTP was determined to be 38 microM, and the apparent K(M) for ATP was 21 microM. Three mutations of basic amino acids to alanine were created in a highly conserved region predicted to be in the vicinity of the nucleotide binding site. A deletion was also constructed to remove the C-terminal structural domain defined by the model; it retained about 1% of wild type enzymatic activity using CTP as co-substrate, confirming that detectable catalytic activity is exhibited by the N-terminal domain, as defined by the model. Our results suggest a mechanism of differential ATP and CTP binding, which explains how the tRNA NTase, having only one catalytic site, utilizes different nucleotide triphosphates depending on the nature of the tRNA substrate.

Adenosine Triphosphate↗

PubMed2005

Protein side-chain packing problem: a maximum edge-weight clique algorithmic approach.

"Protein Side-chain Packing" has an ever-increasing application in the field of bio-informatics, dating from the early methods of homology modeling to protein design and to the protein docking. However, this problem is computationally known to be NP-hard. In this regard, we have developed a novel approach to solve this problem using the notion of a maximum edge-weight clique. Our approach is based on efficient reduction of protein side-chain packing problem to a graph and then solving the reduced graph to find the maximum clique by applying an efficient clique finding algorithm developed by our co-authors. Since our approach is based on deterministic algorithms in contrast to the various existing algorithms based on heuristic approaches, our algorithm guarantees of finding an optimal solution. We have tested this approach to predict the side-chain conformations of a set of proteins and have compared the results with other existing methods. We have found that our results are favorably comparable or better than the results produced by the existing methods. As our test set contains a protein of 494 residues, we have obtained considerable improvement in terms of size of the proteins and in terms of the efficiency and the accuracy of prediction.

Algorithms↗

PubMed2001

Modeling of the outer vestibule and selectivity filter of the L-type Ca2+ channel.

Using the KcsA bacterial K+ channel crystal structure [Doyle, D. A., et al. (1998) Science 280, 69-74] and the model of the outer vestibule of the Na+ channel [Lipkind, G. M., and Fozzard, H. A. (2000) Biochemistry 39, 8161-8170] as structural templates, we propose a structural model of the outer vestibule and selectivity filter of the pore of the Ca2+ channel (alpha1C or Ca(v)1.2). The Ca2+ channel P loops were modeled by alpha-helix-turn-beta-strand motifs, with the glutamate residues of the EEEE motif located in the turns. P loops were docked in the extracellular part of the inverted teepee structure formed by S5 and S6 alpha-helices with backbone coordinates from the M1 and M2 helices of the KcsA crystal structure. This construction results in a conical outer vestibule that tapers to the selectivity filter at the bottom. The modeled selectivity ring forms a wide open pore ( approximately 6 A) in the absence of Ca2+. When Ca2+ is present ( approximately 1 microM), all four glutamate side chains move to the center and form a cage around the dehydrated Ca2+ ion, blocking the pore. In the millimolar concentration range, Ca2+ also interacts with two low-affinity sites located externally and internally, which were modeled by the same carboxylate groups of the selectivity filter. Calculation of the resulting electrostatic potentials show that the single Ca2+ ion is located in an electrostatic trap. Only when three Ca2+ ions are bound simultaneously in the high- and low-affinity sites of the selectivity filter is Ca2+ able to overcome electrostatic attraction, permitting Ca2+ flux.

Amino Acid Motifs↗

PubMed1994

Structural modeling and electrostatic properties of aspartate transcarbamylase from Saccharomyces cerevisiae.

In Saccharomyces cerevisiae the first two reactions of the pyrimidine pathway are catalyzed by a multifunctional protein which possesses carbamylphosphate synthetase and aspartate transcarbamylase activities. Genetic and proteolysis studies suggested that the ATCase activity is carried out by an independently folded domain. In order to provide structural information for ongoing mutagenesis studies, a model of the three-dimensional structure of this domain was generated on the basis of the known X-ray structure of the related catalytic subunit from E. coli ATCase. First, a model of the catalytic monomer was built and refined by energy minimization. In this structure, the conserved residues between the two proteins were found to constitute the hydrophobic core whereas almost all the mutated residues are located at the surface. Then, a trimeric structure was generated in order to build the active site as it lies at the interface between adjacent chains in the E. coli catalytic trimer. After docking a bisubstrate analog into the active site, the whole structure was energy minimized to regularize the interactions at the contact areas between subunits. The resulting model is very similar to that obtained for the E. coli catalytic trimer by X-ray crystallography, with a remarkable conservation of the structure of the active site and its vicinity. Most of the interdomain and intersubunit interactions that are essential for the stability of the E. coli catalytic trimer are maintained in the yeast enzyme even though there is only 42% identity between the two sequences. Free energy calculations indicate that the trimeric assembly is more stable than the monomeric form. Moreover an insertion of four amino acids is localized in a loop which, in E. coli ATCase, is at the surface of the protein. This insertion exposes hydrophobic residues to the solvent. Interestingly, such an insertion is present in all the eukaryotic ATCase genes sequences so far, suggesting that this region is interacting with another domain of the multifunctional protein.

Amino Acid Sequence↗

PubMed2005

Identification of molecular determinants that modulate trafficking of DeltaF508 CFTR, the mutant ABC transporter associated with cystic fibrosis.

Cystic fibrosis is a life-shortening inherited disorder associated primarily with a three-base in frame deletion that eliminates Phe508 in the ABC transporter, cystic fibrosis transmembrane conductance regulator (CFTR). Mutant CFTR, designated deltaF508 CFTR, is misprocessed and retained intracellularly. It is unclear what causes the trafficking impairment despite extensive investigative effort and the disease's prevalence. We hypothesize that the trafficking impairment is mediated by "receptors" of the cellular trafficking machinery that at three sequential "trafficking checkpoints" govern (1) exit from the endoplasmic reticulum (ER), (2) Golgi to the ER retrieval, and (3) targeting from post-Golgi compartments to lysosomes. We propose that, because of the Phe508 deletion and polypeptide misfolding: (1) a forward-directing signal recognized by the sec24 component of the COPII complex that mediates ER exit is eliminated; (2) a basic amino acid signal recognized by the COPI machinery involved in Golgi to ER retrieval becomes activated; and (3) a tyrosine-based sorting signal that targets to the lysosomes likewise becomes activated. We employed recently reported crystal structures of CFTR nucleotide binding domain 1 and sec24 in computational docking models to identify the most plausible CFTR-sec24 recognition domain. Site-directed mutagenesis and heterologous expression were also used to identify amino acid sequences that operate in Golgi to ER and post-Golgi to lysosome targeting. The importance of considering a multiple checkpoint model for trafficking is that rationale design of pharmaceutical interventions would require abrogation of all major checkpoints to deliver deltaF508 CFTR to the cell surface.

ATP-Binding Cassette Transporters↗

PubMed2003

A core-weighted fitting method for docking atomic structures into low-resolution maps: application to cryo-electron microscopy.

Cryo-electron microscopy of "single particles" is a powerful method to analyze structures of large macromolecular assemblies that are not amenable to investigation by traditional X-ray crystallographic methods. A key step in these studies is to obtain atomic interpretations of multiprotein complexes by fitting atomic structures of individual components into maps obtained from electron microscopic data. Here, we report the use of a "core-weighting" method, combined with a grid-threading Monte Carlo (GTMC) approach for this purpose. The "core" of an individual structure is defined to represent the part where the density distribution is least likely to be altered by other components that comprise the macromolecular assembly of interest. The performance of the method has been evaluated by its ability to determine the correct fit of (i) the alpha-chain of the T-cell receptor variable domain into a simulated map of the alphabeta complex at resolutions between 5 and 40 A, and (ii) the E2 catalytic domain of the pyruvate dehydrogenase into an experimentally determined map, at 14 A resolution, of the icosahedral complex formed by 60 copies of this enzyme. Using the X-ray structures of the two test cases as references, we demonstrate that, in contrast to more traditional methods, the combination of the core-weighting method and the grid-threading Monte Carlo approach can identify the correct fit reliably and rapidly from the low-resolution maps that are typical of structures determined with the use of single-particle electron microscopy.

Animals↗

PubMed1995

The three-dimensional solution structure of human stefin A.

The three-dimensional solution structure of recombinant human stefin A has been determined by a simulated annealing protocol using a total of 1113 distance and angle constraints obtained from 1H and 15N HMR spectroscopy. The solution structure is represented by a family of 17 conformers with an average root-mean-square deviation relative to the mean structure of 0.44 A for backbone atoms and 0.94 A for all heavy atoms for the main body of the structure. The protein has a well-defined global fold consisting of five anti-parallel beta-strands wrapped around a central five-turn alpha-helix. There is considerable similarity between the structural features of free stefin A in solution and the X-ray structure of the homologous protein stefin B in its complex with papain, but there are also some important differences in the regions which are fundamental to proteinase binding. The differences consist primarily of two regions of high conformational heterogeneity in free stefin A which correspond in stefin B to two of the components of the tripartite wedge that docks into the active site of the target proteinase. These regions, which are shown to be mobile in solution, are the five N-terminal residues and the second binding loop. In the bound conformation of stefin B they form a turn and a short helix, respectively.

Amino Acid Sequence↗

PubMed2000

Identification of ligands for RNA targets via structure-based virtual screening: HIV-1 TAR.

Binding of the Tat protein to TAR RNA is necessary for viral replication of HIV-1. We screened the Available Chemicals Directory (ACD) to identify ligands to bind to a TAR RNA structure using a four-step docking procedure: rigid docking first, followed by three steps of flexible docking using a pseudobrownian Monte Carlo minimization in torsion angle space with progressively more detailed conformational sampling on a progressively smaller list of top-ranking compounds. To validate the procedure, we successfully docked ligands for five RNA complexes of known structure. For ranking ligands according to binding avidity, an empirical binding free energy function was developed which accounts, in particular, for solvation, isomerization free energy, and changes in conformational entropy. System-specific parameters for the function were derived on a training set of RNA/ligand complexes with known structure and affinity. To validate the free energy function, we screened the entire ACD for ligands for an RNA aptamer which binds L-arginine tightly. The native ligand ranked 17 out of ca. 153,000 compounds screened, i.e., the procedure is able to filter out >99.98% of the database and still retain the native ligand. Screening of the ACD for TAR ligands yielded a high rank for all known TAR ligands contained in the ACD and suggested several other potential TAR ligands. Eight of the highest ranking compounds not previously known to be ligands were assayed for inhibition of the Tat-TAR interaction, and two exhibited a CD50 of ca. 1 microM.

Arginine↗

PubMed2000

Modeling and functional analysis of the interaction between von Willebrand factor A1 domain and glycoprotein Ibalpha.

Binding of the von Willebrand factor (vWF) A1 domain to the glycoprotein (GP) Ib-IX-V complex mediates platelet adhesion to reactive substrates under high shear stress conditions, a key event in hemostasis and thrombosis. We have now used the known three-dimensional structure of the A1 domain to model the interaction with the GP Ibalpha sequence 271-279, which has previously been implicated in ligand binding. Docking procedures suggested that A1 domain residues in strand beta3 and preceding loop (residues 559-566) as well as in helix alpha3 (residues 594-603) interact with Asp residues 272, 274, 277 and sulfated Tyr residues 278 and 279 in GP Ibalpha. To verify this model, 14 mutant A1 domain fragments containing single or multiple side chain substitutions were tested for their ability to mediate platelet adhesion under flow. Each of the vWF residues Tyr(565), Glu(596), and Lys(599) proved to be strictly required for A1 domain function, which, in agreement with previous findings, was also dependent on Gly(561). Moreover, an accessory functional role was apparent for a group of positively charged residues, including Arg at positions 629, 632, 636 and Lys at positions 643 and 645, possibly acting in concert. There was, however, no evidence from the model that these residues directly participate in forming the complex with GP Ibalpha. These results provide a partial model of the vWF-GP Ibalpha interaction linked to the manifestation of functional activity in platelet adhesion.

Amino Acids↗

PubMedDate not supplied

How does consensus scoring work for virtual library screening? An idealized computer experiment.

It has been reported recently that consensus scoring, which combines multiple scoring functions in binding affinity estimation, leads to higher hit-rates in virtual library screening studies. This method seems quite independent to the target receptor, the docking program, or even the scoring functions under investigation. Here we present an idealized computer experiment to explore how consensus scoring works. A hypothetical set of 5000 compounds is used to represent a chemical library under screening. The binding affinities of all its member compounds are assigned by mimicking a real situation. Based on the assumption that the error of a scoring function is a random number in a normal distribution, the predicted binding affinities were generated by adding such a random number to the "observed" binding affinities. The relationship between the hit-rates and the number of scoring functions employed in scoring was then investigated. The performance of several typical ranking strategies for a consensus scoring procedure was also explored. Our results demonstrate that consensus scoring outperforms any single scoring for a simple statistical reason: the mean value of repeated samplings tends to be closer to the true value. Our results also suggest that a moderate number of scoring functions, three or four, are sufficient for the purpose of consensus scoring. As for the ranking strategy, both the rank-by-number and the rank-by-rank strategy work more effectively than the rank-by-vote strategy.

Computer Simulation↗