Biomedical subjects

Richard D Taylor

Publications and source records attributed to Richard D Taylor.

8 recordsLinked to original sources

PubMed2005

Modeling water molecules in protein-ligand docking using GOLD.

We implemented a novel approach to score water mediation and displacement in the protein-ligand docking program GOLD. The method allows water molecules to switch on and off and to rotate around their three principal axes. A constant penalty, sigma(p), representing the loss of rigid-body entropy, is added for water molecules that are switched on, hence rewarding water displacement. We tested the methodology in an extensive validation study. First, sigma(p) is optimized against a training set of 58 protein-ligand complexes. For this training set, our algorithm correctly predicts water mediation/displacement in approximately 92% of the cases. We observed small improvements in the quality of the predicted binding modes for water-mediated complexes. In the second part of this work, an entirely independent set of 225 complexes is used. For this test set, our algorithm correctly predicts water mediation/displacement in approximately 93% of the cases. Improvements in binding mode quality were observed for individual water-mediated complexes.

Algorithms↗

PubMed2005

Principal atrial fibrillation discharges by the new ACC/AHA/ESC classification.

BACKGROUND: The American College of Cardiology, American Heart Association, and European Society of Cardiology Board (ACC/AHA/ESC) 2001 guidelines for management of patients with atrial fibrillation (AF) include a new classification system consisting of 4 categories: first-detected episode; recurrent paroxysmal (self-terminating); recurrent persistent (requiring cardioversion); and permanent. The frequency of hospital discharges within these categories has not been reported. METHODS: The new classification system was applied to 135 consecutive hospital discharges with a principal diagnosis of AF. RESULTS: Classification of AF in these discharged patients included 74 (55%) with first-detected episode; 28 (21%) with recurrent paroxysmal AF; 17 (13%) with recurrent persistent AF; and 16 (12%) with permanent AF. Hypertension (n = 48; 35%) was the most common primary cause of AF, followed by alcohol related (n = 23; 17%), coronary artery disease (n = 20; 15%), and valvular heart disease (n = 17; 12%). For the 102 patients with first-detected and recurrent paroxysmal AF, 71 (69%) converted spontaneously to normal sinus rhythm within 48 hours of admission. Of the 48 patients with a discharge diagnosis of AF, 32 (67%) were receiving anticoagulation therapy. CONCLUSIONS: Most hospital discharges with a principal diagnosis of AF represent the first-detected episode. Diverse causes contribute to AF, and to examine them would help direct therapy. Importantly, in our analysis, 69% of those patients with first-detected or recurrent paroxysmal AF converted spontaneously to normal sinus rhythm within 48 hours of admission.

Adult↗

PubMed2005

Comparing protein-ligand docking programs is difficult.

There is currently great interest in comparing protein-ligand docking programs. A review of recent comparisons shows that it is difficult to draw conclusions of general applicability. Statistical hypothesis testing is required to ensure that differences in pose-prediction success rates and enrichment rates are significant. Numerical measures such as root-mean-square deviation need careful interpretation and may profitably be supplemented by interaction-based measures and visual inspection of dockings. Test sets must be of appropriate diversity and of good experimental reliability. The effects of crystal-packing interactions may be important. The method used for generating starting ligand geometries and positions may have an appreciable effect on docking results. For fair comparison, programs must be given search problems of equal complexity (e.g. binding-site regions of the same size) and approximately equal time in which to solve them. Comparisons based on rescoring require local optimization of the ligand in the space of the new objective function. Re-implementations of published scoring functions may give significantly different results from the originals. Ostensibly minor details in methodology may have a profound influence on headline success rates.

Algorithms↗

PubMed2005

Prediction of binding modes for ligands in the cytochromes P450 and other heme-containing proteins.

The cytochromes P450 (P450s) are a family of heme-containing monooxygenase enzymes involved in a variety of functions, including the metabolism of endogenous and exogenous substances in the human body. During lead optimization, and in drug development, many potential drug candidates are rejected because of the affinity they display for drug-metabolising P450s. Recently, crystal structures of human enzymes involved in drug metabolism have been determined, significantly augmenting the prospect of using structure-based design to modulate the binding and metabolizing properties of compounds against P450 proteins. An important step in the application of structure-based metabolic optimization is the accurate prediction of docking modes in heme binding proteins. In this paper we assess the performance of the docking program GOLD at predicting the binding mode of 45 heme-containing complexes. We achieved success rates of 64% and 57% for Chemscore and Goldscore respectively; these success rates are significantly lower than the value of 79% observed with both scoring functions for the full GOLD validation set. Re-parameterization of metal-acceptor interactions and lipophilicity of planar nitrogen atoms in the scoring functions resulted in a significant increase in the percentage of successful dockings against the heme binding proteins (Chemscore 73%, Goldscore 65%). The modified scoring functions will be useful in docking applications on P450 enzymes and other heme binding proteins.

Algorithms↗

PubMed2004

The parameterization and validation of generalized born models using the pairwise descreening approximation.

Generalized Born Surface Area (GBSA) models for water using the Pairwise Descreening Approximation (PDA) have been parameterized by two different methods. The first method, similar to that used in previously reported parameterizations, optimizes all parameters against the experimental free energies of hydration of organic molecules. The second method optimizes the PDA parameters to compensate only for systematic errors of the PDA. The best models are compared to Poisson-Boltzmann calculations and applied to the computation of potentials of mean force (PMFs) for the association of various molecules. PMFs present a more rigorous test of the ability of a solvation model to correctly reproduce the screening of intermolecular interactions by the solvent, than its accuracy at predicting free energies of hydration of small molecules. Models derived with the first method are sometimes shown to fail to compute accurate potentials of mean force because of large errors in the computation of Born radii, while no such difficulties are observed with the second method. Furthermore, accurate computation of the Born radii appears to be more important than good agreement with experimental free energies of solvation. We discuss the source of errors in the potentials of mean force and suggest means to reduce them. Our findings suggest that Generalized Born models that use the Pairwise Descreening Approximation and that are derived solely by unconstrained optimization of parameters against free energies of hydration should be applied to the modeling of intermolecular interactions with caution.

Journal Article↗

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↗

PubMed2003

FDS: flexible ligand and receptor docking with a continuum solvent model and soft-core energy function.

The docking of flexible small molecule ligands to large flexible protein targets is addressed in this article using a two-stage simulation-based method. The methodology presented is a hybrid approach where the first component is a dock of the ligand to the protein binding site, based on deriving sets of simultaneously satisfied intermolecular hydrogen bonds using graph theory and a recursive distance geometry algorithm. The output structures are reduced in number by cluster analysis based on distance similarities. These structures are submitted to a modified Monte Carlo algorithm using the AMBER-AA molecular mechanics force field with the Generalized Born/Surface Area (GB/SA) continuum model. This solvent model is not only less expensive than an explicit representation, but also yields increased sampling. Sampling is also increased using a rotamer library to direct some of the protein side-chain movements along with large dihedral moves. Finally, a softening function for the nonbonded force field terms is used, enabling the potential energy function to be slowly turned on throughout the course of the simulation. The docking procedure is optimized, and the results are presented for a single complex of the arabinose binding protein. It was found that for a rigid receptor model, the X-ray binding geometry was reproduced and uniquely identified based on the associated potential energy. However, when side-chain flexibility was included, although the X-ray structure was identified, it was one of three possible binding geometries that were energetically indistinguishable. These results suggest that on relaxing the constraint on receptor flexibility, the docking energy hypersurface changes from being funnel-like to rugged. A further 14 complexes were then examined using the optimized protocol. For each complex the docking methodology was tested for a fully flexible ligand, both with and without protein side-chain flexibility. For the rigid protein docking, 13 out of the 15 test cases were able to find the experimental binding mode; this number was reduced to 11 for the flexible protein docking. However, of these 11, in the majority of cases the experimental binding mode was not uniquely identified, but was present in a cluster of low energy structures that were energetically indistinguishable. These results not only support the presence of a rugged docking energy hypersurface, but also suggest that it may be necessary to consider the possibility of more than one binding conformation during ligand optimization.

Algorithms↗

PubMedDate not supplied

Virtual screening using protein-ligand docking: avoiding artificial enrichment.

This study addresses a number of topical issues around the use of protein-ligand docking in virtual screening. We show that, for the validation of such methods, it is key to use focused libraries (containing compounds with one-dimensional properties, similar to the actives), rather than "random" or "drug-like" libraries to test the actives against. We also show that, to obtain good enrichments, the docking program needs to produce reliable binding modes. We demonstrate how pharmacophores can be used to guide the dockings and improve enrichments, and we compare the performance of three consensus-ranking protocols against ranking based on individual scoring functions. Finally, we show that protein-ligand docking can be an effective aid in the screening for weak, fragment-like binders, which has rapidly become a popular strategy for hit identification. All results presented are based on carefully constructed virtual screening experiments against four targets, using the protein-ligand docking program GOLD.

Ligands↗