Analysis of Compounds Isolated from Gnetum gnemon L. Seeds as Potential ACE Inhibitors through Molecular Docking and Molecular Dynamics Simulations

Triputra and Yanuar: Analysis of Compounds Isolated from Gnetum gnemon L. Seeds as Potential ACE Inhibitors through Molecular Docking and Molecular Dynamics Simulations

Authors

INTRODUCTION

Indonesia has the opportunity to develop medicines or herbs that are essential for health. One of the plants having medicinal properties is Gnetum gnemon L. (melinjo). Melinjo is widely available and consumed in large quantities in Indonesia.1 Hydrolyzed-proteins of melinjo seed had been proven to lower blood pressure of hypertensive rats.2

Hypertension is a condition in which arterial blood pressure is continuously elevated. To be diagnosed with hypertension when one’s systolic pressure is higher than 140 mm Hg, and the diastolic pressure reaches greater than or equal to 90 mm Hg.3 Hypertension affects about 40% of people aged 25 years or older, both in developing countries, and developed countries.4 Hypertension is a leading problem in Indonesia, because of the high prevalence and the long term effect.5 In 2013, Indonesia has 25.8% prevalence of adult over 18 years with hypertension based on blood pressure measurements using the sphygmomanometer.6 The majority of hypertensive patients rely on ACE (angiotensin-converting enzyme) inhibitors in managing their conditions.7 Angiotensin-converting enzyme (ACE) plays an inportant role in the conversion of angiotensin I to angiotensin II, a known potent vasoconstrictor.8 Furthermore, angiotensin II increases blood pressure which is the significant risk factor for hypertension.3

Gnetum gnemon L. (melinjo) is a plant that spreads in Southeast Asia, Melanesia, and the Pacific Islands. This tree has a height of 10-15 meters, with a slender-shaped trunk with a diameter of approximately 40 cm. This plant has one oval-shaped seed on each fruit.9 Melinjo belongs to Gnetum (Gnetaceae) genus plant, which is known containing many stilbene derivatives.10-12

Gnetum gnemon L. seed extracts have been known to have some biological activities.13 The previous research showed that the ethyl acetate extract of melinjo seed had the highest inhibitory activity towards ACE compared to other extraction solvents in vivo.14 Furthermore, there are compounds found in the ethyl acetate extract of melinjo seed such as resveratrol, gnetol, isorhapontigenin, trans-piceid, gnetin C, trans-ε-viniferin, gnemonol M, gnemonol K, gnemonoside A, gnemonoside B, gnemonoside C, gnemonoside D, and latifolol. A recent research demonstrated that several compounds of melinjo seed extracts have a better affinity than lisinopril to ACE based on molecular docking results using AutoDock Vina.15 Therefore, this study conducted a further analysis of the ligand binding to amino acid residues at ACE active sites and the ligand binding to Zn2+ ions, which both of these mechanisms are important in ACE inhibition.16 These two essential interactions were analyzed through molecular docking using AutoDock4Zn, which generates more accurate free energy binding calculations than AutoDock Vina on ligand-Zn interactions.17 Molecular dynamics simulations were also conducted, involving several compounds in melinjo seed within 20 ns to check the stability of the ligand-enzyme complexes.18

MATERIALS AND METHODS

Materials
Hardware

A single computer with the Linux operating system Ubuntu 12.04 LTS was used to conduct molecular docking and molecular dynamics simulations. Following specifications are used: Intel® Xeon E5620 (Intel® Core ™, American) processor (CPU) processors, Nvidia® GeForce GTX 780 (Nvidia, American) graphics processing unit (GPU) and 32 gigabytes (GB) Random Access Memory (RAM) DDR3.

Software

AutoDock4Zn forcefield was used to carry out the molecular docking, while AMBER PMEMD and VMD were used to conduct molecular dynamics simulations. LigandScout 4.09.2 and PoseView were also used to visualize the ligand-ACE binding from molecular docking and molecular dynamics simulations results.

Three-dimensional structure of ACE

The three-dimensional structure of ACE was obtained from the Research Collaboratory for Structural Bioinformatics Protein Data Bank (RCSB PDB) with identity (ID) 1O86. This protein is a testicular form of ACE (tACE) that contains C domain in the complex with lisinopril. Moreover, the macromolecule was optimized by removing the water molecules and undesired molecules, adding polar hydrogens and Gasteiger charges, and generating the Zn using zinc_pseudo.py.

Three-dimensional structure of ligands

The three-dimensional structure of ligands used in this study were: resveratrol, gnetol, isorhapontigenin, trans-piceid, gnetin C, trans-ε-viniferin, gnemonol K, gnemonol M, gnemonoside A, gnemonoside B, gnemonoside C, gnemonoside D, latifolol and aglycone of gnemonoside B. Lisinopril, captopril, and enalaprilat were also used as reference standards of free energy binding value and ligand-ACE interactions. All these ligands were obtained from http://www.pubchem.com. Charges were added to each ligand using AM1-BCC of Antechamber. Furthermore, all ligands were minimized using Sander to establish the relaxed structures. At last, all ligands were converted into pdbqt formats using python before performing the molecular docking.

Methods

Molecular docking was conducted using the AutoDock4Zn force field to generate the grid parameter files (*.gpf) of ligand-ACE complexes. The grid center was set up with the grid center coordinate of x=40.553, y=32.798, and z=47.286. The grid box size used was configured into 40x40x40 Å with 0,375 Å spacing. The molecular docking was performed for 100 times (ga_run = 100). Then, a scoring of free energy binding (ΔG) values of each ligand was performed based on the best energy and best cluster criteria. Furthermore, the visualization of ligand-ACE interactions was performed using LigandScout 4.09.2 and PoseView. The interactions of ligand-ACE complexes were analyzed based on pharmacophores which were compared with lisinopril, captopril, and enalaprilat.

Meanwhile, the ligand-ACE interactions based on visualization using LigandScout (Figure 1) 4.09.2 and PoseView (Figure 2) indicated that each compound isolated from melinjo seeds was able to bind the active site of ACE via hydrogen bonds or hydrophobic interactions. Subsequently, based on Table 2, several contents were able to bind Zn2+ ion via cation-pi interactions.

Figure 1

The visualization interactions of ligand-ACE complexes by using LigandScout 4.09.2.

https://s3-us-west-2.amazonaws.com/jourdata/jyp/JYoungPharm_10_2_s32_g000.jpg

Molecular dynamics simulations were carried out at temperatures of 300 K and 310 K within 20 ns for each ligand-ACE complex. The simulation showed that resveratrol, gnetol, isorhapontigenin, trans-piceid, gnetin C, trans-ε-viniferin, gnemonol K, gnemonol M, gnemonoside C, and aglycone of gnemonoside B exhibited the free energy binding (∆G) lower and close to lisinopril, captopril, and enalaprilat. Gnemonoside D and trans-piceid were not selected since these ligands were the glucoside of gnetin C and resveratrol respectively that might be hydrolyzed by the hydrolase enzyme. The temperature of 300 K was chosen since it has been the default temperature of molecular dynamics simulations, while the temperature of 310 K was selected in an attempt to follow the heat of a healthy human body. Afterward, each ligand-ACE complex was analyzed on RMSD, RMSF, the hydrogen bond occupancy, free energy binding (∆G) calculations using MM-GBSA and MM-PBSA methods and cation-pi interactions between ligand and Zn2+ ion.

RESULTS

The results of molecular docking method by using AutoDock4Zn (Table 1) showed that resveratrol, gnetol, isorhapontigenin, trans-piceid, gnetin C, trans-ε-viniferin, gnemonol K, gnemonol M, gnemonoside C, and aglycone of gnemonoside B exhibited free energy binding (∆G) values which were lower than or close to the ∆G values of lisinopril, captopril, and enalaprilat. Meanwhile, the ligand-ACE interactions based on visualization using LigandScout (Figure 1) 4.09.2 and PoseView (Figure 2) indicated that each compound isolated from melinjo seeds was able to bind the active site of ACE via hydrogen bonds or hydrophobic interactions. Subsequently, based on Table 2, several compounds were able to bind Zn2+ ion via cation-pi interactions.

Figure 2

The visualization of ligand-ACE complexes interactions by using PoseView.

https://s3-us-west-2.amazonaws.com/jourdata/jyp/JYoungPharm_10_2_s32_g001.jpg
Table 1

The results of free energy binding (∆G) calculations by using AutoDock4Zn.

No.LigandFree energy binding/∆G (kcal/mol)
1Resveratrol-7.05
2Gnetol-7.04
3Isorhapontigenin-7.19
4trans-Piceid-8.11
5Gnetin C-8.51
6trans-ε-Viniferin-8.13
7Gnemonol K-11.07
8Gnemonol M-8.26
9Gnemonoside A-6.1
10Gnemonoside B-1.55
11Gnemonoside C-6.14
12Gnemonoside D-7.29
13Latifolol-6.19
14Aglycone of gnemonoside B-10.68
15Lisinopril-10.05
16Captopril-7.1
17Enalaprilat-9.98

Molecular dynamics simulations at a temperature of 300 K (Figure 3) showed that all ligand-ACE complexes except resveratrol-ACE and trans-ε-viniferin-ACE complexes exhibited the similar pattern of RMSD fluctuation compared to lisinopril-ACE, captopril-ACE, and enalaprilat-ACE. Consequently, these ligands were able to maintain the conformational stability of ligand-ACE complexes to the same extent as lisinopril-ACE, captopril-ACE, and enalaprilat-ACE.

The RMSF values of molecular dynamics simulations on both temperature of 300 K and 310 K (Figure 4) indicated that the temperature difference did not affect the RMSF value of each ligand-ACE complex.

Figure 4

RMSF of ligand-ACE complexes based on molecular dynamics simulations within 20 ns at a temperature of 300 K (top) and 310 K down (bottom).

https://s3-us-west-2.amazonaws.com/jourdata/jyp/JYoungPharm_10_2_s32_g003.jpg

As shown in Table 3, almost all ligands were able to form hydrogen bonds with amino acids in the active site of ACE whose occupancy values were higher than 25% within 20 ns of molecular dynamics simulations at temperatures of 300 K and 310 K. Only resveratrol which did not exhibit the hydrogen bond occupancy greater or equal to 25% at a temperature of 300 K which demonstrated that the stability and the strength of resveratrol-ACE complex were poor.

Based on the MM-GBSA method at a temperature of 300 K shown in Table 4, the lowest ΔG value was exhibited by lisinopril indicating that there was no compound isolated from melinjo seeds which were able to match the affinity of lisinopril towards ACE. However, there were still ligands that had ΔG values ​​smaller than captopril and enalaprilat, such as gnetin C, trans-ε-viniferin, gnemonol K, and aglycone of gnemonoside B. Meanwhile, at a temperature of 310 K, the smallest ΔG was showed by the gnemonol K(-27.15 kcal/mol).

Table 4

Free energy binding (∆G) values of ligands to ACE based on calculations by using MM-GBSA, MM-PBSA, and AutoDock4Zn forcefield.

Ligand∆G (kcal/mol)
MM-GBSAMM-PBSAAutoDock 4Zn
300 K310 K300 K310 K
Resveratrol-6.49-8.64-3.34-5.15-7.05
Gnetol-17.08-14.27-11.62-13.49-7.04
Isorhapontigenin-15.14-13.29-4.00-6.04-7.19
Gnetin C-30.43-26.87-20.83-14.06-8.51
trans-ε-Viniferin-23.08-12.12-14.59-7.14-8.13
Gnemonol K-25.60-27.15-13.64-18.02-11.07
Gnemonol M-18.94-23.27-6.260.00-8.26
Aglycone of gnemonoside B-28.11-22.83-18.53-16.90-10.68
Lisinopril-34.62-18.25-25.22-11.39-10.15
Captopril-14.46-19.40-9.86-14.89-7.1
Enalaprilat-22.63-17.12-12.23-16.36-9.98

Based on the calculation of ΔG by MM-PBSA method at 300 K shown in Table 4, it can be concluded that lisinopril was the ligand with the smallest ΔG value. Moreover, in the calculation of ΔG with the MM-PBSA method at a temperature of 310 K, gnemonol K and the aglycone of gnemonoside B exhibited lower ΔG value than lisinopril, captopril, and enalaprilat.

Based on Figure 5, all ligands except resveratrol were able to form relatively stable distance to the Zn2+ ion. Resveratrol exhibited a very fluctuating distance to the Zn2+ ion during molecular dynamics simulations within 20 ns at a temperature of 300 K. Meanwhile, at a temperature of 310 K, all ligands were able to form relatively stable distance to the Zn2+ ion, including resveratrol.

Figure 5

The distances between the aromatic ring of ligands with Zn2+ ion within 20 ns of molecular dynamics simulations at a temperature of 300 K (top) and 310 K (down).

https://s3-us-west-2.amazonaws.com/jourdata/jyp/JYoungPharm_10_2_s32_g005.jpg
Figure 6

The 3D visualization of cation-pi interactions of ligand-ACE complexes at temperatures of 300 K (top) and 310 K (bottom) using LigandScout 4.09.

https://s3-us-west-2.amazonaws.com/jourdata/jyp/JYoungPharm_10_2_s32_g004.jpg

DISCUSSION

Molecular docking method using AutoDock4Zn demonstrated that resveratrol, gnetol, isorhapontigenin, trans-piceid, gnetin C, trans-ε-viniferin, gnemonol K, gnemonol M, gnemonoside C, and aglycone of gnemonoside B exhibited ∆G values which were lower than or close to positive controls (Table 1). The lower of free energy binding (∆G) values, the better affinity of ligand towards its receptor. Ligands can inhibit ACE catalytic mechanism via ligand-Zn2+ ion or ligand-amino acids interactions in the active site of ACE. Interestingly, resveratrol, gnetol, isorhapontigenin, trans-piceid, gnetin C, trans-ε-viniferin, gnemonoside C were able to bind the Zn2+ ion (Table 2). Therefore, based on the ∆G and visualization interactions result, it suggested that resveratrol, gnetol, isorhapontigenin, trans-piceid, gnetin C, trans-ε-viniferin, gnemonol K, gnemonol M, gnemonoside C, and aglycone of gnemonoside B potentially have ACE inhibitory activity.

Table 2

The visualization interactions summary of ligand-ACE interactions by using LigandScout 4.09.2 and PoseView compared to lisinopril, captopril, and enalaprilat.

The Active Site of ACE LigandGlu162Gln281His353Ala354Glu384Glu411Asp415Asp453Lys454Lys511His513Tyr520Tyr523Zn2+ ion
Resveratrol        
Gnetol        
Isorhapontigenin         
trans-Piceid      
Gnetin C          
trans-ɛ-Viniferin      
Gnemonol K          
Gnemonol M        
Gnemonoside A         
Gnemonoside B           
Gnemonoside C         
Gnemonoside D           
Latifolol          
Aglycone of gnemonoside B          
Lisinopril    
Captopril       
Enalaprilat     

Analysis of RMSD values of ligand-ACE complexes from molecular dynamics simulations was performed to observe changes in protein conformation within a time range of molecular dynamics simulations.19

Based on Figure 3, there were several ACE-ligand complexes which showed considerable differences in RMSD values between the molecular dynamics simulations at temperatures of 300 K and 310 K, such as the resveratrol-ACE, trans-ε-viniferin-ACE, aglycone of gnemonoside B-ACE, and enalaprilat-ACE. Some of them exhibited the lower RMSD at a temperature of 300 K, whereas some others showed the lower RMSD at a temperature of 310 K. Therefore, it can be indicated that the temperature affected to the RMSD of the ACE-ligand complexes.

Figure 3

RMSD of ligand-ACE complexes based on molecular dynamics simulations within 20 ns at a temperature of 300 K (top) and 310 K (bottom).

https://s3-us-west-2.amazonaws.com/jourdata/jyp/JYoungPharm_10_2_s32_g002.jpg

Analysis of RMSF values aims to observe the flexibility of proteins occurring in each amino acid residue during molecular dynamics simulations.20 We focused on the amino acid residues at the active site of ACE, including Gln281, His353, Ala354, Glu384, Glu411, Asp415, Asp453, Lys454, Lys511, His513, Tyr520, and Tyr523.

Based on the graphs in Figure 4, it demonstrated that at temperatures of 300 K and 310 K all ligands that bonded with amino acid residues at the active sites of ACE did not exhibit a high RMSF value, so ligands were able to form stable conformation at the active site of ACE. The temperature differences did not affect the RMSF value between ligands and the active site of ACE.

The hydrogen bond occupancy calculations were performed to assess the strength and stability of hydrogen bonds within 20 ns of molecular dynamics simulations. Based on the hydrogen bond strength, the occupancy percentage is divided into three ranges: very weak hydrogen bond (25-50%), likely hydrogen bond (50-75%), and strong hydrogen bond (75-100%).21 A hydrogen bond is considered to be stable when it has an occupancy higher than 50%.22 The hydrogen bond occupancy rate was calculated by using VMD with 3.0 Å of the bond cut off distance, and 60° of the bond cut off angle.

Based on Table 3, it can be considered that the hydrogen bond occupancy was likely to change by the 10 K temperature differences. Most ligands exhibited higher occupancy value towards amino acids residues in the active site of ACE at a temperature of 310 K rather than 300 K. It can be concluded that the ACE inhibitory activities of these ligands were better at an average human body temperature. Otherwise, only gnemonol K-ACE and captopril-ACE complexes which had better occupancy at a temperature of 300 K. The hydrogen bond occupancy between enalaprilat and Glu411 at temperatures of 300 K and 310 K were higher than 100%, which indicated that there were two hydrogen bonds occurred between them at the same time.

Table 3

The comparison of hydrogen bond occupancy of ligand-ACE complexes between the molecular dynamics simulations at temperatures of 300 K and 310 K.

LigandHydrogen Bond Donor 300 KHydrogen Bond Acceptor 300 KOccupancy (%) 300 KHydrogen Bond Donor 310 KHydrogen Bond Acceptor 310 KOccupancy (%) 310 K
Resveratrol---LIG580-Side-O1ALA315-Main-O29.40
GnetolLIG-Side-O3GLU411-Side-OE280.60LIG-Side-O2ALA354-Main-O97.30
LIG-Side-O3GLU411-Side-OE293.30
LIG-Side-O1GLU62-Side-OE272.10
TYR523-Side-OHLIG-Side-O343.65
LIG-Main-OGLU376-Side-OE243.35
 LYS511-Side-NZLIG-Main-O29.95
IsorhapontigeninLYS511-Side-NZLIG-Side-O333.40HIS513-Side-NE2LIG-Side-O232.95
   LIG-Side-O3GLU384-Side-OE229.50
Gnetin CLIG-Side-O3ASP377-Side-OD289.50LIG-Side-O3ASP377-Side-OD159.90
ALA354-Main-NLIG-Main-O26.90LIG-Side-O2HIS513-Side-ND151.80
ALA354-Main-NLIG-Main-O48.55
LIG-Side-O3ASP377-Side-OD240.70
LIG-Side-O1GLU162-Side-OE234.15
    TYR146-Side-OHLIG-Side-O125.40
trans-ɛ-ViniferinLIG-Main-OALA354-Main-O93.80LIG-Side-O5ASP377-Side-OD178.45
LIG-Side-O4GLU162-Side-OE255.25LIG-Side-O1ASP453-Side-OD239.30
LIG-Side-O1ASP453-Side-OD254.10LIG-Side-O1ASP453-Side-OD131.50
LIG-Side-O4GLU162-Side-OE147.40LIG-Side-O4GLU162-Side-OE226.20
LIG-Side-O1ASP453-Side-OD133.15
THR166-Side-OG1LIG-Side-O530.30
 LYS511-Side-NZLIG-Side-O325.50   
LIG-Side-O4GLU384-Side-OE149.05LIG-Side-O6GLU411-Side-OE274.90
ALA356-Main-NLIG-Side-O448.20LIG-Side-O4GLU384-Side-OE128.20
GLN281-Side-NE2LIG-Side-O732.45LYS511-Side-NZLIG-Side-O728.15
 LIG-Side-O8GLU162-Side-OE127.60   
Gnemonol MLIG-Side-O4ASP453-Side-OD281.10LIG-Side-O4ASP453-Side-OD286.50
LIG-Side-O3ASP415-Side-OD162.55LIG-Side-O5ASP453-Main-O67.25
LIG-Side-O5ASP453-Main-O47.85LIG-Side-O1ALA354-Main-O59.20
LIG-Side-O5THR282-Side-OG134.55LIG-Side-O3ASP415-Side-OD148.00
    LIG-Side-O2ASP415-Side-OD229.60
Gnemonol KLIG-Side-O1ASP415-Side-OD195.95LIG-Side-O7CYS352-Main-O79.00
Aglycone of gnemonoside BLIG-Side-O4GLU384-Side-OE186.35LIG-Side-O5ALA354-Main-O98.00
LIG-Side-O3GLN369-Side-OE180.95LIG-Side-O4GLU411-Side-OE287.05
LIG-Side-O2GLU162-Side-OE133.85
LIG-Side-O2GLU162-Side-OE231.55
    LIG-Side-O3VAL380-Main-O25.90
LisinoprilTYR523-Side-OHLIG-Side-O156.40LIG-Side-O3GLU411-Side-OE298.65
LIG-Side-O3ALA354-Main-O50.35TYR520-Side-OHLIG-Side-O464.60
HIS513-Side-NE2LIG-Side-O132.55HIS513-Side-NE2LIG-Side-O147.60
 LYS511-Side-NZLIG-Side-O428.75HIS353-Side-NE2LIG-Side-O127.05
CaptoprilTYR520-Side-OHLIG-Side-O261.65TYR520-Side-OHLIG-Side-O234.70
    LYS511-Side-NZLIG-Side-O232.45
EnalaprilatLIG-Side-O2GLU411-Side-OE2132.25LIG-Side-O2GLU411-Side-OE2168.95
 TYR523-Side-OHLIG-Side-O227.50HIS353-Side-NE2LIG-Side-N140.00

MM-GBSA and MM-PBSA methods were capable of calculating the free energy binding (∆G) of a ligand towards a macromolecule, either in a solvated state or a vacuum. Based on Table 4, both MM-GBSA and MM-PBSA methods exhibited different free energy binding (∆G). Also, the 10 K temperature differences led the significant changes in free energy binding (∆G) of some ligands both in the MM-GBSA and MM-PBSA methods. Both methods were able to generate more accurate free energy binding (∆G) values than docking and scoring methods.23

Fluctuations in interaction distance of aromatic ring in ligands with Zn2+ ion were observed to assess the stability of the cation-pi interactions that occurred during the simulation of 20 ns molecular dynamics, both at 300 K and 310 K. Then, visualization process was carried out at the closest distance between aromatic ring of ligands with Zn2+ by using LigandScout 4.09.2 to identify the presence of cation-pi interactions. A cation-pi interaction can occur within a distance of 6.0 Å.24

Based on the graph in Figure 5, resveratrol showed a very fluctuating distance to Zn2+ ion after passing a three ns simulation time which indicated that the cation-phi bond between aromatic ring in resveratrol to Zn2+ ion only lasted for a short time. Meanwhile, aromatic rings in other ligands were stable at a distance range of 4-7 Å. Based on the visualization process results by using LigandScout 4.09.2 at the closest distance of aromatic ring of each ligand to Zn2+ ion, resveratrol, gnetol, isorhapontigenin, gnetin C, gnemonol K, and aglycone of gnemonoside B were able to form cation-pi interactions with the Zn2+ ion at a temperature of 300 K. However, the aromatic ring of resveratrol showed a closer distance with Zn2+ at a temperature of 310 K compared to 300 K. Meanwhile, the other ligands did not show significant differences in fluctuation of the aromatic ring distance to Zn2+ ion compared to the 300 K molecular dynamics simulations. Visualization process by using LigandScout 4.09.2 at the aromatic rings of each ligand to Zn2+ at the closest distance between them, indicated that resveratrol, isorhapontigenin, gnetin C, gnemonol K, gnemonol M, and aglycone gnemonoside B were able to form cation-phi interactions at 310 K.

It is necessary to conduct the molecular docking and molecular dynamics analysis of other substances contained in melinjo seed extracts that have not been carried out in this study. Subsequently, it is needed to lengthen the simulation times of molecular dynamics to see further about ligand-ACE conformational changes. Preclinical and clinical tests are also required to perform for proving the ACE inhibitory activity of each substance. Furthermore, toxicity testing of each compound in melinjo seed extract is needed to evaluate the safety of each material when itis consumed.

CONCLUSION

Based on the in silico study through molecular docking and molecular dynamics simulations of several compounds isolated from Gnetum gnemon L., it may be concluded that resveratrol, gnetol, isorhapontigenin, gnetin C, trans-ε-viniferin, gnemonol K, gnemonol M, and aglycone of gnemonoside B potentially have ACE inhibitory activity due to their free energy binding (∆G), RMSD, RMSF, and the hydrogen bond occupancy compared to lisinopril, captopril, and enalaprilat. It was convinced by the visualization process by using LigandScout 4.09.2 and PoseView on each ligand-ACE complex which showed the binding pattern similarities between the compounds and positive controls, including lisinopril, captopril, and enalaprilat.

ACKNOWLEDGEMENT

The authors thank to Publikasi Internasional Terindeks Untuk Tugas Akhir Mahasiswa (PITTA) 2017 provided by the Directorate of Research and Community Engagement, Universitas Indonesia for the financial support to this study.

CONFLICT OF INTEREST

The authors declare no conflict of interest.

ABBREVIATIONS

ACE

Angiotensin Converting Enzyme

MM-GBSA

Molecular Mechanics Generelized Born Surface Area

MM-PBSA

Molecular Mechanics Poisson-Boltzmann Surface Area

PMEMD

Patial Mesh Ewal Molecular Dynamics

VMD

Visual Molecular Dynamics

RMSD

Root Mean Square Deviation; RMSF: Root Mean Square Fluctuation.

REFERENCES

1. 

Narayanan NK , author. et al. Antitumor activity of melinjo (Gnetum gnemon L.) seed extract in human and murine tumor models in vitro and in a colon-26 tumor-bearing mouse model in vivo. Cancer Med. 2015;4(11):1767–80

2. 

Puspitaningrum YT, Efendi E, Siswoyo TA , authors. Analisis in vivo aktivitas antihipertensi dari protein biji melinjo (Gnetum gnemon) terhidrolisis. Artikel Ilmiah Hasil Penelitian Mahasiswa. 2014;2(2):327–31

3. 

Wells BG, DiPiro JT, Schwinghammer TL, DiPiro CV , authors. Pharmacotherapy Handbook. 9. New York: McGraw-Hill; 2015

4. 

World Health Organization , author. World health statistics. 2012 March 15–2017. http://www.who.int/gho/publications/world_health_statistics/EN_WHS2012_Full.pdf

5. 

Lestari P, Katrin ; Rissyelly ; Maya PG , authors. Inhibition Activity of Angiotensin Converting Enzyme (ACE) and Determination of Total Phenolic and Flavonoid Compound from Bitter Melon Leaves (Momordica charantia L.). Pharmacogn J. 2017;9(2):252–6

6. 

Badan Penelitian dan Pengembangan Kementerian Kesehatan RI , author. Riset Kesehatan Dasar Tahun 2013. Jakarta: Riskesdas. 2013

7. 

Sweileh WM, Sawalha AF, Zyoud SH, Al-Jabil SW, Tameem EJ, Shraim NY , authors. Evaluation of antihypertensive therapy in diabetic hypertensive patients: Impact of ischemic heart disease. J Pharm Pract. 2013;7(1):40–6

8. 

Tutor JT, Chichioco-Hernandez CL , authors. Angiotensin-Converting Enzyme Inhibition of Fractions from Eleusine indica Leaf Extracts. Pharmacogn J. 2018;10(1):25–8

9. 

Manner HI, Elevitch CR , authors. Traditional trees of Pacific Island: Their culture, environment, and use. Hawaii: Permanent Agriculture Resources. 2006

10. 

Huang KS, Li RL, Wang YH, Lin M , authors. Three new stilbene trimers from the lianas of Gnetum hainanense. Planta Medica. 2001;67:61–4

11. 

Huang KS, Wang YH, Li RL, Lin M , authors. Five new stilbene dimers from the lianas of Gnetum hainanense. Journal of National Product. 2000;63(1):86–9

12. 

Huang KS, Wang YH, Li RL, Lin M , authors. Stilbene dimers from the lianas of Gnetum hainanense. Phytochemistry. 2000;63(1):86–9

13. 

Hafidz KA , author. et al. HMG-CoA Reductase Inhibitory Activity of Gnetum gnemon Seed Extract and Identification of Potential Inhibitors for Lowering Cholesterol Level. J Young Pharm. 2017;9(4):559–65

14. 

Munadhil MA , author. Uji aktivitas penghambatan angiotensin converting enzyme (ACE) dan penetapan kadar fenolik total dari biji melinjo (Gnetum gnemon L.). Depok: Fakultas Farmasi Universitas Indonesia. 2016

15. 

Mun’im A, Munadhil MA, Puspitasari N, Azminah ; Yanuar A , authors. Angiotensin-converting enzyme inhibitory activity of melinjo (Gnetum gnemon L.) Seed extracts and molecular docking of its stilbene constituents. Asian J Pharm Clin Res. 2017;10(3):243–8

16. 

Shukor NA , author. et al. Angiotensin-converting enzyme inhibitory effects by plant phenolic 2 compounds: A study of structure-activity relationships. J Agric Food Chem. 2013;61(48):11832–9

17. 

Santos-Martins D, Forli S, Ramos MJ, Olson AJ , authors. AutoDock4Zn: An improved AutoDock force field for small-molecule docking to zinc metalloproteins. J Chem Inf Model. 2014;54(8):2371–9

18. 

Priya R , author. et al. Molecular docking and molecular dynamics to identify a novel human immunodeficiency virus inhibitor from alkaloids of Toddalia asiatica. Pharmacogn Mag. 2015;11(4):414–22

19. 

Salsbury FR , author. Molecular dynamics simulations of protein dynamics and their relevance to drug discovery. Curr Opin Pharmacol. 2010;10(6):738–44

20. 

Fraccalvieri D , author. et al. Conformational and functional analysis of molecular dynamics trajectories by self-organising maps. BMC Bioinformatics. 2011;12(1):158

21. 

Kästner J, Loeffler HH, Roberts SK, Fernandez MI, Winn MD , authors. Ectodomain orientation, conformational plasticity and oligomerization of erbB1 receptors investigates by molecular dynamics. J Struct Biol. 2009;167(2):117–28

22. 

Desheng L, Jian G, Yuanhua C, Wei C, Huai Z, Mingjuan J , authors. Molecular dynamics simulations and MMP/GBSA methods to investigate binding mechanism of aminomethylpyrimidine inhibitor with DPP-IV. Bioorg Med Chem Lett. 2011;21(22):6630–5

23. 

Genheden S, Ryde U , authors. The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities. Expert Opin Drug Discov. 2015;10(5):449–61

24. 

Gallivan JP, Dougherty DA , authors. Cation-pi interactions in structural biology. Proc Natl Acad Sci USA. 1999;96(17):9459–64