Molecular docking and molecular dynamics simulation studies of Trypanosoma cruzi triosephosphate isomerase inhibitors. Insights into the inhibition mechanism and selectivity

A comprehensive review on potential candidates for the treatment of chagas disease

Twenty different infectious disorders induced by bacteria, viruses, and parasites are categorized as neglected tropical diseases (NTDs) by WHO. The severity of chagas disease remains a major concern in endemic areas and an emerging public health hazard in nonendemic countries. Trypanosoma cruzi, the etiological agent of this NTD, is mostly transmitted by triatomine vectors and comprises a range of epidemiologically significant variants. Current chemotherapeutics are obsolete, and one of the primary reasons for treatment cessation is their poor safety and effectiveness. Due to the aforementioned challenges, researchers are now focusing on discovering alternative novel safe, and economically reachable therapies for the treatment of trypanosomiasis. Certain target-based drugs that target specific biochemical processes of the causative parasites have been described as potential antichagasic agents that possesses various types of heterocyclic scaffolds. These flexible molecules have a wide range of biological actions, and various synthesized compounds with strong activity have been documented. This review aims to discuss the available literature on synthetic anti-T. cruzi drugs that will give a food for thought to medicinal chemists thriving to design and develop such drugs. Furthermore, some of the studies discussed herein are concerned with the potential of novel drugs to block new viable sites in T. cruzi.

Triose Phosphate Isomerase Structure-Based Virtual Screening and In Vitro Biological Activity of Natural Products as Leishmania mexicana Inhibitors

Cutaneous leishmaniasis (CL) is a public health problem affecting more than 98 countries worldwide. No vaccine is available to prevent the disease, and available medical treatments cause serious side effects. Additionally, treatment failure and parasite resistance have made the development of new drugs against CL necessary. In this work, a virtual screening of natural products from the BIOFACQUIM and Selleckchem databases was performed using the method of molecular docking at the triosephosphate isomerase (TIM) enzyme interface of Leishmania mexicana (L. mexicana). Finally, the in vitro leishmanicidal activity of selected compounds against two strains of L. mexicana, their cytotoxicity, and selectivity index were determined. The top ten compounds were obtained based on the docking results. Four were selected for further in silico analysis. The ADME-Tox analysis of the selected compounds predicted favorable physicochemical and toxicological properties. Among these four compounds, S-8 (IC50 = 55 µM) demonstrated a two-fold higher activity against the promastigote of both L. mexicana strains than the reference drug glucantime (IC50 = 133 µM). This finding encourages the screening of natural products as new anti-leishmania agents.

Ligand-Based Virtual Screening and Molecular Docking of Benzimidazoles as Potential Inhibitors of Triosephosphate Isomerase Identified New Trypanocidal Agents

Trypanosoma cruzi (T. cruzi) is a parasite that affects humans and other mammals. T. cruzi depends on glycolysis as a source of adenosine triphosphate (ATP) supply, and triosephosphate isomerase (TIM) plays a key role in this metabolic pathway. This enzyme is an attractive target for the design of new trypanocidal drugs. In this study, a ligand-based virtual screening (LBVS) from the ZINC15 database using benzimidazole as a scaffold was accomplished. Later, a molecular docking on the interface of T. cruzi TIM (TcTIM) was performed and the compounds were grouped by interaction profiles. Subsequently, a selection of compounds was made based on cost and availability for in vitro evaluation against blood trypomastigotes. Finally, the compounds were analyzed by molecular dynamics simulation, and physicochemical and pharmacokinetic properties were determined using SwissADME software. A total of 1604 molecules were obtained as potential TcTIM inhibitors. BP2 and BP5 showed trypanocidal activity with half-maximal lytic concentration (LC₅₀) values of 155.86 and 226.30 µM, respectively. Molecular docking and molecular dynamics simulation analyzes showed a favorable docking score of BP5 compound on TcTIM. Additionally, BP5 showed a low docking score (−5.9 Kcal/mol) on human TIM compared to the control ligand (−7.2 Kcal/mol). Both compounds BP2 and BP5 showed good physicochemical and pharmacokinetic properties as new anti-T. cruzi agents.

Synergistic effect of compounds directed to triosephosphate isomerase, a combination to develop drug against trichomoniasis

The development of new drugs is continuous in the world; currently, saving resources (both economic ones and time) and preventing secondary effects have become a necessity for drug developers. Trichomoniasis is the most common nonviral sexually transmitted infection affecting more than 270 million people around the world. In our research group, we focussed on developing a selective and more effective drug against Trichomonas vaginalis, and we previously reported on a compound, called A4, which had a trichomonacidal effect. Later, we determined another compound, called D4, which also had a trichomonacidal effect together with favorable toxicity results. Both A4 and D4 are directed at the enzyme triosephosphate isomerase. Thus, we made combinations between the two compounds, in which we determined a synergistic effect against T. vaginalis, determining the IC₅₀ and the toxicity of the best relationship to obtain the trichomonacidal effect. With these results, we can propose a combination of compounds that represents a promising alternative for the development of a new therapeutic strategy against trichomoniasis.

Advances in Phenazines over the Past Decade: Review of Their Pharmacological Activities, Mechanisms of Action, Biosynthetic Pathways and Synthetic Strategies

Phenazines are a large group of nitrogen-containing heterocycles, providing diverse chemical structures and various biological activities. Natural phenazines are mainly isolated from marine and terrestrial microorganisms. So far, more than 100 different natural compounds and over 6000 synthetic derivatives have been found and investigated. Many phenazines show great pharmacological activity in various fields, such as antimicrobial, antiparasitic, neuroprotective, insecticidal, anti-inflammatory and anticancer activity. Researchers continued to investigate these compounds and hope to develop them as medicines. Cimmino et al. published a significant review about anticancer activity of phenazines, containing articles from 2000 to 2011. Here, we mainly summarize articles from 2012 to 2021. According to sources of compounds, phenazines were categorized into natural phenazines and synthetic phenazine derivatives in this review. Their pharmacological activities, mechanisms of action, biosynthetic pathways and synthetic strategies were summarized. These may provide guidance for the investigation on phenazines in the future.

Recent Advances in the Development of Triose Phosphate Isomerase Inhibitors as Antiprotozoal Agents

BACKGROUND

Parasitic diseases caused by protozoa such as Chagas disease, leishmaniasis, malaria, African trypanosomiasis, amebiasis, trichomoniasis, and giardiasis are considered serious public health problems in developing countries. Drug-resistance among parasites justifies the search for new therapeutic drugs and the identification of new targets becomes a valuable approach. In this scenario, glycolysis pathway which consists of the conversion of glucose into pyruvate plays an important role in the protozoa energy supply and it is therefore considered as a promising target. In this pathway, triose phosphate isomerase (TIM) plays an essential role in efficient energy production. Furthermore, protozoa TIM show structural differences with human enzyme counterparts suggesting the possibility of obtaining selective inhibitors. Therefore, TIM is considered a valid approach to develop new antiprotozoal agents, inhibiting the glycolysis in the parasite.

OBJECTIVE

In this review, we discuss the drug design strategies, structure-activity relationship, and binding modes of outstanding TIM inhibitors against Trypanosoma cruzi, Trypanosoma brucei, Plasmodium falciparum, Giardia lamblia, Leishmania mexicana, Trichomonas vaginalis, and Entamoeba histolytica.

RESULTS

TIM inhibitors showed mainly aromatic systems and symmetrical structure, where the size and type of heteroatom are important for enzyme inhibition. This inhibition is mainly based on the interaction with i) the interfacial region of TIM inducing changes on the quaternary and tertiary structure or ii) with the TIM catalytic region were the main pathways that disabled the catalytic activity of the enzyme.

CONCLUSION

Benzothiazole, benzoxazole, benzimidazole, and sulfhydryl derivatives stand out as TIM inhibitors. In silico and in vitro studies demonstrate that the inhibitors bind mainly at the TIM dimer interface. In this review, the development of new TIM inhibitors as antiprotozoal drugs is demonstrated as an important pharmaceutical strategy that may lead to new therapies for these ancient parasitic diseases.

A Multi-target Drug Designing for BTK, MMP9, Proteasome and TAK1 for the Clinical Treatment of Mantle Cell Lymphoma

BACKGROUND

Mantle cell lymphoma (MCL) is a type of non-Hodgkin lymphoma characterized by the mutation and overexpression of the cyclin D1 protein by the reciprocal chromosomal translocation t(11;14)(q13:q32).

AIM

The present study aims to identify potential inhibition of MMP9, Proteasome, BTK, and TAK1 and determine the most suitable and effective protein target for the MCL.

METHODOLOGY

Nine known inhibitors for MMP9, 24 for proteasome, 15 for BTK and 14 for TAK1 were screened. SB-3CT (PubChem ID: 9883002), oprozomib (PubChem ID: 25067547), zanubrutinib (PubChem ID: 135565884) and TAK1 inhibitor (PubChem ID: 66760355) were recognized as drugs with high binding capacity with their respective protein receptors. 41, 72, 102 and 3 virtual screened compounds were obtained after the similarity search with compound (PubChem ID:102173753), PubChem compound SCHEMBL15569297 (PubChem ID:72374403), PubChem compound SCHEMBL17075298 (PubChem ID:136970120) and compound CID: 71814473 with best virtual screened compounds.

RESULT

MMP9 inhibitors show commendable affinity and good interaction profile of compound holding PubChem ID:102173753 over the most effective established inhibitor SB-3CT. The pharmacophore study of the best virtual screened compound reveals its high efficacy based on various interactions. The virtual screened compound's better affinity with the target MMP9 protein was deduced using toxicity and integration profile studies.

CONCLUSION

Based on the ADMET profile, the compound (PubChem ID: 102173753) could be a potent drug for MCL treatment. Similar to the established SB-3CT, the compound was non-toxic with LD50 values for both the compounds lying in the same range.

Amoebicidal effect of 5,5'-[(4-nitrophenyl)methylene]bis-6-hydroxy-2-mercapto-3-methyl-4(3H)-pyrimidinone), a new drug against Entamoeba histolytica

Entamoeba histolytica is a cosmopolitan protozoan parasite that can produce infections in the intestine and some organs (liver, lungs, and brain), with worldwide prevalence. There are treatments against E. histolytica (antiparasitics), but as the drugs used in these treatments have presented some type of resistance and/or side effects, there are cases with complications of this disease. Therefore, it is necessary to develop new drugs aimed at a specific therapeutic target against this parasite. Here, we used the compound 5,5'-[(4-nitrophenyl)methylene]bis(6-hydroxy-2-mercapto-3-methyl-4(3H)-pyrimidinone) in the patenting process (called D4). D4 has a reported specific use against a glycolytic enzyme, the triosephosphate isomerase of Trichomonas vaginalis (TvTIM). We determined that D4 has an amoebicidal effect in in vitro cultures, with an IC₅₀ value of 18.5 µM, and we proposed a specific site of interaction (Lys77, His110, Gln115, and Glu118) in the triosephosphate isomerase of E. histolytica (EhTIM). Furthermore, compound D4 has favorable experimental and theoretical toxicity results. Therefore, D4 should be further investigated as a potential drug against E. histolytica.

Theoretical evaluation of corrosion inhibition performance of six thiadiazole derivatives

Corrosion inhibition mechanism of six 2-amino-5-alkyl-1,3,4-thiadiazole compounds, for metal surface corrosion was studied by combining quantum chemistry, molecular mechanics and molecular dynamics simulation methods. Molecular reactivity parameters such as [Formula: see text], [Formula: see text], Fukui index were obtained and revealed that the change in alkyl chain length has little influence on the reactivity of thiadiazole inhibitor molecules. Molecular mechanics calculation results show that the molecule with elongated alkyl chain could form self-assembled membrane with higher stability and coverage rate to prevent the diffusion of corrosive substances to metal surface. Molecular dynamics simulation revealed the -equilibrium adsorption behavior of these thiadiazole molecules on metal surface and the calculated results were in agreement with the experimentally determined inhibition efficiencies.

Novel and selective inactivators of Triosephosphate isomerase with anti-trematode activity

Trematode infections such as schistosomiasis and fascioliasis cause significant morbidity in an estimated 250 million people worldwide and the associated agricultural losses are estimated at more than US$ 6 billion per year. Current chemotherapy is limited. Triosephosphate isomerase (TIM), an enzyme of the glycolytic pathway, has emerged as a useful drug target in many parasites, including Fasciola hepatica TIM (FhTIM). We identified 21 novel compounds that selectively inhibit this enzyme. Using microscale thermophoresis we explored the interaction between target and compounds and identified a potent interaction between the sulfonyl-1,2,4-thiadiazole (compound 187) and FhTIM, which showed an IC₅₀ of 5 µM and a K_d of 66 nM. In only 4 hours, this compound killed the juvenile form of F. hepatica with an IC₅₀ of 3 µM, better than the reference drug triclabendazole (TCZ). Interestingly, we discovered in vitro inhibition of FhTIM by TCZ, with an IC₅₀ of 7 µM suggesting a previously uncharacterized role of FhTIM in the mechanism of action of this drug. Compound 187 was also active against various developmental stages of Schistosoma mansoni. The low toxicity in vitro in different cell types and lack of acute toxicity in mice was demonstrated for this compound, as was demonstrated the efficacy of 187in vivo in F. hepatica infected mice. Finally, we obtained the first crystal structure of FhTIM at 1.9 Å resolution which allows us using docking to suggest a mechanism of interaction between compound 187 and TIM. In conclusion, we describe a promising drug candidate to control neglected trematode infections in human and animal health.

Uncovering the Role of Key Active-Site Side Chains in Catalysis: An Extended Brønsted Relationship for Substrate Deprotonation Catalyzed by Wild-Type and Variants of Triosephosphate Isomerase

We report results of detailed empirical valence bond simulations that model the effect of several amino acid substitutions on the thermodynamic (ΔG°) and kinetic activation (ΔG^⧧) barriers to deprotonation of dihydroxyacetone phosphate (DHAP) and d-glyceraldehyde 3-phosphate (GAP) bound to wild-type triosephosphate isomerase (TIM), as well as to the K12G, E97A, E97D, E97Q, K12G/E97A, I170A, L230A, I170A/L230A, and P166A variants of this enzyme. The EVB simulations model the observed effect of the P166A mutation on protein structure. The E97A, E97Q, and E97D mutations of the conserved E97 side chain result in ≤1.0 kcal mol^-1 decreases in the activation barrier for substrate deprotonation. The agreement between experimental and computed activation barriers is within ±1 kcal mol^-1, with a strong linear correlation between ΔG^⧧ and ΔG° for all 11 variants, with slopes β = 0.73 (R² = 0.994) and β = 0.74 (R² = 0.995) for the deprotonation of DHAP and GAP, respectively. These Brønsted-type correlations show that the amino acid side chains examined in this study function to reduce the standard-state Gibbs free energy of reaction for deprotonation of the weak α-carbonyl carbon acid substrate to form the enediolate phosphate reaction intermediate. TIM utilizes the cationic side chain of K12 to provide direct electrostatic stabilization of the enolate oxyanion, and the nonpolar side chains of P166, I170, and L230 are utilized for the construction of an active-site cavity that provides optimal stabilization of the enediolate phosphate intermediate relative to the carbon acid substrate.

The inhibitory mechanism of aurintricarboxylic acid targeting serine/threonine phosphatase Stp1 in Staphylococcus aureus: insights from molecular dynamics simulations

Serine/threonine phosphatase (Stp1) is a member of the bacterial Mg2+- or Mn2+- dependent protein phosphatase/protein phosphatase 2C family, which is involved in the regulation of Staphylococcus aureus virulence. Aurintricarboxylic acid (ATA) is a known Stp1 inhibitor with an IC50 of 1.03 μM, but its inhibitory mechanism has not been elucidated in detail because the Stp1-ATA cocrystal structure has not been determined thus far. In this study, we performed 400 ns molecular dynamics (MD) simulations of the apo-Stp1 and Stp1-ATA complex models. During MD simulations, the flap subdomain of the Stp1-ATA complex experienced a clear conformational transition from an open state to a closed state, whereas the flap domain of apo-Stp1 changed from an open state to a semi-open state. In the Stp1-ATA complex model, the hydrogen bond (H-bond) between D137 and N142 disappeared, whereas critical H-bond interactions were formed between Q160 and H13, Q160/R161 and ATA, as well as N162 and D198. Finally, four residues (D137, N142, Q160, and R161) in Stp1 were mutated to alanine and the mutant enzymes were assessed using phosphate enzyme activity assays, which confirmed their important roles in maintaining Stp1 activity. This study indicated the inhibitory mechanism of ATA targeting Stp1 using MD simulations and sheds light on the future design of allosteric Stp1 inhibitors.

Novel and Selective Rhipicephalus microplus Triosephosphate Isomerase Inhibitors with Acaricidal Activity

The cattle tick Rhipicephalus microplus is one of the most important ectoparasites causing significant economic losses for the cattle industry. The major tool of control is reducing the number of ticks, applying acaricides in cattle. However, overuse has led to selection of resistant populations of R. microplus to most of these products, some even to more than one active principle. Thus, exploration for new molecules with acaricidal activity in R. microplus has become necessary. Triosephosphate isomerase (TIM) is an essential enzyme in R. microplus metabolism and could be an interesting target for the development of new methods for tick control. In this work, we screened 227 compounds, from our in-house chemo-library, against TIM from R. microplus. Four compounds (50, 98, 14, and 161) selectively inhibited this enzyme with IC₅₀ values between 25 and 50 μM. They were also able to diminish cellular viability of BME26 embryonic cells by more than 50% at 50 μM. A molecular docking study showed that the compounds bind in different regions of the protein; compound 14 interacts with the dimer interface. Furthermore, compound 14 affected the survival of partially engorged females, fed artificially, using the capillary technique. This molecule is simple, easy to produce, and important biological data—including toxicological information—are available for it. Our results imply a promising role for compound 14 as a prototype for development of a new acaricidal involving selective TIM inhibition.

Species-Specific Inactivation of Triosephosphate Isomerase from Trypanosoma brucei: Kinetic and Molecular Dynamics Studies

Human African Trypanosomiasis (HAT), a disease that provokes 2184 new cases a year in Sub-Saharan Africa, is caused by Trypanosoma brucei. Current treatments are limited, highly toxic, and parasite strains resistant to them are emerging. Therefore, there is an urgency to find new drugs against HAT. In this context, T. brucei depends on glycolysis as the unique source for ATP supply; therefore, the enzyme triosephosphate isomerase (TIM) is an attractive target for drug design. In the present work, three new benzimidazole derivatives were found as TbTIM inactivators (compounds 1, 2 and 3) with an I₅₀ value of 84, 82 and 73 µM, respectively. Kinetic analyses indicated that the three molecules were selective when tested against human TIM (HsTIM) activity. Additionally, to study their binding mode in TbTIM, we performed a 100 ns molecular dynamics simulation of TbTIM-inactivator complexes. Simulations showed that the binding of compounds disturbs the structure of the protein, affecting the conformations of important domains such as loop 6 and loop 8. In addition, the physicochemical and drug-like parameters showed by the three compounds suggest a good oral absorption. In conclusion, these molecules will serve as a guide to design more potent inactivators that could be used to obtain new drugs against HAT.

Difference FTIR Studies of Substrate Distribution in Triosephosphate Isomerase

Triosephosphate isomerase (TIM) catalyzes the interconversion between dihydroxyacetone phosphate (DHAP) and d-glyceraldehyde 3-phosphate (GAP), via an enediol(ate) intermediate. Determination of substrate population distribution in the TIM/substrate reaction mixture at equilibrium and characterization of the substrate-enzyme interactions in the Michaelis complex are ongoing efforts toward the understanding of the TIM reaction mechanism. By using isotope-edited difference Fourier transform infrared studies with unlabeled and 13C-labeled substrates at specific carbon(s), we are able to show that in the reaction mixture at equilibrium the keto DHAP is the dominant species and the populations of aldehyde GAP and enediol(ate) are very low, consistent with the results from previous X-ray structural and 13C NMR studies. Furthermore, within the DHAP side of the Michaelis complex, there is a set of conformational substates that can be characterized by the different C2═O stretch frequencies. The C2═O frequency differences reflect the different degree of the C2═O bond polarization due to hydrogen bonding from active site residues. The C2═O bond polarization has been considered as an important component for substrate activation within the Michaelis complex. We have found that in the enzyme-substrate reaction mixture with TIM from different organisms the number of substates and their population distribution within the DHAP side of the Michaelis complex may be different. These discoveries provide a rare opportunity to probe the interconversion dynamics of these DHAP substates and form the bases for the future studies to determine if the TIM-catalyzed reaction follows a simple linear reaction pathway, as previously believed, or follows parallel reaction pathways, as suggested in another enzyme system that also shows a set of substates in the Michaelis complex.

Protein-Ligand Docking Methodologies and Its Application in Drug Discovery

Molecular docking is a key tool in structural biology and computer-assisted drug design. Molecular docking is a method which predicts the preferred orientation of a ligand when bound in an active site to form a stable complex. It is the most common method used as a structure-based drug design. Here, the authors intend to discuss the various types of docking methods and their development and applications in modern drug discovery. The important basic theories such as sampling algorithm and scoring functions have been discussed briefly. The performances of the different available docking software have also been discussed. This chapter also includes some application examples of docking studies in modern drug discovery such as targeted drug delivery using carbon nanotubes, docking of nucleic acids to find the binding modes and a comparative study between high-throughput screening and structure-based virtual screening.

Structure–activity relationship and binding mode studies for a series of diketo-acids as HIV integrase inhibitors by 3D-QSAR, molecular docking and molecular dynamics simulations

We explored the main factors affecting the activity of compounds by different statistical and computational methods.

Potent and Selective Inhibitors of Trypanosoma cruzi Triosephosphate Isomerase with Concomitant Inhibition of Cruzipain: Inhibition of Parasite Growth through Multitarget Activity

Triosephosphate isomerase (TIM) is an essential Trypanosoma cruzi enzyme and one of the few validated drug targets for Chagas disease. The known inhibitors of this enzyme behave poorly or have low activity in the parasite. In this work, we used symmetrical diarylideneketones derived from structures with trypanosomicidal activity. We obtained an enzymatic inhibitor with an IC50 value of 86 nm without inhibition effects on the mammalian enzyme. These molecules also affected cruzipain, another essential proteolytic enzyme of the parasite. This dual activity is important to avoid resistance problems. The compounds were studied in vitro against the epimastigote form of the parasite, and nonspecific toxicity to mammalian cells was also evaluated. As a proof of concept, three of the best derivatives were also assayed in vivo. Some of these derivatives showed higher in vitro trypanosomicidal activity than the reference drugs and were effective in protecting infected mice. In addition, these molecules could be obtained by a simple and economic green synthetic route, which is an important feature in the research and development of future drugs for neglected diseases.

3-H-[1,2]Dithiole as a New Anti-Trypanosoma cruzi Chemotype: Biological and Mechanism of Action Studies

The current pharmacological Chagas disease treatments, using Nifurtimox or Benznidazole, show limited therapeutic results and are associated with potential side effects, like mutagenicity. Using random screening we have identified new chemotypes that were able to inhibit relevant targets of the Trypanosoma cruzi. We found 3H-[1,2]dithioles with the ability to inhibit Trypanosoma cruzi triosephosphate isomerase (TcTIM). Herein, we studied the structural modifications of this chemotype to analyze the influence of volume, lipophilicity and electronic properties in the anti-T. cruzi activity. Their selectivity to parasites vs. mammalian cells was also examined. To get insights into a possible mechanism of action, the inhibition of the enzymatic activity of TcTIM and cruzipain, using the isolated enzymes, and the inhibition of membrane sterol biosynthesis and excreted metabolites, using the whole parasite, were achieved. We found that this structural framework is interesting for the generation of innovative drugs for the treatment of Chagas disease.

Exploring the role of quantum chemical descriptors in modeling acute toxicity of diverse chemicals to Daphnia magna

Various quantum-mechanically computed molecular and thermodynamic descriptors along with physico-chemical, electrostatic and topological descriptors are compared while developing quantitative structure-activity relationships (QSARs) for the acute toxicity of 252 diverse organic chemicals towards Daphnia magna. QSAR models based on the quantum-chemical descriptors, computed with routinely employed advanced semi-empirical and ab-initio methods, along with the electron-correlation contribution (CORR) of the descriptors, are analyzed for the external predictivity of the acute toxicity. The models with reliable internal stability and external predictivity are found to be based on the HOMO energy along with the physico-chemical, electrostatic and topological descriptors. Besides this, the total energy and electron-correlation energy are also observed as highly reliable descriptors, suggesting that the intra-molecular interactions between the electrons play an important role in the origin of the acute toxicity, which is in fact an unexplored phenomenon. The models based on quantum-chemical descriptors such as chemical hardness, absolute electronegativity, standard Gibbs free energy and enthalpy are also observed to be reliable. A comparison of the robust models based on the quantum-chemical descriptors computed with various quantum-mechanical methods suggests that the advanced semi-empirical methods such as PM7 can be more reliable than the ab-initio methods which are computationally more expensive.

The unbinding studies of vascular endothelial growth factor receptor-2 protein tyrosine kinase type II inhibitors

Vascular endothelial growth factor receptor-2 (VEGFR-2) tyrosine kinase has two conformations, active and inactive conformations. Type II inhibitors bind to inactive conformation. It has two possible binding/unbinding paths. To explore the unbinding path of inhibitor 01-435 that was generated by fragment build in the binding pocket of VEGFR-2, molecular dynamics (MD) simulation was performed on the crystal structure of VEGFR-2 in complex with 01-435, then steered molecular dynamics (SMD) simulation was executed on the crystal structure of VEGFR-2 in complex with 01-435. Pull force, van der Waals and electrostatic interaction along the two paths were calculated by using SMD simulation. The SMD simulation results indicate that the more favorable path for inhibitor dissociation is along with the traditional ATP-channel rather than the allosteric-pocket-channel, which is mainly due to the less electrostatic interaction that the ligand suffers during dissociation process along the traditional ATP-channel.