key: cord-0870335-el3ilypg
authors: Umar, Abd. Kakhar
title: Flavonoid compounds of buah merah (Pandanus conoideus Lamk) as a potent SARS-CoV-2 main protease inhibitor: in silico approach
date: 2021-08-10
journal: Futur J Pharm Sci
DOI: 10.1186/s43094-021-00309-0
sha: c58d2ea68f2591291c8e43e8d68fb88aa05dce6c
doc_id: 870335
cord_uid: el3ilypg

BACKGROUND: COVID19 is a global pandemic that threatens all nations. As there is no effective antiviral drug for COVID19, we examined the potency of natural ingredients against the SARS-CoV-2 main protease (PDB ID 6YNQ). Buah merah is a typical fruit from Papua, Indonesia, which is known to contain high levels of carotenoids and flavonoids. The contents have been proven to be effective as antiparasitic and anti-HIV. An in silico approach to 16 metabolites of buah merah (Pandanus conoideus Lamk) was carried out using AutoDock Vina. Furthermore, the study of the dynamics of ligand–protein interactions was carried out using CABS Flex 2.0 server to determine the test ligand and receptor complexes' stability. ADMET prediction was also carried out to study the pharmacokinetic profile of potential antiviral candidates. RESULT: The docking results showed that 3 of the 16 buah merah metabolites were potent inhibitors against the SARS-CoV-2 main protease. The flavonoid compounds are quercetin 3′-glucoside, quercetin 3-O-glucose, and taxifolin 3-O-α-arabinopyranose with a binding affinity of − 9.7, − 9.3, and − 8.8, respectively, with stable ligand–protein complex. ADMET study shows that the three compounds are easily dissolved, easily absorbed orally and topically, have a high unbound fraction, low toxicity, and non-irritant. CONCLUSION: We conclude that quercetin 3′-glucoside, quercetin 3-O-glucose, and taxifolin 3-O-α-arabinopyranose can be used and improved as potential anti-SARS-CoV-2 agents in further study.

The 2019 coronavirus disease (COVID19) has had a significant impact on all countries in the world. Based on a WHO report accessed on Apr 12, 2021, there were 4 million new cases in 1 week, with an increase of 11% compared to last week, with over 71.000 new deaths reported [1] . The virus that causes COVID19 is severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which belongs to the β-coronavirus group. Inhibiting virus replication is one method of reducing the severity of infection. The major proteins known to be responsible for the life cycle of viruses have been reported. This protein plays a role in the maturation process of viral proteins to become a target in developing new antivirals [2] [3] [4] .

Several secondary metabolites from plants have been reported to have good inhibitory activity against the SARS-CoV-2 main protease. The terpenoid and triterpenoid compounds are the strongest inhibitor [2, 3] . Previous studies reported that rutin had the best inhibitory effect, even compared to the FDA-approved COVID19 antiviral (Remdesivir) [2, 5] . This terpenoid is known to bind stably to the SARS-CoV-2 main protease receptor (PDB ID 6YNQ). Since the 6YNQ protein has never been identified to have mutations, it has been an appealing focus in developing antivirals using the in silico approach [2] .

Buah merah is a typical Indonesian plant that has been widely used by Indonesian people, especially Papua, as medicine or daily food. Buah merah contains high levels of carotenoid and flavonoid metabolites. Buah merah (Pandanus conoideus Lamk) has been reported to have antioxidant, antitumor, immunomodulatory, antiparasitic, and anti-HIV effects [6] [7] [8] [9] . As buah merah has good antiparasitic and anti-HIV potential, suggesting potential activity against the SARS-CoV-2 main protease [3] . Therefore, this study was conducted to see whether the flavonoid and carotenoid content of buah merah inhibited SARS-CoV-2 main protease compared to rutin and remdesivir. Molecular dynamics were studied to see how the relationship and stability of the ligand-protein complexes. ADMET prediction was also performed to assess the pharmacokinetic profile of potential drug candidates.

Pandanus conoideus Lamk was reported to contain carotenoids and flavonoids. (The list can be seen in Table 1 .) [10, 11] The ligand structure was drawn using ChemDraw Pro 12.0. The ligand structure was then trimmed using ChemDraw's "clean structure" feature, and their energy was minimized (MM2) using Chem3D.

The ligand structure was then saved into PDB format. The ligands were optimized again using AutoDock-Tools 1.5.6 (ADT) (TheScripps Research Institute, the USA) to add Gasteiger charges, set rotatable bonds, and 

The docking process was carried out using AutoDock Vina. The operating system used was Windows 10 Home Single Language 64 bit with AMD Ryzen 5 3500U, Radeon Vega Mobile Gfx 2.10 GHz, and RAM of 8 GB. The energy range was set to 4 and exhaustiveness to 8. The output file was made in PDBQT format used for visualization of docking results. 2D and 3D visualization was done using DS.

The ligand-protein interaction dynamics study was carried out to determine the most active amino acid residues at the binding site of the SARS-CoV-2 main protease. The output file from the docking process produces 9 ligandprotein interaction models for each ligand. All the active amino acid residues bind to the ligand, and their number of occurrences has been observed and recorded. The protein's stable structure was studied using the CABS Flex 2.0 server, which is based on coarse-grained simulations of protein motion [12] . The number of cycles and trajectory frames was set to 50, with a global weight of 1.0 and a temperature of 1.4. The distance restraints generator was set to default values. This test aims to see whether the ligand-protein interaction remains stable during attachment [2] .

The pharmacokinetics profile of the selected potential ligands was studied using pkCSM ADMET to determine the sterol compound's quality and safety. The SMILES string for each ligand is obtained from a PDB ligand file converted to SMI format using DS.

The total active amino acid residues at the binding site of the SARS-CoV-2 main protease were 23, including Asn142, Arg188, Cys44, Cys145, Gln189, Gln192, Glu166, Gly143, His43, His164, Leu141, Leu167, Met49, Met165, Phe140, Pro168, Ser46, Ser144, Thr24, Thr45, Thr45, and Thr190. The protein used (6YNQ) binding site appears to have approximately the same amino acid residue as 6LU7 [2, 3] . Most of the amino acid residues in the 6YNQ binding pocket are hydrogen donors or acceptors (see Fig. 1 ).

Based on the docking results, the flavonoid component of Pandanus conoideus Lamk had the strongest binding affinity, while the carotenoid component only ranged from − 6.8 to − 7.8. In Table 1 , it can be seen that quercetin 3′-glucoside has the highest binding affinity, which is − 9.7, then followed by quercetin 3-O-glucose (− 9.3) and taxifolin 3-O-α-arabinopyranose (− 8.8). This value was significant than the reference flavonoid and antiviral. Rutin only gets a binding affinity score of − 8.4, then followed by astragalin (− 8.2), trifolin (− 8.7), and remdesivir (− 7.5). Previously, rutin, astragalin, trifolin, and remdesivir have been investigated through an in silico approach as potential inhibitors of the SARS-CoV-2 main protease at the same protein code tested in this study (6YNQ) [2] . Compared to the rutin structure, quercetin 3′-glucoside has a structure that fits perfectly with the binding pocket of the SARS-CoV-2 main protease (see Fig. 1 ). Quercetin 3′-glucoside also binds to more amino acid residues than the reference ligands. Hydrogen bonds formed on quercetin 3-O-glucose are two times more than rutin with double hydrogen bonds on the Asn142 and Cys145 residues (see Fig. 2 ). Asn142 and Cys145 are known to be the catalytic active site residue of the SARS-CoV-2 main protease, so that forming bonds to these residues will produce strong inhibition [2, 12] . Apart from Asn142 and Cys145, the test ligands with the strongest binding affinity have similar interactions with several amino acid residues, including Gln189, Glu166, and His41. The reference ligands (astragalin and remdesivir) have a lower binding affinity than rutin. This result is consistent with the previous study [2] . It has also been shown that astragalin and remdesivir form unfavorable donor bonds in the residues of Thr190 and Glu166 (see Fig. 3 ).

The dynamics of ligand-protein interactions were studied to determine which amino acid residues interact with the ligands most frequently and rank the most active amino acid residues in all ligands. This test was only performed on a potent compound obtained from previous molecular docking studies, including quercetin 3′-glucoside (− 9.7), quercetin 3-O-glucose (− 9.3), dan taxifolin 3-O-α-arabinopyranose (− 8.8) as well as with reference flavonoid (rutin, − 8.4) and FDA approved antiviral (remdesivir, − 7.5). In Fig. 4 When interacting with SARS-CoV-2 main protease, all potent test ligands and reference ligands showed more than 90% amino acid residue yielding RMSD < 2 Å. The fluctuation of the root means square can be seen in Fig. 5 . From these results, the interaction of each test ligand and reference ligand with protein forms a stable complex. The stable structure of each of the ligand-protein complexes can be seen in Fig. 6 .

The ADMET prediction results showed that the three potent test ligands showed a good pharmacokinetic profile. ADMET prediction data can be seen in Table 2 .

The SARS-CoV-2 main protease does have many hydrogen donors and acceptors in its binding pocket. This can be seen in the interaction of quercetin 3′-glucoside with the receptor, where nine hydrogen bonds are formed, likewise, for other ligands where the hydrogen bond is dominant. This can be utilized for more optimal ligand development by targeting the hydrogen bonds in the amino acids Asn142 and Cys145, which are crucial amino acids [13] [14] [15] . Among the potent test ligands, quercetin 3′-glucoside most frequently interacts with Cys145 on all modes, acting as a catalytic active site residue. The ranking of amino acid residues' occurrence showed that Cys145, Glu166, Asn142, and His41 were the residues that played the most significant role in interacting with ligands.

Drugs can be classified based on their solubility. Drugs with a LogS value > − 2 show high solubility, the range − 2 to − 4 is slightly soluble, and < − 4 is insoluble. Based on the results of the ADMET prediction study, it can be seen that all potent ligands have good solubility. Value of HIA > 30% and LogKp < − 2.5 demonstrated that all potent ligands have good oral absorption and skin penetration. All potent ligands are also not included as substrates or inhibitors of P-glycoprotein I/II. This shows that P-glycoprotein does not assist the absorption of all potent ligands. All potent ligands also have no contraindication with other drugs whose absorption is assisted by P-glycoprotein [16] .

All potent ligands' distribution is also excellent where the log Vdss value is > − 0.15, and the free fraction in plasma is > 20%. The higher the logVdss value, the more drug fraction distributed to the tissue than in plasma. The more free fraction, the more efficient and the smaller the dose of drug needed. All the test ligands also showed low blood barrier penetration (logBB < − 1 and logPS < − 3), so that it can be said that the ligands would not directly affect the central nervous system. In terms of metabolism, all potent ligands are not substrates or inhibitors of cytochrome P450, so it can be said that all test ligands are not metabolized by cytochrome P450 and do not interfere with the metabolism of other drugs [16] . Quercetin 3′-glucoside, quercetin 3-O-glucose, and taxifolin 3-O-α-arabinopyranose have a total clearance of 0.437, 0.568, and − 0.007, respectively, and are not a substrate of OCT2.

The maximum human tolerable dose of quercetin 3′-glucoside, quercetin 3-O-glucose, and taxifolin 3-O-αarabinopyranose is 4.15, 5.67, and 8.57 mg/KgBB/day, respectively. All test ligands are not hERG I and II inhibitors and therefore do not potentially cause fatal ventricular arrhythmia. The oral rat acute and chronic toxicity of each potent ligand can be seen in Table 2 . All potent ligands are not hepatotoxic and non-irritant.

The terpenoid compounds of buah merah (Pandanus conoideus Lamk) have potent inhibitory activity against the SARS-CoV-2 main protease. Quercetin 3′-glucoside, quercetin 3-O-glucose, and taxifolin 3-O-α-arabinopyranose are potent inhibitors with a binding affinity − 9.7, − 9.3, and − 8.8, respectively. The three compounds that have an excellent pharmacokinetic profile are non-hepatotoxic and non-irritant. Based on this in silico study, we conclude that quercetin 3′-glucoside, quercetin 3-O-glucose, and taxifolin 3-O-α-arabinopyranose can be used and improved as potential anti-SARS-CoV-2 agents in further study.

Abbreviations COVID-19: Coronavirus disease 2019; SARS-CoV-2: Severe acute respiratory syndrome coronavirus-2; PDB: Protein data bank; ADMET: Absorption, distribution, metabolism, excretion, and toxicity; HIV: Human immunodeficiency virus; ADT: AutoDockTools 1.5.6; DS: Discovery studio 2021 client; RMSD: Root mean square distance; RMSF: Root mean square fluctuation.

World Health Institute (WHO) (2021) Weekly epidemiological update on COVID-19-6

Identification of potential flavonoid inhibitors of the SARS-CoV-2 main protease 6YNQ: a molecular docking study

Mushroom-derived bioactive compounds potentially serve as the inhibitors of SARS-CoV-2 main protease: an in silico approach

Molecular mechanisms of the novel coronavirus SARS-CoV-2 and potential anti-COVID19 pharmacological targets since the outbreak of the pandemic

Discovery of potential flavonoid inhibitors against COVID-19 3CL proteinase based on virtual screening strategy

The effect of pandanus conoideus lamk extract to the serum level of TNF-α, IL-10 and parasitemia of plasmodium berghei infected in mice

Tumor necrosis factor-alpha (TNF-Alpha) and intercellular adhesion molecule-1

7:158 expression of plasmodium berghei infected swiss mice treated with red fruit (Pandanus Conoideus Lam) ethanol extract

Immunomodulatory effect of orally red fruit (Pandanus conoideus) extract on the expression of CC chemokine receptor 5 mRNA in HIV patients with antiretroviral therapy

Carotenoid content of five accessions red fruit

Carotenoid composition in buah merah (Pandanus conoideus Lam.), an indigenous red fruit of the Papua Islands

Computational study of antioxidant activity and bioavailability of papua red fruit (Pandanus conoideus Lam.) flavonoids through docking toward human serum albumin

CABS-flex standalone: a simulation environment for fast modeling of protein flexibility

In silico screening of FDA approved drugs reveals ergotamine and dihydroergotamine as potential coronavirus main protease enzyme inhibitors

Remdesivir (GS-5734) as a therapeutic option of 2019-nCOV main protease-in silico approach

In-silico analysis of the inhibition of the SARS-CoV-2 main protease by some active compounds from selected African plants

pkCSM: predicting smallmolecule pharmacokinetic and toxicity properties using graph-based signatures

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

The authors are thankful to Chusnul Nur Ramadhani for her assistance in finding the buah merah content.

The study has only one author. All of the study and data were managed by the author. The author read and approved the final manuscript.

Not applicable.

All data and materials are available upon request.

Ethics approval and consent to participate Not applicable.

Not applicable.

The author declares no conflict of interest.Received: 20 April 2021 Accepted: 1 August 2021