Inhibitory Reactivity of Capsaicin with α-Amylase and α-Glucosidase Related to Antidiabetes using Molecular Docking and Quantum Calculation Methods

This work aims to investigate the inhibitory activity of capsaicin, which is one of capsaicinoid compounds, on these enzymes using a molecular docking and quantum calculation. Acarbose, a commercial diabetes drug, was also investigated for comparison. The docking results revealed that acarbose yields better inhibition efficiency with binding free energy (DGbinding) of about -8.2 to -11.9 kcal/mol, and inhibition constant (Ki) of about 0.0002 to 0.4 μM, where as capsaicin provided the DGbinding of -5.8 to -6.1 kcal/mol and Ki of 23.7 to 45.9 μM. The total binding energy (DEbinding) between each inhibitor and amino acids in active site of enzyme obtained from quantum calculation with MP2/6-31G(d,p) level is in agreement with the DGbinding, i.e. the DEbinding of acarbose was larger negative than that of capsaicin. The amino acids interacting with inhibitor as hydrogen bond mainly contribute to the total binding energy. Nevertheless, it could be concluded that capsaicinoids have high potential to be developed as an alternative drug for diabetes disease.


INTRODUCTION
Nowadays, diabetes is known as the most dangerous disease killing more people than cancer and AIDS combined.Actually, it is said to be a modern disease, and has traditionally been viewed as a disease of rich countries 1 .Around 3.2 million individuals in Thailand have diagnosed currently being diabetes.Thai people are expected to fall in this disease over than 4.3 million in 2035 2 .Diabetes is caused by a metabolic disorder, in which the body cannot properly store and expense the energy found in food 3,4 .Two main types of diabetes are type I and type II diabetes 5,6 .Type I diabetes is dominant when the body's immune system attacks the cells in the pancreas that generate insulin.Type II diabetes is the other when the body does not produce enough insulin, or the cells ignore the insulin 7,8 .One of therapeutic approach for diabetes is to control the blood sugar level by the process of an inhibition of carbohydrate hydrolyzing enzymes such as α-amylase and α-glucosidase [9][10][11][12][13] .
α-Amylase and α-glucosidase inhibitors can help individual case with type II diabetes by preventing the digestion of carbohydrates such as starch 14 .Carbohydrates are normally converted into simple sugars or monosaccharides, which can be absorbed through the intestine 15 .The known natural inhibitors of digestive enzymes include phenolic compounds, anthocyanins, terpenoids and others 16,17 .They operate through various mechanisms.A acarbose, miglitol and voglibose are the commercially available drugs that inhibit pancreatic α-amylase and α-glucosidase in the small intestine (Fig. 1a).The side effects of both enzymatic inhibitors are flatulence, low blood sugar and hepatitis.In recent years, a variety of research has been done on plants.In addition, several herbs and fruits have the property of inhibitors [18][19][20][21] .
Chili peppers are very well-known economic fruit routinely consumed in many tropical countries.One of the plant's products is chili peppers of genus Capsicum and belongs to the Solanaceae family.The five domesticated species of chili peppers include Capsicum annuum, Capsicum frutescens, Capsicum chinense, Capsicum pubescens, and Capsicum baccatum 22 .Capsaicinoids are the chemicals, which give the rise in heat level of the chilli.They present in almost chilli fruits.The two majorities of capsaicinoids presenting in most varieties of the spice from hot chilli peppers are capsaicin (8-Methyl-N-vanillyl-6-nonenamide) and dihydro-capsaicin (8-Methyl-N-vanillylnonanamide) as shown in Fig. 1b & Fig. 1c, in which capsaicin consists of about 69% and dihydrocapsaicin consists of about 22% of the fruits 23 .The distinctive properties of capsaicinoids in food are such pungent additives.In pharmaceutical industries, they are used for the treatment of pain and inflammation in different kinds of diseases.The Protein Data Bank (PDB) (RCSB PDB) is the depository for crystal structures of biological macromolecules from the experiment.Currently, the PDB provides X-ray crystallography of more than 112,968 crystal structures from 1971 to 2015.In addition, the crystal structures of α-amylase and α-glucosidase are available in the PDB of more than 334 and 183 crystal structures, respectively.The PDB codes of the crystalline structure of the wild-type of α-amylase complex and α-glucosidase complex are 3OLD and 3TOP, respectively, which are introduced for this study.
AutoDock program has been widely used and is the most cited docking software.It is very fast, provides high quality predictions of ligand conformations, and good correlations between predicted inhibition constants and experimental ones.In addition, AutoDockTools is a software used to prepare the input files for AutoDock 24 .In order to understand deeper in the molecular interactions between the inhibitor and the residues in any binding pocket which is not observed by an experiment, theoretical investigation has thus been an alternative method to investigate the enzyme and inhibitor interaction in details for large molecular system.Recently, quantum chemical calculations were used to study the effect of mutation related to inhibitor resistant.It was found that the calculations are significant for guiding the cause of resistant 25 .
The present study aims to evaluate capsaicin whether it can inhibit both α-amylase and α-glucosidase in human compared with acarbose.
Then, the particular interactions between these inhibitors with individual residues in binding site of both enzymes have been investigated.The scope and limitation of this study composes of (1) investigating the binding free energy and inhibition constant between inhibitors (capsaicin and acarbose) and enzymes (α-amylase and α-glucosidase) by using a molecular docking method, and (2) calculating the interaction between the inhibitors and amino acids in active site of the enzymes by using quantum mechanics calculation.
The objectives of the study are (1) to evaluate the inhibitory activity of capsaicin and acarbose (for comparison) on α-amylase and α-glucosidase; (2) to investigate the optimum conditions for docking of α-amylase and α-glucosidase; (3) to calculate the interaction between the inhibitors (acarbose and capsaicin) and amino acid in the binding pocket of the enzymes (α-amylase and α-glucosidase) based on the quantum chemical calculations at MP2/6-31G(d,p) level; and (4) to investigate the main interactions between the amino acid residues in the binding site of both enzymes and the inhibitors contributing to the inhibitory efficiency.
α-Glucosidase and α-amylase inhibitors are found in microorganisms, plants and animals 26,27 .Capsaicin, an alkaloid or capsaicinoids, is the principal pungent and irritating constituent of hot chilli peppers that are widely used as food additives and possess antimicrobial properties 28 .Several instrumental techniques have been used for the structural insights of capsaiciniods including molecular docking 29,30 , molecular dynamics simulation 31,32 , homology modeling 33,34 , virtual screening 35,36 , and bioassay validations 37 .
Recently, molecular docking has mostly been used for capsaicinoid analysis because of its rapidness and reliability.Molecular docking is one of the modeling techniques used to determine the position docking of capsaicinoids in the binding pocket of enzymes 38 .Amylase, maltase and sucrase are the enzymes that digest carbohydrate found in the alimentary canal.Amylases are one of the main enzymes used in industry.Amylase has been reported to occur in microorganisms, although they are also found in plants and animals.Two major classes of amylases have been identified in microorganisms, namely α-amylase and glucoamylase.α-Amylases (endo-1,4-α-D-glucan glucohydrolase, E.C. 3.2.1.1)are extracellular enzymes that randomly cleave the 1,4-α-D-glucosidic linkages between adjacent glucose units in the linear amylose chain 39 .α-Amylase is the enzyme responding for the hydrolysis of starch and carbohydrates to monosaccharides 40 .Amylase is found in saliva and small intestine by synthesizing from pancreas.The functional activity of amylase is to break down α-glycosidic bond in starch to be maltose, glucose and dextrin.The inhibition function of α-amylase is beneficial in reducing sugar that your body will absorb it.Currently, medical trend is interested in using amylase inhibitor from natural plant used in treatment of diabetes.
Maltase or also known as α-glucosidase and sucrase or known as invertase work in the area of microvilli in small intestine; glucosidase reaction is well known.Glucosidase inhibitor such as acarbose inhibits both of maltase and sucrase.Thus, in carbohydrate catabolism, disaccharide sugars cannot digest to monosaccharide ones.As a result, it reduces the sugar absorption to the blood system.Glucosidases are glycosidic hydrolase enzymes categorized under the EC number 3.2.1.20which have also been found in microorganisms.
α-Amylases are universally distributed throughout the animal, plant and microbial kingdoms.Two kinds of α-amylase are produced by many mammals, saliva α-amylase from the parotid gland and the pancreatic α-amylase from the pancreas 41 .α-Amylase (AIs) and α-glucosidase (AGIs) inhibitors are widely used in the treatment of patients with type II diabetes.Diet and exercise are the first step in its treatment.But if these measures alone fail to sufficiently control blood glucose levels, starting oral drug therapy is recommended 42 .α-Glucosidase and α-amylase inhibitors are classified into the two major groups as proteinaceous and non-proteinaceous inhibitors.
One benefit of hot chili peppers is therapeutic approach for treating diabetes which decreases the postprandial hyperglycemia.This is done by retarding the absorption of glucose through the inhibition of carbohydrate hydrolyzing enzymes, α-glucosidase and α-amylase, in the digestive tract 43 .The control of postprandial hyperglycemia is an important strategy in the management of diabetes mellitus, especially type II diabetes mellitus (2DM), and reducing chronic complications associated with the disease [44][45][46] .Hence, the inhibition of enzymes (α-amylase and α-glucosidase) involved in the digestion of carbohydrates can significantly decrease the postprandial increase of the blood glucose after a mixed carbohydrate diet, by delaying the process of carbohydrate hydrolysis, absorption, phenolic phytochemicals from spices and have shown promising potentials 47,48 .Presently, the computational studies have not been reported yet about the inhibitory effect of the enzymes by using capsaicinoids as inhibitors.
The non-proteinaceous α-amylase and α-glucosidase inhibitors contains diverse types of organic compounds such as acarbose, acarbose analogues, hibiscus acid, tannins, flavonoids and glucopyranosylidene-spiro-thiohydantoin.α-Glucosidase inhibitors (AGIs) reversibly inhibits a number of α-glucosidase enzymes (e.g.maltase), consequently delaying the absorption of sugars from the gut 49 .In a recent study among healthy subjects it was suggested that the therapeutic effects of AGIs are not only based on a delayed digestion of complex carbohydrates, but also on metabolic effects of colonic starch fermentation50.Acarbose (Glucobay®) is the most widely prescribed AGI.The other AGIs are miglitol (Glyset®) and voglibose (Volix®, Basen®).AGIs might be a reasonable option as first-line drug in the treatment of patients with DM2 as it specifically targets postprandial hyperglycemia, a possible independent risk factor for cardiovascular complications 51 .Although rare cases of hepatic injury are described, AGIs are expected to cause no hypoglycemic events or other life-threatening events, even at overdoses, and cause no weight gain 52 .The α-amylase inhibitor contains diverse types of organic compounds as mentioned above.The inhibitory activity of these compounds against α-amylase is due in part to their cyclic structures, which resemble substrates at catalytic sites of α-amylase 53 .
A number of docking programs have been developed during the last two decades and made available to academic institutions at little or no charge.Basic characteristics such as supported platforms, license terms, as well as applied docking algorithms and scoring functions are presented.DOCK 1.0 54 was the first automated receptor-ligand docking program.It was designed in 1982 by Irwin Kuntz in Department of Pharmacology at the University of California at San Francisco.At present, there are at least a dozen docking tools in the market, the most commonly used being: AutoDock 55 , DOCK 56 , FlexX 56 , GOLD 57 , LigandFit 58 and the relatively new tools: Glide 59 , FRED 60 and the youngest Surflex 54 .

Molecular docking
Molecular docking is a computational method used to study the formation of intermolecular complexes of one smaller molecule (called ligand) with a larger molecule (called receptor), which usually is a protein of known three-dimensional structure.Different types of interactions between the molecules can be distinguished: protein-protein, protein-DNA, DNA-ligand and protein-ligand 61 .For an enzyme and inhibitor, docking aims at correct prediction of the structure of the complex under equilibrium conditions, expressed as equation ( 1): (1) The particular interest is the free energy of binding (DG binding ) expressed as Where, G cpx , G ligand and G receptor are free energy of complex, ligand and receptor, respectively.
The free energy of binding is related to binding affinity as shown in equations ( 3) and (4); Where, R is 1.987 cal•K -1 •mol -1 , T is absolute temperature (K), and Ki is inhibitor constant.
Molecular docking uses the Monte Carlo simulation (MC) to determine the receptor molecule in place and add ligand into it then calculate energy between ligand and receptor molecule (E 1 ).Then, change in ligand's position randomly and calculate for new energy (E 2 ), if E 2 is less than E 1 , a ligand will move to this position after that the ligand will change its position continuously until it reaches specified number of times (for this study, it is 2,500,000 times) or get the lowest energy.The last position of ligand with the lowest energy is called 1 run of docking.In general, docking is performed for 50 runs.
The energy calculation between ligand and receptor is as followed.Firstly, the program makes up grid map for probe atom of ligand.The grid box must cover an active site of receptor molecule as shown in Fig. 2. Each probe atom of ligand is placed in the upper left corner of the box.Calculation of energy between probe atom and receptor is done, and then this data is stored.The program will calculate continuously through every grid position.Therefore, one gets grid map for each probe atom.
The calculation of energy is done while the stimulation is performed by considering of both probe atom and grid point, and reading the energy value from grid map that relates with the kind and the position of atom in each box, and then combining them together for total energy of the system that consists of ligand and receptor molecules.
due to polar atoms.To quantify the contribution of each, these interactions can be modeled by a potential function that gives the energy of the interaction as a function of distance, angle, or charge62,63.The total steric energy of a molecule can be written as a sum of the energies of the interactions in equation ( 5): Where, V is potential energy, V bond is bond stretching energy, V angle is angle bending energy, V torsion is torsional angle energy, V vdW is van der Waals force, and V Coulomb is Coulomb force.

Quantum mechanics
Quantum mechanics or quantum theory is the theoretical basis of modern physics that explains the nature and behavior of matter and energy in the atomic and subatomic level.The nature and behavior of matter and energy at that level is sometimes referred to as quantum physics and quantum mechanics.The development of quantum mechanical techniques is more generally applicable and that can be implemented on a computer.When the properties of such particle (e.g.mass, charge, etc.) are expressed in macroscopic units then the values must usually be multiplied, which are too complication.Several mathematical techniques are applied to these problems for simplification 64 .The starting point for any discussion of quantum mechanics is the Schrödinger equation.The equation is the basis for most of the computational scientist.It can be solved exactly for only a few problems such as the particle in the box, the harmonic oscillator, the particle on a sphere and the hydrogen atom.One of them has been widely used in the quantum chemical calculation is the Hartree-Fock (HF) model.

Ab initio
The most common type of ab initio calculation is called Hartree-Fock (HF) calculation, in which the primary approximation is called the mean field approximation.This means that the Coulombic electron-electron repulsion is not explicitly taken into account, however, its average effect is included in the calculation.This is a variation calculation, which implies that the approximate energies calculated are all equal to or greater than the exact energy.The accuracy of the calculation depends on the

Molecular mechanics
Molecular mechanics assumes the steric energy of a molecule to arise from a few specific interactions within a molecule.These interactions include the stretching or compressing of bonds beyond their equilibrium lengths, bending of bond angles, torsional effects of twisting about single bonds, the van der Waals attractions or repulsions of atoms that come close together, and the Coulomb interactions between partial charges in a molecule size of the basis set used, however because of the mean field approximation, the energies from HF calculations are always greater than the exact energy and tend, with increasing basis size, to a limiting value called the Hartree-Fock limit.
A n additional issue that affects the accuracy of the computed results is the form chosen for the basis functions.The actual form of the single electronic molecular wave function (molecular orbital) is of course not known.The forms, used for the basis functions, can provide a better or worse approximation to the exact numerical single electron solution of the HF equation.The basis functions used most often are combinations of either Slater type orbitals (exp(-ax)) namely as STO or Gaussian type orbitals (exp (-ax 2 )), GTO.Molecular orbital is formed from linear combinations of atomic orbitals, which are nothing more than linear combinations of the basis functions with coefficients found from the appropriate atomic HF calculations.Because of this approximation, most HF calculations give a computed energy greater than the Hartree-Fock limit.The exact set of the basis functions used is often specified by an abbreviation, such as STO-3G or 6-311++g**.From this aspect, it can be read about structures and features of some popular basis sets.
An alternative ab initio method is the second Møller-Plesset perturbation theory (MP2).MP2 is perhaps the simplest model to take reasonable account of electron correlation, and generally provides accurate descriptions of equilibrium structure, conformation and energetics of a variety of chemical reactions, including reactions where chemical bonds are broken.MP methods are suppor ted for the same basis sets and pseudopotentials available for Hartree-Fock and density functional models.

Semi-empirical
Semi-empirical calculations are set up with the same general structure as a HF calculation.Within this framework, certain pieces of information, such as two electron integrals, are approximated or completely omitted.In order to correct the errors introduced by omitting these parts of the calculation, the method is parameterized, by curve fitting in a few parameters or numbers, in order to give the best possible agreement with experimental data.The good side of semi-empirical calculations is that they are much faster than the ab initio calculations.The bad side of semi-empirical calculations is that the results can be erratic.If the molecule under study is similar to molecules in the data base used to parameterize the method, the results may then be very good.If this molecule is significantly different from anything in the parameterization set, the answers may be poor.Semi-empirical calculations have been very successful in computational organic chemistry, where there are only a few elements used extensively and the molecules are of moderate size.However, semi-empirical methods have been devised specifically for the description of inorganic chemistry as well.

Program tools and instruments
There are five program tools including AutoDockTools-1.5.6 (ADT), AutoDock 4.2, HyperChem 8.0, Weblab Viewer Pro 4.0 and Gaussian 09.In this study, all calculations were performed using computer cluster located at Department of Chemistry, Faculty of Science, Khon Kaen University, and National Electronics and Computer Technology Center (NECTEC), Bangkok (supported by National e-Science Infrastructure Consortium).

Molecular docking
Molecular docking is a computer-based method which predicts the preferred orientation of one molecule to a second when bound to each other to form a stable complex.The aim of molecular docking is to predict the interaction between two molecules.Its ultimate goal is to find the most stable conformation between ligand and protein.The procedure of docking consists of three important steps: binding site identification and characterization, orientation of the ligand within the binding site, and evaluation of the orientation for appropriateness of fit.A stepwise process for molecular docking is shown in Figure 3.

Ligand set up
For capsaicin, it was generated by using Hyper Chem-8.0 program (Hypercube, 2007), and then optimized with PM3 semi-empirical level, and saved as hin file format.For acarbose, the X-ray structure was downloaded from the Protein Data Bank (PDB code: 3TOP) complexed with α-glucosidase.Acarbose molecule was separated Fig. 3.A stepwise process for molecular docking and then missing hydrogen atoms were added and their geometries were optimized with MM+ forcefield using HyperChem 8.0.The final structure was also saved as hin file.The hin file was then converted to pdb file format through Babel in Linux.After loading the pdb file into ADT, the atomic charges were assigned with Gasteiger charge and then the rotatable bonds were defined and saved as pdbqt file.
For acarbose, all hydroxyl bonds were assigned as rotatable bond, totally 13 bonds, whereas for capsaicin, the rotatable bonds were assigned for 4 sets consisting of one, two, three and four rotatable bonds.These bonds were listed in Table 1 and shown with arrow.

Enzyme set up
The X-ray structures of enzymes were retrieved from the Protein Data Bank (www.rcsb.org).The PDB code of human pancreatic α-amylase is 3OLD, and the PDB code of human α-glucosidase is 3TOP.As the human α-glucosidase is a symmetric dimmer, the chain A of the enzyme was deleted and only chain B was used in the docking.The water molecules in the file were then deleted using Weblab ViewerPro 4.0 program65 and the protein file was saved as protein.pdb.Each pdb file was loaded into ADT program, then all missing hydrogen atoms were added to the protein structure, all non-polar hydrogens were merged with carbon atoms, and Gasteiger charge were assigned.When design on molecular calculation, only selected elements in the active site vicinity of the enzyme are fixed to make any interactions within the configuration of its counter ions of the inhibitor used, thus each monomer (A or B) of alpha-glucosidase dimer can be represented for their substrate-enzyme interactions basis.It is done by molecular docking and quantum calculation with each monomer, but does not work for its dimer.Certainly, it would be affect to the calculated results as a whole, if concerning on that particular aspect.However, it is aimed only to point out what really specific interactions can occur between the active of the enzyme and the molecular configuration of the capsaicin.In order to performed flexible protein docking, the rotatable bonds of protein were defined as listed in Table 2.The files were saved as flexible.pdbqt and rigid.pdbqt.

Grid parameter files
The grid parameter files were prepared for molecular docking by ADT.For re-docking procedure, the grid box was studied at 60×60×60, 66×66×66 and 70×70×70 Å for both 3OLD and 3TOP.The center of grid box was defined as the center of ligand complexed with its enzyme (x = 6.0, y = -16.9,z = -24.4for 3OLD and x = -51.5,y = 9.0, z = -64.8for 3TOP), the spacing of grid box was set at 0.375 Å and saved as protein.gpf.Then, these grid parameter files were used to generate the grip map files (with command autogrid4.2-p protein.gpf-l protein.glg&).The grid map files contain information about the energy between probe atoms and protein at each grid point and were used for docking calculation.

Docking parameter files
The docking parameter files were prepared for molecular docking by ADT.The main parameters were energy evalution of 2,500,000 and 50 runs with Lamarckian genetic algorithm search.The default values were used for other parameters.Docking calculations were performed by using AutoDock4.2program (with command autodock4.2-p ligand.dpf-l ligand.dlg&).

Quantum chemical calculation
Quantum chemical calculations were used to calculate particular interactions of all complex structures.These calculated results can be related to the inhibitory efficiency obtained from molecular docking.Therefore, the quantum chemical calculations are important for investigating molecular interactions in various types of enzymes (α-amylase and α-glucosidase) to inhibitors (capsaicin and acarbose).The complex structures obtained from docking result were employed by selecting amino acids surrounding inhibitors within the inter-atomic distance of 4 Å.The particular interaction energy was focused on the couple of inhibitor and each amino acid.Then, hydrogen atoms were added into the cutting regions using HyperChem 8.0 program and was investigated using Gaussian 09 program (Gaussian, 2013) with MP2/6-31G(d,p) level of calculations.

Conformational and energetic properties of the enzyme and the inhibitor
The optimal conditions for molecular docking of α-amylase and α-glucosidase enzymes were determined with the re-docking step based on ligand with its enzyme.The docked structures of enzymes and inhibitors were generated using AutoDock4.2program.Also the binding free energy and inhibition constant were evaluated.The energies such as van der Waals, hydrogen bonding, electrostatic and hydrophobic between enzymes and inhibitors were investigated.Quantum calculations with MP2 method were performed in order to determine the main interactions between the inhibitors and selected amino acids surrounding the inhibitors within the inner-atomic distance of 4 Å.

RESULTS AND DISCUSSION
This research aims to evaluate the inhibitory activity of capsaicin on α-amylase and α-glucosidase enzymes in order to search for an alternative diabetes drug.Acarbose, a commercial drug, was also evaluated for comparison.The procedure of the study consists of three main steps: (1) re-docking, (2) docking calculation and (3) quantum calculation.The results and discussion of these steps are as followed.

Re-docking
The re-docking process aims to investigate the optimal conditions used for docking calculations of α-amylase and α-glucosidase.The main conditions or parameters for re-docking are grid box size and number of rotatable bonds of enzymes.As mentioned previously in CHAPTER III, for both enzymes, the grid box size was defined as 60x60x60, 66x66x66 and 70x70x70 Å 3 .For α-amylase, the grid box center was x = 6.0, y = -16.9,z = -24.4which was a center of ligand.A number of rotatable bonds were assigned for three sets (0, 3 and 7 bonds).Similarly, for α-glucosidase, the grid center was x = -51.5,y = 9.0, z = -64.8and rotatable bonds were assigned for two sets as 0 and 4 bonds.The main criteria used to decide the re-docking result is root mean square deviation (RMSD) which is a comparison of ligand position of docked structure and X-ray structure.The acceptable RMSD value must be less than 2.0 Å.The results of re-docking are shown in Table 3 and Table 4.As shown in the tables, the docking results (both DG binding and RMSD values) depend on the number of rotatable bonds and grid box size.
For α-amylase, the rotatable bonds of 3 bonds yielded the RMSD values less than 2.0 Å which is an acceptable value.Considering the grid box size, one found that the box size of 60x60x60 and 66x66x66 Å 3 provided better docking result than the box size of 70x70x70 Å 3 , as seen from the number in cluster value (Table 3).For α-glucosidase (Table 4), both number of rotatable bonds (0 and 4 bonds) resulted the RMSD values less than 2.0 Å, and the number in cluster values were not different.However, the grid box size of 60x60x60 Å3 and 66x66x66 Å 3 yielded the lowest free energy of binding (about -11.9 kcal/mol) compared with grid box size of 70x70x70 Å 3 .Therefore, one can conclude that 3 flexible residues (for α-amylase) and 4 flexible residues (for α-glucosidase) with grid box size of 60x60x60 and 66x66x66 Å 3 are the suitable conditions for docking calculation.Such condition is shown in Table 5.The re-docked positions of acarbose for both enzymes are shown in Fig. 4-7.One can see that the re-docked positions obtained from grid box size of 60x60x60 Å 3 and that obtained from grid box size of 66x66x66 Å 3 are not significant different and are in agreement with the X-ray structure.

Docking of capsaicin compared to acarbose
Using the suitable parameters as listed in Table 5, the docking calculations of capsaicin with both α-amylase and α-glucosidase enzymes were obtained and shown in Tables 6 and 7, respectively.For α-amylase (Table 6), the docking results depend on the number of rotatable bonds of acarbose.The binding free energy varies from about -5.2 to -6.3 kcal/mol and the number in cluster varies from 21 to 50 runs.Clearly, the best docking result was obtained from 2 rotatable bonds of acarbose yielding the number is cluster of 50 runs, although it provided the DG binding of about 0.5 kcal/mol higher than that of 3 rotatable bonds for both grid box sizes.
Similarly, the docking result of α-glucosidase with capsaicin revealed that 2 rotatable bonds of capsaicin yielded the best docking calculation with the number in cluster of 47 and 34 runs for grid box size 60x60x60 and 66x66x66 Å 3 , respectively, and the DG binding of about -6.1 kcal/mol (Table 7).The docking results of capsaicin compared to acarbose were summarized in Table 8.One found that the binding free energies of acarbose were about 1.4 times (for α-amylase) and about 2 times (for α-glucosidase) larger negative than capsaicin.These results yielded the inhibition constants of less than 1 µM for acarbose and of about 20-50 µM for capsaicin.
In conclusion, the inhibitory activity of acarbose is higher than that of capsaicin.In addition, acarbose and capsaicin can inhibit α-glucosidase better than α-amylase.Considering the docking of inhibitors of 60x60x60 Å 3 compared to grid box size of 66x66x66 Å 3 , one can see that the results (DG binding and Ki) are not significant different.This is not surprise, because the positions of inhibitors obtained from both grid box sizes were not different.Therefore, the docked structures obtained from the grid box size of 60x60x60 Å 3 were employed for further analysis.Figures 9-12 show the orientation of inhibitors surrounded by amino acids in active site of enzymes.The interactions between the inhibitors and amino acids were classified as (1) hydrogen bond displayed as dash line, (2) π-π interaction and (3) electrostatic interaction (ELEC).Note that, hydrogen bond is the strongest interaction compared to others.These interactions were listed in Table 9 higher inhibition efficiency as mentioned before.This is not surprised, because acarbose composes of many -OH groups that can form hydrogen bond.
Acarbose can form more hydrogen bonds with amino acids of two enzymes than capsaicin leading to

Quantum chemical calculations
In order to investigate the main inhibitoramino acid interactions contributing to the inhibition efficiency of inhibitor on enzyme, quantum chemical calculations were employed.All amino acids in active side of enzyme surrounding the inhibitor within a distance of 4 Å were taken into account.The MP2/6-31G(d,p) level was used for all calculations.The interactions were investigated in term of binding energy (DE binding ) expressed as DE binding = E inhibitor-amino acid -E inhibitor -E amino acid (6)   Where E inhibitor-amino acid , E inhibitor and E amino acid are the energy of inhibitor-amino acid complex, inhibitor and amino acid, respectively.The inhibitors with molecular surface and amino acid taken into account for the calculations were displayed in Fig. 13-16 for α-amylase and Fig. 17  For α-amylase (Table 11), both inhibitors interact with various amino acids.In the case of acarbose, it can interact with amino acids as both attraction and repulsion.The repulsive interactions are caused by ILE235, HIS299, HIS305 and ALA 307 and the interactions vary from 0.2 to 20.3 kcal/ mol in which ALA307 provides the largest repulsion.Contrastly, the attractions vary from about -0.8 to -50.2 kcal/mol.The main attractions are caused by the amino acids forming hydrogen bonds with acarbose; ARG195, ASP197, LYS200, GLU233 and ASP300.Although GLY308 forms hydrogen bond with acarbose, the interaction is quite small (-2.3 kcal/mol).Other amino acids yielding attraction interact with acarbose as electrostatic interaction.
In the case of capsaicin, it can interact with α-amylasewith less number of amino acids.The attractions vary from about -0.2 to -28.1 kcal/ mol.The largest attraction is obtained from LYS200 caused by hydrogen bond.It can be seen that LYS200 and GLU233 provide very high attraction with both acarbose and capsaicin, where as ILE235 yields repulsion with both inhibitors.
For α-glucosidase (Table 12), acarbose yields strong attraction with amino acids, which are LYS1164, ASP1279, ARG1510, ASP1526 and HIS1584 with the binding energy of -11 to -38 kcal/ mol.Some electrostatic interactions were also large, e.g.LYS 1460.These interactions mainly contribute to the total energy of -165 kcal/mol.In the case of capsaicin, there is only one hydrogen bond but strong (about -29 kcal/mol for ASP1526).Long hydrocarbon chain of capsaicin can interact with amino acids as electrostatic interaction.The important interactions were provided by ASP1157, ASP1420, LYS1460 and ASP1555.In addition, another interaction of capsaicin and α-glucosidase is π-π interaction which occurs with PHE1559 with the energy of -6.2 kcal/ mol.Note that, PRO1159 results the repulsion for both inhibitors with the energy of about 6 and 4 kcal/ mol for acarbose and capsaicin, respectively.
The binding free energy and total binding energy obtained from molecular docking and quantum chemical calculation, respectively, of each system were summarized in Table 13.Clearly, the energy obtained from quantum calculation is consistent with that obtained from molecular docking, in which the total binding energy between acarbose and enzymes is larger than capsaicin.In fact, water salvation on protein interactions is very crucial role concerning their enzyme activity, in particular for in vivo test.However, in this case, insights into the molecular docking and/or quantum calculation aspects, water molecule in the active site vicinity is deleted due to aiming of specific atomic coupling interactions between the enzyme substrate and its counter ions of the molecular configuration of the inhibitor studied.H-bond interactions would be actually pronounced, if we concern on their molecular water solvation among the active sites of the enzyme and the inhibitor used.Therefore, it would be affected to the results obtained.However, we have tested this molecular docking and quantum calculation with Acarbose as drug reference.Thus, we can suggest that these data would be lined with physiolocal functions of both enzymes related to antidiabetes activity as well.In addition, since capsaicin contains in chili paper which is a fruit routinely consumed in many countries, one may conclude that it has a potential to be developed as an alternative drug for diabetes.

CONCLUSION
This work aims to investigate the inhibitory activity of capsaicin on α-amylase and α-glucosidase enzymes which are carbohydrate hydrolyzing enzymes involving diabetes disease.This is to search for a new diabetes drug.The study employed theoretical approach including molecular docking and quantum chemical calculation.Acarbose, a commercial drug, was also studied for comparison.The first step of this work was re-docking procedure aiming to determine the optimal condition for docking calculation.The re-docking result revealed that the grid box size of 60x60x60 Å 3 was suitable for both enzymes.In addition, 3 rotatable bonds of flexible protein (ASP197, GLU233 and ASP300) and 4 rotatable bonds (ASP1157, ASP1279, ASP1526 and HIS1584) were suitable for α-amylase and α-glucosidase, respectively.These conditions were then used for molecular docking.The docking resulted that, for α-amylase, the binding free energy (DG binding ) and inhibition constant (Ki) were about -5.8 kcal/mol and 45.9 µM, respectively.This evaluation was worse compared to acarbose (DG binding = -8.2kcal/mol and Ki = 0.4 µM).For α-glucosidase, similarly, acarbose yielded better inhibition efficiency with DG binding of -11.9 kcal/mol and Ki of 0.2 nM, whereas capsaicin provided the DG binding of -6.1 kcal/mol and Ki of 23.7 µM.Note that, the inhibitory activity of both inhibitors on α-glucosidase is better than that on α-amylase.For in vitro test, anti-α-amylase activities of Thai varieties of hot chilli extracts were ranged of 43.0-147.0mg ACE / g DW or 21-51% inhibition 66 , while anti-α-glucosidase activities of some chilli pepper extracts were also found to be 23-66% 67 .
The docked structures were later used to investigate the inhibitor-amino acid interactions in which only amino acids surrounding the inhibitor within 4 Å were taken into account.This investigation employed the quantum chemical calculation with MP2/6-31G(d,p) level.The total interaction energy is in agreement with DG binding obtained from molecular docking, and the results revealed that amino acids interacting with inhibitor as hydrogen bond mainly contribute to the total energy.Because acarbose consists of several -OH groups which can form hydrogen bond with amino acids in active site of enzyme, this is a reason that the inhibitory activity of acarbose is higher than capsaicin.Nevertheless, one may say that capsaicin has a potential to be developed as an alternative drug for diabetes disease.

Fig. 2 .
Fig. 2. Making of a grid map for probe atom

Fig. 5 .
Fig. 5. Orientation of acarbose in active site of α-amylase obtained from X-ray structure (blue) and re-docking (purple) with grid box size of 66x66x66 Å 3

Table 3 : The results of re-docking for
α-amylase

Table 4 : The results of re-docking for
α-glucosidase