key: cord-0718490-631u9uov authors: Jindal, Divya; Rani, Vibha title: In Silico Studies of Phytoconstituents from Piper longum and Ocimum sanctum as ACE2 and TMRSS2 Inhibitors: Strategies to Combat COVID-19 date: 2022-02-14 journal: Appl Biochem Biotechnol DOI: 10.1007/s12010-022-03827-6 sha: 0a221d79392d7222319af44bf16b54fb9375d833 doc_id: 718490 cord_uid: 631u9uov The recent pandemic due to the COVID-19 virus has caused a global catastrophe. ACE2 and TMPRSS2 are recognized as key targets for viral entry into the host cells. The pandemic has led to the utilization of many synthetic drugs; however, due to various side effects, there is still no effective drug available against the virus. Several natural approaches have been devised, including herbal and ayurvedic medicines, that have proven to be effective against the COVID-19 virus. In the present study, the effect of phytocompounds of Piper longum and Ocimum sanctum on ACE2 and TRMPSS2 proteins has been studied. The in silico study is done using computational tools of networks of protein–protein interaction, molecular docking, and drug assessment in terms of physicochemical properties, drug-likeness, lipophilicity, water solubility, and pharmacokinetics. Out of selected phytoconstituents, vicenin 2, rosmarinic acid, and orientin were found to have the highest efficacy in terms of molecular interaction and drug-likeness properties against ACE2 and TMPRSS2 host receptor proteins. Our in silico study proposes the therapeutic potential of phytocompounds from Piper longum and Ocimum sanctum in modulating ACE2 and TMPRSS2 expression. Targeting ACE2 and TMPRSS2 against the SARS-CoV2 by phytomolecules can serve as a rational approach for designing future anti-COVID drugs. The witnessed outbreak of COVID-19 in 2020, caused by coronavirus familyCoronavirdae, specifically SARS-CoV2, leads to a life crisis of 100,000,000 people with more than 2,200,000 deaths worldwide, lasting for around 2 years now. A major concern is their mortality rate which is increasingly ranging from 1 to > 5% [1] . There are various forms of respiratory syndromes such as MERS and SARS, caused by the family of coronavirus [2] . Diabetes and cancer are the major risk factors for the progression and prognosis of COVID-19 due to the abundance of expression of ACE2 in different organs [26, 27] . In addition, the anecdotal models of evidence evaluated that coronavirus unveils the variations in morbidity and mortality among sexes. Males are almost 3 times more prone to the virus than females [28] . The diversity of plants in the treatment of COVID-19 has shown interaction with other active enzymes involved in its pathogenic pathway. Though certain medications have been used for COVID-19 patients, these synthetic drugs exhibit toxicity which may lead to hyperstimulation of cells and reduction in function [29] . For example, hydroxy-chloroquininone and 2-DG (2-deoxy-D-glucose) have been recently used to suppress COVID infection but they have been reported to have increased melanin, skin rashes, dizziness, ocular, and cardiac toxicity. Hence, there is an immediate and utmost requirement for an alternative that exhibits phytotherapeutic assistance to COVID-19 patients and synergistically improves the activity of an organism. Natural herbs are being used for a century, with their development and production in the form of different formulations. In India, almost 70% population relies on using traditional medicines to cure many prevalent diseases [30] . They produce a wide range of bioactive compounds, with probable anti-microbial and anti-viral activities, etc., against a lot of agents. In our study, we have used two plants, Piper longum and Ocimum sanctum, for their prominent role as ACE2 and TRMPSS2 modulators. Piper longum or Pipali is generally used in houses for cooking purposes as a spice, and its medicinal role has been discussed in several studies. The plant has been a widely used remedy for several diseases like tuberculosis, gonorrhea, respiratory infections, immune-modulatory and anti-tumor activities, CNS depressant, anti-inflammatory, and anti-oxidative properties [31, 32] . Piper longum is a deciduous aromatic climber plant that has perennial roots and shrubs. These are monoecious where male and female flowers are present in separate structures on the same plant. This is mostly grown in India, Nepal, Sri Lanka, and the Philippines. P. longum, in general, comprises piperine, anthraglycosides, sterols, alkaloids, vicenin, and several other compounds [33] . Ocimum sanctum, Tulsi, is holy basil encompassing prominent importance in Indian Hindu Culture. Tulsi is one of the primogenital herbs of the family Lamiaceae and is found in almost every household because of its medicinal and nutritional properties. The ecological supremacy of Tulsi indicates their photosynthetic and pharmaceutical efficiency, as seen in various studies [34] . O. sanctum e is anti-viral, anti-bacterial, and anti-cancerous and treats respiratory problems and many other diseases [35] . They consist of extensive varieties of secondary metabolites like tannins, phenolic, alkaloids, and flavonoids, which support enhancing growth and immunity responses [36] . The present study is designed to observe the effect and efficacy of phytocompounds from selected herbs as ACE2 and TRMPSS2 inhibitors. To the best of our knowledge, there has been no study conducted so far. The proposed study is novel and has huge potential in a current pandemic. 3D structures of ACE2 and TMPRSS2 were downloaded from Protein Data Bank (ACE2 PDB ID: 1R4L and TMPRSS2 PDB ID: 7MEQ). These proteins were initially associated with one of their known inhibitors. The receptor proteins were prepared using the Prime Module of Schrodinger Suite. The ligands were removed from both the receptors, and the final structure was cleaned and saved in PDB format. All the proteins were prepared by removing water, hydrogen, and existing lead components like ions. Chain A from the proteins was deleted, and respected force fields and charges were added. The grid was formed around the active sites of the known inhibitors using the Receptor Grid Generation tool with a dimension of 10 Armstrong [37] . For the discovery of any drug or effective compound, in silico studies provide crucial clues. To combat cost and time, in silico approach becomes the foremost choice in research. For a molecule to be a drug, it has to follow Lipinski's rule which states five categories to be verified, also known as a "rule of thumb." These include that a chemical entity must have a molecular weight less than or equal to 500 Da, 5 or fewer hydrogen bonds donors, 10 or fewer hydrogen acceptors, MLOGP value less than and equal to 4.15, and molecular refractivity has to be in between 30 and 140 [38] . The ligands were downloaded from IMPPAT (Indian Medicinal Plants, Phytochemistry, and Therapeutics) in SDF format. 2-DG and nafamostat were taken as control drugs. The ligands were prepared using the Ligprep module of the Schrodinger suite. The ligands were neutralized and desalted, and Epik ionization was provided. Molecular docking predicts the binding energies based on the affinity of the ligand with the receptor. The prepared proteins and ligands were used for molecular docking studies. Both of them were processed and docked in the ligand docking tab. The zip document of the grid file of protein was uploaded against the ligand output of the file. Extra precision (XP) docking was done to eliminate respective penalties and errors. After docking, the structures were visualized in Maestro software. Top scoring outputs processed were analyzed in the ligand interaction diagram to check interaction and amino acid involved [31, 39] . The SwissADME is a freely accessible web tool (http:// www. swiss adme. ch/) and provides an easy way to analyze results in a computer-aided drug designing platform. The webbased tool provides pharmacokinetics data, physiochemical properties, lipophilicity, water solubility, and drug-likeness, and illustrates the molecule in a boiled egg, showing whether the ligand can cross the blood-brain barrier and gastrointestinal tract [32, 40] . The selected phytomolecules were analyzed for their ADME profiling. To classify cellular pathways which were enriching host factors interacting with coronavirus spike protein, KEGG (https:// www. genome. jp/ kegg/ pathw ay. html) annotations were retrieved [41] . To evaluate combining interactive information from both receptor proteins, i.e., ACE2 and TMPRSS2, pathways involved in different functions were being implied. Search Tool for Retrieval of Interacting Genes (STRING) (https:// string-db. org/) Database predicts and integrates protein-protein interaction to identify functional relationships and interactions between proteins. To seek potential interactions between genes involved in SARS-CoV, the STRING tool was employed. STRING provided a platform that postulated several nodes, edges, average node degrees, protein-protein interaction (PPI) enrichment p-values, and average and local clustering coefficient. STRING provided possible biological processes, molecular function, and cellular components of candidate genes studied [42] . Molinspiration (http:// www. molin spira tion. com/ cgi-bin/ prope rties) predicts the drug resemblance properties of the compound dependent on various descriptors. Drugs entering the body should tie to an organic molecule to communicate its movement. Along their path, the bioactivity of compounds was anticipated by utilizing the Molinspiration tool which gave a bioactivity score of the phytocompound against the human receptors like GPCRs, ionic channels, kinases, various receptors, proteases, and proteins. A complex is considered to be dynamic if the bioactivity score is more than 0.0, modestly dynamic if somewhere in the range of − 5.0 and 0.0, and idle if under − 5.0 [43] . We first conducted molecular docking to predict the probable protein-ligand interactions by minimizing the energy of the ligands and calculating their binding energies. Binding of ligand to respective receptor protein causes an inhibition of ligand toward enzyme, which results in a demonstration of the feasibility of a biochemical reaction. Docking algorithms utilize inhibitory and activator properties of the ligand with receptor protein and form a relationship between the drug's structure and cytotoxicity. We selected approximately a hundred phytocompounds from Piper longum and Ocimum sanctum using the IMPPAT database. ACE2 and TMPRSS2 were docked with all phytostructures. The docking score of significant interactions (~ 15) is shown in Table 1 . However, vicenin 2, rosmarinic acid, and orientin were found to be most efficient. For ACE2, we obtained the highest binding energy of − 11.755 kcal/mol with vicenin 2. Also, for TMPRSS2, vicenin has shown the binding energy of − 7.913 kcal/mol. 2-DG (an ACE2 inhibitor) and nafamostat (a TMPRSS2 inhibitor) were taken as controls for both the targets as it has been currently used and have shown good efficacy as synthetic drugs against SARS-CoV infection. The binding energy of 2-DG with ACE2 was found to be − 6.05 kcal/ mol, while for TMPRSS2, nafamostat binding energy was − 2.188 kcal/mol. The interaction poses of vicenin 2, rosmarinic acid, and orientin with ACE2 and TMPRSS2 are shown in Figs. 1 and 2 , respectively. The amino acids shared by the selected phytocompounds and controls were found to be common suggesting the involvement of a similar biological mechanism. For ACE2, the binding activities were found to be in the following order: vicenin 2 > rosmarinic acid > orientin > L-ascorbic acid > 2,6-diaminohexanoic acid > luteolin 7-O-glucuronide > lignans machilin F > 2,4-dihydroxycinnamic acid > 7,3′,4′,5′-tetrahydroxyflavone > luteolin > oleic acid > 2-coumarinate > 1-(4-hydroxy-3-methoxyphenyl)-1,2,3-tris(4-allyl-2-methoxyphenoxy) propane > palmitic acid > 2-deoxy-glucose, whereas for TMPRSS2, the binding activities were as follows: vicenin2 > rosmarinic acid > luteolin 7-O-glucuronide > orientin > luteolin > eupatin > MLS000877024 > esculetin > 1-(4-hydroxy-3-methoxyphenyl)-1,2,3-tris(4-allyl-2-methoxyphenoxy) propane > molludistin > lignans machilin F > nevadensin > apigenin > nafamostat. Our in silico molecular docking studies suggested the significantly higher binding activities of selected phytocompounds with ACE2 and TMPRSS2 as compared to the presentday synthetic inhibitors used in the study which shows the therapeutic relevance of Piper longum and Ocimum sanctum in combating COVID complications. To further determine the drug-likeliness properties of selected phytocompounds, ADME analysis was done. Lipinski et al. formulated certain parameters by defining their physiochemical ranges to be an oral drug, called their drug-likeness. Rosmarinic acid showed promising results with control 2-DG and nafamostat against ACE2 and TMPRSS2. The physiochemical properties like molecular weight, number of hydrogen bonds acceptors and donors, and number of rotatable bonds show a great similarity with 2-DG and nafamostat. The lipophilicity results of rosmarinic acid showed positive values which demonstrates the greater affinity toward the lipid environment. In comparison to controls, 2-DG has less affinity and nafamostat confers better affinity in a lipid environment. In addition, rosmarinic acid has 0 violations and almost the same bioavailability score of 0.55 in comparison with 2-DG and nafamostat. They all had features to be a candidate drug molecule; however, rosmarinic acid is found to be a better option. The absorption and bioavailability of the selected compounds can be improved in the future with various formulations ( Table 2) . We next predicted the target network to observe the interactors of ACE2 and TMPRSS2 using STRING. ACE2 was found to be interacting with several proteins as shown in Fig. 3 . These proteins were found to engage in various biological processes such as in renin-angiotensin regulation of aldosterone production, GPCR-signaling pathways, cell development, brain renin-angiotensin system, angiotensin maturation, cell growth involved in cardiac muscle, amyloid-beta metabolic process, and vasodilation. Angiotensin type II receptor activity, dipeptidyl-peptidase activity, aminopeptidase, carboxypeptidase, and exopeptidase activity were their majorly concerned molecular functions. On the other hand, TMPRSS2 was observed to be interacting with proteins majorly engaged in the regulation of epithelial cell proliferation, prostate gland development and their growth, protein kinase B signaling, regulation of insulin-like growth factor receptor signaling pathway, and androgen receptor signaling pathway. The foremost KEGG analysis was conducted to further analyze the biological relevance of the identified target and the possible involvement of ACE2 and TMPRSS2. A total of 20 key targets including AR, PTEN, and ERG were identified through PPI network analysis. These targets were mainly focused on the biological processes with ACE2 and TMPRSS2. The KEGG enrichment manifested signaling pathways that were closely related to COVID-19 and renin-angiotensin system, as shown in Figs. 4 and 5. Furthermore, in the COVID-19 pathway, TMPRSS2 and NRP1 are involved in the entry of spike protein via membrane fusion or endocytosis. Another pathway is involved in the renin-angiotensin system wherein Angiotensin (1-7) and Angiotensin II are entered in the host cell by MAS1 and AT1R and further facilitated by ADAM17, TNFR, EGFR, and TLR2/4. Lastly, we predicted bioactivity scores to identify the potency of phytomolecules studied (Table 3 ). Our result showed that rosmarinic acid and orientin might essentially serve as an interface with ACE2 and TMPRSS2 as they acted as enzyme inhibitors. Bioactivity scores of herbal compounds and controls were comparable indicating that these phytocompounds might be utilized as an alternative and in improving treatment to COVID. The prediction of possible targets leads to mechanistic pathways involved in coronavirus integration into the human body and associated chronic diseases. ACE2 and TMPRSS2 facilitate the virus entry into the human genome. The ligands have shown explicit findings in recent times and have been proposed to show major impacts on coronavirus. Molecular docking of the compounds with ligands has provided good binding energy, and the study is projected to assess ligands in vitro analysis. Out of all compounds, rosmarinic acid and orientin have revealed effective energies and follow Lipinski's rule to have drug-likeness properties. Downstream to molecular targets, these proteins are also associated with other molecular functions like cell signaling, apoptosis, DNA-binding transcription activator activity, and amyloid-beta metabolic process, which eventually are targeted when dysregulated. We propose phytoconstituents of Piper longum and Ocimum sanctum as ACE2 and TMRSS2 inhibitors: strategies to combat COVID-19 (Fig. 6) . However, such understandings need further approval utilizing distinctive in vitro experiments and molecular dynamics simulation. The traditional use of Indian medicinal plants is believed to be a huge foundation for the treatment of diseases. For this study, we have investigated the potential of some phytocompounds which are being extracted from the most common medicinal plants: Piper longum and Ocimum sanctum. The results are projected to seize the attention of scientists in the area of drug discovery against COVID-19 for which no specific and successful drug has been discovered with natural phytocompounds. The approach is simply using natural bioactive components, and the same investigation can also be done by various available medicinal herbs. In addition, our study is extended to see the effects of other diseases contributing to the progression of the coronavirus. Acknowledgements The authors acknowledge Jaypee Institute of Information Technology, Noida, for providing the entire infrastructure to complete this project. Funding This study did not require any specific funding support. This work was supported by the Institutional fund of the host organization. In silico data and material are present in the manuscript. Code Availability No code applicable. Fahd Al-Mulla. The two faces of ACE2: The role of ACE2 receptor and its polymorphisms in hypertension and COVID-19 How do viral infections predispose patients to bacterial infections? The origin, transmission and clinical therapies on coronavirus disease 2019 (COVID-19) outbreak -An update on the status. Military A decade after SARS: Strategies for controlling emerging coronaviruses Vitamin D regulation of the renin-angiotensin system Effect of Captopril, an angiotensin-converting enzyme inhibitor, in patients with angina pectoris and heart failure The discovery of angiotensin-converting enzyme 2 and its role in acute lung injury in mice A liquid chromatography-assisted assay for angiotensin-converting enzyme (peptidyl dipeptidase) in serum A novel angiotensin-converting enzyme-related carboxypeptidase (ACE2) converts angiotensin I to angiotensin 1-9 Angiotensin converting enzyme inhibitors, angiotensin receptor blockers, mineralocorticoid receptor antagonists and SARS-CoV-2 infection Expression of the COVID-19 receptor ACE2 in the humanconjunctiva A crucial role of angiotensin converting enzyme 2 (ACE2) in SARS coronavirus-induced lung injury The COVID-19 epidemic Evidence that TMPRSS2 activates the severe acute respiratory syndrome coronavirus spike protein for membrane fusion and reduces viral control by the humoral immune response The role of the renin-angiotensin-aldosterone system in heart failure ACE2: From vasopeptidase to SARS virus receptor Multiple functions of angiotensin-converting enzyme 2 and its relevance in cardiovascular diseases Angiotensin-converting enzyme 2 (ACE2) as a SARS-CoV-2 receptor: Molecular mechanisms and potential therapeutic target The pivotal link between ACE2 deficiency and SARS-CoV-2 infection Transmembrane protease serine 2 (TMPRSS2). Science-Business eXchange The structural biology of type I viral membrane fusion Viral membrane fusion at single pore resolution Mechanisms of coronavirus cell entry mediated by the viral spike protein Proteolytic activation of the SARS-coronavirus spike protein: Cutting enzymes at the cutting edge of antiviral research SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor COVID-19 and androgen-targeted therapy for prostate cancer patients Association of hypertension, diabetes, stroke, cancer, kidney disease, and high-cholesterol with COVID-19 disease severity and fatality: A systematic review Male sex identified by global COVID-19 meta-analysis as a risk factor for death and ITU admission The new classes of synthetic illicit drugs can significantly harm the brain: A neuro imaging perspective with full review of MRI findings The growing use of herbal medicines: Issues relating to adverse reactions and challenges in monitoring safety Extra precision glide: Docking and scoring incorporating a model of hydrophobic enclosure for protein− ligand complexes A BOILED-Egg to predict gastrointestinal absorption and brain penetration of small molecules Potential antiamoebic property of the roots of Piper longum Linn Photosynthetic performance of Ocimum sanctum morphotypes in a semiarid region Ocimum sanctum), excellent source of phytochemicals Phytochemical constituents of some nigerian plants A hierarchical approach to allatom protein loop prediction Lead-and drug-like compounds: The rule-of-five revolution Challenges and advances in computational docking: 2009 in review SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules Database for annotation, visualization, and integrated discovery STRING v10: Protein-protein interaction networks, integrated over the tree of life Pharmacological investigations and Petra/Osiris/Molinspiration (POM) analyses of newly synthesized potentially bioactive organotin(IV) carboxylates