Next Article in Journal
Ultrasound Promoted Synthesis of Bis(substituted pyrazol-4-ylcarbonyl)-Substituted Thioureas
Next Article in Special Issue
Phenolics: From Chemistry to Biology
Previous Article in Journal
A New Insecticidal Sesquiterpene Ester from Celastrus Angulatus
Previous Article in Special Issue
Polyphenols with Antiulcerogenic Action from Aqueous Decoction of Mango Leaves (Mangifera indica L.)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Docking Studies and Anti-enzymatic Activities of Thai Mango Seed Kernel Extract Against Snake Venoms

by
Jiraporn Leanpolchareanchai
1,
Pimolpan Pithayanukul
1,*,
Rapepol Bavovada
2 and
Patchreenart Saparpakorn
3
1
Department of Pharmacy, Faculty of Pharmacy, Mahidol University, Bangkok 10400, Thailand
2
Department of Pharmaceutical Botany, Faculty of Pharmaceutical Sciences, Chulalongkorn University, Bangkok 10330, Thailand
3
Department of Chemistry, Faculty of Science, Kasetsart University, Bangkok 10900, Thailand
*
Author to whom correspondence should be addressed.
Molecules 2009, 14(4), 1404-1422; https://doi.org/10.3390/molecules14041404
Submission received: 16 January 2009 / Revised: 26 March 2009 / Accepted: 30 March 2009 / Published: 31 March 2009
(This article belongs to the Special Issue Phenolics and Polyphenolics)

Abstract

:
The ethanolic extract from seed kernels of Thai mango (MSKE) (Mangifera indica L. cv. ‘Fahlun’) (Anacardiaceae) and its major phenolic principle (pentagalloyl glucopyranose) exhibited dose-dependent inhibitory effects on enzymatic activities of phospholipase A2 (PLA2), hyaluronidase and L-amino acid oxidase (LAAO) of Calloselasma rhodostoma (CR) and Naja naja kaouthia (NK) venoms by in vitro tests. The anti-hemorrhagic and anti-dermonecrotic activities of MSKE against both venoms were clearly supported by in vivo tests. Molecular docking studies indicated that the phenolic molecules of the MSKE could selectively bind to the active sites or their proximity, or modify conserved residues that are critical for the catalysis of PLA2, and selectively bind to the LAAO binding pocket of both CR and NK venoms and thereby inhibit their enzymatic activities. The results imply a potential use of MSKE against snake venoms.

1. Introduction

Venomous snakebite is an important public health problem in remote areas of Thailand since the economic activities in these areas are mainly agricultural. Naja naja kaouthia Lesson (NK) (Thai cobra, Elapidae) causes the highest mortality due to snake envenomation in Thailand [1] and Calloselasma rhodostoma Kuhl (CR) (Malayan pit viper, Viperidae) causes the greatest incidence of envenomation in the country [2]. Among the many venomous snakes in Thailand, CR and NK venoms cause the most severe necrosis on the bite local tissue.
Snake venoms are complex mixtures comprised mainly of proteins and peptides possessing a variety of biological activities. Venom proteins have many diverse enzymatic activities [3]. Enzymes with hydrolytic activity such as phospholipase A2 (PLA2), hyaluronidase and L-amino acid oxidase (LAAO) are found in most snake venoms and are deeply involved in snake envenomation by inducing local effects such as severe inflammatory reactions, hemorrhage and necrosis of local tissues [4]. PLA2s (EC 3.1.1.4) are enzymes that specifically hydrolyzes the sn-2 ester bond of 3-sn membrane glycerophospholipids, generating lysophospholipids and free fatty acids, which themselves cause considerable erythrocyte membrane lytic and damage [5]. In addition, PLA2s exhibit wide varieties of pharmacological effects such as neurotoxicity, cardiotoxicity, myotoxicity, as well as necrotic, anticoagulant, hypotensive, hemolytic, hemorrhagic and edematogenic [6]. PLA2s from snake venoms are classified into groups I or II, based on their sequences and modes of disulphide pairings. Group I PLA2s are found in the venom of Elapidae snakes, whereas group II PLA2s are found in the venom of Viperidae snakes. PLA2s from different snake venoms share a high degree of homology in the amino acid sequence and enzymatic active sites. They hydrolyze phospholipids using a His-Asp doublet plus a conserved water molecule as a nucleophile, and a Ca2+ ion as a cofactor [7]. Hyaluronidase (EC 3.2.1.35), known as endoglycosidase, is frequently referred to a “spreading factor” owing to its ability to cause the distortion in the structural integrity of extracellular matrix in local tissues. This eventually results in the dissemination of target specific toxins. The process described herein is presumed to be the critical step in the enzyme-mediated spreading process [8]. LAAO (EC 1.4.3.2) is a flavoenzyme which catalyses the stereospecific oxidative deamination of an L-amino acid substrate to an α-keto acid along with the production of ammonia and hydrogen peroxide (H2O2) [9]. This enzyme induces apoptosis in mammalian endothelial cells possibly through the production of highly localized concentrations of H2O2[10]. In addition, Sakurai et al. [11] suggested that LAAO from NK venom also functions as a platelet aggregation inhibitor due to H2O2 formation. Recently, Pawelek et al. [12] demonstrated the high-resolution X-ray crystallographic structure of LAAO from CR venom, which indicated that snake venom LAAO is functionally a dimer. Each subunit consists of three domains: a FAD-binding domain, a substrate-binding domain and a helical domain. A funnel is formed by the interface between the substrate-binding and helical domains and provides substrate access to the active site.
Plant extracts constitute a rich source of bioactive compounds with a variety of pharmacological activities. Tannins from plants have been shown to interact with enzymes from snake venoms and act as an antidote [13,14,15,16]. The ethanolic extract of Thai mango seed kernels (MSKE) (Mangifera indica L. cv. ‘Fahlun’, Anacardiaceae) has a relatively high phenolic content and is composed of 1,2,3,4,6-penta-O-galloyl-β-D-glucopyranose (PGG) (61.28%), methyl gallate (MG) (0.68%) and gallic acid (GA) (0.44%) [17]. The extract and its phenolic principles have been shown to exhibit potent antioxidant, anti-tyrosinase, anti-inflammatory and hepatoprotective activities [17,18]. The aim of this study was to investigate the inhibitory effect of MSKE and its isolated phenolic principles on the enzymatic activities of PLA2, hyaluronidase and LAAO of CR and NK venoms by in vitro tests. The inhibitory effects of MSKE on the activities of these enzymes were also observed by following in vivo tests on the anti-hemorrhagic and anti-dermonecrotic activities of CR and NK venoms, respectively. Molecular docking studies using the Gold v3.2 program were carried out to investigate the binding mode of the extract’s principles to the enzymes.

2. Results and Discussion

2.1. In vitro tests for the inhibition of enzymatic activities

2.1.1. PLA2

CR (17.1 µg) and NK (3.4 µg) venoms produced clear zone diameters of 21.1 ± 0.3 mm and 20.6 ± 0.3 mm, respectively, on agarose-erythrocyte egg yolk gel plate and these clear zones were selected as the minimum indirect hemolytic dose for each venom. From the indirect PLA2 assay, the results revealed that the enzyme from both venoms exhibited indirect hemolytic activity. This was probably due to the activity of PLA2s in the hydrolysis of lecithins to lysolecithins and lysis of the cell membranes of red blood cells [19]. Figure 1 shows the log dose-response curves of MSKE and its phenolic principles (GA, MG and PGG) on the inhibition of PLA2s of CR and NK venoms. It can be seen that MSKE and its phenolic principles clearly showed a dose-dependent inhibitory activity against PLA2s of both venoms. The order of potency as judged by the half-effective dose (ED50) was PGG (25.6 ± 1.3 µg) > MG (42.9 ± 2.3 µg) > MSKE (44.4 ± 0.4 µg) > GA (66.0 ± 2.5 µg) for CR venom, and PGG (30.2 ± 2.4 µg) > MSKE (57.9 ± 0.5 µg) > MG (108.7 ± 2.0 µg) > GA (123.5 ± 2.5 µg) for NK venom (Table 1). The extract, GA, MG and PGG alone and solvent systems had no effect on this study. Of these results implied that the inhibition of PLA2s of both venoms by MSKE may result in the prevention of red blood cell rupture.
Table 1. The constituents and contents of Thai MSKE; the effective dose giving 50% inhibition (ED50) of venom enzymatic activities of MSKE and its three isolated pure compounds. The values represent means ± S.E.M. (n = 2).
Table 1. The constituents and contents of Thai MSKE; the effective dose giving 50% inhibition (ED50) of venom enzymatic activities of MSKE and its three isolated pure compounds. The values represent means ± S.E.M. (n = 2).
MSKE and its isolated principles ContentsPLA2 (ED50, µg)Hyaluronidase (ED50, µg)LAAO (ED50, µg)
(mg/g dry weight)CR venomNK venomCR venomNK venomCR venomNK venom
MSKE-44.4 ± 0.457.9 ± 0.56.0 ± 1.95.0 ± 2.61.5 ± 3.24.0 ± 1.7
GA4.40 ± 0.0566.0 ± 2.5123.5 ± 2.516.0 ± 2.514.8 ± 1.41.7 ± 3.44.5 ± 2.8
MG6.80 ± 0.0242.9 ± 2.3108.7 ± 2.05.9 ± 2.48.6 ± 2.12.3 ± 2.44.4 ± 2.0
PGG612.80 ± 34.8025.6 ± 1.330.2 ± 2.43.4 ± 2.42.6 ± 1.70.9 ± 2.42.2 ± 2.0
Figure 1. Semi-log plot of inhibition of PLA2 activity of CR (17.1 µg) (A) and NK (3.4 µg) (B) venoms by MSKE and its phenolic principles (GA, MG and PGG). Values are mean ± S.E.M. of duplicate experiments.
Figure 1. Semi-log plot of inhibition of PLA2 activity of CR (17.1 µg) (A) and NK (3.4 µg) (B) venoms by MSKE and its phenolic principles (GA, MG and PGG). Values are mean ± S.E.M. of duplicate experiments.
Molecules 14 01404 g001

2.1.2. Hyaluronidase

MSKE and its phenolic principles (GA, MG and PGG) inhibited the activity of hyaluronidase in CR and NK venoms in a dose-dependent manner (Figure 2). The order of potency as judged by the half-effective dose (ED50) was PGG (3.4 ± 2.4 µg) > MG (5.9 ± 2.4 µg) > MSKE (6.0 ± 1.9 µg) > GA (16.0 ± 2.5 µg) for CR venom, and PGG (2.6 ± 1.7 µg) > MSKE (5.0 ± 2.6 µg) > MG (8.6 ± 2.1 µg) > GA (14.8 ± 1.4 µg) for NK venom (Table 1). The inhibition of this enzymatic activity by MSKE may contribute to the prevention of extracellular matrix degradation, and thereby decrease the diffusion of toxins through tissues.
Figure 2. Semi-log plot of inhibition of hyaluronidase activity of CR (250 µg) (A) and NK (250 µg) (B) venoms by MSKE and its phenolic principles (GA, MG and PGG). Values are mean ± S.E.M. of duplicate experiments.
Figure 2. Semi-log plot of inhibition of hyaluronidase activity of CR (250 µg) (A) and NK (250 µg) (B) venoms by MSKE and its phenolic principles (GA, MG and PGG). Values are mean ± S.E.M. of duplicate experiments.
Molecules 14 01404 g002

2.1.3. LAAO

MSKE and its phenolic principles (GA, MG and PGG) inhibited the LAAO activity of both venoms in a dose-dependent manner (Figure 3). The order of potency as judged by the half-effective dose (ED50) was PGG (0.9 ± 2.4 µg) > MSKE (1.5 ± 3.2 µg) > GA (1.7 ± 3.4 µg) > MG (2.3 ± 2.4 µg) for CR venom, and PGG (2.2 ± 2.0 µg) > MSKE (4.0 ± 1.7 µg) > MG (4.4 ± 2.0 µg) > GA (4.5 ± 2.8 µg) for NK venom (Table 1). In this assay, H2O2 generated by the enzyme of both venoms was used by horseradish peroxidase to oxidize o-dianisidine to the radical cation, which was monitored at 436 nm using a spectrophotometer. Therefore, the inhibition of LAAO activity of CR and NK venoms by the MSKE might presumably play a role in the prevention of highly localized concentrations of generated H2O2 in the victim’s tissues, since this snake venom enzyme can induce necrosis and apoptosis in mammalian endothelial cells and consequently cause hemorrhage, possibly through the production of highly localized concentrations of H2O2[10]. The inhibition of LAAO activity by the MSKE may contribute to the prevention of hemorrhage and necrosis.
Figure 3. Semi-log plot of inhibition of LAAO activity of CR (50 µg) (A) and NK (50 µg) (B) venoms by MSKE and its phenolic principles (GA, MG and PGG). Values are mean ± S.E.M. of duplicate experiments.
Figure 3. Semi-log plot of inhibition of LAAO activity of CR (50 µg) (A) and NK (50 µg) (B) venoms by MSKE and its phenolic principles (GA, MG and PGG). Values are mean ± S.E.M. of duplicate experiments.
Molecules 14 01404 g003
It is seen from Figure 1, Figure 2 and Figure 3 and Table 1 that MSKE and its phenolic principles (PGG, MG and GA) exhibited dose-dependent inhibitory effects on PLA2, hyaluronidase and LAAO activities of both venoms. These results indicated that the anti-enzymatic potency of MSKE may be attributed to its major principle (PGG) and other unidentified constituents, since PGG exerted its effect at the lowest ED50 value and had the highest percentage content (61.28%) within the MSKE compared with GA (0.44%) and MG (0.68%). It could be noted that although GA and MG inhibited these enzymatic activities of both venoms, their percentage content within the MSKE were very low. Therefore, GA and MG may only possess a negligible effect compared with PGG. It has been reported that other plant extracts also contain compounds that capable of neutralizing PLA2, hyaluronidase and LAAO activities of snake venoms [20,21,22]. The neutralization of these enzymes not only prevents the early effects of envenomation (local tissue damage, paralysis and hypotension), but also limits the systemic effects by reducing the distribution of toxin.

2.2. In vivo tests

2.2.1. Inhibition of hemorrhagic and dermonecrotic activities

The minimum hemorrhagic dose (MHD) of CR venom was determined to be 20.4 µg/mouse, while the venom at three MHDs-induced hemorrhagic lesions with a diameter of 13.4 ± 0.4 mm. The MSKE significantly (P < 0.05) diminished the diameters of hemorrhagic lesions induced by three MHDs of CR venom in a dose-dependent manner, and the hemorrhagic activity could be completely inhibited at the extract amount of 67 mg/kg mouse (Figure 4A). The extract and solvent systems used did not induce hemorrhagic lesion in mice. The minimum necrotizing dose (MND) of NK venom was determined to be 63 µg/rat, which induced a necrotic lesion with a diameter of 4.7 ± 0.1 mm. Figure 4B shows the percentage inhibition of dermonecrotic activity of NK venom by MSKE which was dose-dependent. The extract with an amount of 0.60 mg/kg rat was effective in completely protecting rats from the dermonecrotic activity of one MND of NK venom (P < 0.05). Solvent and the extract alone did not induce necrotic lesion in rats.
Figure 4. Inhibition of hemorrhagic activity of CR venom (61.2 µg/mouse) in mice (n = 8) (A) and inhibition of dermonecrotic activity of NK venom (63 µg/rat) in rats (n = 4) (B) by MSKE. Diameters of hemorrhagic and necrotic lesions (×). Percentage inhibition of hemorrhagic and necrotic activities (●).
Figure 4. Inhibition of hemorrhagic activity of CR venom (61.2 µg/mouse) in mice (n = 8) (A) and inhibition of dermonecrotic activity of NK venom (63 µg/rat) in rats (n = 4) (B) by MSKE. Diameters of hemorrhagic and necrotic lesions (×). Percentage inhibition of hemorrhagic and necrotic activities (●).
Molecules 14 01404 g004

2.3. Molecular modeling

2.3.1. Docking into PLA2-CR and PLA2-NK structures

The results of molecular modeling using a molecular docking method revealed the possible molecular orientation of MSKE constituents (GA, MG and PGG) in the PLA2 binding pocket of CR and NK venoms (Figure 5 and and Table 2). GoldScores revealed that PGG was bound tighter than GA and MG in their binding pockets of CR and NK venoms. The binding orientations of GA, MG and PGG in the PLA2-CR structure are shown in Figure 5A. Strong H-bonding interactions from GA, MG and PGG were found between a hydrogen atom of the hydroxyl group of ligands and an oxygen atom of the carboxylate group of Asp49. Moreover, common H-bonding interactions with Phe6, Gly30 and Lys61 were shown in the docked GA, MG and PGG. GA had a GoldScore of 37.55 and other H-bonding interactions were also formed with Ser2, Leu3 and Tyr52 ( Figure 5B). In Figure 5C, docked MG, with a GoldScore of 37.49, revealed other H-bonding interactions with Leu3, Glu7, Ile10, Tyr28, Cys29 and Trp31. Docked PGG had a GoldScore of 77.26 and revealed H-bonding interactions with Phe6, Glu7, Ile10, Gly18, Phe19, Tyr22, Ser23, Phe24, Tyr28, Gly30, Trp31, Gly32, Cys45, His48, Asp49, Lys61 and Phe97 (Figure 5D). The docked conformations of GA, MG and PGG in PLA2-NK structure are shown in Figure 6. Their GoldScores were 35.52, 34.75 and 64.92, respectively. The docked conformation of GA, MG and PGG showed common H-bonding interactions with Tyr28, Gly30, Arg31, Gly32, Gly33, Cys45, His48, Asp49 and Tyr64. Additional H-bonding with Ser34, Glu56, Trp62, Pro63 and Phe65 were also formed with the docked conformation of PGG. It can be seen from Figure 5 and Figure 6 that the tannin principles of MSKE could form hydrogen bonds with active site residue His48 Nδ1 and Asp49 Oδ1 of the PLA2 from CR and NK venoms. His48 and Asp49 amino acids are essential for enzymatic activity in the calcium-dependent PLA2s; if the inhibitors are bonded to these amino acids the PLA2 catalytic efficiency can be changed [7,23]. Furthermore, the tannins of MSKE also form important hydrogen bonds with the carbonyl oxygen atoms of the residues involved in calcium coordination (Tyr28, Gly30, Gly32 and Asp49) [23]. Therefore, the results of molecular docking indicated that the phenolic molecules of MSKE could selectively bind to the active sites or their proximity, or to modify conserved residues that are critical to the catalysis of PLA2 enzymes. In addition, a previous study found that MSKE was a natural chelator [17] and it is well known that PLA2s depend on Ca2+ to exert its action. This confirmed the presence of a selective mechanism for MSKE against PLA2 activity in CR and NK venoms.
Table 2. GoldScores of ligands docked into LAAO-CR, PLA2-CR and PLA2-NK structures.
Table 2. GoldScores of ligands docked into LAAO-CR, PLA2-CR and PLA2-NK structures.
GoldScore
MSKE constituentsLAAOPLA2
CR venomCR venomNK venom
GA38.9537.5535.52
MG39.7137.4934.75
PGG69.5877.2664.92
Figure 5. (A) Docked conformation of ligand structures in PLA2-CR structure (GA: Red, MG: Margenta and PGG: Green). (B–D) Distances (in Å) between residues in PLA2-CR binding pocket and ligands: GA (B), MG (C) and PGG (D).
Figure 5. (A) Docked conformation of ligand structures in PLA2-CR structure (GA: Red, MG: Margenta and PGG: Green). (B–D) Distances (in Å) between residues in PLA2-CR binding pocket and ligands: GA (B), MG (C) and PGG (D).
Molecules 14 01404 g005
Figure 6. (A) Docked conformation of ligand structures in PLA2-NK structure (GA: Red, MG: Margenta and PGG: Green). (B-D) Distances (in Å) between residues in PLA2-NK binding pocket and ligands: GA (B), MG (C) and PGG (D).
Figure 6. (A) Docked conformation of ligand structures in PLA2-NK structure (GA: Red, MG: Margenta and PGG: Green). (B-D) Distances (in Å) between residues in PLA2-NK binding pocket and ligands: GA (B), MG (C) and PGG (D).
Molecules 14 01404 g006

2.3.2. Docking into LAAO-CR structure

The possible molecular orientation of the MSKE constituents (GA, MG and PGG) in the LAAO binding pocket of CR venom is shown in Figure 7A. Their GoldScores (Table 2) were 38.95, 39.71 and 69.58, respectively. It was proposed that PGG was bound in the LAAO binding pocket of CR venom tighter than GA and MG. From their docked conformation in LAAO-CR structure, common H-bonding interaction of GA, MG and PGG were found with Arg90 and His463. GA also formed H-bonding interaction with Tyr372, Ile430, Gly464 and Trp465, as shown in Figure 7B. Moreover, a H-π interaction was found between GA and phenyl ring of Trp465 with a distance about 3 Å. In the case of docked MG, it revealed H-bonding interactions with Gly464, Trp465, and FAD and a π-π interaction with Trp465 (Figure 7C). The docked PGG showed H-bonding interactions with Leu207, Asn208, Glu209, Glu219, His223, Asp224, Phe227, Ala228, Arg322, Lys345, Phe354, Tyr356, Tyr372, Ile374, Thr431, Thr432 and FAD (Figure 7D). The interactions found between the amino acid residues of LAAO and tannins in this study agreed well with the crystal structures of LAAO complexed with citrate and with o-aminobenzoate [12]. Pawelek et al. [12] reported the amino acid residues of LAAO were involved in the binding of citrate and o-aminobenzoate; i.e. Leu51, Arg90, Asn208, His223, Asp224, Arg283, Gly313, Tyr372, Ile430 and Gly464. It can be seen from Figure 7 that the tannins of MSKE could form hydrogen bonds with amino acid residues in the substrate-binding domain (Arg90, His223, Asp224, Phe227, Ala228, Lys345, Phe354, Tyr356, Tyr372 and Ile374) of the LAAO from CR venom. These results indicated the presence of a selective mechanism of extract against LAAO activity in CR venom, since a protein database search revealed that the N-terminal sequence of LAAO from NK venom had the highest similarity score to LAAO from CR venom and had similarities to several other snake venom LAAOs [11]. Similar to CR venom, the phenolic compounds of the MSKE could therefore possibly bind to the LAAO binding pocket of NK venom and thereby inhibited the LAAO activity.
Figure 7. (A) Docked conformation of ligand structures in LAAO-CR structure (GA: Red, MG: Margenta and PGG: Green). (B–D) Distances (in Å) between residues in LAAO-CR binding pocket and ligands: GA (B), MG (C) and PGG (D).
Figure 7. (A) Docked conformation of ligand structures in LAAO-CR structure (GA: Red, MG: Margenta and PGG: Green). (B–D) Distances (in Å) between residues in LAAO-CR binding pocket and ligands: GA (B), MG (C) and PGG (D).
Molecules 14 01404 g007
The in vivo inhibition of viper venom-induced hemorrhagic and cobra venom-induced dermonecrotic actions by MSKE is proposed. The inhibition was due firstly to a non-selective mechanism as a result of complexation between the venom proteins and tannins in the extract at the injection site [24] and to a selective mechanism where the active sites of the enzymes (PLA2 and LAAO) of CR and NK venoms were blocked by the phenolic principles of MSKE as demonstrated in the molecular modeling study (Figure 5, Figure 6, and Figure 7). In addition, the MSKE has been reported to demonstrate chelating property with ferrous metal ions [17]; therefore, the authors will investigate further on the potential of MSKE in chelating the venom’s hemorrhagic Zn2+-metalloproteinases since these enzymes presumably act by destroying the collagenous basement membrane and other connective tissue collagen with consequent weakening of blood vessel stability causing the hemorrhagic effect [25]. The molecular docking method of the phenolic principles of the MSKE with the crystal structure of venom metalloproteinases will also be evaluated further.

3. Conclusions

This study shows a beneficial relationship between the in vitro anti-enzymatic activities and in vivo anti-hemorrhagic and anti-dermonecrotic activities of MSKE and its isolated principles against CR and NK venoms. The results were clearly supported by molecular docking studies which indicated that the phenolic principles of MSKE could selectively block the active sites of PLA2 and LAAO of both venoms and PGG were bound tighter than GA and MG. These results imply that the anti-enzymatic potency of MSKE may be attributed by its major principle (PGG) since PGG exerted its effect at the lowest value of ED50 with the highest GoldScores and had the highest percentage content within the MSKE compared with GA and MG.

4. Experimental

4.1. Animals and venoms

Male Swiss albino mice (25–30 g) and male Sprague-Dawley rats (250–300 g) were used for the assays of anti-hemorrhagic and anti-dermonecrotic activities, respectively. Experiments were carried out according to the guidelines set by Mahidol University Animal Care and Use Committee (PY-ACUC). Lyophilized CR and NK venoms were obtained from the Queen Saovabha Memorial Institute, Thai Red Cross Society, Bangkok, Thailand.

4.2. Chemicals

Gallic acid (GA; ≥ 98%) and methyl gallate (MG; ≥ 98%) were purchased from Fluka (Buchs, Switzerland). Pentagalloylglucopyranose (PGG; ≥ 95%) was obtained from Endotherm GmbH (Germany). o-Dianisidine (≥ 97%), horseradish peroxidase (987 U/mg solid), hyaluronic acid (≥ 95%) and L-leucine (≥ 98%) were purchased from Sigma Chemical (St. Louis, MO, USA). All reagents used in the in vitro and in vivo experiments were obtained from commercial sources and were of analytical grade.

4.3. Plant material and extraction

Fully grown unripened Thai mango fruits (Mangifera indica L. cv. ‘Fahlun’) were purchased from a local market. A voucher specimen (R.B. 20007) was deposited at the Museum of Natural Medicine, Faculty of Pharmaceutical Sciences, Chulalongkorn University, Bangkok, Thailand. The ethanolic extract from fresh mango seed kernels was obtained by chopping and homogenizing the seed kernels in a blender, using hot ethanol (80 °C) as an extracting solvent at a ratio of seed kernels/solvent = 1:2 (w/v) for 10 min at room temperature (30 °C) prior to centrifugation at 2000 rpm for 15 min in a Hattich Roto magna®. The extraction was performed three times. The ethanolic extracts were filtered, combined and concentrated in a rotary evaporator (Büchi Rotavapor R-200) at 40 °C. The extract was defatted with hexane, evaporated under reduced pressure, and then freeze-dried to afford a crude mango seed kernel extract (MSKE) with a yield of 9.36% (w/w).

4.4. Standardization

GA, MG and PGG were used as chemical markers and the calibration curves for each compound were obtained by densitometric scanning of different quantities of the chemical marker bands on developed chromatographic plates. An aliquot of the crude extract (MSKE) (8 μL, 25 mg/mL) was applied along with serial amounts of the chemical marker stock solution. The thin-layer chromatographic (TLC) plates were developed in a pre-saturated twin trough glass tank using CHCl3/MeOH/EtOAc/ethyl methyl ketone (6:1.6:2:2) with five drops of formic acid as the mobile phase for GA and MG and CHCl3/EtOH/formic acid (3:5:1) for PGG. The developed TLC plates were scanned at 286 nm and the amount of each compound (GA 4.4, MG 6.8 and PGG 612.8 mg/g dry weight) in MSKE was calculated from the calibration curves [17].

4.5. In vitro tests

4.5.1. Inhibition of PLA2 activity

PLA2 activity in CR and NK venoms was measured using an indirect hemolytic assay on agarose-erythrocyte-egg yolk gel plates to define the minimum indirect hemolytic dose (MIHD) [26], the dose of venom that induced a hemolysis halo of 20 mm diameter after incubation for 20 h at 37 °C. The inhibition of PLA2 activity by MSKE and its isolated phenolic principles was determined against one MIHD of CR and NK venoms. Test solutions and the venoms (0.05 mL each) were pre-incubated for 1 h at 37 °C. After centrifugation at 10,000 rpm for 10 min, the supernatant (20 µL) was tested for PLA2 activity. The anti-PLA2 potential of the extract and its phenolic principles was presented as percent inhibition of the enzymatic activity versus dose, in which 100% inhibition should produce no clear zone. ED50, the effective dose that PLA2 activity of the venoms was reduced by 50%, was calculated.

4.5.2. Inhibition of hyaluronidase activity

Hyaluronidase activity was determined according to the methods of Ferrante [27] and Yingprasertchai et al. [28] with slight modification. CR or NK venom (25-1,000 µg) in PBS (pH 7.4) was incubated with hyaluronic acid (50 μg) in 0.2 M sodium acetate buffer (1 mL, pH 5) that contained 0.15 M NaCl for 25 min at 37 °C. The reaction was stopped by adding 2.5% cetyltrimethylammonium bromide in 2% NaOH (2 mL). The optical density at 400 nm of each sample was read after 30 min with a Shimadzu UV-VIS spectrophotometer UV-160A against blank. Turbidity-reducing activity was expressed as the percentage of the remaining hyaluronic acid, taking the absorbance of a tube to which no enzyme was added as 100%. The venom concentration that produced a reduction in turbidity of ~50% was used in inhibition experiments. The inhibition of venom hyaluronidase activity by MSKE and its isolated phenolic principles was determined after pre-incubating an equal volume (100 µL) of test sample (with different amounts) in acetate buffer with CR or NK venom in PBS (2,500 µg/mL) for 30 min at 37 °C, after which hyaluronidase activity was measured. Blanks containing only the buffer and test sample were run in parallel. The anti-hyaluronidase potential of the extract and its phenolic principles was presented as percent inhibition of the enzymatic activity versus dose. ED50, the effective dose that hyaluronidase activity of the venoms was reduced by 50%, was calculated.

4.5.3. Inhibition of LAAO activity

LAAO activity of CR and NK venoms was determined by using an enzyme-coupled assay [29]. The assay mixture contained 0.025 mL 0.014% peroxidase and 0.2 M triethanolamine buffer (1 mL, pH 7.6) that contained 0.1% L-leucine and 0.0065% o-dianisidine. After incubation at 25 °C for 3 min, 0.1 mL CR or NK venom in PBS (1 mg/mL) was added, and an increase in absorbance at 436 nm (Shimadzu UV-VIS spectrophotometer UV-160A) was measured against blank. Enzymatic activity was calculated from the linear portion of the time–absorbance plot, where one unit of enzymatic activity was defined as the amount of enzyme that caused an increase of 0.001 absorbance units per minute. The percent inhibition of LAAO activity by MSKE and its isolated phenolic principles was evaluated by pre-incubating CR or NK venom (0.1 mL, 1 mg/mL) with test sample in sterile water for injection (SWI) (0.1 mL) for 1 h at 37 °C. After pre-incubation, the mixture was centrifuged at 10,000 rpm for 5 min, and supernatant (0.1 mL) was evaluated for LAAO activity. Venoms, test sample, and solvents alone were used as controls. The anti-LAAO potential of the extract and its phenolic principles was presented as percent inhibition of the enzymatic activity versus dose. ED50, the effective dose that LAAO activity of the venoms was reduced by 50%, was calculated.

4.6. In vivo tests

4.6.1. Inhibition of hemorrhagic activity

Hemorrhagic activity of CR venom was determined according to the method of Kondo et al. [30]. The minimum hemorrhagic dose (MHD) of CR venom was established by injecting 40 µL of different amounts of venom in sterile PBS (pH 7.2) intradermally (i.d.), into two separately marked positions on the shaved dorsal skin of unanesthesized mice (n = 8). After 90 minutes, mice were killed with an overdose of ether, the dorsal skin was removed and the diameter of the hemorrhagic lesion on the inner surface of the skin was measured using calipers and background illumination. The mean diameter of the hemorrhagic lesion was calculated for each venom dose. A dose-response curve for venom dose and hemorrhagic lesion diameter was plotted, and the MHD was recorded as the dose that caused a lesion diameter of 10 mm. The anti-hemorrhagic potential of MSKE was determined against three MHDs of CR venom in mice (n = 8) according to the modified methods of Yingprasertchai et al. [28] and Gutiérrez et al. [31]. Twenty microliters of venom in sterile PBS (3.06 μg/μL) was injected (i.d.) into the shaved dorsal skin of mice, immediately followed by test sample in SWI (20 µL) at the same injection site. After 90 minutes, mice were killed by ether overdose. Dorsal skin was removed and the diameter of the hemorrhagic lesion on the inner surface of the skin was measured. Sterile PBS, SWI, venom and MSKE alone were used as controls.

4.6.2. Inhibition of dermonecrotic activity

The methods of Kondo et al. [30] and Theakston and Reid [32] were used to determine the minimum necrotizing dose (MND) of NK venom in rats (n = 4) which is defined as the smallest amount of venom that could cause a necrotic lesion of 5 mm in diameter on the inner dorsal skin of the rat 72 h after intradermal injection. The anti-dermonecrotic potential of MSKE was determined against one MND of NK venom. The venom in sterile PBS (40 µL, 1.57 μg/μL) was injected i.d. into the shaved dorsal skin of rats (n = 4), immediately followed by test sample (40 µL) at the same injection site. Seventy-two hours after injection, rats were killed by ether overdose. The dorsal skin was removed, and the necrotic lesion on the inner surface of the skin was measured. Sterile PBS, venom, and MSKE alone were used as controls.

4.7. Molecular modeling

Structures of gallic acid (GA), methyl gallate (MG) and 1,2,3,4,6-penta-O-galloyl-β-d-glucopyranose (PGG) (Figure 8) were constructed and optimized at the HF/3-21G level of theory using Gaussian 03 program [33]. For LAAO structure, the X-ray structure of LAAO from Calloselasma rhodostoma complexed with citrate was obtained from the Protein Data Bank (PDB code 1f8r) [12]. Protein structure and Flavin Adenine Dinucleotide (FAD) were selected and added hydrogen atoms by using SYBYL version 7.2 program (TRIPOS, Assoc., Inc.: St. Louis, MO). The three-dimensional structures of PLA2 from CR and NK venoms were built by homology modeling technique using Geno3D web server [34]. Geno3D web server used distance geometry, simulated annealing and energy minimization algorithms to build the protein 3D model. After that, the structures were checked the quality of geometry by using PROCHECK [35]. For homology modeling of PLA2 structures from CR venom, the sequence of CRV-S1E6a [36] was based on the structure of PLA­2­ from venom of Agkistridon piscivorus piscivorus (PDB code 1vap) which has the 65% sequence identity. Homology modeling of amino acid sequences of PLA­2­ (CM-III) from NK venom [37] was based on the structure of PLA­2­ from Naja naja naja venom (PDB code 1a3d), which had a high sequence identity of 90%. Hydrogen atoms were added to all structures by using the SYBYL 7.2 program.
Figure 8. Chemical structures of three major constituents of Thai mango seed kernel extract: MG (1), GA (2) and PGG (3).
Figure 8. Chemical structures of three major constituents of Thai mango seed kernel extract: MG (1), GA (2) and PGG (3).
Molecules 14 01404 g008
The molecular docking method was done using Gold version 3.2 program [38] to investigate the binding orientation of GA, MG and PGG in LAAO from CR venom, and PLA2 from CR and NK venoms. The radius of the binding site was set to 10 Å. For LAAO structure, the binding site was defined by using the center of mass of citrate which is in the active site. For both PLA2 structures, the possible binding sites were predicted by using PASS (Putative Active Sites with Spheres) program [39]. PASS identified the positions likely to represent binding sites based upon the size, shape and burial extent of these volumes. The most possible position representing the binding site from PASS, as called calcium-binding loop in the PLA2 structure [40], was used to define the binding site in Gold program. For setting the genetic algorithm parameters, the default parameters of automatic settings were used. The GoldScore fitness function was used to determine the fitness score. The GoldScore fitness function was used to determine the fitness score, as shown in GoldScore.

4.8. Statistical analysis

Results were expressed as mean ± S.E.M. for two determinations. Statistical analysis was carried out using SPSS 13.0 for Windows. Significant differences (P < 0.05) between means were assessed by one-way ANOVA, followed by Tukey’s honesty significant difference test or Dunnett’s T3 test for multiple comparisons.

Acknowledgements

Financial support from the Thailand Research Fund (TRF) through the Royal Golden Jubilee Ph.D. Program (Grant No. PHD/0049/2547) is gratefully acknowledged. We thank Ms. Narumol Pakmanee (Department of Research and Development, the Queen Saovabha Memorial Institute, Bangkok, Thailand) for providing lyophilized CR and NK venoms.

References

  1. Viravan, C.; Veeravat, U.; Warrell, M.J.; Theakston, R.D.; Warrell, D.A. ELISA confirmation of acute and past envenoming by the monocellate Thai cobra (Naja kaouthia). Am. J. Trop. Med. Hyg. 1986, 35, 173–181. [Google Scholar]
  2. Songsumard, S. Snakebite in Thailand. Annual Epidemiological Surveillance Report; Ministry of Public Health: Bangkok, 1995; pp. 324–333. [Google Scholar]
  3. Markland, F.S. Snake venoms. Drugs 1997, 54, 1–10. [Google Scholar] [CrossRef]
  4. Aird, S.D. Ophidian envenomation strategies and the role of purines. Toxicon 2002, 40, 335–393. [Google Scholar] [CrossRef]
  5. Mukherjee, A.K.; Ghosal, S.K.; Maity, C.R. Effect of oral supplementation of vitamin E on the hemolysis and erythrocyte phospholipid-splitting action of cobra and viper venoms. Toxicon 1998, 36, 657–664. [Google Scholar] [CrossRef]
  6. Kini, R.M.; Evans, H.J. Structure-function relationships of phospholipases. The anticoagulant region of phospholipases A2. J. Biol. Chem. 1987, 262, 14402–14407. [Google Scholar]
  7. Berg, O.G.; Gelb, M.H.; Tsai, M.D.; Jain, M.K. Interfacial Enzymology: The Secreted Phospholipase A2-Paradigm. Chem. Rev. 2001, 101, 2613–2653. [Google Scholar] [CrossRef]
  8. Girish, K.S.; Shashidharamurthy, R.; Nagaraju, S.; Gowda, T.V.; Kemparaju, K. Isolation and characterization of hyaluronidase a “Spreading factor” from Indian cobra (Naja naja) venom. Biochimie 2004, 86, 193–202. [Google Scholar] [CrossRef]
  9. Du, X.Y.; Clemetson, K.J. Snake venom L-amino acid oxidases. Toxicon 2002, 40, 659–665. [Google Scholar] [CrossRef]
  10. Suhr, S.M.; Kim, D.S. Identification of the snake venom substance that induces apoptosis. Biochem. Biophys. Res. Commun. 1996, 224, 134–139. [Google Scholar] [CrossRef]
  11. Sakurai, Y.; Takatsuka, H.; Yoshioka, A.; Matsui, T.; Suzuki, M.; Titani, K.; Fujimura, Y. Inhibition of human platelet aggregation by L-amino acid oxidase purified from Naja naja kaouthia venom. Toxicon 2001, 39, 1827–1833. [Google Scholar] [CrossRef]
  12. Pawelek, P.D.; Cheah, J.; Coulombe, R.; Macheroux, P.; Ghisla, S.; Vrielink, A. The structure of L-amino acid oxidase reveals the substrate trajectory into an enantiomerically conserved active site. EMBO J. 2000, 19, 4204–4215. [Google Scholar] [CrossRef]
  13. Okonogi, T.; Hattori, Z.; Ogiso, A.; Mitsui, S. Detoxification by persimmon tannin of snake venoms and bacterial toxins. Toxicon 1979, 17, 524–527. [Google Scholar] [CrossRef]
  14. Abubakar, M.S.; Sule, M.I.; Pateh, U.U.; Abdurahman, E.M.; Haruna, A.K.; Jahun, B.M. In vitro snake venom detoxifying action of the leaf extract of Guiera senegalensis. J. Ethnopharmacol. 2000, 69, 253–257. [Google Scholar] [CrossRef]
  15. Borges, M.H.; Alves, D.L.F.; Raslan, D.S.; Piló-Veloso, D.; Rodrigues, V.M.; Homsi-Brandeburgo, M.I.; De Lima, M.E. Neutralizing properties of Musa paradisiaca L. (Musaceae) juice on phospholipase A2, myotoxic, hemorrhagic and lethal activities of crotalidae venoms. J. Ethnopharmacol. 98, 21–29.
  16. Pithayanukul, P.; Ruenraroengsak, P.; Bavovada, R.; Pakmanee, N.; Suttisri, R.; Saen-oon, S. Inhibition of Naja kaouthia venom activities by plant polyphenols. J. Ethnopharmacol. 2005, 97, 527–533. [Google Scholar] [CrossRef]
  17. Nithitanakool, S.; Pithayanukul, P.; Bavovada, R. Antioxidant and hepatoprotective activities of Thai mango seed kernel extract. Planta Med. (online March 26, 2009).
  18. Nithitanakool, S.; Pithayanukul, P.; Bavovada, R.; Saparpakorn, P. Molecular docking studies and anti-tyrosinase activity of Thai mango seed kernel extract. Molecules 2009, 14, 257–265. [Google Scholar] [CrossRef]
  19. Ketelhut, D.F.J.; Mello, M.H.; Veronese, E.L.G.; Esmeraldino, L.E.; Murakami, M.T.; Arni, R.K.; Giglio, J.R.; Cintra, A.C.O.; Sampaio, S.V. Isolation, characterization and biological activity of acidic phospholipase A2 isoforms from Bothrops jararacussu snake venom. Biochimie 2003, 85, 983–991. [Google Scholar] [CrossRef]
  20. Soares, A.M.; Ticli, F.K.; Marcussi, S.; Lourenco, M.V.; Januario, A.H.; Sampaio, S.V.; Giglio, J.R.; Lomonte, B.; Pereira, P.S. Medicinal plants with inhibitory properties against snake venoms. Curr. Med. Chem. 2005, 12, 2625–2641. [Google Scholar] [CrossRef]
  21. Ushanandini, S.; Nagaraju, S.; Kumar, K.H.; Vedavathi, M.; Machiah, D.K.; Kemparaju, K.; Vishwanath, B.S.; Gowda, T.V.; Girish, K.S. The anti-snake venom properties of Tamarindus indica (leguminosae) seed extract. Phytother. Res. 2006, 20, 851–858. [Google Scholar] [CrossRef]
  22. Leanpolchareanchai, J.; Pithayanukul, P.; Bavovada, R. Anti-necrosis potential of polyphenols against snake venoms. Immunopharmacol. Immunotoxicol. 2009, (in press). [Google Scholar]
  23. Arni, R.K.; Ward, R.J. Phospholipase A2—a structural review. Toxicon 1996, 34, 827–841. [Google Scholar] [CrossRef]
  24. Haslam, E. Natural polyphenols (Vegetable tannins) as drugs: Possible modes of action. J. Nat. Prod. 1996, 59, 205–215. [Google Scholar] [CrossRef]
  25. Bjarnason, J.B.; Fox, J.W. Hemorrhagic toxins from snake venoms. J. Toxicol. Toxin Rev. 1988, 7, 121–209. [Google Scholar] [CrossRef]
  26. Gutiérrez, J.M.; Avila, C.; Rojas, E.; Cerdas, L. An alternative in vitro method for testing the potency of the polyvalent antivenom produced in Costa Rica. Toxicon 1988, 26, 411–413. [Google Scholar] [CrossRef]
  27. Ferrante, N.D. Turbidimetric measurement of acid mucopolysaccharides and hyaluronidase activity. J. Biol. Chem. 1956, 220, 303–306. [Google Scholar]
  28. Yingprasertchai, S.; Bunyasrisawat, S.; Ratanabanangkoon, K. Hyaluronidase inhibitors (sodium cromoglycate and sodium auro-thiomalate) reduce the local tissue damage and prolong the survival time of mice injected with Naja kaouthia and Calloselasma rhodostoma venoms. Toxicon 2003, 42, 635–646. [Google Scholar] [CrossRef]
  29. Tan, N.H.; Tan, C.S. A comparative study of cobra (Naja) venom enzymes. Comp. Biochem. Physiol. B 1988, 90, 745–750. [Google Scholar]
  30. Kondo, H.; Kondo, S.; Ikezawa, H.; Murata, R.; Ohsaka, A. Studies on the quantitative method for determination of hemorrhagic activity of Habu snake venom. Jpn. J. Med. Sci. Biol. 1960, 13, 43–51. [Google Scholar]
  31. Gutiérrez, J.M.; Gené, J.A.; Rojas, G.; Cerdas, L. Neutralization of proteolytic and hemorrhagic activities of Costa Rican snake venoms by a polyvalent antivenom. Toxicon 1985, 23, 887–893. [Google Scholar] [CrossRef]
  32. Theakston, R.D.G.; Reid, H.A. Development of simple standard assay procedures for the characterization of snake venoms. Bull. World Health Organ. 1983, 61, 949–956. [Google Scholar]
  33. Frisch, M.J.; Trucks, G.W.; Schlegel, H.R.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Montgomery, Jr. J.A.; Vreven, T.; Kudin, K.N.; Burant, J.C.; Millam, J.M.; Iyengar, S.S.; Tomasi, J.; Barone, V.; Mennucci, B.; Cossi, M.; Scalmani, G.; Rega, N.; Petersson, G.A.; Nakatsuji, H.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Klene, M.; Li, X.; Knox, J.E.; Hratchian, H.P.; Cross, J.B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R.E.; Yazyev, O.; Austin, A.J.; Cammi, R.; Pomelli, C.; Ochterski, J.W.; Ayala, P.Y.; Morokuma, K.; Voth, G.A.; Salvador, P.; Dannenberg, J.J.; Zakrzewski, V.G.; Dapprich, S.; Daniels, A.D.; Strain, M.C.; Farkas, O.; Malick, D.K.; Rabuck, A.D.; Raghavachari, K.; Foresman, J.B.; Ortiz, J.V.; Cui, Q.; Baboul, A.G.; Clifford, S.; Cioslowski, J.; Stefanov, B.B.; Liu, G.; Liashenko, A.; Piskorz, P.; Komaromi, I.; Martin, R.L.; Fox, D.J.; Keith, T.; Al-Laham, M.A.; Peng, C.Y.; Nanayakkara, A.; Challacombe, M.; Gill, P.M.W.; Johnson, B.; Chen, W.; Wong, M.W.; Gonzalez, C.; Pople, J.A. Gaussian03; Gaussian, Inc.: Pittsburgh, PA, USA, 2003. [Google Scholar]
  34. Combet, C.; Jambon, M.; Deléage, G.; Geourjon, C. Geno3D: Automatic comparative molecular modelling of protein. Bioinformatics 2002, 18, 213–214. [Google Scholar] [CrossRef]
  35. Laskowski, R.A.; MacArthur, M.W.; Moss, D.S.; Thornton, J.M. PROCHECK: A program to check the stereochemical quality of protein structures. J. Appl. Crystallogr. 1993, 26, 283–291. [Google Scholar] [CrossRef]
  36. Tsai, I.H.; Wang, Y.M.; Au, L.C.; Ko, T.P.; Chen, Y.H.; Chu, Y.F. Phospholipases A2 from Callosellasma rhodostoma venom gland cloning and sequencing of 10 of the cDNAs, three-dimensional modelling and chemical modification of the major isozyme. Eur. J. Biochem. 2000, 267, 6684–6691. [Google Scholar] [CrossRef]
  37. Joubert, F.J.; Taljaard, N. Purification, some properties and amino-acid sequences of two phospholipases A (CM-II and CM-III) from Naja naja kaouthia venom. Eur. J. Biochem. 1980, 112, 493–499. [Google Scholar] [CrossRef]
  38. Jones, G.; Willett, P.; Glen, R.C.; Leach, A.R.; Taylor, R. Development and validation of a genetic algorithm for flexible docking. J. Mol. Biol. 1997, 267, 727–748. [Google Scholar]
  39. Brady, G.P.; Stouten, P.F.W. Fast prediction and visualization of protein binding pockets with PASS. J. Comput. Aided Mol. Des. 2000, 14, 383–401. [Google Scholar] [CrossRef]
  40. Gu, L.; Wang, Z.; Song, S.; Shu, Y.; Lin, Z. Crystal structures of an acidic phospholipase A2 from the venom of Naja kaouthia. Toxicon 2002, 40, 917–922. [Google Scholar] [CrossRef]
  • Sample Availability: Samples are available from authors.

Share and Cite

MDPI and ACS Style

Leanpolchareanchai, J.; Pithayanukul, P.; Bavovada, R.; Saparpakorn, P. Molecular Docking Studies and Anti-enzymatic Activities of Thai Mango Seed Kernel Extract Against Snake Venoms. Molecules 2009, 14, 1404-1422. https://doi.org/10.3390/molecules14041404

AMA Style

Leanpolchareanchai J, Pithayanukul P, Bavovada R, Saparpakorn P. Molecular Docking Studies and Anti-enzymatic Activities of Thai Mango Seed Kernel Extract Against Snake Venoms. Molecules. 2009; 14(4):1404-1422. https://doi.org/10.3390/molecules14041404

Chicago/Turabian Style

Leanpolchareanchai, Jiraporn, Pimolpan Pithayanukul, Rapepol Bavovada, and Patchreenart Saparpakorn. 2009. "Molecular Docking Studies and Anti-enzymatic Activities of Thai Mango Seed Kernel Extract Against Snake Venoms" Molecules 14, no. 4: 1404-1422. https://doi.org/10.3390/molecules14041404

APA Style

Leanpolchareanchai, J., Pithayanukul, P., Bavovada, R., & Saparpakorn, P. (2009). Molecular Docking Studies and Anti-enzymatic Activities of Thai Mango Seed Kernel Extract Against Snake Venoms. Molecules, 14(4), 1404-1422. https://doi.org/10.3390/molecules14041404

Article Metrics

Back to TopTop