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.     

Febrifugine analogues as Leishmania donovani trypanothione reductase inhibitors: binding energy analysis assisted by molecular docking,ADMET and molecular dynamics simulation

Visceral leishmaniasis affects people from 70 countries worldwide, mostly from Indian, African and south American continent. The increasing resistance to antimonial, miltefosine and frequent toxicity of amphotericin B drives an urgent need to develop an antileishmanial drug with excellent efficacy and safety profile. In this study we have docked series of febrifugine analogues (n = 8813) against trypanothione reductase in three sequential docking modes. Extra precision docking resulted into 108 ligands showing better docking score as compared to two reference ligand. Furthermore, 108 febrifugine analogues and reference inhibitor clomipramine were subjected to ADMET, QikProp and molecular mechanics, the generalized born model and solvent accessibility study to ensure the toxicity caused by compounds and binding-free energy, respectively. Two best ligands (FFG7 and FFG2) qualifying above screening parameters were further subjected to molecular dynamics simulation. Conducting these studies, here we confirmed that 6-chloro-3-[3-(3-hydroxy-2-piperidyl)-2-oxo-propyl]-7-(4-pyridyl) quinazolin-4-one can be potential drug candidate to fight against Leishmania donovani parasites.     

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     

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.     

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.     

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.     

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.     

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     

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.     

Structure based virtual screening, 3D-QSAR,molecular dynamics and ADMET studies for selection of natural inhibitors against structural and non-structural targets of Chikungunya

The transmission of mosquito-borne Chikungunya virus (CHIKV) has large epidemics worldwide. Till date, there are neither anti-viral drugs nor vaccines available for the treatment of Chikungunya. Accumulated evidences suggest that some natural compounds i.e., Epigallocatechin gallate, Harringtonine, Apigenin, Chrysin, Silybin, etc. have the capability to inhibit CHIKV replication in vitro. Natural compounds are known to possess less or no side effects. Therefore, natural compound in its purified or crude extracts form could be the preeminent and safe mode of therapies for Chikungunya. Wet lab screening and identification of natural compounds against Chikungunya targets is a time consuming and expensive exercise. In the present study, we used in silico techniques like receptor-ligand docking, Molecular dynamic (MD), Three Dimensional Quantitative Structure Activity Relation (3D-QSAR) and ADME properties to screen out potential compounds. Aim of the study is to identify potential lead/s from natural sources using in silico techniques that can be developed as a drug like molecule against Chikungunya infection and replication. Three softwares were used for molecular docking studies. Potential ligands selected by docking studies were subsequently subjected 3D-QSAR studies to predict biological activity. Based on docking scores and pIC₅₀ value, potential anti-Chikungunya compounds were identified. Best docked receptor-ligands were also subjected to MD for more accurate estimation. Lipinski’s rule and ADME studies of the identified compounds were also studied to assess their drug likeness properties. Results of in silico findings, led to identification of few best fit compounds of natural origin against targets of Chikungunya virus which may lead to discovery of new drugs for Chikungunya.

Communicated by Ramaswamy H. Sarma     

Homology modeling,docking and structure-based virtual screening for new inhibitor identification of Klebsiella pneumoniae heptosyltransferase-III

Abstract

Klebsiella pneumoniae (K. pneumoniae) is a Gram-negative opportunistic pathogen commonly associated with hospital-acquired infections that are often resistant even to antibiotics. Heptosyltransferase (HEP) belongs to the family of glycosyltransferase-B (GT-B) and plays an important in the synthesis of lipopolysaccharides (LPS) essential for the formation of bacterial cell membrane. HEP-III participates in the transfer of heptose sugar to the outer surface of bacteria to synthesize LPS. LPS truncation increases the bacterial sensitivity to hydrophobic antibiotics and detergents, making the HEP as a novel drug target. In the present study, we report the 3D homology model of K. pneumoniae HEP-III and its structure validation. Active site was identified based on similarities with known structures using Dali server, and structure-based pharmacophore model was developed for the active site substrate ADP. The generated pharmacophore model was used as a 3D search query for virtual screening of the ASINEX database. The hit compounds were further filtered based on fit value, molecular docking, docking scores, molecular dynamics (MD) simulations of HEP-III complexed with hit molecules, followed by binding free energy calculations using Molecular Mechanics-Poisson–Boltzmann Surface Area (MM-PBSA). The insights obtained in this work provide the rationale for design of novel inhibitors targeting K. pneumoniae HEP-III and the mechanistic aspects of their binding.

Communicated by Ramaswamy H. Sarma     

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.     

Benzoxazinone-thiosemicarbazones as antidiabetic leads via aldose reductase inhibition: Synthesis,biological screening and molecular docking study

Aldose reductase is an important enzyme in the polyol pathway, where glucose is converted to fructose, and sorbitol is released. Aldose reductase activity increases in diabetes as the glucose levels increase, resulting in increased sorbitol production. Sorbitol, being less cell permeable tends to accumulate in tissues such as eye lenses, peripheral nerves and glomerulus that are not insulin sensitive. This excessive build-up of sorbitol is responsible for diabetes associated complications such as retinopathy and neuropathy. In continuation of our interest to design and discover potent inhibitors of aldo-keto reductases (AKRs; aldehyde reductase ALR1 or AKR1A, and aldose reductase ALR2 or AKR1B), herein we designed and investigated a series of new benzoxazinone-thiosemicarbazones (3a-r) as ALR2 and ALR1 inhibitors. Most compounds exhibited excellent inhibitory activities with IC₅₀ values in lower micro-molar range. Compounds 3b and 3l were found to be most active ALR2 inhibitors with IC₅₀ values of 0.52 ± 0.04 and 0.19 ± 0.03 μM, respectively, both compounds were more effective inhibitors as compared to the standard ALR2 inhibitor (sorbinil, with IC₅₀ value of 3.14 ± 0.02 μM).     

Virtual screening,docking, and dynamics of potential new inhibitors of dihydrofolate reductase from Yersinia pestis

In the present work, we propose to design drugs that target the enzyme dihydrofolate redutase (DHFR) as a means of a novel drug therapy against plague. Potential inhibitors of DHFR from Yersinia pestis (YpDHFR) were selected by virtual screening and subjected to docking, molecular dynamics (MD) simulations, and Poisson–Boltzmann surface area method, in order to evaluate their interactions in the active sites of YpDHFR and human DHFR (HssDHFR). The results suggested selectivity for three compounds that were further used to propose the structures of six new potential selective inhibitors for YpDHFR.     

Shape-based virtual screening,docking, and molecular dynamics simulations to identify Mtb-ASADH inhibitors

Aspartate β-semialdehyde dehydrogenase (ASADH) is a key enzyme for the biosynthesis of essential amino acids and several important metabolites in microbes. Inhibition of ASADH enzyme is a promising drug target strategy against Mycobacterium tuberculosis (Mtb). In this work, in silico approach was used to identify potent inhibitors of Mtb-ASADH. Aspartyl β-difluorophosphonate (β-AFP), a known lead compound, was used to understand the molecular recognition interactions (using molecular docking and molecular dynamics analysis). This analysis helped in validating the computational protocol and established the participation of Arg99, Glu224, Cys130, Arg249, and His256 amino acids as the key amino acids in stabilizing ligand–enzyme interactions for effective binding, an essential feature is H-bonding interactions with the two arginyl residues at the two ends of the ligand. Best binding conformation of β-AFP was selected as a template for shape-based virtual screening (ZINC and NCI databases) to identify compounds that competitively inhibit the Mtb-ASADH. The top rank hits were further subjected to ADME and toxicity filters. Final filter was based on molecular docking analysis. Each screened molecule carries the characteristics of the highly electronegative groups on both sides separated by an average distance of 6?Å. Finally, the best predicted 20 compounds exhibited minimum three H-bonding interactions with Arg99 and Arg249. These identified hits can be further used for designing the more potent inhibitors against ASADH family. MD simulations were also performed on two selected compounds (NSC4862 and ZINC02534243) for further validation. During the MD simulations, both compounds showed same H-bonding interactions and remained bound to key active residues of Mtb-ASADH.     

Multiple e-pharmacophore modelling pooled with high-throughput virtual screening,docking and molecular dynamics simulations to discover potential inhibitors of Plasmodium falciparum lactate dehydrogenase (PfLDH)

Development of new antimalarial drugs continues to be of huge importance because of the resistance of malarial parasite towards currently used drugs. Due to the reliance of parasite on glycolysis for energy generation, glycolytic enzymes have played important role as potential targets for the development of new drugs. Plasmodium falciparum lactate dehydrogenase (PfLDH) is a key enzyme for energy generation of malarial parasites and is considered to be a potential antimalarial target. Presently, there are nearly 15 crystal structures bound with inhibitors and substrate that are available in the protein data bank (PDB). In the present work, we attempted to consider multiple crystal structures with bound inhibitors showing affinity in the range of 1.4 × 10²–1.3 × 10⁶ nM efficacy and optimized the pharmacophore based on the energy involved in binding termed as e-pharmacophore mapping. A high throughput virtual screening (HTVS) combined with molecular docking, ADME predictions and molecular dynamics simulation led to the identification of 20 potential compounds which could be further developed as novel inhibitors for PfLDH.     

Capecitabine as a minor groove binder of DNA: molecular docking,molecular dynamics,and multi-spectroscopic studies

The interaction mechanism and binding mode of capecitabine with ctDNA was extensively investigated using docking and molecular dynamics simulations, fluorescence and circular dichroism (CD) spectroscopy, DNA thermal denaturation studies, and viscosity measurements. The possible binding mode and acting forces on the combination between capecitabine and DNA had been predicted through molecular simulation. Results indicated that capecitabine could relatively locate stably in the G-C base-pairs-rich DNA minor groove by hydrogen bond and several weaker nonbonding forces. Fluorescence spectroscopy and fluorescence lifetime measurements confirmed that the quenching was static caused by ground state complex formation. This phenomenon indicated the formation of a complex between capecitabine and ctDNA. Fluorescence data showed that the binding constants of the complex were approximately 2 × 10⁴ M^?1. Calculated thermodynamic parameters suggested that hydrogen bond was the main force during binding, which were consistent with theoretical results. Moreover, CD spectroscopy, DNA melting studies, and viscosity measurements corroborated a groove binding mode of capecitabine with ctDNA. This binding had no effect on B-DNA conformation.     

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.