|Year : 2021 | Volume
| Issue : 3 | Page : 250-253
Encapsulation mechanism of α-mangostin by β-cyclodextrin: Methods of molecular docking and molecular dynamics
Dian Triwahyuningtyas1, Sandra Megantara1, Tai Tze Hong1, Muhammad Yusuf2, Muchtaridi Muchtaridi1
1 Department of Pharmaceutical Analysis and Medicinal Chemistry, Faculty of Pharmacy, Universitas Padjadjaran, Sumedang, Indonesia
2 Department of Chemistry, Faculty of Natural and Matrh Sciences, Universitas Padjadjaran, Sumedang, Indonesia
|Date of Submission||09-Dec-2020|
|Date of Decision||16-Feb-2021|
|Date of Acceptance||04-May-2021|
|Date of Web Publication||16-Jul-2021|
Prof. Muchtaridi Muchtaridi
Department of Pharmaceutical Analysis and Medicinal Chemistry, Faculty of Pharmacy, Universitas Padjadjaran, Jl. Raya, Jatinangor KM 21.5 Bandung-Sumedang, 45363
Source of Support: None, Conflict of Interest: None
The study aimed to investigate the interaction of host-guest between α-mangostin and β-cyclodextrin (βCD) and also to calculate the energy of the complex system between α-mangostin with βCD for drug delivery using methods of 15 molecular dynamics and molecular docking. Simulation of molecular docking and molecular dynamics was utilized to determine molecular interactions and the complex system's bond energy. The docking simulation results showed that α-mangostin-βCD complex has a Gibbs energy value (ΔG) of −6.69 kcal/mol. The Gibbs energy value (ΔG) of molecular dynamics simulation from MMGBSA calculation showed the binding energy of α-mangostin-βCD – 11.73 kcal/mol.
Keywords: Inclusive complex, molecular docking, molecular dynamics, α-mangostin, β-cyclodextrin
|How to cite this article:|
Triwahyuningtyas D, Megantara S, Hong TT, Yusuf M, Muchtaridi M. Encapsulation mechanism of α-mangostin by β-cyclodextrin: Methods of molecular docking and molecular dynamics. J Adv Pharm Technol Res 2021;12:250-3
|How to cite this URL:|
Triwahyuningtyas D, Megantara S, Hong TT, Yusuf M, Muchtaridi M. Encapsulation mechanism of α-mangostin by β-cyclodextrin: Methods of molecular docking and molecular dynamics. J Adv Pharm Technol Res [serial online] 2021 [cited 2021 Jul 23];12:250-3. Available from: https://www.japtr.org/text.asp?2021/12/3/250/321515
| Introduction|| |
The fundamental xanthone compound in the pericarp of the mangosteen is α-mangostin. The mangosteen fruit has attracted a lot of attentions from researchers due to its various pharmacological activities such as antimicrobial, anti-inflammatory, antiviral, anticancer, and antifungal.,,, However, the solubility of α-mangostin in water is only 0.2 μg/ml and it has largely restricted its bioavailability.
Cyclodextrins are cyclic oligosaccharide which contains 6, 7, and 8 glucopyranose units within and are also related to α-cyclodextrin, β-cyclodextrin (βCD), and γ-cyclodextrin, respectively. Cyclodextrins are commonly used in the pharmaceutical field because it has a unique hollow structure., The inclusion complexation cyclodextrin is a technique that is common in enhancing the solubility of poorly water-soluble drugs. Compared to other natural cyclodextrins such as α-cyclodextrin and γ-cyclodextrin, βCD is often used because it can be easily synthesized and the price is cheaper. Cyclodextrins have the ability to form inclusive complexes with various compounds and thus it can help to improve the physicochemical properties of the complex compound. By forming an inclusive complex, the complex will dissolve in the solute and achieve dynamic equilibrium quickly.
βCD is one of the practical cyclodextrins. However, the powerful intramolecular hydrogen bonds of βCD cause its low solubility level in the water, which is only 1.85 g/100 mL at 25°C. According to the studies conducted by Das et al. in 2011 and Zhao et al. in 2012, the formation of inclusive complex between resveratrol and oxyresveratrol with cyclodextrin can increase the solubility of the complex in the water, especially when it is complexed with 2-hydroxypropyl-βCD (HPβCD), which is a derivative of βCD.,
This study aimed to inspect molecular interactions, as well as to calculate the complex system energy between α-mangostin with βCD for drug delivery using methods of molecular dynamics and molecular docking.
| Materials and Methods|| |
Molecular structures constructions
This study used AutoDock 4.2 (The Scripps Research Institute, La Jolla, CA, USA) to estimate the possible arrangement of complex α-mangostin-βCD and α-mangostin-HPβCD. From the Protein Data Bank (ID: 1z0n), the formation of βCD was acquired. α-mangostin structure as a guest molecule was obtained from PubChem and optimized by ChemBio3D Ultra 12.0 (PerkinElmer Inc.) The nonpolar hydrogen atoms from both receptors and guest molecule were emerged, and the Gasteiger charges were added. The grid box used was 60 × 60 × 60 points, featuring grid space of 0.375 Å. The docking parameter used was default setting from AutoDock 4.2, except docking runs, energy evaluations, and number of generations that were adjusted to 100, 2,500,000, and 250, respectively. Conformation results of the docking were grouped by root-mean-square deviation (RMSD) with a tolerance of 2.0 Å. The best docking result was visualized with BIOVIA Discovery Studio Visualizer 2017.
Molecular dynamics simulation
Particle Mesh Ewald Molecular Dynamics from package AMBER 16 was used to perform the simulation of molecular dynamic. The study used the general AMBER force field to acquire the parameters of force field for α-mangostin and βCD molecules using semi-empirical quantum calculation AM1-BCC through an antechamber program. The complexation of α-mangostin and βCD s was performed with tleap and the system was immersed in water using TIP3P water model with a periodic box size of 10 Å. In the preparation system phase, the system was minimized for 9000 steps, including the steepest descent (7000 steps) and conjugate gradient (2000 steps) with cutoff value of 9Å in constant volume periodic boundaries. Next, the system was kept hot under steady mass for 60 ps to 310 K by utilizing a Langevin thermostat with restraints of 5 kcal mol − 1 Å−2. Further on, the system was equilibrated under constant pressure for 500 ps with constant pressure periodic boundary. SHAKE algorithm was applied to constraint hydrogen atoms at their equilibrium distance. In the production phase, a 30ns simulation was done with constant pressure periodic boundary. Trajectory analysis was done with cpptraj module from AMBER 16. The binding energy between the βCD and α-mangostin was measured by Molecular Mechanics-Generalized Born Surface Area (MMGBSA).
| Results|| |
The structure of guest molecule α-mangostin with host molecules βCD obtained is shown in [Figure 1] and [Figure 2].
Molecular docking was implemented to find the best position of α-mangostin inside the cavity of βCD based on the lowest binding energy in the largest cluster.
The results of molecular docking are shown in [Table 1]. From the table, the free binding energy of α-mangostin-βCD was −6.69.
|Table 1: The results of molecular docking of α-mangostin and β-cyclodextrins at 298.15 K in docking 176th|
Click here to view
Molecular dynamic simulation
The best docking complex from molecular docking was further studied with molecular dynamics. MMGBSA was used to measure the inclusion complexes of α-mangostin-cyclodextrin in molecular dynamic simulation.The results of molecular dynamics are shown in [Table 2]. From the table, the free binding energy of α-mangostin-βCD was −11.73 kcal/mol.
| Discussion|| |
The docking process was repeated for 250 times. The grid coordinates for βCD was −6.686; 37.14; −7,932 (x; y; z) and HPβCD −7,214; 37.154; −6.91 (x; y; z) with grid box size of 60 × 60 × 60, to cover the entire surface of the cyclodextrin. From the docking results, no case of α-mangostin coming out of the cyclodextrin cavity during the docking process was found. This indicated that α-mangostin can establish a stable inclusive complex with βCD theoretically. The 176th docking was chosen in the analysis.
The docking results of α-mangostin and βCD [Figure 3] showed that there were 2 hydrogen bonds in the complex with the distance of 2.07Å and 2.10 Å, respectively. There was also a hydrophobic effect that occurred at a distance of 3.90 Å between the aromatic group α-mangostin and C-H from βCD by Pi-Sigma. The average value of binding energy was −6.69 kcal/mol. The negative binding energy values suggested that the formation of all these inclusive complexes was a spontaneous reaction.
|Figure 3: Interaction of α-mangostin (a) from top (b) from front (c) insight look with β-cyclodextrins. Hydrogen bond (green dotted line) and hydrophobic effect (pink dotted line) are shown above|
Click here to view
The snapshots during 30 ns molecular dynamics simulation of α-mangostin and βCD [Figure 4]. The stability of a complex system can be known by calculating RMSD where RMSD is the average difference between atomic positions in a simulation. In the complex system of α-mangostin-βCD, βCD and α-mangostin were stable throughout 30 ns molecular dynamics simulation [Figure 5], while in the complex system of α-mangostin-HPβCD, α-mangostin has stable RMSD graph line throughout the simulation. The system was said to be stable during this simulation.
|Figure 4: α-mangostin- β-cyclodextrins from front view at (a) 1 ns (b) 10 ns (c) 20 ns (d) 30 ns structure|
Click here to view
|Figure 5: The root-mean-square deviation plot of backbone atoms for the simulation of α-mangostin- β-cyclodextrins|
Click here to view
Based on a snapshot of βCD-α-mangostin simulation system, it appeared that only methyl group entered the cavity of cyclodextrin [Figure 4]. The possible reason for this was the large methyl group caused steric hindrance in the βCD cavity. This prevented the other part of α-mangostin from entering into the cavity during the simulation. In addition, the βCD structure was more or less the same throughout the simulation which was in agreement with the RMSD graph which showed that there was not much fluctuation throughout the simulation [Figure 5].
Table 2 shows the results of molecular dynamic simulation of α-mangostin and βCD at 310K in the system water model TIP3P. The study implemented MMGBSA approach to measure binding free energy. Every 100 frames out of 3000 total frames were implemented to measure binding free energy. Van der Waals force (ΔEVDW) made the key contribution in the formation of inclusive complex. In the βCD-α-mangostin complex, the VDW value was −33.74 kcal/mol. From the table, the low value of VDW force indicated the cavity of cyclodextrin was hydrophobic. The total value of binding energy (ΔG) in the βCD-α-mangostin system was − 11.73 kcal/mol.
| Conclusion|| |
In the study, two modeling methods were used to study the complexation of α-mangostin and cyclodextrins. The docking simulation results showed that α-mangostin-βCD complex has a Gibbs energy value (ΔG) of −6.69 kcal/mol. The Gibbs energy value (ΔG) of molecular dynamics simulation from MMGBSA calculation showed the binding energy of α-mangostin-βCD −11.73 kcal/mol. The results showed that α-mangostin-βCD inclusion was a stable and spontaneous process.
Researchers would sincerely and gratefully acknowledge the Rector of Padjadjaran University for funding this project through Academic Leadership Grants 2020 no. 1427/UN6.3.1/LT/2020 from Universitas Padjadjaran.
Financial support and sponsorship
This study was financially supported by the Grant of Thesis Research 2020 (PDD) from Ministry of Research and Technology/BRIN 1827/UN6.3.1/LT/2020 and Academic Leadership Grants 2020 no. 1427/UN6.3.1/LT/2020 from Universitas Padjadjaran.
Conflicts of interest
The author proclaims no conflict of interest.
| References|| |
Nguyen PT, Marquis RE. Antimicrobial actions of α-mangostin against oral streptococci. Can J Microbiol 2011;57:217-8.
Gutierrez OF, Chitchumroonchokchai C, Lesinski GB, Suksamrarn S, Failla ML. α-Mangostin: Anti-Inflammatory activity and metabolism by human cells. J Agric Food Chem 2013;61:3891-9.
Shan T, Ma Q, Guo K, Liu J, Li W, Wang F, et al.
Xanthones from mangosteen extracts as natural chemopreventive agents: Potential anticancer drugs. Curr Mol Med 2011;11:666-77.
Kaomongkolgit R, Jamdee K, Chaisomboon N. Antifungal activity of alpha-mangostin against Candida albicans
. J Oral Sci 2009;51:401-5.
Jadhav P, Petkar B, Pore Y, Kulkarni A, Burade K. Physiochemical and molecular modeling studies of cefixime-L-arginine-cyclodextrin ternary inclusion compounds. Carbohydr Polym 2013;98:1317-25.
Szente L, Szejtli J. Cyclodextrins as food ingredients. Trends Food Sci Technol 2004;15:137-5.
Del Valle EM. Cyclodextrins and their uses: A review. Process Biochem 2004;39:1033-13.
Matencio A. García-Carmona F. López-Nicolás JM. Encapsulation of piceatannol, a naturally occurring hydroxylated analogue of resveratrol, by natural and modified cyclodextrins. Food Funct 2016;7:2367-6.
Liu B, Li Y, Xiao H, Liu Y, Mo H, Ma H, et al.
Characterization of the supermolecular structure of polydatin/6-O-α-Maltosyl-β-cyclodextrin inclusion complex. J Food Sci 2015;80:C1156-61.
Das S, Ng KY. Physicochemical characterization, molecular modeling, and stability of the resveratrolcyclodextrin inclusion complexes. Lat Am J Pharm 2011;30:874-7.
Zhao WJ, Liu GY, Xia Q. Preparation and characterization of resveratrol-loaded hydroxylpropyl-βcyclodextrin. Adv Mater Res 2012; 490-495:3454-8.
Dermawan D, Wathoni N, Muchtaridi M. Host-guest interactions of α− Mangostin with (α, β, γ)−Cyclodextrins: Semi-empirical quantum mechanical methods of PM6 and PM7. J Young Pharm 2019;11:31-8.
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5]
[Table 1], [Table 2]