A Practical Perspective : The Effect of Ligand Conformers on the Negative Image-Based Screening

Negative image-based (NIB) screening is a rigid molecular docking methodology that can also be employed in docking rescoring. During the NIB screening, a negative image is generated based on the target protein’s ligand-binding cavity by inverting its shape and electrostatics. The resulting NIB model is a drug-like entity or pseudo-ligand that is compared directly against ligand 3D conformers, as is done with a template compound in the ligand-based screening. This cavity-based rigid docking has been demonstrated to work with genuine drug targets in both benchmark testing and drug candidate/lead discovery. Firstly, the study explores in-depth the applicability of different ligand 3D conformer…

Comparison of virtual high-throughput screening methods for the identification of phosphodiesterase-5 inhibitors.

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-…

Suitability of MMGBSA for the selection of correct ligand binding modes from docking results

The estimation of the correct binding mode and affinity of a ligand into a target protein using computational methods is challenging. However, docking can introduce poses from which the correct binding mode could be identified using other methods. Here, we analyzed the reliability of binding energy estimation using the molecular mechanics‐generalized Born surface area (MMGBSA) method without and with energy minimization to identify the likely ligand binding modes within docking results. MMGBSA workflow (a) outperformed docking in recognizing the correct binding modes of androgen receptor ligands and (b) improved the correlation coefficient of computational and experimental results of rescor…

Structure-Activity Relationship Analysis of 3-Phenylcoumarin-Based Monoamine Oxidase B Inhibitors

Monoamine oxidase B (MAO-B) catalyzes deamination of monoamines such as neurotransmitters dopamine and norepinephrine. Accordingly, small-molecule MAO-B inhibitors potentially alleviate the symptoms of dopamine-linked neuropathologies such as depression or Parkinson's disease. Coumarin with a functionalized 3-phenyl ring system is a promising scaffold for building potent MAO-B inhibitors. Here, a vast set of 3-phenylcoumarin derivatives was designed using virtual combinatorial chemistry or rationally de novo and synthesized using microwave chemistry. The derivatives inhibited the MAO-B at 100 nM−1 μM. The IC50 value of the most potent derivative 1 was 56 nM. A docking-based structure-activi…

Blocking oestradiol synthesis pathways with potent and selective coumarin derivatives

A comprehensive set of 3-phenylcoumarin analogues with polar substituents was synthesised for blocking oestradiol synthesis by 17-b-hydroxysteroid dehydrogenase 1 (HSD1) in the latter part of the sulphatase pathway. Five analogues produced 62% HSD1 inhibition at 5 mM and, furthermore, three of them produced 68% inhibition at 1 mM. A docking-based structure-activity relationship analysis was done to determine the molecular basis of the inhibition and the cross-reactivity of the analogues was tested against oestrogen receptor, aromatase, cytochrome P450 1A2, and monoamine oxidases. Most of the analogues are only modestly active with 17-b-hydroxysteroid dehydrogenase 2 – a requirement for lowe…

Fragment- and negative image-based screening of phosphodiesterase 10A inhibitors.

A novel virtual screening methodology called fragment- and negative image-based (F-NiB) screening is introduced and tested experimentally using phosphodiesterase 10A (PDE10A) as a case study. Potent PDE10A-specific small-molecule inhibitors are actively sought after for their antipsychotic and neuroprotective effects. The F-NiB combines features from both fragment-based drug discovery and negative image-based (NIB) screening methodologies to facilitate rational drug discovery. The selected structural parts of protein-bound ligand(s) are seamlessly combined with the negative image of the target's ligand-binding cavity. This cavity- and fragment-based hybrid model, namely its shape and electr…

Case-specific performance of MM-PBSA, MM-GBSA, and SIE in virtual screening.

In drug discovery the reliable prediction of binding free energies is of crucial importance. Methods that combine molecular mechanics force fields with continuum solvent models have become popular because of their high accuracy and relatively good computational efficiency. In this research we studied the performance of molecular mechanics generalized Born surface area (MM-GBSA), molecular mechanics Poisson-Boltzmann surface area (MM-PBSA), and solvated interaction energy (SIE) both in their virtual screening efficiency and their ability to predict experimentally determined binding affinities for five different protein targets. The protein-ligand complexes were derived with two different app…

Identification of estrogen receptor α ligands with virtual screening techniques.

Utilization of computer-aided molecular discovery methods in virtual screening (VS) is a cost-effective approach to identify novel bioactive small molecules. Unfortunately, no universal VS strategy can guarantee high hit rates for all biological targets, but each target requires distinct, fine-tuned solutions. Here, we have studied in retrospective manner the effectiveness and usefulness of common pharmacophore hypothesis, molecular docking and negative image-based screening as potential VS tools for a widely applied drug discovery target, estrogen receptor α (ERα). The comparison of the methods helps to demonstrate the differences in their ability to identify active molecules. For example,…

Molecular mechanism of T-cell protein tyrosine phosphatase (TCPTP) activation by mitoxantrone.

T-cell protein tyrosine phosphatase (TCPTP) is a ubiquitously expressed non-receptor protein tyrosine phosphatase. It is involved in the negative regulation of many cellular signaling pathways. Thus, activation of TCPTP could have important therapeutic applications in diseases such as cancer and inflammation. We have previously shown that the α-cytoplasmic tail of integrin α1β1 directly binds and activates TCPTP. In addition, we have identified in a large-scale high-throughput screen six small molecules that activate TCPTP. These small molecule activators include mitoxantrone and spermidine. In this study, we have investigated the molecular mechanism behind agonist-induced TCPTP activation.…

Computational studies of biomolecular screening and interactions

Improving Docking Performance Using Negative Image-Based Rescoring

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 th…

Molecular docking-based design and development of a highly selective probe substrate for UDP-glucuronosyltransferase 1A10

Intestinal and hepatic glucuronidation by the UDP-glucuronosyltransferases (UGTs) greatly affect the bioavailability of phenolic compounds. UGT1A10 catalyzes glucuronidation reactions in the intestine, but not in the liver. Here, our aim was to develop selective, fluorescent substrates to easily elucidate UGT1A10 function. To this end, homology models were constructed and used to design new substrates, and subsequently, six novel C3-substituted (4-fluorophenyl, 4-hydroxyphenyl, 4-methoxyphenyl, 4-(dimethylamino)phenyl, 4-methylphenyl, or triazole) 7-hydroxycoumarin derivatives were synthesized from inexpensive starting materials. All tested compounds could be glucuronidated to nonfluorescen…

Identification of the privileged position in the imidazo[1,2-a]pyridine ring of phosphonocarboxylates for development of Rab geranylgeranyl transferase (RGGT) inhibitors

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 r…

Suitability ofMMGBSAfor the selection of correct ligand binding modes from docking results

The estimation of the correct binding mode and affinity of a ligand into a target protein using computational methods is challenging. However, docking can introduce poses from which the correct binding mode could be identified using other methods. Here, we analyzed the reliability of binding energy estimation using the molecular mechanics-generalized Born surface area (MMGBSA) method without and with energy minimization to identify the likely ligand binding modes within docking results. MMGBSA workflow (a) outperformed docking in recognizing the correct binding modes of androgen receptor ligands and (b) improved the correlation coefficient of computational and experimental results of rescor…

Rocker: Open source, easy-to-use tool for AUC and enrichment calculations and ROC visualization

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…

Reliability of Virtual Screening Methods in Prediction of PDE4Binhibitor Activity

Identification of active ligands using computational methods is a challenging task. For example, molecular docking, pharmacophore modeling, and three dimensional quantitative structure-activity relationship models (3D-QSAR) are widely used methods to identify novel small molecules. However, all these methods have, in addition to advantages, also significant pitfalls. The aim of this study was to compare some commonly used computational methods to estimate their ability to separate highly active PDE4B-inhibitors from less active and inactive ones. Here, 152 molecules with pIC 50 -range of 3.4-10.5, originating from six original studies were used. High correlation coefficients by using dockin…

Fragment‐ and Negative Image‐Based Screening of Phosphodiesterase 10A Inhibitors

A novel virtual screening methodology called fragment‐ and negative image‐based (F‐NiB) screening is introduced and tested experimentally using phosphodiesterase 10A (PDE10A) as a case study. Potent PDE10A‐specific small‐molecule inhibitors are actively sought after for their antipsychotic and neuroprotective effects. The F‐NiB combines features from both fragment‐based drug discovery and negative image‐based (NIB) screening methodologies to facilitate rational drug discovery. The selected structural parts of protein‐bound ligand(s) are seamlessly combined with the negative image of the target's ligand‐binding cavity. This cavity‐ and fragment‐based hybrid model, namely its shape and electr…

2NH and 3OH are crucial structural requirements in sphingomyelin for sticholysin II binding and pore formation in bilayer membranes.

AbstractSticholysin 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 bind…

Discovery of Retinoic Acid-Related Orphan Receptor γt Inverse Agonists via Docking and Negative Image-Based Screening

Retinoic acid-related orphan receptor γt (RORγt) has a vital role in the differentiation of T-helper 17 (TH17) cells. Potent and specific RORγt inverse agonists are sought for treating TH17-related diseases such as psoriasis, rheumatoid arthritis, and type 1 diabetes. Here, the aim was to discover novel RORγt ligands using both standard molecular docking and negative image-based screening. Interestingly, both of these in silico techniques put forward mostly the same compounds for experimental testing. In total, 11 of the 34 molecules purchased for testing were verified as RORγt inverse agonists, thus making the effective hit rate 32%. The pIC50 values for the compounds varied from 4.9 (11 μ…

Identification of the Privileged Position in the Imidazo[1,2-a]pyridine Ring of Phosphonocarboxylates for Development of Rab Geranylgeranyl Transferase (RGGT) Inhibitors

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 r…

Structure–activity relationship of sphingomyelin analogs with sphingomyelinase from Bacillus cereus

AbstractThe aim of this study was to examine how structural properties of different sphingomyelin (SM) analogs affected their substrate properties with sphingomyelinase (SMase) from Bacillus cereus. Using molecular docking and dynamics simulations (for SMase–SM complex), we then attempted to explain the relationship between SM structure and enzyme activity. With both micellar and monolayer substrates, 3O-methylated SM was found not to be degraded by the SMase. 2N-methylated SM was a substrate, but was degraded at about half the rate of its 2NH–SM control. PhytoPSM was readily hydrolyzed by the enzyme. PSM lacking one methyl in the phosphocholine head group was a good substrate, but PSM lack…

Improving Docking Performance Using Negative Image-Based Rescoring

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 th…