What do we know and when do we know it?

Two essential aspects of virtual screening are considered: experimental design and performance metrics. In the design of any retrospective virtual screen, choices have to be made as to the purpose of the exercise. Is the goal to compare methods? Is the interest in a particular type of target or all targets? Are we simulating a ‘real-world’ setting, or teasing out distinguishing features of a method? What are the confidence limits for the results? What should be reported in a publication? In particular, what criteria should be used to decide between different performance metrics? Comparing the field of molecular modeling to other endeavors, such as medical statistics, criminology, or computer hardware evaluation indicates some clear directions. Taken together these suggest the modeling field has a long way to go to provide effective assessment of its approaches, either to itself or to a broader audience, but that there are no technical reasons why progress cannot be made

Application of the PM6 semi-empirical method to modeling proteins enhances docking accuracy of AutoDock

<p>Abstract</p> <p>Background</p> <p>Molecular docking methods are commonly used for predicting binding modes and energies of ligands to proteins. For accurate complex geometry and binding energy estimation, an appropriate method for calculating partial charges is essential. AutoDockTools software, the interface for preparing input files for one of the most widely used docking programs AutoDock 4, utilizes the Gasteiger partial charge calculation method for both protein and ligand charge calculation. However, it has already been shown that more accurate partial charge calculation - and as a consequence, more accurate docking- can be achieved by using quantum chemical methods. For docking calculations quantum chemical partial charge calculation as a routine was only used for ligands so far. The newly developed Mozyme function of MOPAC2009 allows fast partial charge calculation of proteins by quantum mechanical semi-empirical methods. Thus, in the current study, the effect of semi-empirical quantum-mechanical partial charge calculation on docking accuracy could be investigated.</p> <p>Results</p> <p>The docking accuracy of AutoDock 4 using the original AutoDock scoring function was investigated on a set of 53 protein ligand complexes using Gasteiger and PM6 partial charge calculation methods. This has enabled us to compare the effect of the partial charge calculation method on docking accuracy utilizing AutoDock 4 software. Our results showed that the docking accuracy in regard to complex geometry (docking result defined as accurate when the RMSD of the first rank docking result complex is within 2 Å of the experimentally determined X-ray structure) significantly increased when partial charges of the ligands and proteins were calculated with the semi-empirical PM6 method.</p> <p>Out of the 53 complexes analyzed in the course of our study, the geometry of 42 complexes were accurately calculated using PM6 partial charges, while the use of Gasteiger charges resulted in only 28 accurate geometries. The binding affinity estimation was not influenced by the partial charge calculation method - for more accurate binding affinity prediction development of a new scoring function for AutoDock is needed.</p> <p>Conclusion</p> <p>Our results demonstrate that the accuracy of determination of complex geometry using AutoDock 4 for docking calculation greatly increases with the use of quantum chemical partial charge calculation on both the ligands and proteins.</p

Rapid and Accurate Prediction and Scoring of Water Molecules in Protein Binding Sites

Water plays a critical role in ligand-protein interactions. However, it is still challenging to predict accurately not only where water molecules prefer to bind, but also which of those water molecules might be displaceable. The latter is often seen as a route to optimizing affinity of potential drug candidates. Using a protocol we call WaterDock, we show that the freely available AutoDock Vina tool can be used to predict accurately the binding sites of water molecules. WaterDock was validated using data from X-ray crystallography, neutron diffraction and molecular dynamics simulations and correctly predicted 97% of the water molecules in the test set. In addition, we combined data-mining, heuristic and machine learning techniques to develop probabilistic water molecule classifiers. When applied to WaterDock predictions in the Astex Diverse Set of protein ligand complexes, we could identify whether a water molecule was conserved or displaced to an accuracy of 75%. A second model predicted whether water molecules were displaced by polar groups or by non-polar groups to an accuracy of 80%. These results should prove useful for anyone wishing to undertake rational design of new compounds where the displacement of water molecules is being considered as a route to improved affinity

Estrogen Receptor Alpha as a Key Target of Red Wine Polyphenols Action on the Endothelium

BACKGROUND: A greater reduction in cardiovascular risk and vascular protection associated with diet rich in polyphenols are generally accepted; however, the molecular targets for polyphenols effects remain unknown. Meanwhile evidences in the literature have enlightened, not only structural similarities between estrogens and polyphenols known as phytoestrogens, but also in their vascular effects. We hypothesized that alpha isoform of estrogen receptor (ERalpha) could be involved in the transduction of the vascular benefits of polyphenols. METHODOLOGY/PRINCIPAL FINDINGS: Here, we used ERalpha deficient mice to show that endothelium-dependent vasorelaxation induced either by red wine polyphenol extract, Provinols, or delphinidin, an anthocyanin that possesses similar pharmacological profile, is mediated by ERalpha. Indeed, Provinols, delphinidin and ERalpha agonists, 17-beta-estradiol and PPT, are able to induce endothelial vasodilatation in aorta from ERalpha Wild-Type but not from Knock-Out mice, by activation of nitric oxide (NO) pathway in endothelial cells. Besides, silencing the effects of ERalpha completely prevented the effects of Provinols and delphinidin to activate NO pathway (Src, ERK 1/2, eNOS, caveolin-1) leading to NO production. Furthermore, direct interaction between delphinidin and ERalpha activator site is demonstrated using both binding assay and docking. Most interestingly, the ability of short term oral administration of Provinols to decrease response to serotonin and to enhance sensitivity of the endothelium-dependent relaxation to acetylcholine, associated with concomitant increased NO production and decreased superoxide anions, was completely blunted in ERalpha deficient mice. CONCLUSIONS/SIGNIFICANCE: This study provides evidence that red wine polyphenols, especially delphinidin, exert their endothelial benefits via ERalpha activation. It is a major breakthrough bringing new insights of the potential therapeutic of polyphenols against cardiovascular pathologies

A Ligand Channel through the G Protein Coupled Receptor Opsin

The G protein coupled receptor rhodopsin contains a pocket within its seven-transmembrane helix (TM) structure, which bears the inactivating 11-cis-retinal bound by a protonated Schiff-base to Lys296 in TM7. Light-induced 11-cis-/all-trans-isomerization leads to the Schiff-base deprotonated active Meta II intermediate. With Meta II decay, the Schiff-base bond is hydrolyzed, all-trans-retinal is released from the pocket, and the apoprotein opsin reloaded with new 11-cis-retinal. The crystal structure of opsin in its active Ops* conformation provides the basis for computational modeling of retinal release and uptake. The ligand-free 7TM bundle of opsin opens into the hydrophobic membrane layer through openings A (between TM1 and 7), and B (between TM5 and 6), respectively. Using skeleton search and molecular docking, we find a continuous channel through the protein that connects these two openings and comprises in its central part the retinal binding pocket. The channel traverses the receptor over a distance of ca. 70 Å and is between 11.6 and 3.2 Å wide. Both openings are lined with aromatic residues, while the central part is highly polar. Four constrictions within the channel are so narrow that they must stretch to allow passage of the retinal β-ionone-ring. Constrictions are at openings A and B, respectively, and at Trp265 and Lys296 within the retinal pocket. The lysine enforces a 90° elbow-like kink in the channel which limits retinal passage. With a favorable Lys side chain conformation, 11-cis-retinal can take the turn, whereas passage of the all-trans isomer would require more global conformational changes. We discuss possible scenarios for the uptake of 11-cis- and release of all-trans-retinal. If the uptake gate of 11-cis-retinal is assigned to opening B, all-trans is likely to leave through the same gate. The unidirectional passage proposed previously requires uptake of 11-cis-retinal through A and release of photolyzed all-trans-retinal through B

An Innovative Strategy for Dual Inhibitor Design and Its Application in Dual Inhibition of Human Thymidylate Synthase and Dihydrofolate Reductase Enzymes

Due to the diligence of inherent redundancy and robustness in many biological networks and pathways, multitarget inhibitors present a new prospect in the pharmaceutical industry for treatment of complex diseases. Nevertheless, to design multitarget inhibitors is concurrently a great challenge for medicinal chemists. We have developed a novel computational approach by integrating the affinity predictions from structure-based virtual screening with dual ligand-based pharmacophore to discover potential dual inhibitors of human Thymidylate synthase (hTS) and human dihydrofolate reductase (hDHFR). These are the key enzymes in folate metabolic pathway that is necessary for the biosynthesis of RNA,DNA, and protein. Their inhibition has found clinical utility as antitumor, antimicrobial, and antiprotozoal agents. A druglike database was utilized to perform dual-target docking studies. Hits identified through docking experiments were mapped over a dual pharmacophore which was developed from experimentally known dual inhibitors of hTS and hDHFR. Pharmacophore mapping procedure helped us in eliminating the compounds which do not possess basic chemical features necessary for dual inhibition. Finally, three structurally diverse hit compounds that showed key interactions at both activesites, mapped well upon the dual pharmacophore, and exhibited lowest binding energies were regarded as possible dual inhibitors of hTS and hDHFR. Furthermore, optimization studies were performed for final dual hit compound and eight optimized dual hits demonstrating excellent binding features at target systems were also regarded as possible dual inhibitors of hTS and hDHFR. In general, the strategy used in the current study could be a promising computational approach and may be generally applicable to other dual target drug designs

Development and validation of an improved algorithm for overlaying flexible molecules

A program for overlaying multiple flexible molecules has been developed. Candidate overlays are generated by a novel fingerprint algorithm, scored on three objective functions (union volume, hydrogen-bond match, and hydrophobic match), and ranked by constrained Pareto ranking. A diverse subset of the best ranked solutions is chosen using an overlay-dissimilarity metric. If necessary, the solutions can be optimised. A multi-objective genetic algorithm can be used to find additional overlays with a given mapping of chemical features but different ligand conformations. The fingerprint algorithm may also be used to produce constrained overlays, in which user-specified chemical groups are forced to be superimposed. The program has been tested on several sets of ligands, for each of which the true overlay is known from protein–ligand crystal structures. Both objective and subjective success criteria indicate that good results are obtained on the majority of these sets