Structure-based virtual screening,molecular docking,ADMET and molecular simulations to develop benzoxaborole analogs as potential inhibitor against Leishmania donovani trypanothione reductase

Visceral leishmaniasis (VL) is the most fatal form of leishmaniasis and it affects 70 countries worldwide. Increasing drug resistant for antileishmanial drugs such as miltefosine, sodium stibogluconate and pentamidine has been reported in the VL endemic region. Amphotericin B has shown potential antileishmanial activity in different formulations but its cost of treatment and associated nephrotoxicity have limited its use by affected people living in the endemic zone. To control the VL infection in the affected countries, it is necessary to develop new antileishmanial compounds with high efficacy and negligible toxicity. Computer aided programs such as binding free energy estimation; ADMET prediction and molecular dynamics simulation can be used to investigate novel antileishmanial molecules in shorter duration. To develop antileishmanial lead molecule, we performed standard precision (SP) docking for 1160 benzoxaborole analogs along with reference inhibitors against trypanothione reductase of Leishmania parasite. Furthermore, extra precision (XP) docking, ADMET prediction, prime MM-GBSA was conducted over 115 ligands, showing better docking score than reference inhibitors to get potential antileishmanial compounds. Simultaneously, area under the curve (AUC) was estimated using ROC plot to validate the SP and XP docking protocol. Later on, two benzoxaborole analogs with best MM-GBSA ΔG-bind were subjected to molecular simulation and docking confirmation to ensure the ligand interaction with TR. The presented drug discovery based on computational study confirms that BOB27 can be used as a potential drug candidate and warrants further experimental investigation to fight against VL in endemic areas.     

Developing imidazole analogues as potential inhibitor for Leishmania donovani trypanothione reductase: virtual screening,molecular docking,dynamics and ADMET approach

Visceral leishmaniasis (VL) affects Indian subcontinent, African and South American continent, and it covers 70 countries worldwide. Visceral form of leishmaniasis is caused by Leishmania donovani in Indian subcontinent which is lethal if left untreated. Extensive resistance to antileishmanial drugs such as sodium stibogluconate, pentamidine and miltefosine and their decreased efficacy has been reported in the endemic region. Amphotericin B drug has shown good antileishmanial activity with significant toxicity, but its cost of treatment has limited the outreach of this treatment to affected people living in endemic zone. So, there is an urgent need to identify new antileishmanial drugs with excellent activity and minimal toxicity issues. Trypanothione reductase, a component of antioxidant system, is necessary for parasite growth and survival to raise infection. To develop potential inhibitor, we docked nine hundred and eighty-four 5-nitroimidazole analogues along with clomipramine which is a well-known inhibitor for TR. Total one hundred and forty-seven 5-nitroimidazole analogues with better docking score than clomipramine were chosen for ADMET and QikProp studies. Among these imidazole analogues, total twenty-four imidazole analogues and clomipramine were chosen on the basis of their ADMET, QikProp, and prime MM-GBSA study. Later on, two analogues with best MM-GBSA dG bind were undergone molecular dynamic simulation to ensure protein–ligand interactions. Using above approach, we confirm that ethyl 2-acetyl-5-[4-butyl-2-(3-hydroxypentyl)-5-nitro-1H-imidazol-1-yl]pent-2-enoate can be a drug candidate against L. donovani for the treatment of VL in the Indian subcontinent.     

Structural analysis and molecular docking of trypanocidal aryloxy-quinones in trypanothione and glutathione reductases: a comparison with biochemical data

A set of aryloxy-quinones, previously synthesized and evaluated against Trypanosoma cruzi epimastigotes cultures, were found more potent and selective than nifurtimox. One of the possible mechanisms of the trypanocidal activity of these quinones could be inhibition of trypanothione reductase (TR). Considering that glutathione reductase (GR) is the equivalent of TR in humans, biochemical, kinetic, and molecular docking studies in TR and GR were envisaged and compared with the trypanocidal and cytotoxic data of a set of aryloxy-quinones. Biochemical assays indicated that three naphthoquinones (Nq-h, Nq-g, and Nq-d) selectively inhibit TR and the TR kinetic analyses indicated that Nq-h inhibit TR in a noncompetitive mechanism. Molecular dockings were performed in TR and GR in the following three putative binding sites: the catalytic site, the dimer interface, and the nicotinamide adenine dinucleotide phosphate-binding site. In TR and GR, the aryloxy-quinones were found to exhibit high affinity for a site near it cognate-binding site in a place in which the noncompetitive kinetics could be justified. Taking as examples the three compounds with TR specificity (TRS) (Nq-h, Nq-g, and Nq-d), the presence of a network of contacts with the quinonic ring sustained by the triad of Lys62, Met400′, Ser464′ residues, seems to contribute hardly to the TRS. Compound Nq-b, a naphthoquinone with nitrophenoxy substituent, proved to be the best scaffold for the design of trypanocidal compounds with low toxicity. However, the compound displayed only a poor and non-selective effect toward TR indicating that TR inhibition is not the main reason for the antiparasitic activity of the aryloxy-quinones.     

Analysis of Coxiela burnetti dihydrofolate reductase via in silico docking with inhibitors and molecular dynamics simulation

Coxiella burnetii is a gram-negative bacterium able to infect several eukaryotic cells, mainly monocytes and macrophages. It is found widely in nature with ticks, birds, and mammals as major hosts. C. burnetii is also the biological warfare agent that causes Q fever, a disease that has no vaccine or proven chemotherapy available. Considering the current geopolitical context, this fact reinforces the need for discovering new treatments and molecular targets for drug design against C. burnetii. Among the main molecular targets against bacterial diseases reported, the enzyme dihydrofolate reductase (DHFR) has been investigated for several infectious diseases. In the present work, we applied molecular modeling techniques to evaluate the interactions of known DHFR inhibitors in the active sites of human and C. burnetii DHFR (HssDHFR and CbDHFR) in order to investigate their potential as selective inhibitors of CbDHFR. Results showed that most of the ligands studied compete for the binding site of the substrate more effectively than the reference drug trimethoprim. Also the most promising compounds were proposed as leads for the drug design of potential CbDHFR inhibitors.     

Insights about resveratrol analogs against trypanothione reductase of Leishmania braziliensis: Molecular modeling,computational docking and in vitro antileishmanial studies

In this work, we combined molecular modeling, computational docking and in vitro analysis to explore the antileishmanial effect of some resveratrol analogs (ResAn), focusing on their pro-oxidant effect. The molecular target was the trypanothione reductase of Leishmania braziliensis (LbTryR), an essential component of the antioxidant defenses in trypanosomatid parasites. Three-dimensional structures of LbTryR were modeled and molecular docking studies of ResAn1-5 compounds showed the following affinity: ResAn1?>?ResAn2?>?ResAn4?>?ResAn5?>?ResAn3. Positive correlation was observed between these compounds’ affinity to the LbTryR and the IC₅₀ values against Leishmania sp (ResAn1?<?ResAn2?<?ResAn4), which allows for TryR being considered an important target for them. As the compound ResAn1 showed the best antileishmanial activity, and docking studies showed its high affinity for NADP binding site (NS) of TryR, plus having been able to induce ROS production in L. braziliensis promastigotes treated, ResAn1 probably occupies NS interfering in the electron transfer processes responsible for the catalytic reaction. The in silico prediction of ADMET properties suggests that ResAn1 may be a promising drug candidate with properties to cross biological membranes and high gastrointestinal absorption, not violating Lipinski’s rules. Ultimately, the antileishmanial effect of ResAn can be associated with a pro-oxidant effect which, in turn, can be exploited as an antimicrobial agent.

Communicated by Ramaswamy H. Sarma     

Synthesis,molecular docking,binding free energy calculation and molecular dynamics simulation studies of benzothiazol-2-ylcarbamodithioates as Staphylococcus aureus MurD inhibitors

Abstract

A new series of benzothiazol-2-ylcarbamodithioate functional compounds 5a-f has been designed, synthesized and characterized by spectral data. These compounds were screened for their in vitro antibacterial activity against strains of Staphylococcus aureus (NCIM 5021, NCIM 5022 and methicillin-resistant isolate 43300), Bacillus subtilis (NCIM 2545), Escherichia coli (NCIM 2567), Klebsiella pneumoniae (NCIM 2706) and Psudomonas aeruginosa (NCIM 2036). Compounds 5a and 5d exhibited significant activity against all the tested bacterial strains. Specifically, compounds 5a and 5d showed potent activity against K. pneumoniae (NCIM 2706), while compound 5a also displayed potent activity against S. aureus (NCIM 5021). Compound 5d showed minimum IC₅₀ value of 13.37?μM against S. aureus MurD enzyme. Further, the binding interactions of compounds 5a-f in the catalytic pocket have been investigated using the extra-precision molecular docking and binding free energy calculation by MM-GBSA approach. A 30?ns molecular dynamics simulation of 5d/modeled S. aureus MurD enzyme was performed to determine the stability of the predicted binding conformation.     

Pharmacophore generation,atom-based 3D-QSAR,molecular docking and molecular dynamics simulation studies on benzamide analogues as FtsZ inhibitors

FtsZ is an appealing target for the design of antimicrobial agent that can be used to defeat the multidrug-resistant bacterial pathogens. Pharmacophore modelling, molecular docking and molecular dynamics (MD) simulation studies were performed on a series of three-substituted benzamide derivatives. In the present study a five-featured pharmacophore model with one hydrogen bond acceptors, one hydrogen bond donors, one hydrophobic and two aromatic rings was developed using 97 molecules having MIC values ranging from .07 to 957 μM. A statistically significant 3D-QSAR model was obtained using this pharmacophore hypothesis with a good correlation coefficient (R² = .8319), cross validated coefficient (Q² = .6213) and a high Fisher ratio (F = 103.9) with three component PLS factor. A good correlation between experimental and predicted activity of the training (R² = .83) and test set (R² = .67) molecules were displayed by ADHRR.1682 model. The generated model was further validated by enrichment studies using the decoy test and MAE-based criteria to measure the efficiency of the model. The docking studies of all selected inhibitors in the active site of FtsZ protein showed crucial hydrogen bond interactions with Val 207, Asn 263, Leu 209, Gly 205 and Asn-299 residues. The binding free energies of these inhibitors were calculated by the molecular mechanics/generalized born surface area VSGB 2.0 method. Finally, a 15 ns MD simulation was done to confirm the stability of the 4DXD–ligand complex. On a wider scope, the prospect of present work provides insight in designing molecules with better selective FtsZ inhibitory potential.     

Molecular dynamics simulation,binding free energy calculation and molecular docking of human D-amino acid oxidase (DAAO) with its inhibitors

Schizophrenia is a mental illness; most affected people live in developing countries, and neither appropriate treatment nor commercial drugs are currently available. One possibility is to inhibit human-d-amino acid oxidase (h-DAAO). In this study, molecular dynamic simulations of the monomer, dimer and tetramer forms of h-DAAO complexed with the inhibitor 3-hydroxyquinolin-2(1H)-one(2) were performed. Seven residues, Leu51, Gln53, Leu215, Tyr228, Ile230, Arg283 and Gly313, were identified as essential for interacting with the inhibitor. Molecular docking of h-DAAO with pyrrole, quinoline and kojic acid derivatives, representing 69 known or potential h-DAAO inhibitors, was also performed. The results indicated that the activity of the inhibitor can be improved by modifying the compounds to have a substituent group capable of interacting with the side chain of Tyr228. Van der Waals interactions of the inhibitor with the hydrophobic pocket of h-DAAO and electrostatic interactions or H-bonds with Arg283 and Gly313 were important elements in determining the efficiency of the inhibitor. These results provide information on the interaction between h-DAAO and its inhibitors at the molecular level and can aid in the design of novel inhibitors against h-DAAO for new drug development in the treatment of schizophrenia.     

Determination of potential selective inhibitors for ROCKI and ROCKII isoforms with molecular modeling techniques: structure based docking,ADMET and molecular dynamics simulation

Rho-associated protein kinases (ROCKs) are a member of the serine/threonine protein kinase family and potential therapeutic target for various diseases. This enzyme has two isoforms, Rho-associated protein kinase I (ROCKI) and Rho-associated protein kinase II (ROCKII). They share an overall 65% homology in all amino acid sequence and 92% homology in kinase domains. Since, the kinase domains of ROCKI and ROCKII are highly conserved and similar, the discovery and design of isoform-selective inhibitors are more challenging. Thus, most currently available agents that is against ROCKs exhibit low selectivity and severe side effects. Therefore, this study aimed to elucidate the interaction of compounds that indicated high potential in experimental studies against ROCKI and ROCKII enzymes in the molecular level with molecular modeling techniques. Firstly, we determined the interaction property of catalytic sites of the ROCKs by analyzing with molecular docking. Based on these results, the best ligands (50 compounds) corresponding to experimental studies were selected, and then absorption, distribution, metabolism and excretion – toxicity (ADMET) analysis of these compounds were implemented. According to these study results, the compound 40 for ROCKI and the compound 50 for ROCKII were identified as selective and highly potent inhibitors. And finally, molecular dynamics (MD) simulations were performed for the stability of ROCKs with identified compounds. In the light of this study, it will be possible to treat diseases that ROCKs have a role by developing more effective and specific ROCK inhibitors.

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.     

Quercetin derivatives as non-nucleoside inhibitors for dengue polymerase: molecular docking,molecular dynamics simulation,and binding free energy calculation

Dengue is an important public health problem in tropical and subtropical regions of the world. Neither vaccine nor an antiviral medication is available to treat dengue. This insists the need of drug discovery for dengue. In order to find a potent lead molecule, RNA-dependent RNA polymerase which is essential for dengue viral replication is chosen as a drug target. As Quercetin showed antiviral activity against several viruses, quercetin derivatives developed by combinatorial library synthesis and mined from PubChem databases were screened for a potent anti-dengue viral agent. Our study predicted Quercetin 3-(6″-(E)-p-coumaroylsophoroside)-7-rhamnoside as a dengue polymerase inhibitor. The results were validated by molecular dynamics simulation studies which reveal water bridges and hydrogen bonds as major contributors for the stability of the polymerase-lead complex. Interactions formed by this compound with residues Trp795, Arg792 and Glu351 are found to be essential for the stability of the polymerase-lead complex. Our study demonstrates Quercetin 3-(6″-(E)-p-coumaroylsophoroside)-7-rhamnoside as a potent non-nucleoside inhibitor for dengue polymerase.     

Novel urushiol derivatives as HDAC8 inhibitors: rational design,virtual screening,molecular docking and molecular dynamics studies

Three series of novel urushiol derivatives were designed by introducing a hydroxamic acid moiety into the tail of an alkyl side chain and substituents with differing electronic properties or steric bulk onto the benzene ring and alkyl side chain. The compounds’ binding affinity toward HDAC8 was screened by Glide docking. The highest-scoring compounds were processed further with molecular docking, MD simulations, and binding free energy studies to analyze the binding modes and mechanisms. Ten compounds had Glide scores of ?8.2 to ?10.2, which revealed that introducing hydroxy, carbonyl, amino, or methyl ether groups into the alkyl side chain or addition of –F, –Cl, sulfonamide, benzamido, amino, or hydroxy substituents on the benzene ring could significantly increase binding affinity. Molecular docking studies revealed that zinc ion coordination, hydrogen bonding, and hydrophobic interactions contributed to the high calculated binding affinities of these compounds toward HDAC8. MD simulations and binding free energy studies showed that all complexes possessed good stability, as characterized by low RMSDs, low RMSFs of residues, moderate hydrogen bonding and zinc ion coordination and low values of binding free energies. Hie147, Tyr121, Phe175, Hip110, Phe119, Tyr273, Lys21, Gly118, Gln230, Leu122, Gly269, and Gly107 contributed favorably to the binding; and Van der Waals and electrostatic interactions provided major contributions to the stability of these complexes. These results show the potential of urushiol derivatives as HDAC8 binding lead compounds, which have great therapeutic potential in the treatment of various malignancies, neurological disorders, and human parasitic diseases.     

New potential inhibitors of mTOR: a computational investigation integrating molecular docking,virtual screening and molecular dynamics simulation

The mTOR (mammalian or mechanistic Target Of Rapamycin), a complex metabolic pathway that involves multiple steps and regulators, is a major human metabolic pathway responsible for cell growth control in response to multiple factors and that is dysregulated in various types of cancer. The classical inhibition of the mTOR pathway is performed by rapamycin and its analogs (rapalogs). Considering that rapamycin binds to an allosteric site and performs a crucial role in the inhibition of the mTOR complex without causing the deleterious side effects common to ATP-competitive inhibitors, we employ ligand-based drug design strategies, such as virtual screening methodology, computational determination of ADME/Tox properties of selected molecules, and molecular dynamics in order to select molecules with the potential to become non-ATP-competitive inhibitors of the mTOR enzymatic complex. Our findings suggest five novel potential mTOR inhibitors, with similar or better properties than the classic inhibitor complex, rapamycin.     

Atomic insight into designed carbamate-based derivatives as acetylcholine esterase (AChE) inhibitors: a computational study by multiple molecular docking and molecular dynamics simulation

Over 100 variants have been designed and studied, using multiple docking methods such as Autodock Vina, ArgusLab, Molegro Virtual Docker, and Hex-Cuda, to study the effect of alteration in the structure of carbamate-based acetylcholyne esterase (AChE) inhibitors. Sixteen selected systems were then subjected to 14 ns molecular dynamics (MD) simulations. Results from all the docking methods are in agreement. Variants that involved biphenyl substituents possess the most negative binding energies in the ?37.64 to ?39.31 kJ mol^?1 range due to their π–π interactions with AChE aromatic residues. The root mean square deviation values showed that all of these components achieved equilibration after 6 ns. Gyration radius (R_g) and solvent accessibility surface area were calculated to further investigate the AChE conformational changes in the presence of these components. MD simulation results suggested that these components might interact with AChE, possibly with no major changes in AChE secondary and tertiary structures.     

Discover potential inhibitors for PFKFB3 using 3D-QSAR,virtual screening,molecular docking and molecular dynamics simulation

Abstract

The 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase-3 (PFKFB3) is a master regulator of glycolysis in cancer cells by synthesizing fructose-2,6-bisphosphate (F-2,6-BP), a potent allosteric activator of phosphofructokinase-1 (PFK-1), which is a rate-limiting enzyme of glycolysis. PFKFB3 is an attractive target for cancer treatment. It is valuable to discover promising inhibitors by using 3D-QSAR pharmacophore modeling, virtual screening, molecular docking and molecular dynamics simulation. Twenty molecules with known activity were used to build 3D-QSAR pharmacophore models. The best pharmacophore model was ADHR called Hypo1, which had the highest correlation value of 0.98 and the lowest RMSD of 0.82. Then, the Hypo1 was validated by cost value method, test set method and decoy set validation method. Next, the Hypo1 combined with Lipinski's rule of five and ADMET properties were employed to screen databases including Asinex and Specs, total of 1,048,159 molecules. The hits retrieved from screening were docked into protein by different procedures including HTVS, SP and XP. Finally, nine molecules were picked out as potential PFKFB3 inhibitors. The stability of PFKFB3-lead complexes was verified by 40?ns molecular dynamics simulation. The binding free energy and the energy contribution of per residue to the binding energy were calculated by MM-PBSA based on molecular dynamics simulation.     

Exploration of potential RSK2 inhibitors by pharmacophore modelling,structure-based 3D-QSAR,molecular docking study and molecular dynamics simulation

Abstract

The p90 ribosomal s6 kinase 2 (RSK2) is a promising target because of its over expression and activation in human cancer cells and tissues. Over the last few years, significant efforts have been made in order to develop RSK2 inhibitors to treat myeloma, prostatic cancer, skin cancer and etc., but with limited success so far. In this paper, pharmacophore modelling, molecular docking study and molecular dynamics (MD) simulation have been performed to explore the novel inhibitors of RSK2. Pharmacophore models were developed by 95 molecules having pIC₅₀ ranging from 4.577 to 9.000. The pharmacophore model includes one hydrogen bond acceptor (A), one hydrogen bond donor (D), one hydrophobic feature (H) and one aromatic ring (R). It is the best pharmacophore hypothesis that has the highest correlation coefficient (R² = 0.91) and cross validation coefficient (Q² = 0.71) at 5 component PLS factor. It was evaluated using enrichment analysis and the best model was used for virtual screening. The constraints used in this study were docking score, ADME properties, binding free energy estimates and IFD Score to screen the database. Ultimately, 12 hits were identified as potent and novel RSK2 inhibitors. A 15 ns molecular dynamics (MD) simulation was further employed to validate the reliability of the docking results.     

Probing ligand binding modes of human cytochrome P450 2J2 by homology modeling, molecular dynamics simulation, and flexible molecular docking

Cytochrome P450 (P450) 2J2 catalyzes epoxidation of arachidonic acid to eicosatrienoic acids, which are related to a variety of diseases such as coronary artery disease, hypertension, and carcinogenesis. Recent experimental data also suggest that P450 2J2 could be a novel biomarker and a potential target for cancer therapy. However, the active site topology and substrate specificity of this enzyme remain unclear. In this study, a three-dimensional model of human P450 2J2 was first constructed on the basis of the crystal structure of human P450 2C9 in complex with a substrate using homology modeling method, and refined by molecular dynamics simulation. Flexible docking approaches were then employed to dock four ligands into the active site of P450 2J2 in order to probe the ligand-binding modes. By analyzing the results, active site architecture and certain key residues responsible for substrate specificity were identified on the enzyme, which might be very helpful for understanding the enzyme's biological role and providing insights for designing novel inhibitors of P450 2J2.     

A combination of pharmacophore modeling,atom-based 3D-QSAR,molecular docking and molecular dynamics simulation studies on PDE4 enzyme inhibitors

Phosphodiesterases 4 enzyme is an attractive target for the design of anti-inflammatory and bronchodilator agents. In the present study, pharmacophore and atom-based 3D-QSAR studies were carried out for pyrazolopyridine and quinoline derivatives using Schrödinger suite 2014-3. A four-point pharmacophore model was developed using 74 molecules having pIC₅₀ ranging from 10.1 to 4.5. The best four feature model consists of one hydrogen bond acceptor, two aromatic rings, and one hydrophobic group. The pharmacophore hypothesis yielded a statistically significant 3D-QSAR model, with a high correlation coefficient (R^2?=?.9949), cross validation coefficient (Q^2?=?.7291), and Pearson-r (.9107) at six component partial least square factor. The external validation indicated that our QSAR model possessed high predictive power with R² value of .88. The generated model was further validated by enrichment studies using the decoy test. Molecular docking, free energy calculation, and molecular dynamics (MD) simulation studies have been performed to explore the putative binding modes of these ligands. A 10-ns MD simulation confirmed the docking results of both stability of the 1XMU–ligand complex and the presumed active conformation. Outcomes of the present study provide insight in designing novel molecules with better PDE4 inhibitory activity.     

Exploring the effect of inhibitor AKB-9778 on VE-PTP by molecular docking and molecular dynamics simulation

Diabetic macular edema, also known as diabetic eye disease, is mainly caused by the overexpression of vascular endothelial protein tyrosine phosphatase (VE-PTP) at hypoxia/ischemic. AKB-9778 is a known VE-PTP inhibitor that can effectively interact with the active site of VE-PTP to inhibit the activity of VE-PTP. However, the binding pattern of VE-PTP with AKB-9778 and the dynamic implications of AKB-9778 on VE-PTP system at the molecular level are poorly understood. Through molecular docking, it was found that the AKB-9778 was docked well in the binding pocket of VE-PTP by the interactions of hydrogen bond and Van der Waals. Furthermore, after molecular dynamic simulations on VE-PTP system and VE-PTP ^AKB-9778 system, a series of postdynamic analyses found that the flexibility and conformation of the active site undergone an obvious transition after VE-PTP binding with AKB-9778. Moreover, by constructing the RIN, it was found that the different interactions in the active site were the detailed reasons for the conformational differences between these two systems. Thus, the finding here might provide a deeper understanding of AKB-9778 as VE-PTP Inhibitor.