As an important target for the development of novel antibiotics, UDP-3-O-(R-3-hydroxymyristoyl)-N-acetylglucosamine deacetylase (LpxC) has been widely studied. Pyridone methylsulfone hydroxamate (PMH) compounds can effectively inhibit the catalytic activity of LpxC. In this work, the three-dimensional quantitative structure-activity relationships of PMH inhibitors were explored and models with good predictive ability were established using comparative molecular field analysis and comparative molecular similarity index analysis methods. The effect of PMH inhibitors' electrostatic potential on the inhibitory ability of Pseudomonas aeruginosa LpxC (PaLpxC) is revealed at the molecular level via molecular electrostatic potential analyses. Then, two molecular dynamics simulations for the PaLpxC and PaLpxC-inhibitor systems were also performed respectively to investigate the key residues of PaLpxC hydrolase binding to water molecules. The results indicate that orderly alternative water molecules can form stable hydrogen bonds with M62, E77, T190, and H264 in the catalytic center, and tetracoordinate to zinc ion along with H78, H237, and D241. It was found that the conformational transition space of PaLpxC changes after association with PMH inhibitors through free energy landscape and cluster analyses. Finally, a possible inhibitory mechanism of PMH inhibitors was proposed, based on our molecular simulation. This paper will provide a theoretical basis for the molecular design of LpxC-targeting antibiotics.

The increasing resistance to anti-tb drugs has enforced strategies for finding new drug targets against Mycobacterium tuberculosis (Mtb). In recent years enzymes associated with the rhamnose pathway in Mtb have attracted attention as drug targets. The present work is on α-D-glucose-1-phosphate thymidylyltransferase (RmlA), the first enzyme involved in the biosynthesis of L-rhamnose, of Mtb cell wall. This study aims to derive a 3D structure of RmlA by using a comparative modeling approach. Structural refinement and energy minimization of the built model have been done with molecular dynamics. The reliability assessment of the built model was carried out with various protein checking tools such as Procheck, Whatif, ProsA, Errat, and Verify 3D. The obtained model investigates the relation between the structure and function. Molecular docking interactions of Mtb-RmlA with modified EMB (ethambutol) ligands and natural substrate have revealed specific key residues Arg13, Lys23, Asn109, and Thr223 which play an important role in ligand binding and selection. Compared to all EMB ligands, EMB-1 has shown better interaction with Mtb-RmlA model. The information thus discussed above will be useful for the rational design of safe and effective inhibitors specific to RmlA enzyme pertaining to the treatment of tuberculosis.

Primary hyperoxaluria type 1 (PHT1) treatment is mainly focused on inhibiting the enzyme glycolate oxidase, which plays a pivotal role in the production of glyoxylate, which undergoes oxidation to produce oxalate. When the renal secretion capacity exceeds, calcium oxalate forms stones that accumulate in the kidneys. In this respect, detailed QSAR analysis, molecular docking, and dynamics simulations of a series of inhibitors containing glycolic, glyoxylic, and salicylic acid groups have been performed employing different regression machine learning techniques. Three robust models with less than 9 descriptors-based on a tenfold cross (Q2 CV) and external (Q2 EXT) validation-were found i.e., MLR1 (Q2 CV = 0.893, Q2 EXT = 0.897), RF1 (Q2 CV = 0.889, Q2 EXT = 0.907), and IBK1 (Q2 CV = 0.891, Q2 EXT = 0.907). An ensemble model was built by averaging the predicted pIC50 of the three models, obtaining a Q2 EXT = 0.933. Physicochemical properties such as charge, electronegativity, hardness, softness, van der Waals volume, and polarizability were considered as attributes to build the models. To get more insight into the potential biological activity of the compouds studied herein, docking and dynamic analysis were carried out, finding the hydrophobic and polar residues show important interactions with the ligands. A screening of the DrugBank database V.5.1.7 was performed, leading to the proposal of seven commercial drugs within the applicability domain of the models, that can be suggested as possible PHT1 treatment.

We synthesized Schiff base and its complexes derivatives of chitosan (CS) in order to develop antibiotic compounds based on functionalized-chitosan against gram-positive and gram-negative bacteria. IR, UV-Vis, AFM, SEM, Melting point, X-ray diffraction (XRD), elemental analysis, and 1H NMR techniques were employed to characterize the chemical structures and properties of these compounds. XRD, UV-Vis, and 1H NMR techniques confirmed the formation of Schiff base and its functionalized-chitosan to metals. Subsequently, our antibacterial studies revealed that antibacterial activities of [Zn(Schiff base)(CS)] against S. aureus bacteria increased compared to those of their compounds. In addition, hemolysis test of CS-Schiff base-Cu(II) demonstrated better hemolytic activity than vitamin C, CS-Schiff base, CS-Schiff base-Zn(II), and CS-Schiff base-Ni(II). In a computational strategy, we carried out the optimization of compounds with molecular mechanics (MM+), Semi-emprical (AM1), Abinitio (STO-3G), AMBER, BIO+(CHARMM), and OPLS. Frontier orbital density distributions (HOMO and LUMO), and the optimized computational UV of the compounds were assessed. The optimized computational UV-Vis was similar to the experimental UV-Vis. We applied the docking methods to predict the DNA binding affinity, Staphylococcus aureus enoyl-acyl carrier protein reductase (ENRs), and Staphylococcus aureus enoyl-acyl carrier protein reductase (saFabI). Ultimately, the obtained data herein suggested that Schiff base is more selective toward ENRs and saFabI compared to chitosan, its complexes, and metronidazole.

In this study two genistein derivatives (G1 and G2) are reported as inhibitors of acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE), and differences in the inhibition of AChE are described. Although they differ in structure by a single methyl group, the inhibitory effect of G1 (IC50=264 nmol/L) on AChE was 80 times stronger than that of G2 (IC50=21,210 nmol/L). Enzyme-kinetic analysis, molecular docking and molecular dynamics (MD) simulations were conducted to better understand the molecular basis for this difference. The results obtained by kinetic analysis demonstrated that G1 can interact with both the catalytic active site and peripheral anionic site of AChE. The predicted binding free energies of two complexes calculated by the molecular mechanics/generalized born surface area (MM/GBSA) method were consistent with the experimental data. The analysis of the individual energy terms suggested that a difference between the net electrostatic contributions (ΔE ele+ΔG GB) was responsible for the binding affinities of these two inhibitors. Additionally, analysis of the molecular mechanics and MM/GBSA free energy decomposition revealed that the difference between G1 and G2 originated from interactions with Tyr124, Glu292, Val294 and Phe338 of AChE. In conclusion, the results reveal significant differences at the molecular level in the mechanism of inhibition of AChE by these structurally related compounds.

We report molecular characterization of the kisspeptin receptor (kiss1r), an essential gatekeeper for reproduction and onset of puberty in vertebrates. The full-length cDNA sequence of kiss1r is 1786bp which consist of 5' UTR (untranslated region) 261bp, 3' UTR of 424bp and open reading frame of 1101 encoding a putative protein of 366 amino acids. Basal tissue expression pattern of kiss1r mRNA revealed that it is mainly expressed in the brain and testis. We also report the structure of the kiss1r, along with plausible activation mechanism of this receptor by kisspeptin using computational modelling and dynamic simulation approach of multiple 100ns of timescale. A present modelling and simulations studies shed light on the molecular level of interaction, suggesting that direct hydrogen bonds between ASN4, SER5, GLY7, ARG9 and PHE10 of kisspeptin and TRP7, ASN8, GLU11, ILE17, ASN19 and TYR183 of kiss1r could be crucial role players in initial binding of receptor and the kisspeptin towards allosteric modulatory effects of kisspeptin on the receptor. To the best our knowledge, this is the first report on computational modelling and molecular dynamic simulations of kiss1r in animals shedding light on its possible mode of activation.

Plasmodium knowlesi, recognized as the fifth Plasmodium parasite, is the least studied malaria parasite. It is a significant cause of morbidity and mortality in the South-East Asia region. Enzymes of folate synthesis, especially dihydrofolate reductase (DHFR), is a well-approved drug target in other Plasmodium species, but its role in Plasmodium knowlesi is poorly studied. This work characterizes PkDHFR as a drug target and identifies inhibitors that can withstand the upcoming problem of resistance. The 3D structure of the PkDHFR target is modelled using comparative modelling, and further, it is refined and validated using energy minimization and torsional angle analysis methods. We extracted 13 compounds from DrugBank and ZINC databases using the "target similarity search" criteria. These compounds were categorized based on their binding affinity (-4.49 to -10.08 kcal/mol) and pose prediction against the active site of PkDHFR. Later on, the top 5 PkDHFR-compound complexes with high or equivalent binding affinity to its natural ligand (dihydrofolate) have undergone for dynamics. The simulation experiments reveal the higher stability of DB00563-PkDHFR complex and less conformational fluctuations and share a similar degree of compactness throughout the simulation trajectory. The MM/GBSA calculation of free energy of DB00563 is also the least (-72.84 kcal/mol) compared to others. Furthermore, the flexible side chain of DB00563 can bind and block the active site of PkDHFR more efficiently. Thus, the identified drug may be considered as a potential candidate for treating P. knowlesi malaria.

Tankyrase enzymes (TNKS), a core part of the canonical Wnt pathway, are a promising target in the search for potential anti-cancer agents. Although several hundreds of the TNKS inhibitors are currently known, identification of their novel chemotypes attracts considerable interest. In this study, the molecular docking and machine learning-based virtual screening techniques combined with the physico-chemical and ADMET (absorption, distribution, metabolism, excretion, toxicity) profile prediction and molecular dynamics simulations were applied to a subset of the ZINC database containing about 1.7 M commercially available compounds. Out of seven candidate compounds biologically evaluated in vitro for their inhibition of the TNKS2 enzyme using immunochemical assay, two compounds have shown a decent level of inhibitory activity with the IC50 values of less than 10 nM and 10 μM. Relatively simple scores based on molecular docking or MM-PBSA (molecular mechanics, Poisson-Boltzmann, surface area) methods proved unsuitable for predicting the effect of structural modification or for accurate ranking of the compounds based on their binding energies. On the other hand, the molecular dynamics simulations and Free Energy Perturbation (FEP) calculations allowed us to further decipher the structure-activity relationships and retrospectively analyze the docking-based virtual screening performance. This approach can be applied at the subsequent lead optimization stages.

In the recent cancer treatment, B-Raf kinase is one of key targets. Nowadays, a group of imidazopyridines as B-Raf kinase inhibitors have been reported. In order to investigate the interaction between this group of inhibitors and B-Raf kinase, molecular docking, molecular dynamic (MD) simulation and binding free energy (ΔGbind) calculation were performed in this work. Molecular docking was carried out to identify the key residues in the binding site, and MD simulations were performed to determine the detail binding mode. The results obtained from MD simulation reveal that the binding site is stable during the MD simulations, and some hydrogen bonds (H-bonds) in MD simulations are different from H-bonds in the docking mode. Based on the obtained MD trajectories, ΔGbind was computed by using Molecular Mechanics Generalized Born Surface Area (MM-GBSA), and the obtained energies are consistent with the activities. An energetic analysis reveals that both electrostatic and van der Waals contributions are important to ΔGbind, and the unfavorable polar solvation contribution results in the instability of the inhibitor with the lowest activity. These results are expected to understand the binding between B-Raf and imidazopyridines and provide some useful information to design potential B-Raf inhibitors.

To explore the inhibitory mechanism of catechins for digestive enzymes, we investigated the binding mode of catechins to a typical digestive enzyme-trypsin and analyzed the structure-activity relationship of catechins, using an integration of molecular docking, molecular dynamics simulation and binding free energy calculation. We found that catechins with different structures bound to a conservative pocket S1 of trypsin, which is comprised of residues 189-195, 214-220 and 225-228. In the trypsin-catechin complexes, Asp189 by forming strong hydrogen bonding, and Gln192, Trp215 and Gly216 through hydrophobic interactions, all significantly contribute to the binding of catechins. The number and the position of hydroxyl and aromatic groups, the structure of stereoisomers, and the orientation of catechins in the binding pocket S1 of trypsin all affect the binding affinity. The binding affinity is in the order of Epigallocatechin gallate (EGCG) > Epicatechin gallate (ECG) > Epicatechin (EC) > Epigallocatechin (EGC), and 2R-3R EGCG shows the strongest binding affinity out of other stereoisomers. Meanwhile, the synergic conformational changes of residues and catechins were also analyzed. These findings will be helpful in understanding the knowledge of interactions between catechins and trypsin and referable for the design of novel polyphenol based functional food and nutriceutical formulas.

To elucidate interactions between the antifungal cyclic lipopeptides iturin A, fengycin, and surfactin produced by Bacillus bacteria and the microtubular protein β-tubulin in plant pathogenic fungi (Fusarium oxysporum, Colletrotrichum gloeosporioides, Alternaria alternata, and Fusarium solani) in molecular docking and molecular dynamics simulations, we retrieved the structure of tubulin co-crystallized with taxol from the Protein Data Bank (PDB) (ID: 1JFF) and the structure of the cyclic lipopeptides from PubChem (Compound CID: 102287549, 100977820, 10129764). Similarity and homology analyses of the retrieved β-tubulin structure with those of the fungi showed that the conserved domains shared 84% similarity, and the root mean square deviation (RMSD) was less than 2 Å. In the molecular docking studies, within the binding pocket, residues Pro274, Thr276, and Glu27 of β-tubulin were responsible for the interaction with the cyclic lipopeptides. In the molecular dynamics analysis, two groups of ligands were formed based on the number of poses analyzed with respect to the RMSD. Group 1 was made up of 10, 100, and 500 poses with distances 0.080 to 0.092 nm and RMSDs of 0.10 to 0.15 nm. For group 2, consisting of 1000 poses, the initial and final distance was 0.1 nm and the RMSDs were in the range of 0.10 to 0.30 nm. These results suggest that iturin A and fengycin bind with higher affinity than surfactin to β-tubulin. These two lipopeptides may be used as lead compounds to develop new antifungal agents or employed directly as biorational products to control plant pathogenic fungi.

The COVID-19 pandemic outbreak demands the designing of potential drugs as there is no specific treatment available. Thanks to their safety and effectiveness, phytochemicals have been used to treat various diseases, including antiviral therapeutics. Molecular docking is a simple, quick, and effective way to screen a variety of molecules for structure-based drug design. Here, we investigate molecular docking experiments on compounds present in plant species, Cocculus hirsutus and Rhodiola rosea and show their potential for the treatment of COVID-19. Almost all the components showed higher binding affinity than the built-in ligand, and those with significantly higher binding affinity were explored further. Molecular mechanics-based generalized born surface area calculations were used to re-rank the top candidates, rhodionidin and cocsoline, and their stability in association with viral protease was confirmed. Density functional theory was used for detailed investigations of the geometries, and electrical properties of rhodionidin and cocsoline. Using the frontier molecular orbitals method, the charge transfer within the molecule was calculated. Chemical reactivity and intermolecular interactions were studied using molecular electrostatic potential maps. These in silico discoveries will simulate the identification of powerful COVID-19 inhibitors, and similar research is likely to make a significant contribution to antiviral drug discovery.

Treatment of several autoimmune diseases and types of cancer has been an intense area of research over the past two decades. Many signaling pathways that regulate innate and/or adaptive immunity, as well as those that induce overexpression or mutation of protein kinases, have been targeted for drug discovery. One of the serine/threonine kinases, Interleukin-1 Receptor Associated Kinase 4 (IRAK4) regulates signaling through various Toll-like receptors (TLRs) and interleukin-1 receptor (IL1R). It controls diverse cellular processes including inflammation, apoptosis, and cellular differentiation. MyD88 gain-of-function mutations or overexpression of IRAK4 has been implicated in various types of malignancies such as Waldenström macroglobulinemia, B cell lymphoma, colorectal cancer, pancreatic ductal adenocarcinoma, breast cancer, etc. Moreover, over activation of IRAK4 is also associated with several autoimmune diseases. The significant role of IRAK4 makes it an interesting target for the discovery and development of potent small molecule inhibitors. A few potent IRAK4 inhibitors such as PF-06650833, RA9 and BAY1834845 have recently entered phase I/II clinical trial studies. Nevertheless, there is still a need of selective inhibitors for the treatment of cancer and various autoimmune diseases. A great need for the same intrigued us to perform molecular modeling studies on 4,6-diaminonicotinamide derivatives as IRAK4 inhibitors. We performed molecular docking and dynamics simulation of 50 ns for one of the most active compounds of the dataset. We also carried out MM-PBSA binding free energy calculation to identify the active site residues, interactions of which are contributing to the total binding energy. The final 50 ns conformation of the most active compound was selected to perform dataset alignment in a 3D-QSAR study. Generated RF-CoMFA (q2 = 0.751, ONC = 4, r2 = 0.911) model revealed reasonable statistical results. Overall results of molecular dynamics simulation, MM-PBSA binding free energy calculation and RF-CoMFA model revealed important active site residues of IRAK4 and necessary structural properties of ligand to design more potent IRAK4 inhibitors. We designed few IRAK4 inhibitors based on these results, which possessed higher activity (predicted pIC50) than the most active compounds of the dataset selected for this study. Moreover, ADMET properties of these inhibitors revealed promising results and need to be validated using experimental studies.

Nuclear retinoic acid receptors (RARs) are ligand-dependent transcription factors involved in various biological processes, such as embryogenesis, cell proliferation, differentiation, reproduction, and apoptosis. These receptors are regulated by retinoids, i.e., retinoic acid (RA) and its analogs, as receptor agonists. RAR agonists are promising therapeutic agents for the treatment of serious dermatological disorders, including some malignant conditions. By inducing apoptosis, they are able to inhibit the proliferation of diverse cancer cell lines. Also, RAR agonists have recently been identified as therapeutic options for some neurodegenerative diseases. These features make retinoids very attractive molecules for medical purposes. Synthetic selective RAR agonists have several advantages over endogenous ones, but they suffer poor pharmacokinetic properties. These compounds are normally lipophilic acids with unfavorable drug-like features such as poor oral bioavailability. Recently, highly selective, potent, and less toxic RAR agonists with proper lipophilicity, thus, good oral bioavailability have been developed for some therapeutic applications. In the present study, ligand and structure-based virtual screening technique was exploited to introduce some novel RARα agonists. Pharmacokinetic assessment was also performed in silico to suggest those compounds which have optimized drug-like features. Finally, two compounds with the best in silico pharmacological features are proposed as lead molecules for future development of RARα agonists.

Cholinesterase inhibitors (ChE-Is) are the standard for the therapy of AD associated disorders and are the only class of approved drugs by the Food and Drug Administration (FDA). Additionally, acetylcholinesterase (AChE) is the target for many Alzheimer's dementia drugs which block the function of AChE but have some side effects. Therefore, in this paper, an attempt was made to elucidate cholinesterase inhibition potential of secondary metabolite from Cannabis plant which has negligible or no side effect. Molecular docking of 500 herbal compounds, against AChE, was performed using Autodock 4.2 as per the standard protocols. Molecular dynamics simulations have also been carried out to check stability of binding complex in water for 1000 ps. Our molecular docking and simulation have predicted high binding affinity of secondary metabolite (C28H34N2O6) to AChE. Further, molecular dynamics simulations for 1000 ps suggest that ligand interaction with the residues Asp72, Tyr70-121-334, and Phe288 of AChE, all of which fall under active site/subsite or binding pocket, might be critical for the inhibitory activity of AChE. This approach might be helpful to understand the selectivity of the given drug molecule in the treatment of Alzheimer's disease. The study provides evidence for consideration of C28H34N2O6 as a valuable small ligand molecule in treatment and prevention of AD associated disorders and further in vitro and in vivo investigations may prove its therapeutic potential.

Herbal bioactive compounds have captured pronounced attention considering their health-promoting effects as well as their functional properties. In this study, the binding mechanism between milk protein bovine β-lactoglobulin (β-LG), oleuropein (OLE) and safranal (SAF) found in olive leaf extract and saffron, respectively via spectroscopic and in silico studies. Fluorescence quenching information exhibited that interactions with both ligands were spontaneous and hydrophobic interactions were dominant. Also, the CD spectroscopy results demonstrated the increase in β-sheet structure and decrease in the α-helix content for both ligands. Size of β-LG-OLE complex was higher than β-LG-SAF due to the conformation and larger molecular size. Molecular docking and simulation studies revealed that SAF and OLE bind in the central calyx of β-LG and the surface of β-LG next to hydrophobic residues. Lastly, OLE formed a more stabilized complex compared to SAF based on the molecular dynamic simulation results.

This paper provides a feasible model for molecular structure analysis and interaction mechanism of aptamer and micromolecule. In this study, modeling and dynamic simulation of ssDNA aptamer (P-18S2) and target (Palytoxin, PTX) were performed separately. Then, the complex structure between DNA and PTX was predicted, and docking results showed that PTX could combine steadily at the groove's top of DNA model by strong hydrogen-bonds and electrostatic interaction. Thus, we truncated and optimized P-18S2 by simulating. At the same time, we also confirmed the reliability of simulation results by experiments. With the experimental and computational results, the study provided a more reasonable interpretation for the high affinity and specific binding of P-18S2 and PTX, which laid the foundation for further optimization and development of aptamers in molecular diagnostics and therapeutic applications.

Our investigation intended to analyze the chemical composition and the antioxidant activity of Carrichtera annua and to evaluate the antiproliferative effect of C. annua crude and phenolics extracts by MTT assay on a panel of cancerous and non-cancerous breast and liver cell lines. The total flavonoid and phenolic contents of C. annua were 47.3 ± 17.9 mg RE/g and 83.8 ± 5.3 mg respectively. C. annua extract exhibited remarkable antioxidant capacity (50.92 ± 5.64 mg GAE/g) in comparison with BHT (74.86 ± 3.92 mg GAE/g). Moreover, the extract exhibited promising reduction ability (1.17 mMol Fe+2/g) in comparison to the positive control (ascorbic acid with 2.75 ± 0.91) and it displayed some definite radical scavenging effect on DPPH (IC50 values of 211.9 ± 3.7 µg/mL). Chemical profiling of C. annua extract was achieved by LC-ESI-TOF-MS/MS analysis. Forty-nine hits mainly polyphenols were detected. Flavonoid fraction of C. annua was more active than the crude extract. It demonstrated selective cytotoxicity against the MCF-7 and HepG2 cells (IC50 = 13.04 and 19.3 µg/mL respectively), induced cell cycle arrest at pre-G1 and G2/M-phases and displayed apoptotic effect. Molecular docking studies supported our findings and revealed that kaempferol-3,7-O-bis-α-L-rhamnoside and kaempferol-3-rutinoside were the most active inhibitors of Bcl-2. Therefore, C. annua herb seems to be a promising candidate to further advance anticancer research. In extrapolation, the intake of C. annua phenolics might be adventitious for alleviating breast and liver malignancies and tumoral proliferation in humans.

We aim to investigate the mechanism and effective components of Shenghui decoction (SHD), which has been shown to inhibit acetylcholinesterase (AChE) through molecular docking (MD) and molecular dynamics simulation (MDS). The effective ingredients in SHD were collected through the TCMSP database and literature review. All components were docked with AChE using CDOCKER. Receptor ligand interaction analysis was performed for the optimal ligand. Two simulation models (model I and II) containing AChE and acetylcholine (ACh) were constructed, in which model II contained the best-docked ligand. Perform 90ns MDS on 2 models. After the simulation, the distance between ACh and AChE peripheral active sites were calculated in both models. The root mean square deviation (RMSD) curve of ligand and receptor, the radius of gyration (Rog) of the receptor, the distance between ligand center and binding site center, and the binding energy of ligand and receptor were calculated in model II. 98 ingredients of SHD were collected, and the best ligand was Tumulosic acid. The residues that form conventional hydrogen bonds between AChE and Tumulosic acid include Tyr132 and Glu201. MDS revealed that ACh could bind to AChE active site in model I. In model II, ACh cannot bind to the binding cavity because the ligand occupies the active site. The RMSD of AChE and Tumulosic acid tends to be stable, the Rog curve of AChE is relatively stable, and the distance between ligand and binding cavity does not fluctuate greatly, indicating that the structure of receptor and ligand is relatively stable. The binding energy of AChE and Tumulosic acid was -24.14 ± 2.46 kcal/mol. SHD contains many effective ingredients that may inhibit AChE activity. Tumulosic acid can occupy the binding site to prevent ACh from entering the chemical domain, thus exerting AChE inhibitory effect.

Protein prenylation is a posttranslational modification that is indispensable for translocation of membrane GTPases like Ras, Rho, Ras etc. Proteins of Ras family undergo farnesylation by FTase while Rho family goes through geranylgeranylation by GGTase1. There is only an infinitesimal difference in signal recognition between FTase and GGTase1. FTase inhibitors mostly end up selecting the cells with mutated Ras proteins that have acquired affinity towards GGTase1 in cancer microcosms. Therefore, it is of interest to identify GGTase1 and FTase dual inhibitors using the docking tool AutoDock Vina. Docking data show that curcumin (from turmeric) has higher binding affinity to GGTase1 than that of established peptidomimetic GGTase1 inhibitors (GGTI) such as GGTI-297, GGTI-298, CHEMBL525185. Curcumin also interacts with FTase with binding energy comparable to co-crystalized compound 2-[3-(3-ethyl-1-methyl-2-oxo-azepan-3-yl)-phenoxy]-4-[1-amino-1-(1-methyl-1h-imidizol-5-yl)-ethyl]-benzonitrile (BNE). The docked complex was further simulated for 10 ns using molecular dynamics simulation for stability. Thus, the molecular basis for curcumin binding to GGTase1 and FTase is reported.