Kelor Docing
Kelor Docing
JURNAL
BIOTEKNOLOGI & BIOSAINS INDONESIA
ABSTRACT
Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) disease (COVID-19) is a threat
to human health. This infection is determined by the interaction of the spike S1 domain protein
with angiotensin-converting enzyme 2 (ACE2) in the epithelial cells of the respiratory tract,
especially the lungs. ACE2 inhibition is an important target in controlling COVID-19. Therefore,
this study aims to explore the ability of apiin, epicatechin, and hesperetin from Moringa oleifera
in interacting with the ACE2 using MOE 2008.10. The ligand molecules were prepared from
PubChem database. The ACE2 protein was retrieved from Protein Data Bank (ID 1R4L) and
analyzed for the active sites. Analysis of docking scores and hydrogen bonds of ACE2-ligand
complex and active site showed that the affinity of flavonoids can be ranked as hesperetin >
epicatechin > apiin > C19H23Cl2N3O4. The results provided computational information that apiin,
epicatechin, and hesperetin have the potential to prevent COVID-19 infection. The prediction
of activity spectra for substances (PASS) score showed the ligand displays antiviral activity.
ABSTRAK
Infeksi severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) pada pandemi
coronavirus disease 2019 (COVID-19) menjadi ancaman dunia kesehatan saat ini. Infeksi SARS-
CoV-2 ditentukan oleh interaksi protein spike envelope S1 domain dengan reseptor angiotensin-
converting enzyme 2 (ACE2) yang diekspresikan pada sel epitel saluran pernafasan terutama
paru-paru. Mekanisme penghambatan ACE2 menjadi target penting dalam pengendalian COVID-
19. Penelitian ini bertujuan mengeksplorasi kemampuan senyawa apiin, epicatechin, dan
hesperetin dari Moringa oleifera dalam berinteraksi dengan sisi aktif ACE2 menggunakan metode
penambatan molekul. Studi dilakukan dengan preparasi struktur molekul ligan dari PubChem
database dan diolah dengan MOE 2008.10. Selanjutnya, data protein ACE2 (Protein Data Bank
ID 1R4L) dianalisis sisi aktifnya untuk mengetahui lokasi penambatan ligan senyawa. Analisis skor
docking dan ikatan hydrogen komplek ligan dan sisi aktif ACE2 menunjukkan bahwa afinitas
flavonoid dapat diperingkatkan sebagai afinitas hesperetin > epicatechin > apiin > C19H23Cl2N3O4.
Ketiga ligan senyawa yang terkandung dalam M. oleifera secara in silico mampu mengikat sisi
aktif ACE2, sehingga berpotensi mencegah infeksi COVID-19. Skor PASS (prediction of activity
spectra for substances) menunjukkan aktivitas biologis ligan yang menyerupai antiviral.
259
Potensi Senyawa Bioaktif Tanaman Kelor... Koentjoro et al.
260
J Bioteknol Biosains Indones – Vol 7 No 2 Thn 2020
Moringa oleifera atau yang lebih dikenal pada penelitian ini adalah YASARA version
sebagai kelor merupakan tanaman tropis yang 20.4.24 (de Groot et al. 1997), Molecular
dibudidayakan dan menjadi salah satu sumber Operating Environment (MOE) MOE 2008.10
nutrisi masyarakat Indonesia (Riastiwi et al. (Vilar et al. 2008), dan PyMOL versi 2.3.4
2018). Pada penelitian sebelumnya, M. oleifera (Schrödinger, LLC, USA, http://www.pymol.org/)
memiliki aktivitas antioksidan anti-aging
(Imamsari et al. 2018). Kandungan tanaman Evaluasi aktivitas biologis dan preparasi ligan
tersebut diketahui terutama flavonoid mampu Senyawa bioaktif yang dipilih dalam
mendenaturasi dan mengganggu fungsi docking ini adalah apiin, epicatechin, dan
banyak makromolekul termasuk ACE (Khan et hesperetin. Struktur 2D (2 Dimensi) dan 3D (3
al. 2019). Oleh karena itu, tujuan penelitian ini Dimensi) dari apiin, epicatechin, dan
adalah untuk memprediksi apakah senyawa hesperetin diperoleh dari PubChem database
bioaktif flavonoid dari M. oleifera antara lain (https://pubchem.ncbi.nlm.nih.gov/), yang
apiin, epicatechin, dan hesperetin mampu disajikan dalam Tabel 1. Langkah pertama
berinteraksi dengan reseptor ACE2, sehingga yang dilakukan untuk mengetahui apakah
berpeluang sebagai inhibitor interaksi SARS- senyawa yang dipilih memiliki spektrum
CoV-2 dan sel inang. Metode yang digunakan aktivitas biologis adalah analisis PASS
adalah melalui molecular docking. Hasil studi ini (prediction of activity spectra for substances)
diharapkan mampu memberikan pemodelan online (Lagunin et al. 2000). Software ini
molekuler dan pendekatan berbasis informatika memprediksi hampir 900 jenis aktivitas
dalam menjelaskan senyawa-senyawa bioaktif berdasarkan struktur suatu senyawa. Masing-
dari M. oleifera yang berpotensi sebagai masing struktur kemudian diubah menjadi
kandidat pencegahan infeksi SARS-CoV-2 bentuk file SMILES menggunakan software
secara alami. ChemDraw Professional 16.0 dan dianalisis
aktivitas biologisnya menggunakan software
BAHAN DAN METODE PASS online (http://www.way2drug.com/passonline/).
Untuk persiapan docking, maka
Tempat dan waktu penelitian masing-masing struktur selanjutnya
Penelitian ini dilakukan di Laboratorium dipreparasi menggunakan software MOE
Biologi Molekuler, Fakultas Kesehatan, 2008.10. Persiapan meliputi optimasi
Universitas Nahdlatul Ulama Surabaya pada geometri ligand dan pembuangan molekul air
bulan April sampai dengan Mei 2020. (H2O). Molekul ligan selanjutnya disimpan
dalam format sdf.
Bahan
Hardware yang digunakan adalah Persiapan struktur protein
Laptop Intel® Core™ i5–3340 M CPU Struktur protein ACE2 diperoleh dari
@2.700GhZ, 2.70 Ghz dengan sistem PDB (Protein Data Bank) pada website
operasi Microsoft Windows 10 pro dan 64 bits https://www.rcsb.org/, dengan PDB ID 1R4L
operating systems. Software yang digunakan (Tabel 2.) yang berasal dari organisme H.
Keterangan: Pa menunjukkan probalititas aktivitas senyawa, Pi mewakili probabilitas senyawa tidak aktif
261
Potensi Senyawa Bioaktif Tanaman Kelor... Koentjoro et al.
Tabel 2. Data PDB makromolekul protein yang penempatan (triangle matcher) diulang
digunakan dengan pembacaan energi tiap posisi dan
parameter lain yang disesuaikan dengan
Spesifikasi MOE 2008.10. Hasil yang ditampilkan
PDB ID 1R4L adalah yang terbaik dari 100 pengulangan.
Organisme asal H. sapiens
Analisis docking
Metode analisis struktur X-ray crystallography Analisis docking dilakukan
Resolusi 3.00 Å menggunakan beberapa parameter dari
interaksi ligan-reseptor, yaitu energi ikatan
sapiens (Towler et al. 2004). Molekul ACE2 dan konstansta inhibisi (S) terendah, ikatan
disiapkan dengan menggunakan software hydrogen. Ikatan yang terjadi terutama
YASARA versi 20.4.24. Persiapan yang pada sisi aktif yaitu pada asam amino yang
dilakukan adalah pemisahan gugus H2O, sesuai akan dijadikan prediksi sementara
pemisahan ligan alami ((S,S)-2-{1-carboxy-2- serta kemampuan senyawa steroid dan
[3-(3,5-Dichloro-Benzyl)-3H-imidazol-4-YL]- flavonoid yang tertambat pada molekul
ethylamino}-4-methyl-pentanoic acid, N- protein sasaran yang terbesar akan
acetyl-D-glucosamine, ion zinc, dan klorida). dijadikan sebagai dasar penentuan
Struktur protein dan ligan yang telah senyawa dengan aktivitas terbaik (Mukund
dipisahkan selanjutnya disimpan dalam satu et al. 2019). Hasil docking kemudian di-
folder dengan format pdb. superposisikan dengan data asli (PDB ID
1R4L) dan divisualisasikan dengan program
Proses molecular docking PyMOL versi 2.3.4.
Molecular docking (penambatan
molekul) dilakukan menggunakan software HASIL DAN PEMBAHASAN
MOE 2008.10. Proses pertama dari
molecular docking adalah optimasi Interaksi berperan penting dalam
geometri struktur ligan yang akan diuji untuk perancangan obat berbasis struktural.
menghasilkan energi molekul terendah. Dalam penelitian ini, senyawa bioaktif
Senyawa ligan diperbaiki dengan dipilih berasal dari data sekunder Lin et al.
penambahan atom hidrogen dan muatan (2019) yang melakukan analisis M. oleifera
parsial gasteiger charges (Sliwoski et al. menggunakan ultra-performance liquid
2013). Masing-masing senyawa ligan chromatography quadrupole time-of-flight
selanjutnya disimpan dalam format mdb. mass spectrometry (UPLC-QTOF-MS).
Penentuan sisi aktif dari reseptor ACE2 Analisis PASS (Tabel 1) memprediksi apiiin
selanjutnya dilakukan untuk mengetahui memiliki aktivitas antivirus Hepatitis B
lokasi dan menghitung kemungkinan dengan nilai Pa (probablitias untuk aktif)
interaksi sisi aktif dengan ligan senyawa 0,423. Epicatehin memiliki nilai Pa 0,692
bioaktif dari koordinat atom (Rahman et al. sebagai antivirus influenza. Sedangkan
2019). 1R4L yang diunduh dari PDB adalah hesperetin memiliki prediksi dengan nilai Pa
dalam bentuk berikatan dengan ligan dan 0,673 untuk antiviral influenza.
molekul air, sehingga pada tahap persiapan Tabel 3 menyajikan hasil docking ligan
reseptor kedua molekul tersebut senyawa apiin, epicatechin, dan hesperetin
dihilangkan. Molekul air pada analisis ini terhadap struktur protein kristal ACE2. Pada
dihilangkan karena pada umumnya tidak analisis ini, molekul target protein ACE2 (PDB
terlibat pada pengikatan protein dan ligan, 1R4L) dipilih karena merupakan hasil dari
selain itu agar perhitungan docking lebih analisis X-ray kristalografi. Sebagaimana
mudah karena molekul air dapat telah diketahui, X-ray kristalografi merupakan
menyebabkan distorsi pada pencarian sisi teknik visualisasi struktur molekul pada
aktif (Parikh dan Kellogg 2014). resolusi atom dan menunjukkan struktur
Reseptor ACE2 dan ligan yang telah mendekati yang ada di dalam sel (Aitipamula
dipersiapkan selanjutnya disimulasikan dan Vangala 2017). 1R4L yang dipilih ini
dengan menu compute-simulation-dock berasal dari spesies H. sapiens agar hasil
pada software MOE 2008.10. Hasil docking docking yang dilakukan mendekati kondisi
disimpan dalam format mdb. Metode yang sebenarnya.
262
J Bioteknol Biosains Indones – Vol 7 No 2 Thn 2020
Tabel 3. Hasil docking ligan senyawa terhadap struktur protein kristal ACE2
263
Potensi Senyawa Bioaktif Tanaman Kelor... Koentjoro et al.
Ikatan antara residu asam amino pada apigenin, luteolin, dan quercetin dari M.
sisi aktif dan ligan diharapkan mampu oleifera memiliki aktivitas mengikat ACE2 dan
menghilangkan fungsi enzimatik dari ACE2, berpotensi dalam menghambat infeksi SARS-
sehingga menurunkan interaksi SARS-CoV- CoV-2. Oo et al. (2016) melaporkan bahwa
2. Analisis in silico yang dilakukan Laksmiani hesperetin dan apiin memiliki aktivitas
et al. (2020) menunjukkan bahwa senyawa penghambatan infeksi Virus Chikungunya.
Tabel 4. Ikatan hidrogen antara ligan dan residu asam amino, serta panjang interaksinya
Ikatan hidrogen terbentuk jika jarak antara hidrogen dengan atom elektronegatif berada pada rentang 2.5–3.5 Å
264
J Bioteknol Biosains Indones – Vol 7 No 2 Thn 2020
Gambar 2. Konformasi–docking ACE2 dan ligan senyawa apiin. (A) Animasi 2D ligan (warna kuning) dalam sisi aktif
1R4L; (B) Interaksi ligan dengan residu asam amino dari ACE2. Komplek dihasilkan dari interaksi ligan
MOE 2008.10
Gambar 3. Konformasi–docking ACE2 dan ligan senyawa epicatechin. (A) Animasi 2D ligan (warna kuning) dalam sisi
aktif 1R4L; (B) Interaksi ligan dengan residu asam amino dari ACE2. Komplek dihasilkan dari interaksi ligan
MOE 2008.10
265
Potensi Senyawa Bioaktif Tanaman Kelor... Koentjoro et al.
Gambar 4. Konformasi–docking ACE2 dan ligan senyawa hesperetin. (A) Animasi 2D ligan (warna kuning) dalam sisi
aktif 1R4L; (B) Interaksi ligan dengan residu asam amino dari ACE2. Komplek dihasilkan dari interaksi ligan
MOE 2008.10
Gambar 5. Representasi simulasi docking dari komplek ACE2 dengan ligan. Ketiga ligan ditampilkan dalam mode stik
berada di dalam sisi aktif dari 1R4L (warna silver), yaitu 5-[3-[4-(3,5-dichloro-2-methoxyphenyl)piperidin-1-
yl]-3-oxopropyl]-5-methylimidazolidine-2,4-dione (warna hijau), apiin (warna biru), epicatechin (warna
kuning), dan hesperetin (warna merah). Komplek dihasilkan dari PyMOL versi 2.3.4
266
J Bioteknol Biosains Indones – Vol 7 No 2 Thn 2020
267
Potensi Senyawa Bioaktif Tanaman Kelor... Koentjoro et al.
268
J Bioteknol Biosains Indones – Vol 7 No 2 Thn 2020
simple: Including explicit water with AIDS. AIDS 34: 1575–1576. doi:
molecules in protein-protein docking 10.1097/QAD.0000000000002553
simulations improves model quality. Sun J, Chen YT, Fan XD, Wang XY, Han QY,
Proteins 82: 916-932. doi: Liu ZW (2020) Advances in the use of
10.1002/prot.24466 chloroquine and hydroxychloroquine for
Rahman N, Muhammad I, Gul-E-Nayab HK, the treatment of COVID-19. Postgrad
Aschner M, Filosa R, Daglia M (2019) Med 132: 604-613. doi:
Molecular docking of isolated alkaloids 10.1080/00325481.2020.1778982
for possible α-glucosidase inhibition. Towler P, Staker B, Prasad SG, Menon S,
Biomolecules 9: 544. doi: Tang J, Parsons T, Ryan D, Fisher M,
10.3390/biom9100544 Williams D, Dales NA, Patane MA,
Riastiwi I, Damayanto IPGP, Ridwan R, Pantoliano MW (2004) ACE2 X-ray
Handayani T, Leksonowati A (2018) structures reveal a large hinge-bending
Moringa oleifera distribution in Java and motion important for inhibitor binding
Lesser Sunda islands attributed with and catalysis. J Biol Chem 279: 17996-
annual rainfall. Biosaintifika 10: 613-621. 18007. doi: 10.1074/jbc.M311191200
doi: 10.15294/biosaintifika.v10i3.16115 Vergara-Jimenez M, Almatrafi MM,
Saha P, Banerjee AK, Tripathi PP, Srivastava Fernandez ML (2017) Bioactive
AK, Ray U (2020) A virus that has gone components in Moringa oleifera leaves
viral: Amino acid mutation in S protein protect against chronic disease.
of Indian isolate of coronavirus COVID- Antioxidants (Basel) 6: 91. doi:
19 might impact receptor binding, and 10.3390/antiox6040091
thus, infectivity. Biosci Rep 40: Vilar S, Cozza G, Moro S (2008) Medicinal
BSR20201312. doi: chemistry and the molecular operating
10.1042/bsr20201312 environment (MOE): Application of
Schrezenmeier E, Dorner T (2020) QSAR and molecular docking to drug
Mechanisms of action of discovery. Curr Top Med Chem 8:
hydroxychloroquine and chloroquine: 1555-1572. doi:
Implication for rheumatology. Nat Rev 10.2174/156802608786786624
Rheumatol 16: 155-156. doi: Wang HD, Du GH, Liu A (2009) Evaluation for
10.1038/s41584-020-0372-x inhibitory effects of natural flavonoids
Sharifi N, Souri E, Ziai AS, Amin G, Amanlou on neuraminidases. Chinese J New
M (2013) Discovery of new angiotensin Drugs 15: 1435-1439
converting enzyme (ACE) inhibitors WHO (2020) Report of the WHO-China Joint
from medicinal plants to treat Mission on Coronavirus Disease 2019
hypertension using an in vitro assay. (COVID-19). World Health
Daru 21: 74. doi: 10.1186/2008-2231- Organization, Geneva Switzerland
21-74 Wiese O, Zemlin AE, Pillay TS (2020)
Sliwoski G, Kothiwale S, Meiler J, Lowe Jr EW Molecules in pathogenesis: Angiotensin
(2013) Computational methods in drug converting enzyme 2 (ACE2). J Clin
discovery. Pharmacol Rev 66: 334‐395. Pathol 2020: 1-6. doi:10.1136/jclinpath-
doi: 10.1124/pr.112.007336 2020-206954
Su L, Ma X, Yu H, Zhang Z, Bian P, Han Y, Xia S, Liu M, Wang C, Xu W, Lan Q, Feng S,
Sun J, Liu Y, Yang C, Geng J, Zhang Z, Qi F, Bao L, Du L, Liu S, Qin C, Sun F,
Gai Z (2020a) The different clinical Shi Z, Zhu Y, Jiang S, Lu L (2020)
characteristics of corona virus disease Inhibition of SARS-CoV-2 (previously
cases between children and their 2019-nCoV) infection by a highly potent
families in China – the character of pan-coronavirus fusion inhibitor
children with COVID-19. Emerg targeting its spike protein that harbors a
Microbes Infect 9: 707-713. doi: high capacity to mediate membrane
10.1080/22221751.2020.1744483 fusion. Cell Res 30: 343-355. doi:
Su J, Shen X, Ni Q, Zhao H, Cai J, Zhu B, Wu 10.1038/s41422-020-0305-x
W, Lang G, Xu K, Sheng J (2020b) Xu J, Zhao S, Teng T, Abdalla AE, Zhu W, Xie
Infection of severe acute respiratory L, Wang Y, Guo X (2020) Systmetic
syndrome coronavirus 2 in a patient comparison of two animal-to-human
269
Potensi Senyawa Bioaktif Tanaman Kelor... Koentjoro et al.
270