β-Glucoronidase is an important glycosidase enzyme which catalyzes the hydrolysis of complex carbohydrates into simplest monomeric units. Its over-expression relates with, several type of cancers, including breast, colon and prostate cancer. To treat these disorders the available drug brands on the market are silymyrin, and its derivatives. Some other drugs such as Nialamide, Isocarboxazid, and Phenelzine have also been reported to inhibit GUS activity. However, camptothecin, a plant alkaloid, its derivatives hycamptin and camptosar have been reported to have been granted approval for clinical use, but cause severe side effects including cirrhosis of the liver. Therefore, to overcome these adverse effects, there is a strong need to search and identify lead candidates which possesses therapeutic potential against this target receptor.
In searching and identifying inhibitors, we conducted structure-based pharmacophore based virtual screening of an in-house database with large chemical space of a diverse class of compounds. We developed five structure-based-pharmacophore models. Three individual structure-based-pharmacophore models derived from the available PDB I.D, 3LPF, 3LPG, and 3K4D and two structure-based-shared feature and merged feature pharmacophore models derived by using Ligand Scout software 3.0 version.
Pharmacophore-based virtual screening of in-house data-base identified 1,249 hits, along with 66 reported inhibitors dataset these hit candidates (1,315) were subjected for docking studies, by using FRED 3.0.1 version which successfully docked the pharmacophore-based virtually screened hit candidates, FRED docked and score the candidates by using its scoring function Chemgauss-4, which was further rescored by using GOLD software of 5.1 version into Gold-score, Chem-score and ASP score.
Enrichment factor is an essential parameter to evaluate the efficiency of the docking and scoring comparative to a random selection of compounds, therefore enrichment factor was calculated for 5%, 10%, 15% and 20% for hit candidates of in-house data-base. For 5% of data-set enrichment factor of Chemgauss-4 scoring function was found to be as most efficient, while the rest of the 10%, 15% and 20% of data-base Chem-score scoring function of GOLD was found to be as efficient one. Therefore we selected the docked molecules of top ranked 5% enriched data-base and subjected for in-vitro screening. Out of 5% enrichment (68) compounds, 33 compounds were made available for in-vitro screening, Out of these, eleven (11) compounds showed potent inhibitory potential comparative to the standard (D-saccharic acid, 1,4-lactone). These compounds were also evaluated for cytotoxicity assay, and three compounds were found to be completely non-cytotoxic, however the remaining showed moderate cytotoxicity.
In the present application virtual screening hit results, (scaffold hopping) has been successfully performed, and we identified in top 5% enriched data-base, new classes of compound with potent biological activity against β-glucoronidase, which are not cytotoxic against 3T3 mouse fibroblast cell line. Therefore, these compounds will be used to evaluate the β-glucoronidase activity at the in-vivo level, and other further later steps of drug designing and discovery process.
Structure-based pharmacophore mapping was the keen step which took the ligand-receptor information and developed the model, which searched and identified the inhibitors in the large chemical space (8,262) compounds. Structure-based Pharmacophore (SBPs) model was derived from protein-ligand complexes which illustrate the potential interactions exist between ligand and protein. It is a useful tool for medicinal chemists to identify novel ligands which fulfill the pharmacophore requirements and have a high probability of being biologically active. This has been proven and validated from my structure-based pharmacophore mapping and virtual screening. We developed five structure-based-pharmacophore models, three individual structure-based-pharmacophore models derived from the available PDB I.D, 3LPF, 3LPG, and 3K4D and two structure-based shared feature and merged feature pharmacophore models were derived by using Ligand Scout software 3.0 version.
In the models, the most repeated, keen interactions are b/w Glu413, Tyr472 and Phe161 of the active site amino-acid residues with ligands, the similar interactions are also observed in our potent inhibitors. Therefore, an efficient virtual screening defined in terms of new scaffolds hopping (searching of structurally new and novel compounds). Low hit rates of interesting scaffolds are always preferable over high hit rates of already known scaffolds. Usually a series of compound becomes active against a targeted receptor, but here a scaffold hopping results due to structure-based Pharmacophore model. The details of structure-based Pharmacophores along with the interactions are as follows.
In-silico based theoretical step were used at first than experimentally evaluated the identified Hits, which is the rational approach towards drug designing and discovery process. Usually people use in-silico techniques after the experimental work.
During in-silico based screening we first used Lipiniski ROF based filters (Omega filter from Open eye), which filtered the unstable and toxic compounds from data-base, and selected those compounds which followed the drug ability criteria so that our compounds would show non-cytotoxicity.
These compounds were also used to evaluate the cytotoxicity against normal fibroblast 3T3 cell line of mouse. Three compounds were found to be completely non-cytotoxic while, the remaining compounds showed moderate cytotoxicity. The list of the potent inhibitors, compounds with activity data, are as follows:
Following are the bio-assay protocol used to evaluate the biological activities of compounds against the enzyme β-glucoronidase and normal cell line of mouse fibroblast.
β-Glucuronidase inhibition assay protocol: Inhibitory activity of β-Glucuronidase was determined with the help of spectrophotometric method by measuring the absorbance at 405 nm of p-nitro phenol formed from the substrate (p-nitro phenyl-β-D-glucuronide, N1627-250 mg, Sigma Aldrich). The total reaction volume was 250 μL. The compound dissolved in DMSO (100%), which becomes 2% in the ultimate assay (250 μL) and the similar conditions were used for standard (D-saccharin acid 1,4-lactone, Sigma Aldrich). The reaction mixture contained 185 μL of 0.1 M acetate buffer, 5 μL of test compound solution, 10 μL of (1U) enzyme solution (G7396-25KU, Sigma Aldrich) was incubated at 37° C. for 30 min. The plates were read on a multiplate reader (SpectraMax plus 384) at 405 nm after the addition of 50 μL of 0.4 mM p-nitrophenyl-β-D-glucuronide. All assays were performed in triplicate. IC50 Values were calculated by using EZ-Fit software (Perrella Scientific Inc., Amherst, Mass., U.S.A.). These values are the mean of three independent readings.
Cytotoxicity assay Protocol: Cytotoxic activity of compounds was evaluated in 96-well flat-bottomed micro plates by using the standard MTT (3-[4,5-dimethylthiazole-2-yl]-2,5-diphenyl-tetrazolium bromide) colorimetric assay (15). For this purpose, 3T3 (mouse fibroblast) cells were cultured in Dulbecco's Modified Eagle Medium, supplemented with 5% of fetal bovine serum (FBS), 100 IU/ml of penicillin and 100 μg/ml of streptomycin in 75 cm2 flasks, and kept in 5% CO2 incubator at 37° C. Exponentially growing cells were harvested, counted with haemocytometer and diluted with a particular medium. Cell culture with the concentration of 5×104 cells/ml was prepared and introduced (100 μL/well) into 96-well plates.
After overnight incubation, medium was removed and 200 μL of fresh medium was added with different concentrations of compounds (1-30 μM). After 48 hrs, 200 ptL MTT (0.5 mg/ml) was added to each well and incubated further for 4 hrs. Subsequently, 100 μLof DMSO was added to each well. The extent of MTT reduction to formazan within cells was calculated by measuring the absorbance at 540 nm, using a micro plate reader (Spectra Max plus, Molecular Devices, Calif., USA). The cytotoxicity was recorded as concentration causing 50% growth inhibition (IC50) for 3T3 cells. The percent inhibition was calculated by using the following formula
% inhibition=100−((mean of O.D of test compound−mean of O.D of negative control)/(mean of O.D of positive control−mean of O.D of negative control)*100).
The results (% inhibition) were processed by using Soft-Max Pro software (Molecular Device, USA).