Cheminformatics-Based Drug Design Approach for Identification of Inhibitors Targeting the Characteristic Residues of MMP-13 Hemopexin Domain

Background

MMP-13, a zinc dependent protease which catalyses the cleavage of type II collagen, is expressed in osteoarthritis (OA) and rheumatoid arthritis (RA) patients, but not in normal adult tissues. Therefore, the protease has been intensively studied as a target for the inhibition of progression of OA and RA. Recent reports suggest that selective inhibition of MMP-13 may be achieved by targeting the hemopexin (Hpx) domain of the protease, which is critical for substrate specificity. In this study, we applied a cheminformatics-based drug design approach for the identification and characterization of inhibitors targeting the amino acid residues characteristic to Hpx domain of MMP-13; these inhibitors may potentially be employed in the treatment of OA and RA.

Methodology/Principal Findings

Sequence-based mutual information analysis revealed five characteristic (completely conserved and unique), putative functional residues of the Hpx domain of MMP-13 (these residues hereafter are referred to as HCR-13_pf). Binding of a ligand to as many of the HCR-13_pf is postulated to result in an increased selective inhibition of the Hpx domain of MMP-13. Through the in silico structure-based high-throughput virtual screening (HTVS) method of Glide, against a large public library of 16908 molecules from Maybridge, PubChem and Binding, we identified 25 ligands that interact with at least one of the HCR-13_pf. Assessment of cross-reactivity of the 25 ligands with MMP-1 and MMP-8, members of the collagenase family as MMP-13, returned seven lead molecules that did not bind to any one of the putative functional residues of Hpx domain of MMP-1 and any of the catalytic active site residues of MMP-1 and -8, suggesting that the ligands are not likely to interact with the functional or catalytic residues of other MMPs. Further, in silico analysis of physicochemical and pharmacokinetic parameters based on Lipinski''s rule of five and ADMET (absorption, distribution, metabolism, excretion and toxicity) respectively, suggested potential utility of the compounds as drug leads.

Conclusions/Significance

We have identified seven distinct drug-like molecules binding to the HCR-13_pf of MMP-13 with no observable cross-reactivity to MMP-1 and MMP-8. These molecules are potential selective inhibitors of MMP-13 that can be experimentally validated and their backbone structural scaffold could serve as building blocks in designing drug-like molecules for OA, RA and other inflammatory disorders. The systematic cheminformatics-based drug design approach applied herein can be used for rational search of other public/commercial combinatorial libraries for more potent molecules, capable of selectively inhibiting the collagenolytic activity of MMP-13.     

Identification of potential inhibitors for HCV NS3 genotype 4a by combining protein–ligand interaction fingerprint, 3D pharmacophore,docking, and dynamic simulation

HCV NS3 protease domain has been one of the most attractive targets for developing new drugs for HCV infection and many drugs were successfully developed, but all of them were designed for targeting HCV genotype 1 infection. HCV genotype 4a dominant in Egypt has paid less attention. Here, we describe our protocol of virtual screening in identification of novel potential potent inhibitors for HCV NS3 of genotype 4a using homology modeling, PLIF (protein–ligand interaction fingerprint), docking, pharmacophore, and dynamic simulation. A high-quality 3D model of HCV NS3 protease of genotype 4a was constructed using crystal structure of HCV NS3 protease of genotype 1b (PDB ID: 4u01) as a template. PLIF was generated using five crystal structures of HCV NS3 (PDB ID: 4u01, 3kee, 4ktc, 4i33, and 5epn) which revealed the most important residues and their interactions with the co-crystalized ligands. A 3D pharmacophore model consisting of six features was developed from the generated PLIF data and then used as a screening filter for 11,244 compounds. Only 423 compounds passed the pharmacophore filter and entered the docking-based virtual screening stage. The highest ranked five hits from docking result (compound (C1–C5)) were selected for further analysis. They exhibited stronger interaction and higher binding affinity than HCV NS3 protease ligands. Dynamic simulation of the protein–best lead complex was performed to validate and augment the virtual screening results and it showed that these compounds have a strong binding affinity and could be very effective in treating HCV genotype 4a infections.     

Expedient solid-phase synthesis of both symmetric and asymmetric diol libraries targeting aspartic proteases

C₂-symmetric diols have been shown to be highly potent against HIV-1 protease (PR). However, gaining access to these compounds has been hampered by the need of multistep solution-phase reactions which are often tedious and inefficient. In this Letter, we have disclosed a solid-phase strategy for rapid preparation of small molecule-based, symmetric and asymmetric diols as potential HIV-1 protease inhibitors. Upon biological screening, we found one of them, SYM-5, to be a potent and selective inhibitor (K_i = 400 nM) against HIV-1 protease.     

HIV-1 Protease in the Fission Yeast Schizosaccharomyces pombe

Background

HIV-1 protease (PR) is an essential viral enzyme. Its primary function is to proteolyze the viral Gag-Pol polyprotein for production of viral enzymes and structural proteins and for maturation of infectious viral particles. Increasing evidence suggests that PR cleaves host cellular proteins. However, the nature of PR-host cellular protein interactions is elusive. This study aimed to develop a fission yeast (Schizosaccharomyces pombe) model system and to examine the possible interaction of HIV-1 PR with cellular proteins and its potential impact on cell proliferation and viability.

Results

A fission yeast strain RE294 was created that carried a single integrated copy of the PR gene in its chromosome. The PR gene was expressed using an inducible nmt1 promoter so that PR-specific effects could be measured. HIV-1 PR from this system cleaved the same indigenous viral p6/MA protein substrate as it does in natural HIV-1 infections. HIV-1 PR expression in fission yeast cells prevented cell proliferation and induced cellular oxidative stress and changes in mitochondrial morphology that led to cell death. Both these PR activities can be prevented by a PR-specific enzymatic inhibitor, indinavir, suggesting that PR-mediated proteolytic activities and cytotoxic effects resulted from enzymatic activities of HIV-1 PR. Through genome-wide screening, a serine/threonine kinase, Hhp2, was identified that suppresses HIV-1 PR-induced protease cleavage and cell death in fission yeast and in mammalian cells, where it prevented PR-induced apoptosis and cleavage of caspase-3 and caspase-8.

Conclusions

This is the first report to show that HIV-1 protease is functional as an enzyme in fission yeast, and that it behaves in a similar manner as it does in HIV-1 infection. HIV-1 PR-induced cell death in fission yeast could potentially be used as an endpoint for mechanistic studies, and this system could be used for developing a high-throughput system for drug screenings.     

Identification of Novel HIV 1- Protease Inhibitors: Application of Ligand and Structure Based Pharmacophore Mapping and Virtual Screening

A combined ligand and structure-based drug design approach provides a synergistic advantage over either methods performed individually. Present work bestows a good assembly of ligand and structure-based pharmacophore generation concept. Ligand-oriented study was accomplished by employing the HypoGen module of Catalyst in which we have translated the experimental findings into 3-D pharmacophore models by identifying key features (four point pharmacophore) necessary for interaction of the inhibitors with the active site of HIV-1 protease enzyme using a training set of 33 compounds belonging to the cyclic cyanoguanidines and cyclic urea derivatives. The most predictive pharmacophore model (hypothesis 1), consisting of four features, namely, two hydrogen bond acceptors and two hydrophobic, showed a correlation (r) of 0.90 and a root mean square of 0.71 and cost difference of 56.59 bits between null cost and fixed cost. The model was validated using CatScramble technique, internal and external test set prediction. In the second phase of our study, a structure-based five feature pharmacophore hypothesis was generated which signifies the importance of hydrogen bond donor, hydrogen bond acceptors and hydrophobic interaction between the HIV-1 protease enzyme and its inhibitors. This work has taken a significant step towards the full integration of ligand and structure-based drug design methodologies as pharmacophoric features retrieved from structure-based strategy complemented the features from ligand-based study hence proving the accuracy of the developed models. The ligand-based pharmacophore model was used in virtual screening of Maybridge and NCI compound database resulting in the identification of four structurally diverse druggable compounds with nM activities.     

Synthesis of dammarane-type triterpene derivatives and their ability to inhibit HIV and HCV proteases

We synthesized dammarane-type triterpene derivatives and evaluated their ability to inhibit HIV-1 and HCV proteases to understand their structure–activity relationships. All of the mono- and di-succinyl derivatives (5a–5f) were powerful inhibitors of HIV-1 protease (IC₅₀ < 10 μM). However, only di-succinyl (5e) and 2,3-seco-2,3-dioic acid (3b) derivatives similarly inhibited HCV protease (IC₅₀ < 10 μM). A-nor dammarane-type triterpenes (4a and 4b, IC₅₀ 10.0 and 29.9 μM, respectively) inhibited HIV-1 protease moderately or strongly, but were inactive against HCV protease. All compounds that powerfully inhibited HIV-1 or HCV protease did not appreciably inhibit the general human proteases, renin and trypsin (IC₅₀ > 1000 μM). These findings indicated that the mono-succinyl dammarane type derivatives (5a–5d) selectively inhibited HIV-1 protease and that the di-succinyl (5e, 5f) as well as 2,3-seco-2,3-dioic acid (3b) derivatives preferably inhibited both viral proteases.     

Bioassays to monitor Taspase1 function for the identification of pharmacogenetic inhibitors

Background

Threonine Aspartase 1 (Taspase1) mediates cleavage of the mixed lineage leukemia (MLL) protein and leukemia provoking MLL-fusions. In contrast to other proteases, the understanding of Taspase1''s (patho)biological relevance and function is limited, since neither small molecule inhibitors nor cell based functional assays for Taspase1 are currently available.

Methodology/Findings

Efficient cell-based assays to probe Taspase1 function in vivo are presented here. These are composed of glutathione S-transferase, autofluorescent protein variants, Taspase1 cleavage sites and rational combinations of nuclear import and export signals. The biosensors localize predominantly to the cytoplasm, whereas expression of biologically active Taspase1 but not of inactive Taspase1 mutants or of the protease Caspase3 triggers their proteolytic cleavage and nuclear accumulation. Compared to in vitro assays using recombinant components the in vivo assay was highly efficient. Employing an optimized nuclear translocation algorithm, the triple-color assay could be adapted to a high-throughput microscopy platform (Z''factor = 0.63). Automated high-content data analysis was used to screen a focused compound library, selected by an in silico pharmacophor screening approach, as well as a collection of fungal extracts. Screening identified two compounds, N-[2-[(4-amino-6-oxo-3H-pyrimidin-2-yl)sulfanyl]ethyl]benzenesulfonamide and 2-benzyltriazole-4,5-dicarboxylic acid, which partially inhibited Taspase1 cleavage in living cells. Additionally, the assay was exploited to probe endogenous Taspase1 in solid tumor cell models and to identify an improved consensus sequence for efficient Taspase1 cleavage. This allowed the in silico identification of novel putative Taspase1 targets. Those include the FERM Domain-Containing Protein 4B, the Tyrosine-Protein Phosphatase Zeta, and DNA Polymerase Zeta. Cleavage site recognition and proteolytic processing of these substrates were verified in the context of the biosensor.

Conclusions

The assay not only allows to genetically probe Taspase1 structure function in vivo, but is also applicable for high-content screening to identify Taspase1 inhibitors. Such tools will provide novel insights into Taspase1''s function and its potential therapeutic relevance.     

Molecular docking and structure-based virtual screening studies of potential drug target,CAAX prenyl proteases,of Leishmania donovani

Targeting CAAX prenyl proteases of Leishmania donovani can be a good approach towards developing a drug molecule against Leishmaniasis. We have modeled the structure of CAAX prenyl protease I and II of L. donovani, using homology modeling approach. The structures were further validated using Ramachandran plot and ProSA. Active site prediction has shown difference in the amino acid residues present at the active site of CAAX prenyl protease I and CAAX prenyl protease II. The electrostatic potential surface of the CAAX prenyl protease I and II has revealed that CAAX prenyl protease I has more electropositive and electronegative potentials as compared CAAX prenyl protease II suggesting significant difference in their activity. Molecular docking with known bisubstrate analog inhibitors of protein farnesyl transferase and peptidyl (acyloxy) methyl ketones reveals significant binding of these molecules with CAAX prenyl protease I, but comparatively less binding with CAAX prenyl protease II. New and potent inhibitors were also found using structure-based virtual screening. The best docked compounds obtained from virtual screening were subjected to induced fit docking to get best docked configurations. Prediction of drug-like characteristics has revealed that the best docked compounds are in line with Lipinski’s rule. Moreover, best docked protein–ligand complexes of CAAX prenyl protease I and II are found to be stable throughout 20 ns simulation. Overall, the study has identified potent drug molecules targeting CAAX prenyl protease I and II of L. donovani whose drug candidature can be verified further using biochemical and cellular studies.     

Whiskers-less HIV-protease: a possible way for HIV-1 deactivation

Background

Among viral enzymes, the human HIV-1 protease comprises the most interesting target for drug discovery. There are increasing efforts focused on designing more effective inhibitors for HIV-1 protease in order to prevent viral replication in AIDS patients. The frequent and continuous mutation of HIV-1 protease gene creates a formidable obstacle for enzyme inhibition which could not be overcome by the traditional single drug therapy. Nowadays, in vitro and in silico studies of protease inhibition constitute an advanced field in biological researches. In this article, we tried to simulate protease-substrate complexes in different states; a native state and states with whiskers deleted from one and two subunits. Molecular dynamic simulations were carried out in a cubic box filled with explicit water at 37°C and in 1atomsphere of pressure.

Results

Our results showed that whisker truncation of protease subunits causes the dimer structure to decrease in compactness, disrupts substrate-binding site interactions and changes in flap status simultaneously.

Conclusions

Based on our findings we claim that whisker truncation even when applied to a single subunit, threats dimer association which probably leads to enzyme inactivation. We may postulate that inserting a gene to express truncated protease inside infected cells can interfere with protease dimerization. The resulted proteases would presumably have a combination of native and truncated subunits in their structures which exert no enzyme activities as evidenced by the present work. Our finding may create a new field of research in HIV gene therapy for protease inhibition, circumventing problems of drug resistance.     

Striped murrel S1 family serine protease: immune characterization,antibacterial property and enzyme activities

We report a molecular characterization of S1 family serine protease (SP-1) from snakehead murrel (or called striped murrel) Channa striatus (Cs). CsSP-1 polypeptide contained a catalytic core domain (otherwise known as serine protease trypsin domain) between H20 and I237 along with a catalytic triad at H⁶¹, D¹⁰⁴ and S¹⁹⁷. Phylogenetic analysis confirmed that CsSP-1 belongs to serine protease S1 family. The tertiary structure showed that CsSP-1 contains 14 β-sheets as 2 separate β-barrels (the first β-barrel consists of 8 β-sheets in the N-terminal region and the second β-barrel consists of 6 β-sheets in the C-terminal region) and 3 α-helical regions. Significantly (P < 0.05) the highest CsSP-1 mRNA expression was observed in intestine, liver and kidney, moderate expression was seen in spleen, head kidney, skin and blood, and the lowest one in brain, gill, muscle and heart. Further, the expression was induced in intestine with fungus Aphanomyces invadans and bacteria Aeromonas hydrophila. The recombinant CsSP-1 protein showed antibacterial activity against both gram-negative and gram-positive bacteria. The optimum CsSP-1 enzyme activity against the substrate casein was determined at 8 mM casein concentration. Moreover, the activity was highly influenced by 5 mM phenyl-methylsulfonyl fluoride followed by ethylenediaminetetraacetic acid, 4-(2-aminoethyl)benzenesulfonylfluoride hydrochloride and calpain inhibitor I. The CsSP-1 enzyme exhibited the highest activity at pH 7.5 and temperature 35°C. The overall results showed the potential involvement of CsSP-1 in the immune system of murrels. However, further research is necessary to study the mechanism of implicit trypsin association in the defence process.     

Computer aided screening of Accacia nilotica phytochemicals against HCV NS3/4a

Background: HCV has become a leading cause of liver cirrhosis and hepatocellular carcinoma and is a major health concern worldwide. To date, there is no vaccine available in the market to tackle this disease, therefore there is a strong need to develop antiviral compounds that can target all genotypes of HCV with the same efficiency. Medicinal plants have low cost and are less toxic therefore, extracts of medicinal plants can serve as important antiviral agents against HCV. This study was designed to screen phytochemicals of Accacia nilotica to find a potent drug candidate that can inhibit HCV infection effectively.Results: Docking of NS3/4A protease and Flavonoids of Accacia nilotica revealed that most of the flavonoids bound deeply with the active site of NS3/4A protease. Compound 01 showed a high ranking on docking score. All other compounds also showed reliable docking scores and had interactions with the binding cavity of NS3/4A protease, suggesting them as a potent drug candidate to block HCV replication.Conclusion: To recognize binding interactions of Accacia nilotica phytochemicals with NS3/4A protease, molecular docking was performed to find potential inhibitor against NS3/4A protease of HCV. After post docking analysis, important interactions were found between active compounds and active site of NS3/4A protease. It can be concluded from the study that phytochemicals of Accacia nilotica may serve as a potential drug candidate with relatively simple structural changes against HCV NS3/4A protease.     

Antibacterial and antioxidant activities of four kaempferol methyl ethers isolated from Dodonaea viscosa Jacq. var. angustifolia leaf extracts

Fractionation of dichloromethane and acetone fractions obtained by serial extraction from the leaf powder of Dodonaea viscosa Jacq. var. angustifolia (Sapindaceae) resulted in the isolation of four kaempferol methyl ethers. The compounds were identified by spectral data (¹H NMR, ¹³C NMR and MS) as: 3, 5, 7-trihydroxy-4'-methoxyflavone (1); 5, 7, 4'-trihydroxy-3, 6-dimethoxyflavone (2); 5, 7-dihydroxy-3, 6, 4'-trimethoxyflavone (santin) (3); and 5-hydroxy -3, 7, 4'-trimethoxyflavone (4) together with 3,4',5,7-tetrahydroxy flavone (kaempferol) (5). Antioxidant potential of the compounds was evaluated using a DPPH spectrophotometric assay, while antibacterial activity was determined using a serial dilution microplate technique. The isolates demonstrated varying degrees of antioxidant and antibacterial activities. Of all the compounds investigated, compounds 1 and 5 demonstrated some antioxidant activity (EC₅₀ = 75.49 ± 1.76 µM and 35.06 ± 0.85 respectively) but lower than l-ascorbic acid (EC₅₀ = 13.55 ± 0.28 µM) used as a standard antioxidant agent. The minimum inhibitory concentration (MIC) of isolated compounds against Staphylococcus aureus, Enterococcus faecalis, Escherichia coli and Pseudomonas aeruginosa varied from 16 µg/ml to more than 250 µg/ml. Some structure activity relationships could be established for these compounds.     

Inhibitors of Helicobacter pylori protease HtrA found by 'virtual ligand' screening combat bacterial invasion of epithelia

Background

The human pathogen Helicobacter pylori (H. pylori) is a main cause for gastric inflammation and cancer. Increasing bacterial resistance against antibiotics demands for innovative strategies for therapeutic intervention.

Methodology/Principal Findings

We present a method for structure-based virtual screening that is based on the comprehensive prediction of ligand binding sites on a protein model and automated construction of a ligand-receptor interaction map. Pharmacophoric features of the map are clustered and transformed in a correlation vector (‘virtual ligand’) for rapid virtual screening of compound databases. This computer-based technique was validated for 18 different targets of pharmaceutical interest in a retrospective screening experiment. Prospective screening for inhibitory agents was performed for the protease HtrA from the human pathogen H. pylori using a homology model of the target protein. Among 22 tested compounds six block E-cadherin cleavage by HtrA in vitro and result in reduced scattering and wound healing of gastric epithelial cells, thereby preventing bacterial infiltration of the epithelium.

Conclusions/Significance

This study demonstrates that receptor-based virtual screening with a permissive (‘fuzzy’) pharmacophore model can help identify small bioactive agents for combating bacterial infection.     

Evaluation of the amino acid binding site of Mycobacterium tuberculosis glutamine synthetase for drug discovery

A combination of a literature survey, structure-based virtual screening and synthesis of a small library was performed to identify hits to the potential antimycobacterial drug target, glutamine synthetase. The best inhibitor identified from the literature survey was (2S,5R)-2,6-diamino-5-hydroxyhexanoic acid (4, IC(50) of 610+/-15microM). In the virtual screening 46,400 compounds were docked and subjected to a pharmacophore search. Of these compounds, 29 were purchased and tested in a biological assay, allowing three novel inhibitors containing an aromatic scaffold to be identified. Based on one of the hits from the virtual screening a small library of 15 analogues was synthesized producing four compounds that inhibited glutamine synthetase.