|Year : 2021 | Volume
| Issue : 3 | Page : 261-266
In silico identification of natural products from Centella asiatica as severe acute respiratory syndromecoronavirus 2 main protease inhibitor
Putu Gita Maya1, Widyaswari Mahayasih1, Harizal1, Herman2, Islamudin Ahmad3
1 Department of Pharmacy, Faculty of Health Sciences, Universitas Esa Unggul, West Jakarta, DKI Jakarta, Indonesia
2 Laboratory of Pharmaceutical Research and Development of FARMAKA TROPIS, Faculty of Pharmacy, Universitas Mulawarman, Samarinda, Indonesia
3 Department of Pharmaceutical Sciences, Faculty of Pharmacy, Universitas Mulawarman, Samarinda, East Kalimantan, Indonesia
|Date of Submission||08-Dec-2020|
|Date of Decision||08-Feb-2021|
|Date of Acceptance||28-Mar-2021|
|Date of Web Publication||16-Jul-2021|
Dr. Islamudin Ahmad
Department of Pharmaceutical Sciences, Faculty of Pharmacy, Universitas Mulawarman, Samarinda, East Kalimantan
Putu Gita Maya
Department of Pharmacy, Faculty of Health Sciences, Universitas Esa Unggul, West Jakarta, DKI Jakarta
Source of Support: None, Conflict of Interest: None
Severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2) main protease (S-CoV-2 Mpro) is one of the main targets in designing antiviral against SARS-CoV-2. Centella asiatica contains several triterpenoids, polyacetylenes, and benzoic ester derivative with various biological activities including anti-inflammation and antiviral. Triterpenoids from C. asiatica could act as inhibitors of S-CoV-2 Mpro. The main objective of this study was to identify potential natural products from C. asiatica as S-CoV-2 Mpro inhibitor with better pharmacokinetic through in silico molecular docking method. : As much as 11 compounds from C. asiatica were docked with S-CoV-2 Mpro (PDB ID: 6LU7) using AutoDock v4.2.6. Pharmacokinetic parameters of these compounds were assessed using SwissADME (free access webserver). Molecular docking results of 11 natural products indicated that asiatate 6 and asiatate 10 have strong interaction with quite similar binding free energy compared to native ligand (‒9.00 and‒9.58 kcal/mol compared to ‒9.18 kcal/mol, respectively) with proper interaction to the catalytic dyad (His41 and Cys145). Pharmacokinetic analysis revealed that asiatate 4, asiatate 10, and asiatate 11 have poor pharmacokinetic properties. These results indicated that asiatate 6 could be recommended for further study as S-CoV-2 Mpro inhibitor.
Keywords: Asiaticoside, coronavirus disease-2019, isothankunic acid, molecular docking, polyacetylenes, triterpenoids
|How to cite this article:|
Maya PG, Mahayasih W, Harizal, Herman, Ahmad I. In silico identification of natural products from Centella asiatica as severe acute respiratory syndromecoronavirus 2 main protease inhibitor. J Adv Pharm Technol Res 2021;12:261-6
|How to cite this URL:|
Maya PG, Mahayasih W, Harizal, Herman, Ahmad I. In silico identification of natural products from Centella asiatica as severe acute respiratory syndromecoronavirus 2 main protease inhibitor. J Adv Pharm Technol Res [serial online] 2021 [cited 2022 Oct 4];12:261-6. Available from: https://www.japtr.org/text.asp?2021/12/3/261/321514
| Introduction|| |
Coronavirus disease-2019 outbreak is caused by severe acute respiratory syndrome coronavirus-2 (S-CoV-2), a betacoronavirus strain. Development of antiviral drugs has been performed by targeting certain proteins either in the virus including structural or nonstructural proteins in S-CoV-2 or host targets.
Main protease (Mpro) or 3-chymotrypsin-like protease is a nonstructural protein that has a significant role in viral replication and governing the response of the cell host. Together with papain-like protease, Mpro cleaves and transforms two large polyproteins (pp1a and pp1ab) translated from the viral RNA selectively at the Leu-Gln↓(Ser, Ala, Gly) cleavage site. Due to its crucial role and specific cleavage site, inhibiting the activity of Mpro enzyme would stop viral replication with no undesired toxic effect on human.
Several potential synthetic drugs have been evaluated clinically against S-CoV-2 such as ritonavir/lopinavir and cobicistat/darunavir. However, these drugs still have not been approved by FDA for S-CoV-2 patients due to several limitations such as poor pharmacokinetic, low effectivity, and certain adverse effects. On the other side, repurposing of various classes of natural products in several herbs for S-CoV-2 Mpro inhibitor also showed promising results. These repurposings were mainly conducted through in silico approaches such as terpenoids from Cacospongia mycofijiensis, diterpenoids and bioflavonoids from Torreya nucifera, glycosylated flavonoids from various herbs, and saponins and tannins from various herbs.
Centella asiatica is a traditional herb from Southeast Asia that is widely used as medicines, cosmetics, foods, and beverages. This herb is known with different local names in different countries such as tapak kuda (Indonesia), pegaga (Indonesia and Malaysia), and gotu kola (Sri Lanka). C. asiatica contains various classes of active compounds (triterpenes, polyacetylenes, carotenoids, flavonoids, etc.) with different pharmacological activities including antioxidant, anti-inflammation, antidiabetic, and antimicrobial. C. asiatica also has antiviral activity against several viral strain such as human immunodeficiency virus (HIV), hepatitis B virus (HBV), and herpes simplex virus (HSV).,
Among these classes of active compounds, aglycone and glycoside forms of triterpenoids were considered as the most active ingredients in this herb, especially for therapeutic purposes. Several triterpenoids also have been identified in silico as active S-Cov-2 inhibitors in several reports such as asiatic acid derivatives and asiaticoside derivatives. In this report, we tried to identify potential triterpenoids, polyacetylenes, and benzoic ester derivatives from C. asiatica as S-CoV-2 Mpro inhibitor through molecular docking as part of repurposing of active ingredients in C. asiatica.
| Materials and Methods|| |
Materials used in this investigation were 11 known compounds from C. asiatica obtained from various literatures, and receptor protein S-CoV-2 Mpro (PDB ID:6LU7). The in silico molecular docking study was conducted using Autodock v4.2.6 and AutodockTools (http://autodock.scrips.edu/), ChemOffice Pro v15.00 PerkinElmer, Python Molecular Viewer (PMV 1.5.6), Open Babel graphical user interface (GUI), and Discovery Studio Visualizer.
Native ligand and protein receptor preparation
Protein structure of S-CoV-2 Mpro complexed with N3 as native ligand at PDB 6LU7 (2.16 Å resolution) was obtained from the Protein Data Bank via website https://rcsb.org. Native ligand from each Protein Data Bank was separated using PMV 1.5.6. Receptor and native ligand were separated in PDBQT (.pdbqt) format using AutodockTools and Open Babel GUI programs.
Test ligands preparation
The structure of test ligands from C. asiatica was obtained from some literatures [as shown in [Figure 1]]. Preparation of ligands structure was performed using ChemDraw® Pro v15 and Chem3D® Pro v15 then process, as described by Wang et al.
|Figure 1: Structure of test and native ligands (a) centellin, (b) asiaticin, (c) cenetellicin, (d) betulinic acid, (e) madasiatic acid, (f) isothankunic acid, (g) terminolic acid, (h) asiatic acid, (i) madecassic acid, (j) asiaticoside, (k) madecassoside, (l) N3, native ligand|
Click here to view
In silico molecular docking analysis
Autodock v4.2.6 (The Scripps Research Institute) was used for in silico molecular docking of 11 compounds from C. asiatica as S-CoV-2 Mpro inhibitor. Lamarckian parameters were used to perform the docking simulation. The grid box size and position were also validated to analyze the 11 known compounds from C. asiatica and visualized using Accelrys Discovery Studio Visualizer 4.0.
Pharmacokinetic properties analysis
Pharmacokinetic parameters of test and native ligands were predicted using SwissADME (http://swissadme.ch/).
| Results|| |
In silico molecular docking study
Inhibition activity of 11 natural products isolated C. asiatica against S-CoV-2 Mpro has been conducted through in silico molecular docking study. This investigation aimed to identify activity of S-CoV-2 Mpro inhibition based on the bonding form and interaction between the ligands and active site of receptor. The structure of S-CoV-2 Mpro consisted of 3 domains including domain I (residues 8–101), domain II (residues 102–184), and domain III (residues 185–200). Residue histidine 41 (His41) located in domain I and cysteine 145 (Cys145) located in domain II were the active sites of S-CoV-2 Mpro.
Docking results of ligands-receptor were evaluated by verifying re-docking result of native ligand (N3) and S-CoV-2 Mpro. Re-docking verification was performed based on the RSMD value of lowest interaction energy between backbone atoms and active site of S-CoV-2 Mpro of 1.112 Å (with Gibbs free energy of ‒6.91 kcal/mol at 86 clusters) and 0.978 Å (with Gibbs free energy of ‒9.18 kcal/mol at 36 clusters) for 100 runs with 2 Å RMSD tolerance, as described in [Figure 2].
|Figure 2: Native ligand (blue) and native ligand re-docked (red) comparisons|
Click here to view
According to the interactions shown in [Table 1], the interaction pattern of native ligand (N3)-amino acids located inside the active site after the re-docking was still similar to the interaction before the re-docking in which the native ligand (N3) had closely interacted with His41 in domain I dan Cys145 in domain II. Native ligand (N3) was formed pi-alkyl interaction with His41, covalent bonding with Cys145, and also interaction with other residues. Molecular docking results of 11 compounds obtained from C. asiatica indicated that asiatate 6 and asiatate 10 [Figure 3] have strong interaction with quite similar binding free energy to native ligand (‒9.00 and ‒9.58 kcal/mol compared to ‒9.18 kcal/mol, respectively). Asiatate 11 even had stronger interaction with S-CoV-2 Mpro active site (binding free energy of ‒10.98 kcal/mol); however, it did not interact with His41 residue in domain I as one of the residues of catalytic dyad.
|Figure 3: Three dimension of S-CoV-2 Mpro with (a) native ligand and (b) asiatate 10|
Click here to view
Pharmacokinetics properties analysis
Pharmacokinetics properties of native and test ligands were analyzed using free web server SwissADME. Molecular weight of 11 compounds showed a wide range of values from 250 to 1000 g/mol. In this range, 9 compounds (asiatate 1–9) have MW below 505 g/mol which are easy to be transported, diffused, and absorbed in the human body, while 2 compounds (asiatate 10–11) have about 950 g/mol which are hard to be distributed inside the body. These also could be figured out from high topological polar surface area (>140 Å2), low bioavailability score (0.17 for small molecules that fail to pass rule of five), low gastrointestinal absorption, and no blood–brain barrier (BBB) permeation values that indicated that asiatate 10 and asiatate 11 have poor pharmacokinetic properties. The remaining molecules (asiatate 1–asiatate 9) showed better pharmacokinetic parameters except for asiatate 4 with high lipophilicity (log P >5) which would cause rapid metabolic turnover, low solubility, and poor absorption.
| Discussion|| |
C. asiatica is a traditional herb containing various active compounds such as terpenoids and polyacetylenes. with antiviral activity against certain viruses such as HIV, HSV, and HBV., In this report, it was revealed that several active compounds might have anti-S-CoV-2 activity by inhibiting S-CoV-2 Mpro.
In general, all of the test ligands showed high affinity toward S-CoV-2 Mpro with binding free energy range from ‒7.52 to ‒10.98 kcal/mol; however, only asiatate 6 (EBind = ‒9.00 kcal/mol) and asiatate 10 (EBind = ‒9.58 kcal/mol) showed the highest affinity with proper interaction with His41 and Cys145 as catalytic dyad residues. Although both asiatate 6 and asiatate 10 were triterpenoid derivatives, they interacted with amino acid residues in the active site of S-CoV-2 Mpro with different parts of molecules. Asiatate 6 (a nonglycosylated triterpenoid) entered the active site and interacted with the residue through its pentacyclic triterpenoid backbone, while asiatate 10 (a glycosylated triterpenoid) entered the active site through glycosides moiety. This phenomenon was probably due to the bigger size of the asiatate 6 and asiatate 10 compared to the size of active site. Higher affinity of asiatate 10 and asiatate 11 toward S-CoV-2 Mpro was also probably due to the formation of many hydrogen bonds between amino acids residues and glycoside moiety.
From pharmacokinetic viewpoint, asiatate 6 showed better pharmacokinetic properties than asiatate 10 in which the later molecule violated 3 of 5 Lipinski rules, while the former one only violated 1 of 5 Lipinski rules. Other poor pharmacokinetic parameters also showed for asiatate 10 such as lower gastrointestinal absorption and smaller bioavailability score [Table 2]. This high affinity-poor pharmacokinetic properties of asiatate 10 also can be observed from glycyrrhizin which have the same structural features as asiatate 10. From in silico molecular docking and pharmacokinetics analysis, it can be stated that asiatate 6 was the most recommended compound for further in vitro study.
|Table 2: Pharmacokinetics properties of native ligand (N3) and test ligands|
Click here to view
Both asiatate 6 and asiatate 10 were pentacyclic triterpenoid derivatives. This class of compound showed a high affinity toward S-CoV-2 Mpro in some reports including several test ligands used in these investigations such as asiatic acid (asiatate 8) and asiaticoside derivatives (asiatate 10)., Other triterpenoid derivatives either from semisynthetic or natural also showed strong interaction including withaferin, acetylated asiatic acid, soyasaponin I, glycyrrhizin, and quinone-methide triterpenes. From these reports, it was revealed that triterpenoid derivatives need special attention for S-CoV-2 Mpro inhibitor.
| Conclusions|| |
Of 11 natural products isolated from C. asiatica, asiatate 6 indicated high potency as S-CoV-2 Mpro inhibitor. Asiatate 6 has quite similar binding free energy with S-CoV-2 Mpro as compared to native ligand (N3) with proper interaction to the catalytic dyad residues. Pharmacokinetic parameters' analysis also revealed that asiatate 6 has good pharmacokinetic performance. From these results, asiatate 6 could be recommended for further evaluation as S-CoV-2 Mpro inhibitor.
Financial support and sponsorship
This research was supported by Ministry of Research and Technology/National Research and Innovation Agency in funding aspect through Inter-Higher Education Collaborative Research grant with contract number 183/SP2H/AMD/LT/DRPM/2020.
Conflicts of interest
There are no conflicts of interest.
| References|| |
Gil C, Ginex T, Maestro I, Nozal V, Barrado-Gil L, Cuesta-Geijo MÁ, et al
. COVID-19: Drug targets and potential treatments. J Med Chem 2020;63:12359-86.
Hilgenfeld R. From SARS to MERS: Crystallographic studies on coronaviral proteases enable antiviral drug design. FEBS J 2014;281:4085-96.
Cao B, Wang Y, Wen D, Liu W, Wang J, Fan G, et al
. A trial of lopinavir-ritonavir in adults hospitalized with severe Covid-19. N Engl J Med 2020;382:1787-99.
Chen J, Xia L, Liu L, Xu Q, Ling Y, Huang D, et al
. Antiviral activity and safety of darunavir/cobicistat for the treatment of COVID-19. Open Forum Infect Dis 2020;7:ofaa241.
Osborne V, Davies M, Lane S, Evans A, Denyer J, Dhanda S, et al
. Lopinavir-ritonavir in the treatment of COVID-19: A dynamic systematic benefit-risk assessment. Drug Saf 2020;43:809-21.
Sepay N, Sekar A, Halder UC, Alarifi A, Afzal M. Anti-COVID-19 terpenoid from marine sources: A docking, admet and molecular dynamics study. J Mol Struct 2020;1228:129433.
Ghosh R, Chakraborty A, Biswas A, Chowdhuri S. Computer aided identification of potential SARS CoV-2 main protease inhibitors from diterpenoids and biflavonoids of Torreya nucifera leaves. J Biomol Struct Dyn 2020:1-6.
Falade VA, Adelusi TI, Adedotun IO, Abdul-Hammed M, Lawal TA, Agboluaje SA. In silico investigation of saponins and tannins as potential inhibitors of SARS-CoV-2 main protease (Mpro). In Silico Pharmacol. 2021;9:9. doi: 10.1007/s40203-020-00071-w.
Sardrood SG, Saadatmand S, Assareh MH, Satan TN. Chemical composition and biological activity of essential oils of Centella asiatica
(L.). Toxicol Environ Health Sci 2019;11:125-31.
Thanigaivel S, Durgadevi H, Balasubramaniam J, Mythily V, Elanchezhiyan M. Comparative evaluation of the anti-Hepatitis B virus activity of Centella asiatica
and Camellia sinensis
(green tea). BMC Infect Dis 2014;14:P21.
Lamorde M, Tabuti JR, Obua C, Kukunda-Byobona C, Lanyero H, Byakika-Kibwika P, et al
. Medicinal plants used by traditional medicine practitioners for the treatment of HIV/AIDS and related conditions in Uganda. J Ethnopharmacol 2010;130:43-53.
Sun B, Wu L, Wu Y, Zhang C, Qin L, Hayashi M, et al
. Therapeutic potential of Centella asiatica
and its triterpenes: A review. Front Pharmacol 2020;11:1-24.
Azerad R. Chemical structures, production and enzymatic transformations of sapogenins and saponins from Centella asiatica
(L.) Urban. Fitoterapia 2016;114:168-87.
Siddiqui BS, Aslam H, Ali ST, Khan S, Begum S. Chemical constituents of Centella asiatica
. J Asian Nat Prod Res 2007;9:407-14.
Morris GM, Huey R, Lindstrom W, Sanner MF, Belew RK, Goodsell DS, et al
. AutoDock4 and AutodockTools4: Automated docking with selectivity receptor flexibility. J Comput Chem 2009;30:2785-91.
Boyle NM, Banck M, James CA, Morley C, Vandermeersch T, Hutchison GR. Open babel : An open chemical toolbox. J Cheminform 2011;3:1-4.
Wang L, Wu Y, Deng Y, Kim B, Pierce L, Krilov G, et al
. Accurate and reliable prediction of relative ligand binding potency in prospective drug discovery by way of a modern free-energy calculation protocol and force field. J Am Chem Soc 2015;137:2695-703.
Daina A, Michielin O, Zoete V. SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci Rep 2017;7:1-3.
Sk MF, Roy R, Jonniya NA, Poddar S, Kar P. Elucidating biophysical basis of binding of inhibitors to SARS-CoV-2 main protease by using molecular dynamics simulations and free energy calculations. J Biomol Struct Dyn. 2020:1-13.
Matsson P, Kihlberg J. How big is too big for cell permeability? J Med Chem 2017;60:1662-4.
Martin YC. A bioavailability score. J Med Chem 2005;48:3164-70.
Arnott JA, Planey SL. The influence of lipophilicity in drug discovery and design. Expert Opin Drug Discov 2012;7:863-75.
Chen D, Oezguen N, Urvil P, Ferguson C, Dann SM, Savidge TC. Regulation of protein-ligand binding affinity by hydrogen bond pairing. Sci Adv 2016;2:e1501240.
Narkhede RR, Pise AV, Cheke RS, Shinde SD. Recognition of natural products as potential inhibitors of COVID-19 main protease (Mpro): In silico
evidences. Nat Products Bioprospect 2020;10:297-306.
Musfiroh I, Azura AR, Rahayu D. Prediction of Asiatic acid derivatives affinity against SARS-CoV-2 main protease using molecular docking. Pharm Sci Res 2020;7:57-64.
Ryu YB, Park SJ, Kim YM, Lee JY, Seo WD, Chang JS, et al
. SARS-CoV 3CLpro
inhibitory effects of quinone-methide triterpenes from Tripterygium regelii
. Bioorganic Med Chem Lett 2010;20:1873-6.
[Figure 1], [Figure 2], [Figure 3]
[Table 1], [Table 2]