key: cord-0944523-j0u9yzds authors: My, Tran Thi Ai; Loan, Huynh Thi Phuong; Hai, Nguyen Thi Thanh; Hieu, Le Trung; Hoa, Tran Thai; Thuy, Bui Thi Phuong; Quang, Duong Tuan; Triet, Nguyen Thanh; Anh, Tran Thi Van; Dieu, Nguyen Thi Xuan; Trung, Nguyen Tien; Hue, Nguyen Van; Tat, Pham Van; Tung, Vo Thanh; Nhung, Nguyen Thi Ai title: Evaluation of the Inhibitory Activities of COVID‐19 of Melaleuca cajuputi Oil Using Docking Simulation date: 2020-06-02 journal: ChemistrySelect DOI: 10.1002/slct.202000822 sha: d1c4040605c82ec7b32ad685dc41d18173251863 doc_id: 944523 cord_uid: j0u9yzds GC‐MS was applied to identify 24 main substances in Melaleuca cajuputi essential oil (TA) extracted from fresh cajeput leaves through steam distilling. The inhibitory capability of active compounds in the TA from Thua Thien Hue, Vietnam over the Angiotensin‐Converting Enzyme 2 (ACE2) protein in human body ‐ the host receptor for SARS‐CoV‐2 and the main protease (PDB6LU7) of the SARS‐CoV‐2 using docking simulation has been studied herein. The results indicate that the ACE2 and PDB6LU7 proteins were strongly inhibited by 10 out of 24 compounds accounting for 70.9% in the TA. The most powerful anticoronavirus activity is expressed in the order: Terpineol (TA2) ≈ Guaiol (TA5) ≈ Linalool (TA19) > Cineol (TA1) > β‐Selinenol (TA3) > α‐Eudesmol (TA4) > γ‐Eudesmol (TA7). Interestingly, the synergistic interactions of these 10 substances of the TA exhibit excellent inhibition into the ACE2 and PDB6LU7 proteins. The docking results orient that the natural Melaleuca cajuputi essential oil is considered as a valuable resource for preventing SARS‐CoV‐2 invasion into human body. The SARS-CoV-2 pandemic has been the hottest issue all over the world currently. The number of newly infected and died patients has been increasing day by day. [1] Finding solutions to treat and prevent SARS-CoV-2 is concerned by scientists in the world. In particular, medicine derived from natural and safe herbs is one of the research directions that deserve attention. Given the abundant supply in Vietnam and the specific healing properties of Melaleuca cajuputi, [2] we recommend a solution for prevention of SARS-CoV-2 by using cajeput essential oil. Melaleuca cajuputi (Figure 1 ) is one of hundreds species of Melaleuca genus (Myrtaceae family) that is commonly grown in Vietnam, Indonesia, and Thailand. [3] Melaleuca cajuputi can be used as fuelwood, piles and frame poles in construction. Its flowers attract honey bees. [4] Its timber is used for pulp and paper, fiber and particle board. It can produce qualified charcoal and potential sawn timber. [5] Especially, Melaleuca cajuputi essential oil is widely used for medicinal application. [6] Traditionally, cajeput oil is further used in herbal medicinal products to relieve symptoms of coughs, colds, as a laxative, general muscle relaxant and sedative. [7] Many pharmacological studies have shown that the major constituents of the Melaleuca cajuputi essential oil are 1,8cineol and α-terpineo l [4] which have antibacterial, antifungal and antiparasitic activities. Furthermore, there were also reports on the antiviral properties of essential oil from Melaleuca species , [7] [8] that the major monoterpenes and their derivatives of the essential oil from tea tree (Melaleuca alternifolia) including α-terpinene, γ-terpinene, α-pinene, p-cymene, 1,4terpineol, α-terpineol, thymol, citral, and 1,8-cineol had an antiviral activity against herpes simplex virus type 1 (HSV-1) in vitro. [9] The other results showed that tea tree oil and a few of its compounds (e. g. 1,4-terpinenol, terpinolene and αterpineol) possessed an inhibitory effect on influenza A/PR/ 8virus subtype H1 N1 replication at non-cytotoxic concentrations. [10] The literature review indicated that cajeput oil might be a promising medicine in preventing the infection of SARS-CoV-2. The fact is that cajeput essential oil is proved to have antifungal, anti-inflammatory, antioxidant activities. , [2] , [8] [11] According to some folk remedies, people use cajeput oil to help increase resistance, inhibit some bacteria and viruses, freshen air in the room, which are useful for providing postnatal care to women and newborns. [2] In this current SARS-CoV-2 pandemic, there has been a great demand for preventive and curative vaccines as well as measures to prevent and stop the spread of the disease. Based on the traditional usage of using cajeput essential oil in Vietnam, the idea of our research is using the docking simulation to predict the capability of molecular structure in the cajeput oil in inhibiting Angiotensin-converting enzyme 2 (ACE2) protein in human body, causing SARS-CoV-2 to lose its host receptor and destroy its protein (PDB6LU7) at the same time. Noticeably, the ACE2 protein is the host receptor of the SARS-CoV-2 and SARS-CoV, , [12] , [13] [14] therefore, if the ACE2 protein is inhibited, the SARS-CoV-2 could be prevented and treated. Further information of the ACE2 and PDB6LU7 proteins can be found at UniProtKB [15] and Worldwide Protein Data Bank, [16] respectively. In this study, we got the cajeput essential oil derived from Thua Thien Hue, Central Vietnam, and determined the ability to inhibit ACE2 protein in human body and PDB6LU7 protein of SARS-CoV-2. This is important in discovering medicinal herbs and developing treatments to prevent the SARS-CoV-2 pandemic. [16] Results and Discussion Docking simulation of compounds in Melaleuca cajuputi essential oil into ACE2 protein in human body and PBD6LU7 protein in the SARS-CoV-2 So far, to the best of our knowledge, there has been only available result related to coronavirus (SARS-CoV-2) resistance through docking simulation studies of organic compounds in garlic essential oil. [17] We demonstrate successful docking molecules of compounds in the Melaleuca cajuputi essential oil into the ACE2 and PDB6LU7 proteins to contribute to the orientation and encourage the use of natural products for prevention of SARS-COV-2 invasion. The 24 active compounds in the Melaleuca cajuputi essential oil were docked into the ACE2 and PDB6LU7 proteins using MOE 2015.10 program. The results indicate that the 10 compounds (identification details in Table 1 and Figure S1 in the Supporting Information Summary) have the strongest interactions on the selected proteins that are the ones with the most powerful inhibitory activity for SAR-CoV-2. The remaining compounds that have weak or almost no interactions with the ACE2 and PDBLU7 proteins are not presented in the main text. Please find more details in Table S1 in the Supporting Information Summary. Docking simulation results with the interactions between TA1-TA5; TA7; TA10; TA17-TA19 compounds (accounting for 70.9% in Melaleuca cajuputi oil) and 2 proteins (ACE2 and PDB6LU7) are presented in Figures 2 and 4 and Figure S2 (Supporting Information Summary). The results of docking score energy (DS) and root mean square deviation (RMSD) between these 10 compounds in Melaleuca cajuputi oil and proteins with various interactions including hydrogen bonds, cation-π bonds, π-π bonds, and ionic interactions as well as the interaction distance between amino acids and the active sites of compounds are shown in Tables 2 and 3 and Table S2 (Supporting Information Summary). It is found that Terpineol (TA2), Guaiol (TA5), and Linalool (TA19) are the best 3 compounds to have the binding ability toward the ACE2 and PDB6LU7 proteins with lots of convergence points and the inhibitory intensity of these compounds on ACE2 and PDB6LU7 proteins is almost similar (Table 2 and Figure 2) . Specifically, the TA2 has the strongest inhibitory effects on the ACE2 protein in human body and PDB6LU7 protein in the SARS-CoV-2 ( Figure 2 ) with the docking simulation results as follows: i) Docking of TA2-ACE2 has DS of À 11.0 kcal.mol À 1 and RMSD of 1.97 Å, while the site-site bonding interaction exhibiting the length of the binding of compound with amino acid between -OH of TA2 and amino acids Asn 103 is 2.25 Å and 2.08 Å for Gln 101; ii) Docking of TA2-SARS-CoV-2 has DS =-10.9 kcal.mol À 1 and RMSD = 1.16 Å, and site-site bonding interaction between -OH and amino acid His 163 is 2.38 Å. It can be realized that terpineol (TA2) is a monoterpene alcohol easily found in several essential oils such as petitgrain oil, cajuput oil, and pine oil. TA2 has a pleasant scent similar to lilacs and it is a common ingredient in perfumes, cosmetics and fragrances. [18] The monoterpene compounds including αterpinene, γ-terpinene, α-pinene, p-cymene, terpinen-4-ol, αterpineol, thymol, citral and 1,8-cineol in the Melaleuca alternifolia essential oil were reported against HSV-1 in vitro. [9] Besides, α-terpineol attracts great attention from scientists because it has a wide range of biological applications such as an effective inhibition of pathogenic viruses that helps protect our health, treatment of insect bites or wounds, relieving uncomfortable itching and reducing swelling rapidly, and curing respiratory infections. Moreover, terpineol in peppermint essential oil was proved as inhibitory agent for liver microsomes with IC 50 value of 14.8 μM. [19] The use of α-terpineol in [19] Therefore, we recognize that the results will support future research for the use of α-terpineol in the SARS-CoV-2 prevention. The inhibitory effects of Guliol (TA5) on the ACE2 and PDB6LU7 proteins are not significantly different compared to the TA2 based on the docking results as follows: iii) Docking of TA5-ACE2 has DS of À 11.1 kcal.mol À 1 and RMSD of 1.49 Å, and the site-site bonding interaction exhibiting the length of the binding of ligand with amino acid between -OH of TA2 and Asn 103 is 2.22 Å and 1.80 Å for Gln 101; iv) Docking of TA5-SARS-CoV-2 has DS of À 10.9 kcal.mol À 1 and RMSD of 0.84 Å, and site-site bonding interaction between -OH and amino acid His 163 is 2.51 Å. The fact is that TA5 is a sesquiterpenoid alcohol, and its molecule containing -OH group has strong interactions with amino acids Asn 103 and His 163 of ACE2 and PDB6LU7 proteins, respectively. Recent study has indicated that TA5 with strong antibacterial activity is also present in guaiacum and cypress pine oils . [20] This docking simulation gives a clear evidence for this activity of TA5. Lonalool (TA19) is capable of binding to ACE2 and PDB6LU7 proteins and has quite good interactions with amino acid Asn 210 of ACE2 protein, and Gly 143 as well as Cys 145 of PDB6LU7 protein (Table 2 and Figure 2) : v) Docking of TA19-ACE2 has DS =-10.9 kcal.mol À 1 , RMSD = 1.77 Å, and site-site bonding interaction between -OH and amino acid Asn 210 is 1.96 Å. vi) Docking of TA19-SARS-CoV-2 has DS of À 11.1 kcal.mol À 1 , and RMSD of 1.18 Å, and site-site bonding interactions between -OH and amino acid Gly 143 and Cys 145 are 2.89 Å and 2.82 Å, respectively. TA19 is a monoterpenoid having very flexible -OH group at C3 position, which is the active center of the TA19 that interacts with proteins. This compound can be found in oranges, tangerines, mangoes, other herbs such as basil, coriander, and lavender. It helps reduce stress, inhibit cancer cells, and fight inflammation. [21] Moreover, TA19 is also a major ingredient in some essential oils and used as aroma to produce perfumes, pharmaceuticals, and cosmetics. [21] Figures 3 and 4, and Table 3 present the docking results of ACE2 and PDB6LU7 proteins by Cineol (TA1), β-Selinenol (TA3), α-Eudesmol (TA4), and γ-Eudesmol (TA7). The anti-SARS-CoV-2 and anti-ACE2 activities of the cajuput essential oil composition into the PDB6LU7 and ACE2 proteins are shown in the following order: TA1 > TA3 > TA4 > TA7. Although TA1 is the highest content in the cajuput essential oil (31.6%) in cajeput essential oil, its interactions with the ACE2 and PDB6LU7 proteins are weaker compared to TA2, TA5, and TA19 (Figure 4) as follows: vii) TA1-ACE2 has DS =-10.8 kcal.mol À 1 and RMSD of 1.28 Å, the site-site bonding interaction between -Owith amino acid Tyr 196 of 2.32 Å. viii) TA1-SARS-CoV-2 exhibits the DS value of À 10.9 kcal.mol À 1 and RMSD = 1.00 Å, and the site-site bonding interaction between -O-with amino acid Glu 143 is 2.31 Å. TA1 is present in Cinnamomum camphora, Laurus nobilis, Artemisia vulgaris, Ocimum basilicum, Salvia rosmarinus. TA1 is commonly used in dentistry or used as cough remedy at an acceptable amount. [22] Some studies have shown that TA1 has good antimicrobial, antiviral, antiinflammatory and antioxidant properties. , [22] [23] TA3 exhibits the interactions with ACE2 and PDB6LU7 proteins in which the docking results are displayed as follows: ix) TA3-ACE2 has DS value of À 10.1 kcal.mol À 1 and RMSD = 1.80 Å. The site-site bond length between -OH and amino acid Tyr 202 is 2.15 Å. x) TA3-SARS-CoV-2 has DS =-10.8 kcal.mol À 1 , RMSD of 2.89 Å, and the bonding interaction between -OH ligand with amino acid Leu 141 is 1.97 Å. The TA3 compound can be found in fischeriana essential oil, devrra tortuosa oil, which have been recently studied on antioxidant activity and antifungal activity. , [24] [25] Docking simulation of the interactions between TA4 and the ACE2 and PDB6LU7 proteins are presented as TA4-ACE2 and TA4-PDB6LU7. TA4-ACE2 has DS =-9.9 kcal.mol À 1 , RMSD = 1.96 Å and the bonding length interaction between -OH with amino acid Asp 208 is 2.83 Å while TA4-SARS-CoV-2 consists of DS =-9.4 kcal.mol À 1 , RMSD gives 2.18 Å, and site-site distance between -OH and amino acid His 163 of 2.23 Å. TA4 is an eudesmane sesquiterpenoid which exhibits a double bond between C3 and C4 in compound; the flexible -OH group at C11 has strong interactions with amino acids Asp 208 of ACE2 and His 163 of PDB6LU7. The fact is that TA4 can also be found in Asteraceae oil, Centaurea appendicigera oil, and Centaurea helenioides oil and has great antibacterial and anti-oxidant activities. [26] The docking simulation results of TA7 and the two proteins indicate a similar trend but rather weak interaction compared to TA4. TA7-ACE2 has DS =-9.9 kcal.mol À 1 and RMSD = 1.21 Å, the site-site bonding between -OH of compound with Glu 208 of 2.37 Å. TA7-SARS-CoV-2 has the similar DS value of À 9.4 kcal.mol À 1 and RMSD = 1.15 Å, and the bond distance of -OH and amino acid Leu 141 is 1.90 Å. The above results can affirm that the strongest interactions of the most active compounds of cajeput essential oil with the ACE2 and PDB6LU7 proteins are described in this order: TA2 � TA5 � TA19 > TA1 > TA3 > TA4 > TA7. The compounds TA10, TA17, and TA18 indicate the weak inhibitory effects on the ACE2 and PDB6LU7 proteins (Table S2 and Figure S1 in the Supporting Information Summary). The docking results of TA10-ACE2 exhibit DS =-9.7 kcal.mol À 1 , RMSD = 0.95 Å and the site-site distance between -OH and Glu 208 is 1.77 Å. TA10-SARS-CoV-2 has DS =-9.2 kcal.mol À 1 , RMSD = 2.75 Å and site-site distance between -OH with Gln 189 is 2.90 Å. DS and RMSD of TA18-ACE2 are bit larger than those of TA18-SARS-COV-2. Noticeably, TA17 has the weakest inhibitory effects on the ACE2 and PDB6LU7 proteins in the 10 compounds of the cajeput essential oil. TA17 -ACE2 has DS of À 8.3 kcal.mol À 1 and RMSD of 1.30 Å, and site-site bonding between carbon of compound and aroma of amino acid Phe 40 is 3.03 Å. TA17-SARS-CoV-2 gives DS =-7.3 kcal.mol À 1 , RMSD = 2.12 Å, and the site-site bonding interaction between carbon of ligand and aroma of amino acid His 41 is 2.72 Å. Based on the docking simulation results analyzed above, the interaction order of the compounds in the cajeput essential oil is listed as follows: Interestingly, the co-interaction of the above 10 compounds (10TA) with ACE2 and PDB6LU7 proteins was also investigated and shown in Figure 5 . There were 12 and 13 interactions between 10TA in cajuput essential oil with the ACE2 and PDB6LU7 proteins, respectively. The results exhibit the main interactions of -OH and -O-groups in 10TA compounds TA1-TA5, TA7, TA10, TA17-TA19 with amino acids including His 163, His 41, Leu 141, Gly 143, Glu 143, Gln 189 of PDB6LU7 protein, and average DS value is À 9.9 kcal.mol À 1 . The main interactions of 10TA with amino acids in ACE2 are Asn 103, Asn 210, Tyr 196, Tyr 202, Glu 208, Glu 209, Asp 208, Phe 40 of ACE2 protein, and its average DS value is À 10.1 kcal.mol À 1 . The study demonstrates the strong inhibitory effects of active sites of the compounds in cajeput essential oil on the ACE2 and PDB6LU7 proteins. Therefore, they are crucial for the use of the cajeput essential oil in preventing viruses that cause flu, pneumonia in general and the SARS-CoV-2 in specific. Melaleuca cajuputi essential oil is a valuable natural pharmaceutical source that can inhibit the ACE2 protein -the host receptor of PDB6LU7 protein in the SARS-CoV-2 as well as attack the PDB6LU7 protein at the same time. Docking simulation results indicated the active bonding sites of 10 (TA1-TA5, TA7, TA10, TA17 -TA19) out of 24 compounds of the Melaleuca cajuputi essential oil are capable of inhibiting ACE2 and resisting PDB6LU7 protein in SARS-CoV-2. The three substances Terpineol (TA2), Guaiol (TA5), and Linalool (TA19) have the strongest inhibitory effects on ACE2 and PDB6LU7 proteins. The order of the active compounds inhibiting the ACE2 and PDB6LU7 proteins is: TA2 � TA5 � TA19 > TA1 > TA3 > TA4 > TA7 > TA10 > TA18 > TA17. Interestingly, the docking simulation results also indicate the synergistic interactions of 10 substances in the Melaleuca cajuputi essential oil exhibit the significant inhibition into the ACE2 and PDB6LU7 proteins. This prevents protein maturation of the virus and the spread of infection. This study brings out a recommendation to use cajeput essential oil in preventing infections and limiting the spread of SARS-CoV-2. Docking simulation: The docking method requires information of the protein structures and the ligand structures. When the docking is conducted, all of these structures are used to simulate the interactions between the ligands and the proteins and then the bonding and docking score energy results are evaluated. The result is displayed as the ligands configurations and the docking score (DS), root mean square deviation (RMSD), types of interactions, and respective distances between the ligands and proteins. The molecular docking modeling has 3 steps. , [27] , [28] , [29] [30] Structures of ACE2 and PDB6LU7 proteins are available at UniProtKB [14] and Worldwide Protein Data Bank, [15] respectively. Water molecules are deleted in the protein structure with the Sequence Editor of the MOE 2015.10; Quickprep tool is then used to repair the protein structure and 3D protonation: The protein active zone is determined based on the ligand position within a radius of 4.5 Å and the presence of important amino acids, then the protein structures are saved in *.pdb format. The ligands, which are the 24 compounds in the Melaleuca cajuputi essential oil, were optimized via: Conj Grad for minina energy; Termination for energy change = 0.0001 kcal.mol À 1 ; Max interactions = 1000; Modify charge: gasteiger -Huckel. Next, molecular dynamics on the MOE 2015.10 system are performed and the molecule structures in *.sdf format are saved. The docking simulation parameters displayed via the number of poses to keep for further analysis of interaction is 10; the maximum number of solutions per iteration is 1000; maximum number of solutions per fragmentation is 200. Docking score energy (DS) demonstrates the binding between ligands and the site À site distance of ligands and proteins. It is also necessary to analyze the interactions formed between ligands and important amino acids in the site À site distances. Hydrogen bonds, ion bonds, arene -arene (π -π), cationarenes (cations -π), and van der Waals interactions are detected by contacts with hydrophilic, hydrophobic, and solvent interaction. Evaluation of docking score energy (DS) exhibits the analysis of interaction between the compounds of the essential oil and targeted proteins, and performance of interaction on 2D, 3D plane using MOE 2015.10 software. Various interactions such as hydrogen bonds, cation-π, π-π bonds, and ionic interactions, interactions distance between amino acids and the active sites of compounds are plotted. Besides, van der Waals interactions are detected by contact with hydrophilic and hydrophobic surfaces between the compounds and the bonding point. Finally, the ligand-protein interactions are determined. Supporting Information Summary is supplied for Experiment Section, docking simulation results of other compounds in Melaleuca cajuputi essential oil that are not shown in the main text of this paper. Microbes Infect Essential leaf oils from Melaleuca cajuputi Australian trees and shrubs: species for land rehabilitation and farm planting, ACIAR Monograph Proceedings IEEE INFOCOM 2000. Conference on Computer Communications. Nineteenth Annual Joint Conference of the IEEE Computer and Communications Societies (Cat. No. 00CH37064) Plant resources of South-East AsiaI: Essential Oil Plants Proceedings 3rd IEEE Symposium on Application-Specific Systems and Software Engineering Technology This research is supported by the Vietnam National Foundation for Science and Technology Development (NAFOSTED) under grant number 02/2020/ĐX. The authors declare no conflict of interest.Keywords: ACE2 · Melaleuca cajuputi oil · Biological activity · PDB6LU7 · Computational Chemistry · SARS-CoV-2