Molecular docking is by far the most common method used in protein structure-based virtual screening. This paper presents Panther, a novel ultrafast multipurpose docking tool. In Panther, a simple shape-electrostatic model of the ligand-binding area of the protein is created by utilizing the protein crystal structure. The features of the possible ligands are then compared to the model by using a similarity search algorithm. On average, one ligand can be processed in a few minutes by using classical docking methods, whereas using Panther processing takes <1 s. The presented Panther protocol can be used in several applications, such as speeding up the early phases of drug discovery projects, reducing the number of failures in the clinical phase of the drug development process, and estimating the environmental toxicity of chemicals. Panther-code is available in our web pages (http://www.jyu.fi/panther) free of charge after registration.

Despite the large computational costs of molecular docking, the default scoring functions are often unable to recognize the active hits from the inactive molecules in large-scale virtual screening experiments. Thus, even though a correct binding pose might be sampled during the docking, the active compound or its biologically relevant pose is not necessarily given high enough score to arouse the attention. Various rescoring and post-processing approaches have emerged for improving the docking performance. Here, it is shown that the very early enrichment (number of actives scored higher than 1% of the highest ranked decoys) can be improved on average 2.5-fold or even 8.7-fold by comparing the docking-based ligand conformers directly against the target protein's cavity shape and electrostatics. The similarity comparison of the conformers is performed without geometry optimization against the negative image of the target protein's ligand-binding cavity using the negative image-based (NIB) screening protocol. The viability of the NIB rescoring or the R-NiB, pioneered in this study, was tested with 11 target proteins using benchmark libraries. By focusing on the shape/electrostatics complementarity of the ligand-receptor association, the R-NiB is able to improve the early enrichment of docking essentially without adding to the computing cost. By implementing consensus scoring, in which the R-NiB and the original docking scoring are weighted for optimal outcome, the early enrichment is improved to a level that facilitates effective drug discovery. Moreover, the use of equal weight from the original docking scoring and the R-NiB scoring improves the yield in most cases.

Receiver operating characteristics (ROC) curve with the calculation of area under curve (AUC) is a useful tool to evaluate the performance of biomedical and chemoinformatics data. For example, in virtual drug screening ROC curves are very often used to visualize the efficiency of the used application to separate active ligands from inactive molecules. Unfortunately, most of the available tools for ROC analysis are implemented into commercially available software packages, or are plugins in statistical software, which are not always the easiest to use. Here, we present Rocker, a simple ROC curve visualization tool that can be used for the generation of publication quality images. Rocker also includes an automatic calculation of the AUC for the ROC curve and Boltzmann-enhanced discrimination of ROC (BEDROC). Furthermore, in virtual screening campaigns it is often important to understand the early enrichment of active ligand identification, for this Rocker offers automated calculation routine. To enable further development of Rocker, it is freely available (MIT-GPL license) for use and modifications from our web-site (http://www.jyu.fi/rocker).Graphical Abstract

Sticholysin II (StnII) is a pore-forming toxin from the sea anemone Stichodactyla heliantus which belongs to the large actinoporin family. The toxin binds to sphingomyelin (SM) containing membranes, and shows high binding specificity for this lipid. In this study, we have examined the role of the hydrogen bonding groups of the SM long-chain base (i.e., the 2NH and the 3OH) for StnII recognition. We prepared methylated SM-analogs which had reduced hydrogen bonding capability from 2NH and 3OH. Both surface plasmon resonance experiments, and isothermal titration calorimetry measurements indicated that StnII failed to bind to bilayers containing methylated SM-analogs, whereas clear binding was seen to SM-containing bilayers. StnII also failed to induce calcein release (i.e., pore formation) from vesicles made to contain methylated SM-analogs, but readily induced calcein release from SM-containing vesicles. Molecular modeling of SM docked to the phosphocholine binding site of StnII indicated that the 2NH and 3OH groups were likely to form a hydrogen bond with Tyr135. In addition, it appeared that Tyr111 and Tyr136 could donate hydrogen bonds to phosphate oxygen, thus stabilizing SM binding to the toxin. We conclude that the interfacial hydrogen bonding properties of SM, in addition to the phosphocholine head group, are crucial for high-affinity SM/StnII-interaction.

Reliable and effective virtual high-throughput screening (vHTS) methods are desperately needed to minimize the expenses involved in drug discovery projects. Here, we present an improvement to the negative image-based (NIB) screening: the shape, the electrostatics, and the solvation state of the target protein's ligand-binding site are included into the vHTS. Additionally, the initial vHTS results are postprocessed with molecular mechanics/generalized Born surface area (MMGBSA) calculations to estimate the favorability of ligand-protein interactions. The results show that docking produces very good early enrichment for phosphodiesterase-5 (PDE-5); however, in general, the NIB and the ligand-based screening performed better with or without the added electrostatics. Furthermore, the postprocessing of the NIB screening results using MMGBSA calculations improved the early enrichment for the PDE-5 considerably, thus, making hit discovery affordable.

No abstract

No abstract

Members of the Rab GTPase family are master regulators of vesicle trafficking. When disregulated, they are associated with a number of pathological states. The inhibition of RGGT, an enzyme responsible for post-translational geranylgeranylation of Rab GTPases represents one way to control the activity of these proteins. Because the number of molecules modulating RGGT is limited, we combined molecular modeling with biological assays to ascertain how modifications of phosphonocarboxylates, the first reported RGGT inhibitors, rationally improve understanding of their structure-activity relationship. We have identified the privileged position in the core scaffold of the imidazo[1,2-a]pyridine ring, which can be modified without compromising compounds' potency. Thus modified compounds are micromolar inhibitors of Rab11A prenylation, simultaneously being inactive against Rap1A/Rap1B modification, with the ability to inhibit proliferation of the HeLa cancer cell line. These findings were rationalized by molecular docking, which recognized interaction of phosphonic and carboxylic groups as decisive in phosphonocarboxylate localization in the RGGT binding site.