Curcumin specifically binds to the human calcium–calmodulin-dependent protein kinase IV: fluorescence and molecular dynamics simulation studies

Calcium–calmodulin-dependent protein kinase IV (CAMK4) plays significant role in the regulation of calcium-dependent gene expression, and thus, it is involved in varieties of cellular functions such as cell signaling and neuronal survival. On the other hand, curcumin, a naturally occurring yellow bioactive component of turmeric possesses wide spectrum of biological actions, and it is widely used to treat atherosclerosis, diabetes, cancer, and inflammation. It also acts as an antioxidant. Here, we studied the interaction of curcumin with human CAMK4 at pH 7.4 using molecular docking, molecular dynamics (MD) simulations, fluorescence binding, and surface plasmon resonance (SPR) methods. We performed MD simulations for both neutral and anionic forms of CAMK4-curcumin complexes for a reasonably long time (150 ns) to see the overall stability of the protein–ligand complex. Molecular docking studies revealed that the curcumin binds in the large hydrophobic cavity of kinase domain of CAMK4 through several hydrophobic and hydrogen-bonded interactions. Additionally, MD simulations studies contributed in understanding the stability of protein–ligand complex system in aqueous solution and conformational changes in the CAMK4 upon binding of curcumin. A significant increase in the fluorescence intensity at 495 nm was observed (λ_exc = 425 nm), suggesting a strong interaction of curcumin to the CAMK4. A high binding affinity (K_D = 3.7 × 10^?8 ± .03 M) of curcumin for the CAMK4 was measured by SPR further indicating curcumin as a potential ligand for the CAMK4. This study will provide insights into designing a new inspired curcumin derivatives as therapeutic agents against many life-threatening diseases.     

In silico approaches to evaluate the molecular properties of organophosphate compounds to inhibit acetylcholinesterase activity in housefly

Organophosphate compounds (OPC) have become the primary choice as insecticides and are widely used across the world. Additionally, OPCs were also commonly used as a chemical warfare agent that triggers a great challenge to public safety. Exposure of OPCs to human causes immediate excitation of cholinergic neurotransmission through transient elevation of synaptic acetylcholine (ACh) levels and accumulations. Likewise, prolonged exposure of OPCs can affect the processes in immune response, carbohydrate metabolism, cardiovascular toxicity, and several others. Studies revealed that the toxicity of OPCs was provoked by inhibition of acetylcholinesterase (AChE). Therefore, combined in silico approaches – pharmacophore-based 3D-QSAR model; docking and Molecular Dynamics (MD) – were used to assess the precise and comprehensive effects of series of known OP-derived compounds together with its ?log LD₅₀ values. The selected five-featured pharmacophore model – AAHHR.61 – displayed the highest correlation (R² = .9166), cross-validated coefficient (Q² = .8221), F = 63.2, Pearson-R = .9615 with low RMSE = .2621 values obtained using five component PLS factors. Subsequently, the well-validated model was then used as a 3D query to search novel OPCs using a high-throughput virtual screening technique. Simultaneously, the docking studies predicted the binding pose of the most active OPC in the MdAChE binding pocket. Additionally, the stability of docking was verified using MD simulation. The results revealed that OP22 and predicted lead compounds bound tightly to S315 of MdAChE through potential hydrogen bond interaction over time. Overall, this study might provide valuable insight into binding mode of OPCs and hit compounds to inhibit AChE in housefly.     

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.     

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.     

Insights into molecular interactions between CaM and its inhibitors from molecular dynamics simulations and experimental data

In order to contribute to the structural basis for rational design of calmodulin (CaM) inhibitors, we analyzed the interaction of CaM with 14 classic antagonists and two compounds that do not affect CaM, using docking and molecular dynamics (MD) simulations, and the data were compared to available experimental data. The Ca²⁺-CaM-Ligands complexes were simulated 20 ns, with CaM starting in the “open” and “closed” conformations. The analysis of the MD simulations provided insight into the conformational changes undergone by CaM during its interaction with these ligands. These simulations were used to predict the binding free energies (ΔG) from contributions ΔH and ΔS, giving useful information about CaM ligand binding thermodynamics. The ΔG predicted for the CaM’s inhibitors correlated well with available experimental data as the r² obtained was 0.76 and 0.82 for the group of xanthones. Additionally, valuable information is presented here: I) CaM has two preferred ligand binding sites in the open conformation known as site 1 and 4, II) CaM can bind ligands of diverse structural nature, III) the flexibility of CaM is reduced by the union of its ligands, leading to a reduction in the Ca²⁺-CaM entropy, IV) enthalpy dominates the molecular recognition process in the system Ca²⁺-CaM-Ligand, and V) the ligands making more extensive contact with the protein have higher affinity for Ca²⁺-CaM. Despite their limitations, docking and MD simulations in combination with experimental data continue to be excellent tools for research in pharmacology, toward a rational design of new drugs.     

Cytotoxicity and comparative binding mechanism of piperine with human serum albumin and α-1-acid glycoprotein

Human serum albumin (HSA) and α-1-acid glycoprotein (AGP) (acute phase protein) are the plasma proteins in blood system which transports many drugs. To understand the pharmacological importance of piperine molecule, here, we studied the anti-inflammatory activity of piperine on mouse macrophages (RAW 264.7) cell lines, which reveals that piperine caused an increase in inhibition growth of inflammated macrophages. Further, the fluorescence maximum quenching of proteins were observed upon binding of piperine to HSA and AGP through a static quenching mechanism. The binding constants obtained from fluorescence emission were found to be K_piperine?=?5.7 ± .2 × 10⁵ M^?1 and K_piperine = 9.3± .25 × 10⁴ M^?1 which correspond to the free energy of ?7.8 and ?6.71 kcal M^?1at 25 °C for HSA and AGP, respectively. Further, circular dichrosim studies revealed that there is a marginal change in the secondary structural content of HSA due to partial destabilization of HSA–piperine complexes. Consequently, inference drawn from the site-specific markers (phenylbutazone, site I marker) studies to identify the binding site of HSA noticed that piperine binds at site I (IIA), which was further authenticated by molecular docking and molecular dynamic (MD) studies. The binding constants and free energy corresponding to experimental and computational analysis suggest that there are hydrophobic and hydrophilic interactions when piperine binds to HSA. Additionally, the MD studies have showed that HSA–piperine complex reaches equilibration state at around 3 ns, which prove that the HSA–piperine complex is stable in nature.     

Probing the interaction between cHAVc3 peptide and the EC1 domain of E-cadherin using NMR and molecular dynamics simulations

The goal of this work is to probe the interaction between cyclic cHAVc3 peptide and the EC1 domain of human E-cadherin protein. Cyclic cHAVc3 peptide (cyclo(1,6)Ac-CSHAVC-NH₂) binds to the EC1 domain as shown by chemical shift perturbations in the 2D ¹H,-¹⁵N-HSQC NMR spectrum. The molecular dynamics (MD) simulations of the EC1 domain showed folding of the C-terminal tail region into the main head region of the EC1 domain. For cHAVc3 peptide, replica exchange molecular dynamics (REMD) simulations generated five structural clusters of cHAVc3 peptide. Representative structures of cHAVc3 and the EC1 structure from MD simulations were used in molecular docking experiments with NMR constraints to determine the binding site of the peptide on EC1. The results suggest that cHAVc3 binds to EC1 around residues Y36, S37, I38, I53, F77, S78, H79, and I94. The dissociation constants (K_d values) of cHAVc3 peptide to EC1 were estimated using the NMR chemical shifts data and the estimated K_ds are in the range of .5 × 10^?5–7.0 × 10^?5 M.     

Synthesis and molecular docking studies of some 4-phthalimidobenzenesulfonamide derivatives as acetylcholinesterase and butyrylcholinesterase inhibitors

A series of 4-phthalimidobenzenesulfonamide derivatives were designed, synthesized and evaluated for the inhibitory activities against acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE). Structures of the title compounds were confirmed by spectral and elemental analyses. The cholinesterase (ChE) inhibitory activity studies were carried out using Ellman’s colorimetric method. The biological activity results revealed that all of the title compounds (except for compound 8) displayed high selectivity against AChE. Among the tested compounds, compound 7 was found to be the most potent against AChE (IC₅₀=?1.35?±?0.08?μM), while compound 3 exhibited the highest inhibition against BuChE (IC₅₀=?13.41?±?0.62?μM). Molecular docking studies of the most active compound 7 in AChE showed that this compound can interact with both the catalytic active site (CAS) and the peripheral anionic site (PAS) of AChE.     

In silico studies on p21-activated kinase 4 inhibitors: comprehensive application of 3D-QSAR analysis,molecular docking,molecular dynamics simulations,and MM-GBSA calculation

Abstract

P21-activated kinase 4 (PAK4) is a serine/threonine protein kinase, which is associated with many cancer diseases, and thus being considered as a potential drug target. In this study, three-dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulations were performed to explore the structure-activity relationship of a series of pyrropyrazole PAK4 inhibitors. The statistical parameters of comparative molecular field analysis (CoMFA, Q ² = 0.837, R ² = 0.990, and R ² _pred = 0.967) and comparative molecular similarity indices analysis (CoMSIA, Q ² = 0.720, R ² = 0.972, and R ² _pred = 0.946) were obtained from 3D-QSAR model, which exhibited good predictive ability and significant statistical reliability. The binding mode of PAK4 with its inhibitors was obtained through molecular docking study, which indicated that the residues of GLU396, LEU398, LYS350, and ASP458 were important for activity. Molecular mechanics generalized born surface area (MM-GBSA) method was performed to calculate the binding free energy, which indicated that the coulomb, lipophilic and van der Waals (vdW) interactions made major contributions to the binding affinity. Furthermore, through 100?ns MD simulations, we obtained the key amino acid residues and the types of interactions they participated in. Based on the constructed 3D-QSAR model, some novel pyrropyrazole derivatives targeting PAK4 were designed with improved predicted activities. Pharmacokinetic and toxicity predictions of the designed PAK4 inhibitors were obtained by the pkCSM, indicating these compounds had better absorption, distribution, metabolism, excretion and toxicity (ADMET) properties. Above research provided a valuable insight for developing novel and effective pyrropyrazole compounds targeting PAK4.     

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.     

Phase-transition properties of glycerol-monopalmitate lipid bilayers investigated by molecular dynamics simulation: influence of the system size and force-field parameters

Phase-transition properties of glycerol-1-monopalmitate (GMP) bilayers are investigated using explicit-solvent molecular dynamics (MD) simulations, initiated from structures appropriate for the gel (GL) or liquid crystal (LC) phases, and carried out at different hydration levels and temperatures. Building up on a previous study and based on 600 ns simulations, the influence of the system size and of the force field on the equilibrium thermodynamic and dynamic parameters of the bilayers in the GL and LC phases, as well as on the temperature T_m and properties of the GL ? LC phase transition, are analysed. Qualitatively speaking, the results agree with the available experimental data for the area per lipid in the two phases and for the phase-transition temperatures at the three hydration levels irrespective of the selected model parameters. They also suggest that the total number of hydrogen bonds formed between a lipid headgroup and its environment is essentially constant, amounting to about four in both the LC and the GL phases. Quantitatively speaking, the dependence of T_m on the hydration level is found to be non-systematic across the different combinations of model parameters. This results in part from a sensitivity of the results on the system size and force-field parameters but also from the limited accuracy of the bracketing approach employed here to estimate T_m. Finally, a simple kinetic model is proposed to account for the timescales of the transitions. This model involves enthalpy and entropy increases of about 26 kJ mol^? 1 and 83 J mol^? 1 K^? 1 per lipid, upon going from the GL to the LC phase. The transition state is associated with activation parameters corresponding to 13% and 11%, respectively, of these values along the GL → LC transition, resulting in an activation free energy of about 0.3 kJ mol^? 1 per lipid at T_m.     

Differential interactions and structural stability of chitosan oligomers with human serum albumin and α-1-glycoprotein

Chitosan is a naturally occurring deacetylated derivative of chitin with versatile biological activities. Here, we studied the interaction of chitosan oligomers with low degree of polymerization such as chitosan monomer (CM), chitosan dimer (CD), and chitosan trimer (CT) with human serum albumin (HSA) a major blood carrier protein and α-1-glycoprotein (AGP). Since, HSA and AGP are the two important plasma proteins that determine the drug disposition and affect the fate of distribution of drugs. Fluorescence emission spectra indicated that CM, CD, and CT had binding constants of K_CM = 6.2 ± .01 × 10⁵ M^?1, K_CD = 5.0 ± .01 × 10⁴ M^?1, and K_CT = 1.6 ± .01 × 10⁶ M^?1_, respectively, suggesting strong binding with HSA. However, binding of chitooligomers with AGP was insignificant. Thermodynamic and molecular docking analysis indicated that hydrogen bonds and also hydrophobic interaction played an important role in stabilizing the HSA-chitooligomer complexes with free energies of ?7.87, ?6.35, and ?8.4?Kcal/mol for CM, CD, and CT, respectively. Further, circular dichroism studies indicated a minor unfolding of HSA secondary structure, upon interaction with chitooligomers, which are supported with fluctuations of root mean square deviation (RMSD) and radius of gyration (R_g) of HSA. Docking analysis revealed that all three chitooligomers were bound to HSA within subdomain IIA (Site I). In addition, RMSD and R_g analysis depicted that HSA-chitooligomer complexes stabilized at around 4.5 ns. These results suggest that HSA might serve as a carrier in delivering chitooligomers to target tissues than AGP which has pharmacological importance.     

Westalsan: A New Acetylcholine Esterase Inhibitor from the Endophytic Fungus Westerdykella nigra

A new cytochalasan alkaloid, westalsan ( 1 ), along with two known cytochalasan compounds, phomacin B ( 2 ) and 19-hydroxy-19,20-dihydrophomacin C ( 3 ), were isolated from the solid rice culture of Westerdykella nigra, a marine-derived endophytic fungus, isolated from the roots of mangrove Avicennia marina (Forssk.) Vierh. The structures of compounds 1 – 3 were established on the basis of extensive 1D and 2D NMR spectroscopic techniques in combination with HR-ESI-MS. The ability of the isolated compounds to inhibit acetylcholine esterase activity was evaluated. Compound 3 showed the highest acetylcholine esterase inhibitory activity (IC₅₀ 0.056±0.003 μM), followed by compound 1 (IC₅₀ 0.088±0.005 μM) and compound 2 (IC₅₀ 0.140±0.007 μM) compared to donepezil (IC₅₀ 0.035±0.002 μM). This was further confirmed by molecular docking experiment.     

In vitro and in silico evaluation of fucosterol from Sargassum horridum as potential human acetylcholinesterase inhibitor

The fucosterol has been reported numerous biological activities. In this study, the activity in vitro of the fucosterol from Sargassum horridum as potential human acetylcholinesterase inhibitor was evaluated. The structural identification was obtained by nuclear magnetic resonance (NMR) spectroscopy and based on experimental data, we combined docking and molecular dynamics simulations coupled to the molecular-mechanics-generalized-born-surface-area approach to evaluating the structural and energetic basis for the molecular recognition of fucosterol and neostigmine at the binding site of acetylcholinesterase (AChE). In addition, the Lineweaver–Burk plot showed the nature of a non-competitive inhibition. The maximum velocity (V_max) and the constant of Michaelis–Menten (K_m) estimated for fucosterol (0.006 µM) were 0.015 1/V_o (ΔA/h and 6.399 1/[ACh] mM^?1, respectively. While, for neostigmine (0.14 µM), the V_max was 0.022 1/V_o (ΔA/h) and K_m of 6.726 1/[ACh] mM^?1, these results showed a more effective inhibition by fucosterol respect to neostigmine. Structural analysis revealed that neostigmine reaches the AChE binding site reported elsewhere, whereas fucosterol can act as a no-competitive and competitive acetylcholinesterase inhibitor, in agree with kinetic enzymatic experiments. Binding free energy calculations revealed that fucosterol reaches the acetylcholinesterase binding site with higher affinity than neostigmine, which is according to experimental results. Whereas the per-residue decomposition free energy analysis let us identify crucial residues involved in the molecular recognition of ligands by AChE. Results corroborate the ability of theoretical methods to provide crucial information at the atomic level about energetic and structural differences in the binding interaction and affinity from fucosterol with AChE.

Communicated by Ramaswamy H. Sarma     

Exploration of human serum albumin binding sites by docking and molecular dynamics flexible ligand–protein interactions

Five‐nanosecond molecular dynamics (MD) simulations were performed on human serum albumin (HSA) to study the conformational features of its primary ligand binding sites (I and II). Additionally, 11 HSA snapshots were extracted every 0.5 ns to explore the binding affinity (K_d) of 94 known HSA binding drugs using a blind docking procedure. MD simulations indicate that there is considerable flexibility for the protein, including the known sites I and II. Movements at HSA sites I and II were evidenced by structural analyses and docking simulations. The latter enabled the study and analysis of the HSA–ligand interactions of warfarin and ketoprofen (ligands binding to sites I and II, respectively) in greater detail. Our results indicate that the free energy values by docking (K_d observed) depend upon the conformations of both HSA and the ligand. The 94 HSA–ligand binding K_d values, obtained by the docking procedure, were subjected to a quantitative structure‐activity relationship (QSAR) study by multiple regression analysis. The best correlation between the observed and QSAR theoretical (K_d predicted) data was displayed at 2.5 ns. This study provides evidence that HSA binding sites I and II interact specifically with a variety of compounds through conformational adjustments of the protein structure in conjunction with ligand conformational adaptation to these sites. These results serve to explain the high ligand‐promiscuity of HSA. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 161–170, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Unraveling the stability of plasma proteins upon interaction of synthesized uridine products: biophysical and molecular dynamics approach

Abstract

Most of the drugs binding to human serum albumin (HSA) are transported to various parts of the body. Here, we have studied the molecular interaction between HSA and synthesized uridine derivatives, 1-[(3R, 4S, 5?R)-2-methyl-3, 4-dihydroxy-5-(hydroxymethyl)oxolan-2-yl]pyrimidine-2,4-dion.)(C-MU); [(2R,3R,4R,5R)-5-(2,4-dioxo-3,4-dihydropyrimidin-1(2H)-yl)-3,4-dihydroxy-4-methyl-tetrahydrofuran-2-yl] methyl methyl phosphochloridate (CM-MU) and [(2R,3S,4R,5R)-5-(2,4-dioxopyrimidin-1-yl)-2-methyl-3,4-dihydroxyoxolan-2-yl] methyl dihydrogen phosphate (P-MU). Cytotoxic studies of these synthesized compounds with mouse macrophages (RAW 246.7) and HeLa cells (human cervical cancer cells) and binding mechanism of these uridine derivatives with HSA were performed. Subsequently, fluorescence quenching was observed upon titration of uridine derivatives with HSA via static mode of quenching, and the binding constants (K_2-C-MU = 4?±?0.03?×?10⁴M^?1, K_5-CM-MU = 1.95?±?0.03?×?10⁴ M^?1 and K_5-P-MU =1.56?±?0.03?×?10⁴ M^?1) were found to be in sync with the computational results. Further, molecular displacement and molecular docking data revealed that all the derivatives are binding in the subdomain IIA and IIB regions of HSA. The protein secondary structure of complexes was determined by circular dichroism, indicating partial unfolding of the protein upon addition of the uridine derivatives. Furthermore, atomic force microscopy data reveal the change in topology upon binding of 2-C-MU, 5-CM-MU and 5-P-MU with HSA, indicating change in the microenvironment around tryptophan region. Additionally, cytotoxicity studies on HeLa and Raw Cell lines suggested that these molecules have significant anti-proliferative and anti-inflammatory properties. Hence, the study may be of help for development of new drugs based on uridine derivatives which may be helpful for combating various potential diseases.

Communicated by Ramaswamy H. Sarma     

Binding of NDGA and morin with phospholipase A2: experimental and computational evidences

The effects of morin and nordihydroguaiaretic acid (NDGA), two plant secondary metabolites, on porcine pancreatic phospholipase A₂ (PLA₂) were investigated by isothermal titration calorimetry (ITC) and in silico docking analyses. The binding energies obtained for NDGA and morin from the ITC studies are ? 6.36 and ? 5.91 kcal mol^? 1, respectively. Similarly, the glide scores obtained for NDGA and morin towards PLA₂ were ? 7.32 and ? 7.23 kcal mol^? 1, respectively. Further the docked complexes were subjected to MD simulation in the presence of explicit water molecules to check the binding stability of the ligands in the active site of PLA₂. The bound ligands make hydrogen bonds with the active site residues of the enzyme and coordinate bonds with catalytically important Ca²⁺ ion. The binding of ligands at the active site of PLA₂ may also contribute to the reported anti-inflammatory properties of NDGA and morin.