Ligand-based and structure-based virtual screening to identify carbonic anhydrase IX inhibitors

A three-dimensional pharmacophore model of CA IX inhibitors was generated and used to screen the ZINC database of commercially available compounds. The hits were docked in a CA IX homology model. By visualizing the binding mode and score of these compounds, six derivatives were selected and evaluated for their inhibitory potency against CA IX. A highly active CA IX inhibitor was identified which may be used as a lead to design novel such derivatives.     

Discovering novel carbonic anhydrase type IX (CA IX) inhibitors from seven million compounds using virtual screening and in vitro analysis

Carbonic anhydrase type IX (CA IX) enzyme is mostly over expressed in different cancer cell lines and tumor tissues. Potent CA IX inhibitors can be effective for adjusting the pH imbalance in tumor cells. In the present work, we represented the successful application of high throughput virtual screening (HTVS) of large dataset from ZINC database included of ～7 million compounds to discover novel inhibitors of CA IX. HTVS and molecular docking were performed using consequence Glide/standard precision (SP), extra precision (XP) and induced fit docking (IFD) molecular docking protocols. For each compound, docking code calculates a set of low-energy poses and then exhaustively scans the binding pocket of the target with small compounds. Novel CA IX inhibitor candidates were suggested based on molecular modeling studies and a few of them were tested using in vitro analysis. These compounds were determined as good inhibitors against human CA IX target with K_i in the range of 0.85–1.58?μM. In order to predict the pharmaceutical properties of the selected compounds, ADME (absorption, distribution, metabolism and excretion) analysis was also carried out.     

An Insight into the Structural Requirements and Pharmacophore Identification of Carbonic Anhydrase Inhibitors to Combat Oxidative Stress at High Altitudes: An In-Silico Approach

Carbonic anhydrases (CA) inhibitory action could be linked to the treatment of a number of ailments, including cancer, osteoporosis, glaucoma, and several neurological problems. For the development of effective CA inhibitors, a variety of heterocyclic rings have been investigated. Furthermore, at high altitudes, oxygen pressure drops, resulting in the formation of reactive oxygen and nitrogen species, and CA inhibitors having role in combating this oxidative stress. Acetazolamide contains thiadiazole ring, which has aroused researchers’ interest because of its CA inhibitory action. In the present study, we used a number of drug design tools, such as pharmacophore modeling, 3D QSAR, docking, and virtual screening on twenty-seven 1,3,4-thiadiazole derivatives that have been described as potential CA inhibitors in the literature. An atom-based 3D-QSAR analysis was carried out to determine the contribution of individual atoms to model generation, while a pharmacophore mapping investigation was carried out to find the common unique pharmacophoric properties required for biological activity. The coefficient of determination for both the training and test sets were statistically significant in the generated model. The best QSAR model was chosen based on the values of R² (0.8757) and Q² (0.7888). A molecular docking study was also conducted against the most potent analogue 4m, which has the highest SP docking score (−5.217) (PDB ID: 6g3v). The virtual screening revealed a number of promising compounds. The screened compound ZINC77699643 interacted with the amino acid residues, Pro201 and Thr199, in the virtual screening study (PDB ID: 6g3v). These interactions demonstrated the significance of the CA inhibitory activity of the compound. Furthermore, ADME study revealed useful information regarding compound’s drug-like properties. Therefore, the findings of the present investigation could aid in the development of more potent CA inhibitors, which could benefit the treatment of oxidative stress at high altitudes.     

Identification of novel natural inhibitors of Opisthorchis felineus cytochrome P450 using structure-based screening and molecular dynamic simulation

Opisthorchis felineus is the etiological agent of opisthorchiasis in humans. O. felineus cytochrome P450 (OfCYP450) is an important enzyme in the parasite xenobiotic metabolism. To identify the potential anti-opisthorchid compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,65,869) against the OfCYP450. The ligands were screened against OfCYP450 in four sequential docking modes that resulted in 361 ligands having better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 10 compounds were found to fit well with re-docking studies. After refinement by docking and drug-likeness analyses, four potential inhibitors (ZINC2358298, ZINC8790946, ZINC70707116, and ZINC85878789) were identified. These ligands with reference compounds (itraconazole and fluconazole) were further subjected to molecular dynamics simulation (MDS) and binding energy analyses to compare the dynamic structure of protein after ligand binding and the stability of the OfCYP450 and bound complexes. The binding energy analyses were also calculated. The results suggested that the compounds had a negative binding energy with ?259.41, ?110.09, ?188.25, ?163.30, ?202.10, and ?158.79 kJ mol^?1 for itraconazole, fluconazole, and compounds with IDs ZINC2358298, ZINC8790946, ZINC70707116, and ZINC85878789, respectively. These lead compounds displayed significant pharmacological and structural properties to be drug candidates. On the basis of MDS results and binding energy analyses, we concluded that ZINC8790946, ZINC70707116, and ZINC85878789 have excellent potential to inhibit OfCYP450.     

Identification of potential inhibitors of Fasciola gigantica thioredoxin1: computational screening,molecular dynamics simulation,and binding free energy studies

Fasciola gigantica is the causative organism of fascioliasis and is responsible for major economic losses in livestock production globally. F. gigantica thioredoxin1 (FgTrx1) is an important redox-active enzyme involved in maintaining the redox homeostasis in the cell. To identify a potential anti-fasciolid compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,67,740) against the FgTrx1 structure. The ligands were docked against FgTrx1 and 309 ligands were found to have better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 30 compounds were found to fit well for re-docking studies. After refinement by molecular docking and drug-likeness analysis, three potential inhibitors (ZINC15970091, ZINC9312362, and ZINC9312661) were identified. These three ligands were further subjected to molecular dynamics simulation (MDS) to compare the dynamics and stability of the protein structure after binding of the ligands. The binding free energy analyses were calculated to determine the intermolecular interactions. The results suggested that the two compounds had a binding free energy of –82.237, and –109.52 kJ.mol^?1 for compounds with IDs ZINC9312362 and ZINC9312661, respectively. These predicted compounds displayed considerable pharmacological and structural properties to be drug candidates. We concluded that these two compounds could be potential drug candidates to fight against F. gigantica parasites.     

Structure-based screening and molecular dynamics simulations offer novel natural compounds as potential inhibitors of Mycobacterium tuberculosis isocitrate lyase

Mycobacterium tuberculosis is the etiological agent of tuberculosis in humans and is responsible for more than two million deaths annually. M. tuberculosis isocitrate lyase (MtbICL) catalyzes the first step in the glyoxylate cycle, plays a pivotal role in the persistence of M. tuberculosis, which acts as a potential target for an anti-tubercular drug. To identify the potential anti-tuberculosis compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,67,748) against the MtbICL structure. The ligands were docked against MtbICL in three sequential docking modes that resulted in 340 ligands having better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 27 compounds were found to fit well with re-docking studies. After refinement by molecular docking and drug-likeness analyses, three potential inhibitors (ZINC1306071, ZINC2111081, and ZINC2134917) were identified. These three ligands and the reference compounds were further subjected to molecular dynamics simulation and binding energy analyses to compare the dynamic structure of protein after ligand binding and the stability of the MtbICL and bound complexes. The binding free energy analyses were calculated to validate and capture the intermolecular interactions. The results suggested that the three compounds had a negative binding energy with ?96.462, ?143.549, and ?122.526 kJ mol^?1 for compounds with IDs ZINC1306071, ZINC2111081, and ZINC2134917, respectively. These lead compounds displayed substantial pharmacological and structural properties to be drug candidates. We concluded that ZINC2111081 has a great potential to inhibit MtbICL and would add to the drug discovery process against tuberculosis.     

Structural and energetic understanding of novel natural inhibitors of Mycobacterium tuberculosis malate synthase

Persistent infection by Mycobacterium tuberculosis requires the glyoxylate shunt. This is a bypass to the tricarboxylic acid cycle in which isocitrate lyase (ICL) and malate synthase (MS) catalyze the net incorporation of carbon during mycobacterial growth on acetate or fatty acids as the primary carbon source. To identify a potential antitubercular compound, we performed a structure-based screening of natural compounds from the ZINC database (n = 1 67 740) against the M tuberculosis MS (MtbMS) structure. The ligands were screened against MtbMS, and 354 ligands were found to have better docking score. These compounds were assessed for Lipinski and absorption, distribution, metabolism, excretion, and toxicity prediction where 15 compounds were found to fit well for redocking studies. After refinement by molecular docking and drug-likeness analysis, four potential inhibitors (ZINC1483899, ZINC1754310, ZINC2269664, and ZINC15729522) were identified. These four ligands with phenyl-diketo acid were further subjected to molecular dynamics simulation to compare the dynamics and stability of the protein structure after ligand binding. The binding energy analysis was calculated to determine the intermolecular interactions. Our results suggested that the four compounds had a binding free energy of −201.96, −242.02, −187.03, and −169.02 kJ·mol⁻¹, for compounds with IDs ZINC1483899, ZINC1754310, ZINC2269664, and ZINC15729522, respectively. We concluded that two compounds (ZINC1483899 and ZINC1754310) displayed considerable structural and pharmacological properties and could be probable drug candidates to fight against M tuberculosis parasites.     

In silico approaches to identify novel myeloid cell leukemia-1 (Mcl-1) inhibitors for treatment of cancer

Myeloid cell leukemia-1 (Mcl-1) has been a validated and attractive target for cancer therapy. Over-expression of Mcl-1 in many cancers allows cancer cells to evade apoptosis and contributes to the resistance to current chemotherapeutics. Here, we identified new Mcl-1 inhibitors using a multi-step virtual screening approach. First, based on two different ligand-receptor complexes, 20 pharmacophore models were established by simultaneously using ‘Receptor-Ligand Pharmacophore Generation’ method and manual build feature method, and then carefully validated by a test database. Then, pharmacophore-based virtual screening (PB-VS) could be performed by using the 20 pharmacophore models. In addition, docking study was used to predict the possible binding poses of compounds, and the docking parameters were optimized before performing docking-based virtual screening (DB-VS). Moreover, a 3D QSAR model was established by applying the 55 aligned Mcl-1 inhibitors. The 55 inhibitors sharing the same scaffold were docked into the Mcl-1 active site before alignment, then the inhibitors with possible binding conformations were aligned. For the training set, the 3D QSAR model gave a correlation coefficient r² of 0.996; for the test set, the correlation coefficient r² was 0.812. Therefore, the developed 3D QSAR model was a good model, which could be applied for carrying out 3D QSAR-based virtual screening (QSARD-VS). After the above three virtual screening methods orderly filtering, 23 potential inhibitors with novel scaffolds were identified. Furthermore, we have discussed in detail the mapping results of two potent compounds onto pharmacophore models, 3D QSAR model, and the interactions between the compounds and active site residues.     

3D QSAR studies,pharmacophore modeling,and virtual screening of diarylpyrazole–benzenesulfonamide derivatives as a template to obtain new inhibitors,using human carbonic anhydrase II as a model protein

A 3D-QSAR modeling was performed on a series of diarylpyrazole-benzenesulfonamide derivatives acting as inhibitors of the metalloenzyme carbonic anhydrase (CA, EC 4.2.1.1). The compounds were collected from two datasets with the same scaffold, and utilized as a template for a new pharmacophore model to screen the ZINC database of commercially available derivatives. The datasets were divided into training, test, and validation sets. As the first step, comparative molecular field analysis (CoMFA), CoMFA region focusing and comparative molecular similarity indices analysis (CoMSIA) in parallel with docking studies were applied to a set of 41 human (h) CA II inhibitors. The validity and the prediction capacity of the resulting models were evaluated by leave-one-out (LOO) cross-validation approach. The reliability of the model for the prediction of possibly new CA inhibitors was also tested.     

Identification of potential leukocyte antigen-related protein (PTP-LAR) inhibitors through 3D QSAR pharmacophore-based virtual screening and molecular dynamics simulation

Abstract

Owing to its negative regulatory role in insulin signaling, protein tyrosine phosphatase of leukocyte antigen-related protein (PTP-LAR) was widely thought as a potential drug target for diabetes. Now, it was urgent to search for potential LAR inhibitors targeting diabetes. Initially, the pharmacophore models of LAR inhibitors were established with the application of the HypoGen module. The cost analysis, test set validation, as well as Fischer’s test was used to verify the efficiency of pharmacophore model. Then, the best pharmacophore model (Hypo-1-LAR) was applied for the virtual screening of the ZINC database. And 30 compounds met the Lipinski’s rule of five. Among them, 10 compounds with better binding affinity than the known LAR inhibitor (BDBM50296375) were discovered by docking studies. Finally, molecular dynamics simulations and post-analysis experiments (RMSD, RMSF, PCA, DCCM and RIN) were conducted to explore the effect of ligands (ZINC97018474 and Compound 1) on LAR and preliminary understand why ZINC97018474 had better inhibitory activity than Compound 1 (BDBM50296375).

Communicated by Ramaswamy H. Sarma     

Discovery of promising FtsZ inhibitors by E-pharmacophore, 3D-QSAR,molecular docking study,and molecular dynamics simulation

Filamentous temperature-sensitive protein Z (FtsZ), playing a key role in bacterial cell division, is regarded as a promising target for the design of antimicrobial agent. This study is looking for potential high-efficiency FtsZ inhibitors. Ligand-based pharmacophore and E-pharmacophore, virtual screening and molecular docking were used to detect promising FtsZ inhibitors, and molecular dynamics simulation was used to study the stability of protein-ligand complexes in this paper. Sixty-three inhibitors from published literatures with pIC₅₀ ranging from 2.483 to 5.678 were collected to develop ligand-based pharmacophore model. 4DXD bound with 9PC was selected to develop the E-pharmacophore model. The pharmacophore models validated by test set method and decoy set were employed for virtual screening to exclude inactive compounds against ZINC database. After molecular docking, ADME analysis, IFD docking and MM-GBSA, 8 hits were identified as potent FtsZ inhibitors. A 50?ns molecular dynamics simulation was implemented on the compounds to assess the stability between potent inhibitors and FtsZ. The results indicated that the candidate compounds had a high docking score and were strongly combined with FtsZ by forming hydrogen bonding interactions with key amino acid residues, and van der Waals forces and hydrophobic interactions had significant contribution to the stability of the binding. Molecular dynamics simulation results showed that the protein-ligand compounds performed well in both the stability and flexibility of the simulation process.     

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

Abstract

HCV NS5B polymerase 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 Hepatitis C Virus genotype 1 (HCV GT1). Hepatitis C virus genotype 4a (HCV GT4a) dominant in Egypt has paid less attention. Here, we describe our protocol of virtual screening in identification of novel potential potent inhibitors for HCV NS5B polymerase of GT4a using homology modeling, protein–ligand interaction fingerprint (PLIF), docking, pharmacophore, and 3D CoMFA quantitative structure activity relationship (QSAR). Firstly, a high-quality 3D model of HCV NS5B polymerase of GT4a was constructed using crystal structure of HCV NS5B polymerase of GT1 (PDB ID: 3hkw) as a template. Then, both the model and the template were simulated to compare conformational stability. PLIF was generated using five crystal structures of HCV NS5B (PDB ID: 4mia, 4mib, 4mk9, 4mka, and 4mkb), which revealed the most important residues and their interactions with the co-crystalized ligands. After that, a 3D pharmacophore model was developed from the generated PLIF data and then used as a screening filter for 17000328 drug-like zinc database compounds. 900 compounds passed the pharmacophore filter and entered the docking-based virtual screening stage. Finally, a 3D CoMFA QSAR was developed using 42 compounds as a training and 19 compounds as a test set. The 3D CoMFA QSAR was used to design and screen some potential inhibitors, these compounds were further evaluated by the docking stage. The highest ranked five hits from docking result (compounds (p1–p4) and compound q1) were selected for further analysis.

Communicated by Ramaswamy H. Sarma     

3D-QSAR pharmacophore-based virtual screening,molecular docking and molecular dynamics simulation toward identifying lead compounds for NS2B–NS3 protease inhibitors

NS2B–NS3 protease has been identified to serve as lead drug design target due to its significant role in West Nile viral (WNV) and dengue virus (DENV) reproduction and replication. There are currently no approved chemotherapeutic drugs and effective vaccines to inhibit DENV and WNV infections. In this work, 3D-QSAR pharmacophore model has been developed to discover potential inhibitory candidates. Validation through Fischer’s model and decoy test indicate that the developed 3D pharmacophore model is highly predictive for DENV inhibitors, which was then employed to screen ZINC chemical library to obtain reasonable hits. Following ADMET filtering, 15 hits were subjected to further filter through molecular docking and CoMFA modeling. Finally, top three hits were identified as lead compounds or potential inhibitory candidates with IC₅₀ values of ～0.4637?µM and fitness of ～57.73. It is implied from CoMFA modeling that substituents at the side site of benzotriazole such as a p-nitro group (e.g. biphenyl head) and a carbonyl (e.g. carboxylate function) at the side site of furan or amino group may improve bioactivity of ZINC85645245, respectively. Molecular dynamics simulations (MDS) were performed to discover new interactions and reinforce the binding modes from docking for the hits also. The QSAR and MDS results obtained from this work should be useful in determining structural requirements for inhibitor development as well as in designing more potential inhibitors for NS2B–NS3 protease.     

Docking and 3D QSAR study of thiourea analogs as potent inhibitors of influenza virus neuraminidase

Surflex-Dock was applied to study interactions between 30 thiourea analogs and neuraminidase (NA). The docking results showed that hydrogen bonding and electrostatic interactions were highly correlated with the activities of neuraminidase inhibitors (NIs), followed by hydrophobic and steric factors. Moreover, there was a strong correlation between the predicted binding affinity (total score) and experimental pIC₅₀ (correlation coefficient r = 0.870; P < 0.0001). A three dimensional holographic vector of atomic interaction field (3D-HoVAIF) was employed to construct a QSAR model. The r ², q ² and r ² _test values of the optimal QSAR model were 0.849, 0.724 and 0.689, respectively. From the QSAR model, it could be seen that electrostatic, hydrophobic and steric interactions were closely related to inhibitory activity, which was consistent with the docking results. Based on the docking and QSAR results, five new compounds with high predicted activities were designed.     

Screening of hepatitis C NS5B polymerase inhibitors containing benzothiadiazine core: a steered molecular dynamics

Hepatic C virus (HCV) is a global health problem, resulting in liver cirrhosis and inflammation that can develop to hepatocellular carcinoma and fatality. The NS5B polymerase of HCV plays an important role in viral RNA replication process, making it an attractive therapeutic target for design and development of anti-HCV drugs. To search new potent compounds against the HCV NS5B polymerase, the molecular docking and the steered molecular dynamics (SMD) simulation techniques were performed. The potential potent inhibitors of the NS5B polymerase were screened out from the ZINC database using structural similarity search and molecular docking technique. Five top-hit compounds (the ZINC compounds 49888724, 49054741, 49777239, 49793673, and 49780355) were then studied by the SMD simulations based on the hypothesis that a high rupture force relates to a high binding efficiency. The results demonstrated that the ZINC compound 49888724 had a greater maximum rupture force, reflecting a good binding strength and inhibitory potency than known inhibitors and the rest four ZINC compounds. Therefore, our finding indicated that the ZINC compound 49888724 is a potential candidate to be a novel NS5B inhibitor for further design. Besides, the van der Waals interaction could be considered as the main contribution for stabilizing the NS5B-ligand complex.