INTRODUCTION
Spices were known for their medicinal value right after they were established for use in the society (Duke, 2002). India is known for a wide variety of spices due to its diverse climatic conditions. Amongst this pepper is prominent in India owing to its biodiversity. There are different species of pepper found in India such as Piper nigrum Linn. (black pepper), Piper longum Linn. (long pepper), Piper cubeba Linn. (tailed pepper) etc. all belonging to the family Piperaceae. Piper nigrum is native to Malabar coast of India as the vine requires tropical hot and humid climate with optimum temperature of 28oC and high rainfall of about 125-200cm for its growth (Prange, 2010) . Piper longum is native to Indo-Malayan realm which includes different regions of India such as Kerala, Tamil Nadu, Assam, Madhya Pradesh, Meghalaya and many more. It requires a hot and humid climate, porous soil that is well drained and rich in organic matter with partial shade for its growth (Satyavati, 2021) . Piper cubeba on the other hand is a native of Indonesia that also grows wildly in the foothills of Himalayas in India, where it finds well drained soil with high humidity and optimum temperature (30-32oC) and rainfall (150cm).
Pepper has been widely used in AYUSH drug formulations (Ayurveda, Unani and Siddha medicine) due to its potential medicinal properties (Ahmad et al., 2012) . Piper longum is the principal ingredient of some Ayurvedic medicines such as Trikatu churnam, Chavanprash and Amrit kalash. Trikatu churnam, a combination of Piper nigrum L. (black pepper), Piper longum L. (long pepper) and Zingiber officinale (ginger) is the major ingredient in more than 1500 Ayurvedic/Siddha formulations for treating various ailments (Committee, 1978) . Trikatu is known to reduce Vata and Kapha and increase Pitta. Its Unani properties include control of temperament (Mizaj) (Kaushik, Jain, Khan, & Rai, 2018) . A range of biological activities exhibited is due to a wide spectrum of bioactive molecules present in the seeds, roots, stem and leaves of different pepper species (Mgbeahuruike, Yrjönen, Vuorela, & Holm, 2017) . Medicinal use of a plant derived drugs in modern medicine is based on a single active compound in the purified form. Several bioactive compounds isolated and characterised from different varieties of pepper finds its use in various disorders afflicting humans. Herein, in the present study the bioactive molecule piperine in three Piper species commonly available in South India were determined using HPLC analysis. Further, since, inflammation is the major pathology involved in most of the disease conditions and since pepper is traditionally used in managing respiratory inflammation, we ventured to identify its anti-inflammatory potential by molecular docking.
MATERIALS AND METHODS
Chemicals
Analytical grade solvents from Sisco Research Laboratories (SRL), India. were used for phytochemical extraction. All other chemicals used were obtained from Himedia (India) and Sigma-Aldrich (St Louis, MO), unless otherwise specified.
Preparation of the extracts
Black pepper (Piper nigrum), long pepper (Piper longum) and tailed pepper (Piper cubeba) were purchased from Pondicherry market and they were verified by a botanist and voucher specimens were stored in the Natural Products Research Laboratory, Department of Biotechnology, Pondicherry University. Fifty gram of each sample was ground using mortar and pestle and then they were transferred to a conical flask and dissolved in 100% methanol. The conical flask was cotton plugged and covered using an aluminium foil to avoid oxidation by light and extracted by cold percolation method by shaking at 100 rpm for 48 hours. The supernatant was then collected in an amber coloured bottle using a Whatman filter paper. Methanol was removed using a rotary vacuum evaporator and the extracts were stored at 4°C until further use. The supernatant was then concentrated using a rotavapour at 30°C and 350 mbar. Concentrated sample was stored in a glass bottle at 4ºC until further use.
Phytochemical Analysis of the Extracts
Preliminary qualitative phytochemical screening of the three varieties of the methanolic extract of the pepper species were checked for the presence of the primary and secondary metabolites such as: carbohydrates, proteins, steroids, glycosides, alkaloids, phenolics, flavonoids, tannins and saponins (Vasanthi et al., 2009). Further, the phenolic, flavonoid and alkaloid content were quantified spectrophotometrically (Harborne, 1984).
Chromatographic Analysis of the Pepper varieties
High Performance Thin Layer Chromatography (HPTLC) fingerprinting of the pepper extracts
Fifty milligram of each pepper extract was dissolved in 1 mL of methanol by vortexing and filtered using a syringe filter. The precoated silica gel 60F254 aluminium sheets (20 × 10 cm) used as stationary phase was activated in an oven at 50oC for 15 minutes. Toluene: ethyl acetate: formic acid (3:2:0.5 v/v/v) served as the mobile phase. The chromatographic conditions maintained includes: 30 min of plate and chamber saturation time; 70 mm of migration distance allowed; wavelength scanning done at 254 and 366 nm keeping the slit dimension at 10×0.6 mm. Ten microlitres of each of the extracts of Piper longum, Piper nigrum, Piper cubeba (5mg/ml) and the phytoconstituents such as Piperine, Stigmasterol, Cinnamic acid, Eugenol, Isoeugenol, Quercetin and Gallic acid (1mg/mL) dissolved in methanol were spotted on the TLC plate and developed at constant temperature. Photometric measurements were performed at 254 and 366nm in reflectance mode using CAMAG-HPTLC (Switzerland) scanner and WINCATS software.
High performance liquid chromatography (HPLC quantification of piperine in the pepper extracts
The pepper extract was subjected to HPLC analysis by using Shimadzu UFLC equipped with UV detector. Quantification of Piperine in Piper longum, Piper nigrum and Piper cubeba against the standard Piperine was performed by Reverse Phase High Performance Liquid Chromatography (RP-HPLC) following the method of (Upadhyay, Sharma, & Joshi, 2013) . C18-250X4.60 mm 5 μm, 100 Å column was used to separate piperine. Acetonitrile and HPLC grade water were used as the mobile phase in the ratio (60:40). Isocratic elution was performed with 60% of acetonitrile with a run time of 30 minutes. A flow rate of 1 ml/min was preserved throughtout the run and the phytomolecule was detected using an UV–detector t at 254 nm. The standard concentration of piperine was made at 2 mg/mL and its injection volume was kept at 25μl. The sample concentration of Piper nigrum and Piper longum were made at 62.5 ng/mL and for Piper cubeba, it was made at 250 ng/mL. The injection volume for all the pepper varieties were kept at 25μl. Temperature was maintained at 25°C.
Protein Preparation
Four target proteins viz., TNFα, NF-κB, COX-1, and COX-2 that are well known for mediating inflammatory response were selected for the present study. Three-dimensional (3D) conformers of the target proteins of interest TNFα (PDB code: 2AZ5), NF-κB (PDB code: 3GUT), COX-1 (PDB code: 5WBE) and COX-2 (PDB code: 1CX2) were downloaded from www.rcsb.org as Protein Data Bank (PDB) format, one of the primary archival resources in areas of structural biology for further in-silico docking analysis (Burley, Berman, & Bhikadiya, 2019). Using PubMed (http://www.ncbi.nlm.nih.gov/pubmed), all the relevant literatures required for the present study were collected. All protein structures were visualized in PYMOL to check for structural integrity and unwanted amino acid chains were removed to avoid overlapping results derived from homo, di / tri / tetramer proteins. Polar hydrogens were then added to the structure in Auto dock tools along with addition of Kollman charges. Atoms were further assigned as AD4 type and the PDB file was saved in. pdbqt format. The proteins such as TNFα (PDB code: 2AZ5), NF-κB (PDB code: 3GUT), COX-1 (PDB code: 5WBE) and COX-2 (PDB code: 1CX2) were then analysed by checking the intactness of proteins and by removing the extra chains associated with the protein. Simultaneously, polar hydrogen bonds were added to the protein molecules along with Kollman charges. Then the proteins were assigned AD4 type atoms to strengthen the molecule. All the files were then saved in pdbqt format to prepare them for interacting with the ligands of interest.
Ligand Preparation
The bioactive phytocompound in pepper which is piperine was the ligand of interest to carry out this study. Piperine was chosen as it is the main bioactive compound in pepper and known for its pharmacological effects (Derosa, Maffioli, & Sahebkar, 2016; Gorgani, Mohammadi, Najafpour, & Nikzad, 2017) . Accordingly, the structures of the ligand were drawn using ChemBioDraw, and the 2D conformation of the ligand was subjected to 3D-optimisation by ChemSketch, one of the best chemical drawing software and was converted to 3D conformer. Subsequently, Chimera software was used for converting ligand and protein files from .sdf to .pdb format. Ligands (.pdb) were opened in autodock tools and were saved in.pdbqt format after detecting the route.
Docking
Auto dock tools (ADT) v1.5.6 and Autodock v4.2 programs were used to carry out docking analysis (Morris, Ruth, & Lindstrom, 2009). To run Auto dock 1.5.6, a searching grid is generally used which would be extended over the selected amino acids in the receptor protein along with the addition of polar hydrogens to the ligand moieties, atomic salvation parameters and assignment of Kollman charges. Gasteiger-type polar hydrogen charges were assigned and carbons were merged with the nonpolar hydrogens and the internal degrees of freedom and torsions were established. Thereafter, grid parameter files were saved in .gpf format followed by autogrid run. Computational binding or docking took place in between the bioactive molecule piperine and the respective receptor proteins wherein the ligands were considered as flexible and the macromolecule was regarded as a rigid body. Lamarckian Genetic Algorithm was applied to carry out the search and 20 runs were performed for this study. Gromos force field was used to minimize energy. The whole receptor protein was used for blind docking. Affinity maps for all the atom types present, as well as an electrostatic map were computed with a grid spacing of 0.375 Å. Search and docking parameters were selected before running the final auto dock. After accepting docking parameters, output files were saved as .dpf (docking parameter file) format followed by selection of program path name, parameter file and autodock was run. After completion of docking, the resulted files were saved as .dlg (docking log file) format. The final projection of the interaction and docking pose analysis was done by LigPlus, a software used to visualize the environment of a ligand either by downloading from the PDB ligand-protein complex or by extracting a docked pose from a docking run complexed with the target protein. To perform this, .dlg file was opened in auto dock and different docking pose conformation were analysed using conformation play option in auto dock and hydrogen bonds were also built. The docked protein- ligand complex was then saved as .pdb format which was further analysed using LigPlot+ tool. Binding energies of different generated protein-ligand complex were also recorded by LigPlus. Finally, the final .pdb complex was run using LigPlus to show the active amino acid sites and hydrogen bonds between the protein and ligand. The different complexes were then sorted with respect to the predicted binding energy and analysed.
RESULTS AND DISCUSSION
Pepper has been used for centuries for its culinary and medicinal properties. It has been used in Ayurveda, Siddha, and traditional Chinese medicine due to the wide range of bio-activities they exhibit such as anti-microbial, anti-diabetic, anti-inflammatory, gastro-intestinal stimulant, anti-asthmatic, anti-carcinogenic, anti-hyperlipidemic, anti-diarrheal and anti-oxidant activities (Saxena, Saxena, Nema, Singh, & Gupta, 2022). Further, it is well known as a bioenhancer in most AYUSH preparations and is now being used in developing nutraceutical supplements/formulations enhancing the bioactivity of some lead molecules. The reason behind their diverse pharmacological activity is the presence of multiple bioactive compounds in pepper. However, to further refine and understand the intricacies of their mechanism, an in-depth study of their bioactive molecules is required. In the present study the variations in the phytochemical constituents in three different varieties of pepper commonly available in South India were analysed and the bioactive molecule piperine was quantified and checked for its anti-inflammatory potential using in-silico techniques.
Table 1
Phytochemical profiling of the Pepper Extracts
The phytochemicals in various medicinal plants have been studied extensively for their role in preventing disease and promoting health by establishing their efficacy and understanding their probable mode of action (Sangeetha, Priya, & Vasanthi, 2013; Sundaram, Vasanthi, & R, 2019). Accordingly, the preliminary qualitative phytochemical analysis revealed the presence of the primary metabolites such as carbohydrates and reducing sugars in the form of glycosides, amino acids and proteins and steroids. The secondary metabolites include the polyphenols, flavonoids, saponins and alkaloids. It was observed that all the three pepper varieties showed consistent results throughout the test with only slight variation in the final colour intensities in the respective tests; Piper nigrum showing the highest intensity followed by Piper longum and the least intense being Piper cubeba (Table 1). The variation in the intensity owes to the presence of different levels of phytochemicals in the different species of the genus piper. Therefore, we further quantified the total phenolic content, flavonoids and total alkaloid content spectrophotometrically. It is interesting to note that although the polyphenols and flavonoids in general exhibit various pharmacological potential, the alkaloid content was higher than the polyphenols and flavonoids in all the three species (Table 2). Among the three species tested, the alkaloid content was found to be maximum in Piper longum (1.478 mg/g) followed by Piper nigrum (0.601 mg/g) and Piper cubeba (0.316 mg/g). It has been reported that the presence of alkaloids could express many pharmacological activities such as antihypertensive effects, antiarrhythmic effects, and anticancer actions (Alencar, Oldoni, & Castro, 2007). The above results confirms that there is variation in the secondary metabolite content in the three pepper varieties tested. The differences identified is not only because they belong to different species but also due to the difference in the geographical and climatic conditions where they grow (Zou & Zou, 2021).
Table 2
Subsequently, the TLC elution profiles of the three pepper extracts viz. Piper longum, Piper nigrum, and Piper cubeba were visualized under near UV (254 nm) & far UV (366 nm) (Figure 1). At least 15 distinguishable bands were observed for both Piper longum and Piper nigrum whereas Piper cubeba showed only 8 such distinguishable bands. It was also interesting to note that the intensity of the bands was quite high for both Piper longum and Piper nigrum whereas Piper cubeba showed faint bands. This result concluded that Piper nigrum and Piper longum have a variety of secondary metabolites in higher concentration as compared to Piper cubeba which gave few distinguishable bands of low intensity.
Further, as TLC is a crude method used for initial band visualization, High-performance thin-layer chromatography (HPTLC) profiling provided a much more refined and detailed analysis.
The relative front values of the individual bands in each of the extract lanes viz. Piper longum (Pl); 5 mg/ml, Piper nigrum (Pn); 5 mg/ml and Piper cubeba (Pc); 5 mg/ml were matched against phytochemical standards viz., Piperine (P); 1 mg/ml, Stigmasterol (St); 1 mg/ml, Cinnamic acid (Ca); 1 mg/ml, Eugenol (Eu); 5 mg/ml, Isoeugenol (IEu); 1 mg/ml, Quercetin (Q); 1 mg/ml and Gallic acid (Ga); 1 mg/ml were visualized in both near UV (254 nm) and far UV (366 nm) (Figure 2). It was observed that only the standards - Piperine, Eugenol, Isoeugenol, and Quercetin matched against their respective bands present in the extracts. Among the three extracts, Piper cubeba showed the least intense bands when compared with Piper longum and Piper nigrum. The results implied that the current solvent system was able to simultaneously elute and resolve discrete bands of only Piperine, Eugenol, Isoeugenol, and Quercetin along with Piper longum, Piper nigrum, and Piper cubeba. Also, from the intensity of the bands in the three Piper extracts, it was concluded that the individual components present in them were in a higher concentration in Piper longum and Piper nigrum but not for Piper cubeba which was noticed in our TLC analysis also.
It is well known that high pressure liquid chromatography (HPLC) is more sensitive than any other analytical technique for qualitative and quantitative analysis of various compounds (Blum, 2014). However, method development and optimization of sensitivity to detect analytes in various samples is still a challenge (Upadhyay et al., 2013). In the present study, the total run time of piperine was found to be 15 minutes, and the piperine appeared on chromatogram at 7.105 minutes for Piper nigrum extract, at 7.060 minutes for Piper longum extract, at 7.076 minutes for Piper cubeba extract (Figure 3). The retention time of the reference standard, piperine (Figure 3) was observed at 7.099 minutes. It was also noted that apart from the piperine peaks all the pepper extracts showed small peaks; Piper nigrum showed a small peak at 6.398 minutes, Piper longum showed a small peak at 6.357 minutes and Piper cubeba showed a small peak at 5.121 minutes. The amount of piperine present in the three varieties of pepper is computed in (Table 3). The results confirm that the current HPLC method was a stable, convenient, and rapid method for determining the purity and concentration of piperine. Throughout the experiment chromatographic traces were completely clean and none exhibited any type of interference throughout the entire running time confirming the sensitivity and specificity of the method used in quantification.
Table 3
Modern medicine is based upon several parameters such as identification and isolation of a particular molecule and its subsequent synthetic modification to achieve a greater degree of efficacy or potency which includes its dosage and toxicity profiling before it reaches the market. Piperine or piperoylpiperidine constitutes about 98% of the total alkaloid found in pepper, it is hydrophobic and interacts with different cellular proteins and factors to either inhibit them or activate them (Meghwal & Goswami, 2012). It also has bio-enhancing activity wherein, it increases the efficacy of certain drugs by increasing its bioavailability either by altering the membrane dynamics, or inhibiting the gastrointestinal and hepatic metabolism (Peterson et al., 2019). The amount of piperine present in Piper nigrum was 235.05 µg/mL and Piper longum accounted for a piperine content of 268.50 µg/mL whereas, Piper cubeba possessed only 8.56 µg/mL of piperine. Our results corroborate to an earlier report on the piperine content in black pepper (Hamrapurkar, Jadhav, Kundnani, & Zine, 2011). Piperine exists in three more isomeric states namely isopiperine (cis–trans isomer), chavicine (cis–cis isomer), and isochavicine (trans–cis isomer) of which only the former is known to contribute to the pungency of black pepper. Apart from the piperine peaks, all the pepper extracts showed small peaks with different retention times which is probably due to the isomeric forms or other piperine derivatives such as piperamide, piperamine, etc.
In silico Docking Analysis:
The present work aimed to evaluate the role of piperine on target proteins that modulate inflammation as it is one of the key pathologies involved in most chronic disorders and since pepper is traditionally used in managing respiratory inflammation (Rehman, Mehmood, & Haneef, 2015). The inflammatory response is associated with the expression of several proteins such as NF-κB, COX-2, COX-1, and TNF-α (Sundaram, Vasanthi, & R, 2019). Therefore, the influence of piperine on NF-κB, COX-2, COX-1, and TNF-α was evaluated using in-silico analysis. Several scientific studies have shown that COX has two isoforms responsible for the synthesis of prostaglandins among which COX-1 derived prostaglandins mainly mediates cytoprotective actions on the gastrointestinal mucosa whereas Cyclooxygenase-2 (COX 2) is responsible for the production of inflammatory mediators such as prostaglandins and its metabolites which are involved in vasodilation, thermoregulation, etc. (Zarghi & Arfaei, 2011); whereas NF-κB can induce the expression of pro-inflammatory cytokines, chemokines, and genes involved in inflammasome complex formation (Dorrington & Fraser, 2019). On the other hand, TNF-α is a multifunctional cytokine involved in systematic inflammation and is one of the cytokines that make up the acute phase reaction and their respective metabolites (Parameswaran & Patial, 2010) .
Table 4
Receptor | Ligand (piperine) | Interacting amino acid |
NF-κB | -2.68 | (H-Glutamine 132) |
COX-2 | -5.33 | (Hydrophobic interaction) |
COX-1 | -3.0 | (H-Threonine 94) |
TNF-α | -4.07 | (H-Serine 95) |
Docking studies are important for predicting the functionality of the bioactive compounds before conducting wet-lab experiments and most frequently it is an important aspect in structure-based drug designing (Ferreira, Santos, Oliva, & Andricopulo, 2015). Docking simulation of the alkaloid piperine with the ligands (inflammatory markers) viz. NF-κB, COX-2, COX-1, and TNF-α were performed and the binding conformations and binding affinities were obtained by using the docking protocol, Autodock (autodoc tools- 1.5.4 version). All the possible rotatable bonds of the compound piperine were taken into account in the process of docking to identify the best binding confirmation with the above-mentioned inflammatory proteins. NF-κB showed the binding score of -2.68 in the conformation where piperine bonded with Glutamine-132 via hydrogen bond in the binding pocket of the former protein whereas COX-2 showed the best binding score of -5.33 in the conformation mentioned in the figure where piperine hydrophobically interacted with the protein molecule. On the other hand, COX-1 showed the best binding score of -3.0 in the conformation where piperine bonded with Threonine 94 via hydrogen bond whereas TNF- α showed the best binding score of -4.07 in the conformation where piperine bonded with Serine 95 via hydrogen bond (Table 4 ) and (Figure 4 ). This proves that piperine has a significant binding interaction with these inflammatory proteins thereby modulating their expression. In the present study, molecular docking studies of inflammatory marker proteins with piperine provides valuable insights about binding stability and scoring functions with the targeted receptors. Our results are in congruence with the earlier studies on the anti-inflammatory potential of piperine on fibroblast cell lines (Bang, Oh, & Choi, 2009), and in experimental animal models of inflammation wherein, piperine significantly reduced the proinflammatory mediator TNF α in experimental rats (Umar, Sarwar, Umar, & K, 2013). However, further studies are warranted to validate the traditional claims and identify the detailed mode of action of piperine before translating it to clinical use.