<?xml version="1.0" encoding="utf-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Microbiol.</journal-id>
<journal-title>Frontiers in Microbiology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Microbiol.</abbrev-journal-title>
<issn pub-type="epub">1664-302X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmicb.2022.1073419</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Microbiology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Targeting envelope proteins of poxviruses to repurpose phytochemicals against monkeypox: An <italic>in silico</italic> investigation</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Gulati</surname>
<given-names>Pallavi</given-names>
</name>
<xref rid="aff1" ref-type="aff"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Chadha</surname>
<given-names>Jatin</given-names>
</name>
<xref rid="aff2" ref-type="aff"><sup>2</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/1230449/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Harjai</surname>
<given-names>Kusum</given-names>
</name>
<xref rid="aff2" ref-type="aff"><sup>2</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/699103/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Singh</surname>
<given-names>Sandeepa</given-names>
</name>
<xref rid="aff3" ref-type="aff"><sup>3</sup></xref>
<xref rid="c001" ref-type="corresp"><sup>&#x002A;</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/2057559/overview"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Department of Microbiology, University of Delhi</institution>, <addr-line>New Delhi</addr-line>, <country>India</country></aff>
<aff id="aff2"><sup>2</sup><institution>Department of Microbiology, Panjab University</institution>, <addr-line>Chandigarh</addr-line>, <country>India</country></aff>
<aff id="aff3"><sup>3</sup><institution>Department of Botany, Maitreyi College, University of Delhi</institution>, <addr-line>New Delhi</addr-line>, <country>India</country></aff>
<author-notes>
<fn id="fn0001" fn-type="edited-by"><p>Edited by: Naveen Kumar, ICAR-National Institute of High Security Animal Diseases (ICAR-NIHSAD), India</p></fn>
<fn id="fn0002" fn-type="edited-by"><p>Reviewed by: Anand Anbarasu, VIT University, India; Rajesh Kumar Pathak, Chung-Ang University, South Korea</p></fn>
<corresp id="c001">&#x002A;Correspondence: Sandeepa Singh, &#x02709; <email>ssingh1@maitreyi.du.ac.in</email></corresp>
<fn id="fn0003" fn-type="other"><p>This article was submitted to Virology, a section of the journal Frontiers in Microbiology</p></fn>
</author-notes>
<pub-date pub-type="epub">
<day>05</day>
<month>01</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>13</volume>
<elocation-id>1073419</elocation-id>
<history>
<date date-type="received">
<day>18</day>
<month>10</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>16</day>
<month>12</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2023 Gulati, Chadha, Harjai and Singh.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Gulati, Chadha, Harjai and Singh</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>The monkeypox virus (MPXV) has become a major threat due to the increasing global caseload and the ongoing multi-country outbreak in non-endemic territories. Due to limited research in this avenue and the lack of intervention strategies, the present study was aimed to virtually screen bioactive phytochemicals against envelope proteins of MPXV <italic>via</italic> rigorous computational approaches. Molecular docking, molecular dynamic (MD) simulations, and MM/PBSA analysis were used to investigate the binding affinity of 12 phytochemicals against three envelope proteins of MPXV, <italic>viz</italic>., D13, A26, and H3. Silibinin, oleanolic acid, and ursolic acid were computationally identified as potential phytochemicals that showed strong binding affinity toward all the tested structural proteins of MPXV through molecular docking. The stability of the docked complexes was also confirmed by MD simulations and MM/PBSA calculations. Results from the iMODS server also complemented the findings from molecular docking and MD simulations. ADME analysis also computationally confirmed the drug-like properties of the phytochemicals, thereby asserting their suitability for consumption. Hence, this study envisions the candidature of bioactive phytochemicals as promising inhibitors against the envelope proteins of the MPXV, serving as template molecules that could further be experimentally evaluated for their efficacy against monkeypox.</p>
</abstract>
<kwd-group>
<kwd>poxviruses</kwd>
<kwd>monkeypox virus</kwd>
<kwd>envelope proteins</kwd>
<kwd>drug repurposing</kwd>
<kwd>bioactive phytochemicals</kwd>
<kwd>molecular docking</kwd>
<kwd>molecular simulation</kwd>
</kwd-group>
<counts>
<fig-count count="11"/>
<table-count count="4"/>
<equation-count count="0"/>
<ref-count count="52"/>
<page-count count="19"/>
<word-count count="9014"/>
</counts>
</article-meta>
</front>
<body>
<sec id="sec1" sec-type="intro">
<label>1.</label>
<title>Introduction</title>
<p>In recent years, the world has witnessed the resurgence of several viral infections, including the Nipah virus, influenza virus (H1N1, H5N1), Ebola virus, Zika virus, and the recent coronavirus outbreaks (MERS-CoV and SARS-CoV-2; <xref ref-type="bibr" rid="ref34">Mourya et al., 2019</xref>; <xref ref-type="bibr" rid="ref9">Chadha et al., 2022c</xref>). The latest addition to this list is the monkeypox virus (MPXV), one of the four human pathogenic orthopoxvirus species, grouped along with variola, cowpox, and vaccinia virus (<xref ref-type="bibr" rid="ref28">Lam et al., 2022</xref>). Monkeypox (MPX) is an emerging zoonotic disease caused by the MPXV, which was first discovered in 1958 following local outbreaks of a smallpox-like disease in monkey colonies maintained at research facilities in Africa and Denmark (<xref ref-type="bibr" rid="ref8">Chadha et al., 2022b</xref>). However, the first clinical case of MPX (in humans) was reported in 1970 in the Equatorial province of the Democratic Republic of Congo (DRC; <xref ref-type="bibr" rid="ref32">Marennikova et al., 1972</xref>). Since then, sporadic MPX outbreaks have been reported across several African countries, particularly the DRC, Republic of The Congo, Central African Republic, Nigeria, Liberia, Ghana, and Cameroon (<xref ref-type="bibr" rid="ref4">Bunge et al., 2022</xref>). The MPXV has been endemic to sub-Saharan countries for decades, but in mid-2022 (amid the COVID-19 pandemic), an unexpected surge in MPX infections was recorded globally (<xref ref-type="bibr" rid="ref90">WHO, 2022</xref>). Multi-country outbreak of MPX was reported across several non-endemic countries, raising concern among clinicians and epidemiologists (<xref ref-type="bibr" rid="ref8">Chadha et al., 2022b</xref>). With the increasing global caseload, the World Health Organization also declared the ongoing multi-country outbreak of MPX as a Public Health Emergency of International Concern (PHEIC) on 23rd July 2022.</p>
<p>The MPXV is an enveloped virus with a dsDNA genome of nearly 190 kbp and harbors a dumbbell-shaped core with lateral bodies (<xref ref-type="bibr" rid="ref22">Kugelman et al., 2014</xref>). The virus is known to infect various mammalian species, including primates, but its natural animal reservoir remains unknown. However, some investigations suggest that rodents are the primary reservoir of the MPXV (<xref ref-type="bibr" rid="ref38">Reynolds et al., 2012</xref>). To date, two distinguishable clades of MPXV have been identified: West African and Central African (Congo Basin). The latter is highly transmissible among humans and exhibits potentially high virulence and case fatality rate (~ 11%) (<xref ref-type="bibr" rid="ref44">Sklenovsk&#x00E1; and Van Ranst, 2018</xref>). MPX transmission has been strongly associated with animal bites, skin abrasions, prolonged or direct contact with bodily fluids of infected animals or humans, respiratory droplets, and contaminated inanimate objects (<xref ref-type="bibr" rid="ref29">Ligon, 2004</xref>). However, the transmission of presently-circulating MPXV strains through sexual routes is still debated (<xref ref-type="bibr" rid="ref8">Chadha et al., 2022b</xref>). The incubation period of MPXV ranges from 3 to 14 days, and the clinical manifestations are seldom mistaken for other viral infections like herpes simplex virus and chickenpox. Symptoms of MPX are similar to that of smallpox, including fever accompanied by blisters, muscle pain, rashes, chills, and respiratory illness (<xref ref-type="bibr" rid="ref45">Titanji et al., 2022</xref>). The present-day treatment and clinical management of MPX involve the administration of repurposed antiviral drugs (cidofovir, brincidofovir, and tecovirimat) and preventive vaccines (ACAM2000 and JYNNEOS; <xref ref-type="bibr" rid="ref8">Chadha et al., 2022b</xref>). Moreover, around 50 single-nucleotide polymorphisms have been identified in the genome of 2022 MPXV, which remained undetected in the previously circulating MPXV strains between 2018 and 2019 (<xref ref-type="bibr" rid="ref19">Isidro et al., 2022</xref>). This increases the gravity of the ongoing multi-country outbreak, pointing toward the &#x201C;accelerated evolution&#x201D; of MPXV. Hence, there is a pressing need to develop novel therapeutics against the MPXV. In this direction, drug repurposing proves to be an effective alternative since it bypasses the need for drug development, thereby reducing the time frame for drug discovery (<xref ref-type="bibr" rid="ref18">Gulati et al., 2021</xref>).</p>
<p>Due to their natural existence, easy availability, extensive antimicrobial, antiviral, antivirulence, and pharmacological properties, the application of plant-derived products and bioactive phytochemicals in modern medicine has gained momentum among researchers (<xref ref-type="bibr" rid="ref6">Chadha et al., 2021</xref>, <xref ref-type="bibr" rid="ref7">2022a</xref>). Numerous <italic>in silico</italic> and pharmacoinformatics-based reports have predicted the effectiveness of phytochemicals and other drugs against SARS-CoV-2, Ebola, and Nipah virus (<xref ref-type="bibr" rid="ref27">Kwofie et al., 2019</xref>; <xref ref-type="bibr" rid="ref36">Pathania et al., 2019</xref>; <xref ref-type="bibr" rid="ref12">Choudhury et al., 2021</xref>; <xref ref-type="bibr" rid="ref25">Kumar A. et al., 2021</xref>). On similar lines, we employed a comprehensive <italic>in silico</italic> approach to identify phytochemicals that could potentially act as inhibitors against poxviruses for repurposing them against the MPXV. In a recent study, Lam et al. demonstrated a high sequence similarity among the structural proteins of viruses belonging to <italic>Poxviridae</italic>, especially between the MPXV and vaccinia virus (<xref ref-type="bibr" rid="ref28">Lam et al., 2022</xref>). Contingent to this, we utilized the structural proteins of the vaccinia virus, which is a close relative of MPXV, as the drug target(s). The first step toward any successful viral infection is cellular entry, which is facilitated by the envelope proteins that are involved in the attachment of the virus to the host cell. Hence, drugs that target the envelope proteins could potentially inhibit viral infection at an early stage. Considering this fact, the phytochemicals included in this study were targeted against the three structural proteins (envelope) of poxviruses <italic>viz</italic>., D13, A26, and H3 (<xref rid="fig1" ref-type="fig">Figure 1</xref>) using molecular docking coupled with dynamic simulations and also evaluated for their pharmacokinetic properties and druglikeliness.</p>
<fig position="float" id="fig1">
<label>Figure 1</label>
<caption>
<p>Crystal structures of envelope proteins of MPXV. <bold>(A)</bold> D13 protein trimer (PDB ID: 6BED), <bold>(B)</bold> A26 protein (PDB ID: 6A9S), <bold>(C)</bold> H3 protein (PDB ID: 5EJ0). The structures were retrieved from the PDB online server (<ext-link xlink:href="http://www.rcsb.org" ext-link-type="uri">www.rcsb.org</ext-link>).</p>
</caption>
<graphic xlink:href="fmicb-13-1073419-g001.tif"/>
</fig>
</sec>
<sec id="sec2" sec-type="materials|methods">
<label>2.</label>
<title>Materials and methods</title>
<sec id="sec3">
<label>2.1.</label>
<title>ADME analysis</title>
<p>The ligands under consideration (phytochemicals) were initially screened to eliminate the compounds violating Lipinski&#x2019;s rule of five (<xref ref-type="bibr" rid="ref2">Benet et al., 2016</xref>). The molecular and physicochemical properties of the phytochemicals were investigated using the SwissADME tool (<xref ref-type="bibr" rid="ref13">Daina et al., 2017</xref>). The phytochemicals in the PDB format were converted to standard Simplified Molecule Input Line Entry Specification (SMILES) format and uploaded to the SwissADME prediction tool.<xref rid="fn0004" ref-type="fn"><sup>1</sup></xref> Further, to computationally predict the drug absorption, distribution, metabolism, excretion and toxicity (ADMET) profile of the phytochemicals, the algorithm-based pharmacokinetics server &#x2013; pkCSM was used.<xref rid="fn0005" ref-type="fn"><sup>2</sup></xref> The server helps in predicting the physiochemical and pharmacological properties of a small molecule. SMILES of the phytochemicals were retrieved from PubChem and uploaded to the pkCSM &#x2013; pharmacokinetics server. Various <italic>in vivo</italic> absorption parameters like Caco2 permeability, intestinal absorption, P-glycoprotein I and II inhibition were predicted. Drug distribution attributes like penetration of blood&#x2013;brain barrier (BBB) and central nervous system (CNS), volume of distribution at steady-state (VDss) along with <italic>in vivo</italic> metabolic variables such as inhibition of cytochrome P450 superfamily enzymes were computed. Drug excretion was studied using parameters like total renal clearance and renal OCT2 substrate. Lastly, the toxigenic properties of the phytochemicals were also computationally predicted using the AMES test, hERG I and II inhibition, hepatotoxicity, skin sensitization, and toxicity toward <italic>Tetrahymena pyriformis</italic> and minnows, etc.</p>
</sec>
<sec id="sec4">
<label>2.2.</label>
<title>Preparation of receptor</title>
<p>The 3D structures of the D13 protein (PDB ID: 6BED, 2.75&#x2009;&#x00C5; resolution), A26 protein (PDB ID: 6A9S, Chain A, 1.18&#x2009;&#x00C5; resolution), and H3 envelope protein (PDB ID: 5EJ0, Chain A, 1.90&#x2009;&#x00C5; resolution) of vaccinia virus were retrieved from the Protein Data Bank (PDB).<xref rid="fn0006" ref-type="fn"><sup>3</sup></xref> All water molecules, heteroatoms, and ligands were removed for the preparation of protein receptor files, followed by the addition of polar hydrogen and Kollman charges. Grid box was generated covering the whole proteins by setting up grid box values as <italic>x</italic>&#x2009;=&#x2009;27.173, <italic>y</italic>&#x2009;=&#x2009;3.424, <italic>z</italic>&#x2009;=&#x2009;33.695 for A26 protein; <italic>x</italic>&#x2009;=&#x2009;105.514, <italic>y</italic>&#x2009;=&#x2009;105.444, <italic>z</italic>&#x2009;=&#x2009;13.868 for D13 protein and <italic>x</italic>&#x2009;=&#x2009;27.579, <italic>y</italic>&#x2009;=&#x2009;33.921, <italic>z</italic>&#x2009;=&#x2009;43.260 for H3 protein. The grid spacing was set to 1&#x2009;&#x00C5;.</p>
</sec>
<sec id="sec5">
<label>2.3.</label>
<title>Preparation of ligands</title>
<p>3D structures of the phytochemicals were downloaded from PubChem in sdf format.<xref rid="fn0007" ref-type="fn"><sup>4</sup></xref> These structures were converted to pdb format using the Open Babel tool. These phytochemicals are known to exhibit antiviral as well as anti-virulent activities (<xref ref-type="bibr" rid="ref23">Kumar et al., 2020</xref>; <xref ref-type="bibr" rid="ref24">Kumar S. et al., 2021</xref>; <xref ref-type="bibr" rid="ref7">Chadha et al., 2022a</xref>). The phytochemicals included in the study are as follows: curcumin (PubChem CID: 969516), E-guggulsterone (PubChem CID: 6439929), licoflavonol (PubChem CID: 5481964), myricetin (PubChem CID: 5281672), oleanolic acid (PubChem CID: 10494), piperidine (PubChem CID: 8082), quercetin (PubChem CID: 5280343), rosmarinic acid (PubChem CID: 5281792), silibinin (PubChem CID: 31553), stigmasterol (PubChem CID: 5280794), ursolic acid (PubChem CID: 64945), withanone (PubChem CID: 21679027). Rifampicin (PubChem CID: 135398735) was used as a positive control for this study.</p>
</sec>
<sec id="sec6">
<label>2.4.</label>
<title>Molecular docking, analysis, and visualization</title>
<p>AutoDock Vina (Version 1.5.7) was used for the docking studies. The receptors and ligands were prepared in the.pdbqt format. Molecular docking was carried out using instructed command prompt, and the exhaustiveness was set to &#x201C;24.&#x201D; Each docking experiment resulted in nine confirmations, and the best binding pose was selected based on the predicted binding energy (in kcal/mol) and molecular interactions. The protein-ligand interactions and PDB structures were visualized and analyzed in PyMOL molecular visualization tool<xref rid="fn0008" ref-type="fn"><sup>5</sup></xref> and LigPlot (<xref ref-type="bibr" rid="ref48">Wallace et al., 1995</xref>).</p>
</sec>
<sec id="sec7">
<label>2.5.</label>
<title>Molecular dynamics simulation</title>
<p>Molecular dynamics (MD) simulations were performed for complexes which demonstrated convincing docking score and protein-ligand interactions. Simulations were performed using the GROMACS software (version 2022.4) with GROMOS96 54a7 force field as described earlier (<xref ref-type="bibr" rid="ref46">Van Der Spoel et al., 2005</xref>; <xref ref-type="bibr" rid="ref40">Schmid et al., 2011</xref>). Ligand topologies was generated using the PRODRG server (<xref ref-type="bibr" rid="ref41">Sch&#x00FC;ttelkopf and van Aalten, 2004</xref>). The docked complex was emplaced in a dodeahedron box with minimum 1&#x2009;nm distance from the box edge. The system was solvated using SPC water model, followed by neutralization through the addition of 7 sodium ions for the H3 protein, 2 chloride ions for the A26 protein, and 34 sodium ions for the D13 protein. This was followed by energy minimzation to eliminate steric clashes or inappropriate geometry using covergence criteria of maximum force &#x003C;1,000&#x2009;kJ/mol/nm. The energy minimized system was subjected to equilibration under NVT and NPT ensembles successively to optimize the temperature and pressure of the system at 300&#x2009;K and 1&#x2009;bar, respectively, using the V-rescale thermostat and Parrinello-Rahman pressure coupling. PME was used to compute long-range electrostatic interactions. Following equilibration, unrestrained final production run was performed for a time duration of 100&#x2009;ns using the leap frog dynamics integrator with a step size of 2&#x2009;fs. The periodic boundary conditions were considered in three dimensions. The analysis of simulation trajectory was done using modules available in GROMACS. Similar methodology was used for simulation of unliganded protein to serve as control.</p>
</sec>
<sec id="sec8">
<label>2.6.</label>
<title>MM/PBSA free energy calculation</title>
<p>The MM/PBSA (Molecular Mechanics Poisson Boltzmann Surface Area) technique was utilized for determining the binding energy along with per residue energy contribution. In MM-PBSA, the polar part of the solvation energy &#x0394;Gpsolv is calculated by solving the Poisson-Boltzmann equation while &#x0394;Gnpsolv is calculated using the linear relation to the solvent accessible surface area. In the present study, g_mmpbsa tool was utilized for the estimation of different components of the binding free energy of complexes (<xref ref-type="bibr" rid="ref26">Kumari et al., 2014</xref>). The last 10&#x2009;ns of the trajectory were used for analysis.</p>
</sec>
<sec id="sec9">
<label>2.7.</label>
<title>Normal mode analysis with iMODS server</title>
<p>Normal mode analysis (NMA) of strongly docked complexes were carried out using iMODS Server<xref rid="fn0009" ref-type="fn"><sup>6</sup></xref> at 300&#x2009;K constant temperature, 1&#x2009;atm constant pressure at molecular mechanics level (<xref ref-type="bibr" rid="ref20">Jos et al., 2014</xref>). iMODS is a computation-based user-friendly server that examines the protein motion within the internal coordinates <italic>via</italic> normal mode analysis. This online server calculates the dihedral coordinates of C&#x03B1; atoms, structural deformability, and eigenvalues of the docked complex.</p>
</sec>
<sec id="sec10">
<label>2.8.</label>
<title>Examining the multi-target affinity, cross-reactivity, and inhibition potential of phytochemicals</title>
<p>The LigTMap server was used to map the alternate targets of the virtually screened lead phytochemicals selected for study.<xref rid="fn0010" ref-type="fn"><sup>7</sup></xref> To evaluate drug target specificity, the server provides a ligand similarity score ranging from 0 to 0.9 (0&#x2009;=&#x2009;lowest, 0.9&#x2009;=&#x2009;highest) along with the list of target proteins (predicted). The inhibition potential of the lead phytochemcials was also determined in terms of inhibition contant (K<italic><sub>i</sub></italic>) from the binding emergy (&#x2206;G), using the formula: K<italic><sub>i</sub></italic>&#x2009;=&#x2009;(&#x2206;G/RT), where R is the universal gas constant (1.985&#x2009;&#x00D7;&#x2009;10<sup>&#x2013;3</sup>&#x2009;kcal/mol-K) and T is the temperature (298.15&#x2009;K), as previously described (<xref ref-type="bibr" rid="ref35">Ortiz et al., 2019</xref>).</p>
</sec>
</sec>
<sec id="sec11" sec-type="results">
<label>3.</label>
<title>Results</title>
<sec id="sec12">
<label>3.1.</label>
<title>Phytochemical screening and ADME analysis</title>
<p>After undertaking an exhaustive literature survey, we prepared a library of bioactive phytochemicals known to harbor antiviral properties that could possibly be used for inhibiting poxviruses. A total of 12 phytochemicals were shortlisted for the study and further subjected to ADME analysis. All 12 phytochemicals were tested for Lipinski&#x2019;s rule of five to examine their drug likeliness and physiochemical properties. The rule outlines molecular properties that are important for the pharmacokinetic properties of a drug in the human body, including its absorption, distribution, metabolism, and excretion. These parameters require molecular weight to be &#x003C;500&#x2009;Da, the number of hydrogen bond acceptors &#x003C;10, hydrogen bond donors to be &#x2264;5, and log P, i.e., octanol water coefficient to be &#x003C;5 (<xref ref-type="bibr" rid="ref31">Lipinski, 2004</xref>). This rule is important in determining the suitability of the compound to be used as an active drug for oral consumption in humans (<xref ref-type="bibr" rid="ref31">Lipinski, 2004</xref>). It has been well documented that molecules that fail to follow Lipinski&#x2019;s rule of five have poor absorption and bioavailability in the body. All the test phytochemicals adhered to Lipinski&#x2019;s rule of five with no more than 2 violations (at max; <xref rid="tab1" ref-type="table">Table 1</xref>) and displayed ideal drug-like properties, with no toxicity or adverse effects on cellular metabolism (<xref ref-type="supplementary-material" rid="SM10">Supplementary Tables S1</xref>, <xref ref-type="supplementary-material" rid="SM10">S2</xref>), thereby being suitable for consumption.</p>
<table-wrap position="float" id="tab1">
<label>Table 1</label>
<caption>
<p>Lipinski rule of five values for the test phytochemicals.</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="top">S. No.</th>
<th align="left" valign="top">Name of phytochemical</th>
<th align="center" valign="top">PubChem ID</th>
<th align="center" valign="top">Molecular weight (&#x003C;&#x2009;500&#x2009;Da)</th>
<th align="center" valign="top">Lipophilicity (LogP&#x2009;&#x003C;&#x2009;5)</th>
<th align="center" valign="top">#H-bond acceptors (&#x003C;&#x2009;10)</th>
<th align="center" valign="top">#H-bond donors (&#x003C;&#x2009;5)</th>
<th align="center" valign="top">No. of violations</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="top">1</td>
<td align="left" valign="top">Curcumin</td>
<td align="center" valign="top">969,516</td>
<td align="left" valign="top">368.38</td>
<td align="left" valign="top">1.47</td>
<td align="center" valign="top">6</td>
<td align="center" valign="top">2</td>
<td align="center" valign="top">0</td>
</tr>
<tr>
<td align="left" valign="top">2</td>
<td align="left" valign="top">E-Guggulsterone</td>
<td align="center" valign="top">6,439,929</td>
<td align="left" valign="top">312.45</td>
<td align="left" valign="top">3.86</td>
<td align="center" valign="top">2</td>
<td align="center" valign="top">0</td>
<td align="center" valign="top">0</td>
</tr>
<tr>
<td align="left" valign="top">3</td>
<td align="left" valign="top">Licoflavone A</td>
<td align="center" valign="top">5,481,964</td>
<td align="left" valign="top">322.35</td>
<td align="left" valign="top">2.2</td>
<td align="center" valign="top">4</td>
<td align="center" valign="top">2</td>
<td align="center" valign="top">0</td>
</tr>
<tr>
<td align="left" valign="top">4</td>
<td align="left" valign="top">Myricitrin</td>
<td align="center" valign="top">5,281,672</td>
<td align="left" valign="top">464.38</td>
<td align="left" valign="top">&#x2212;2.32</td>
<td align="center" valign="top">12</td>
<td align="center" valign="top">8</td>
<td align="center" valign="top">2</td>
</tr>
<tr>
<td align="left" valign="top">5</td>
<td align="left" valign="top">Oleanolic acid</td>
<td align="center" valign="top">10,494</td>
<td align="left" valign="top">456.7</td>
<td align="left" valign="top">5.82</td>
<td align="center" valign="top">3</td>
<td align="center" valign="top">2</td>
<td align="center" valign="top">1</td>
</tr>
<tr>
<td align="left" valign="top">6</td>
<td align="left" valign="top">Piperine</td>
<td align="center" valign="top">8,082</td>
<td align="left" valign="top">285.34</td>
<td align="left" valign="top">2.39</td>
<td align="center" valign="top">3</td>
<td align="center" valign="top">0</td>
<td align="center" valign="top">0</td>
</tr>
<tr>
<td align="left" valign="top">7</td>
<td align="left" valign="top">Quercetin</td>
<td align="center" valign="top">5,280,343</td>
<td align="left" valign="top">302.24</td>
<td align="left" valign="top">&#x2212;0.56</td>
<td align="center" valign="top">7</td>
<td align="center" valign="top">5</td>
<td align="center" valign="top">0</td>
</tr>
<tr>
<td align="left" valign="top">8</td>
<td align="left" valign="top">Rosmarinic acid</td>
<td align="center" valign="top">5,281,792</td>
<td align="left" valign="top">360.31</td>
<td align="left" valign="top">0.9</td>
<td align="center" valign="top">8</td>
<td align="center" valign="top">5</td>
<td align="center" valign="top">0</td>
</tr>
<tr>
<td align="left" valign="top">9</td>
<td align="left" valign="top">Silibinin</td>
<td align="center" valign="top">31,553</td>
<td align="left" valign="top">482.44</td>
<td align="left" valign="top">&#x2212;0.4</td>
<td align="center" valign="top">10</td>
<td align="center" valign="top">5</td>
<td align="center" valign="top">0</td>
</tr>
<tr>
<td align="left" valign="top">10</td>
<td align="left" valign="top">Stigmasterol</td>
<td align="center" valign="top">5,280,794</td>
<td align="left" valign="top">412.69</td>
<td align="left" valign="top">6.62</td>
<td align="center" valign="top">1</td>
<td align="center" valign="top">1</td>
<td align="center" valign="top">1</td>
</tr>
<tr>
<td align="left" valign="top">11</td>
<td align="left" valign="top">Ursolic acid</td>
<td align="center" valign="top">64,945</td>
<td align="left" valign="top">456.7</td>
<td align="left" valign="top">5.82</td>
<td align="center" valign="top">3</td>
<td align="center" valign="top">2</td>
<td align="center" valign="top">1</td>
</tr>
<tr>
<td align="left" valign="top">12</td>
<td align="left" valign="top">Withanone</td>
<td align="center" valign="top">21,679,027</td>
<td align="left" valign="top">470.6</td>
<td align="left" valign="top">2.75</td>
<td align="center" valign="top">6</td>
<td align="center" valign="top">2</td>
<td align="center" valign="top">0</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p>All the compounds were found to obey Lipinski&#x2019;s rules.</p>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="sec13">
<label>3.2.</label>
<title>Phytochemicals exhibit high binding affinity toward envelope proteins of poxvirus</title>
<p>Molecular docking was performed to computationally predict the interactions between the phytochemicals and envelope proteins of the vaccinia virus. AutoDock Vina tool was employed to analyze the binding affinities, poses, and confirmation of the ligands with respect to the target receptor. Each docking resulted in 9 different poses for each phytochemical, yielding different binding energies. The pose with the least binding energy (indicating the strongest interaction) was taken into consideration for further analysis. Generally, binding energy lower than the upper threshold of &#x2212;6.0&#x2009;kcal/mol is regarded as the cut-off value for determining the robust interaction between a test ligand and protein. All 12 phytochemicals exhibited excellent binding scores with values below &#x2212;6.0&#x2009;kcal/mol (<xref rid="fig2" ref-type="fig">Figure 2</xref>). Some phytochemicals such as oleanolic acid, ursolic acid, withanone, silibinin, and stigmasterol displayed strong affinities with binding energy even lower than &#x2212;7.5&#x2009;kcal/mol against all the 3 protein receptors (<xref rid="fig2" ref-type="fig">Figure 2</xref>). The binding energies for most of the phytochemicals were found to be either similar or lower than that of rifampicin, which was used as a control. The results for each envelope protein have been summarized in detail in <xref ref-type="supplementary-material" rid="SM10">Supplementary Tables S3</xref>&#x2013;<xref ref-type="supplementary-material" rid="SM10">S5</xref> and also discussed in the subsequent sections.</p>
<fig position="float" id="fig2">
<label>Figure 2</label>
<caption>
<p>Bar graph depicting the binding energies of different phytochemicals docked against D13, A26, and H3 envelope proteins of MPXV. Rifampicin was used as a control.</p>
</caption>
<graphic xlink:href="fmicb-13-1073419-g002.tif"/>
</fig>
<sec id="sec14">
<label>3.2.1.</label>
<title>Molecular interactions between D13 protein and phytochemicals</title>
<p>Since D13 is a viral protein that provides rigidity to the membrane of poxviruses and plays a critical role in viral morphogenesis, it was selected as a drug target for the study. The phytochemicals under investigation demonstrated high binding affinity toward D13 through a plethora of interactions, as compared to rifampicin (<xref rid="fig3" ref-type="fig">Figure 3</xref> and <xref ref-type="supplementary-material" rid="SM1">Supplementary Figure S1</xref>). The amino acid residues on the D13 protein found to be interacting with the phytochemicals have been enlisted in <xref ref-type="supplementary-material" rid="SM10">Supplementary Table S3</xref>. Among these, the majority of phytochemicals showed more hydrophobic associations as compared to hydrogen bonds. However, with silibinin and myricetin, the number of hydrogen bonds outnumbered the hydrophobic interactions (<xref rid="fig3" ref-type="fig">Figure 3</xref>). Specific residues of the D13 protein, such as Ser152, Thr153, Ser256, and Asn435, significantly facilitated interactions with the phytochemicals. Out of the 12 phytochemicals, oleanolic acid, ursolic acid, and silibinin showed the highest affinity exhibiting binding energies lower than &#x2212;9.0&#x2009;kcal/mol with multiple hydrogen bonds and hydrophobic interactions. Briefly, oleanolic acid formed hydrogen bonds with His128, Thr153, and Asn472, and hydrophobic interactions with Lys127, Val154, Tyr155, Glu280, Lys436, Lys434, Arg461, and Ile462. In contrast, ursolic acid showed hydrogen bonding with His128 and Thr153 and hydrophobic interactions with Lys127, Val134, Tyr155, Glu280, Lys434, Lys436, Arg461, and Ileu462. Similarly, silibinin also formed numerous hydrogen bonds with Lys127, Ser152, Thr153, Ile440, Ser441, Asn456, Arg461, and Glu465 along with five hydrophobic interactions with Glu74, His128, Asp438, Gly457, and Pro458. Interestingly, the binding positions of these phytochemicals were similar to that of rifampicin as they occupy the central channel of the D13 trimer. Therefore, it can be speculated that the binding of these phytochemicals at this position may interrupt the assembly of membrane proteins mediated by D13, thereby outlining the potential of phytochemicals as perspective candidates against poxviruses.</p>
<fig position="float" id="fig3">
<label>Figure 3</label>
<caption>
<p>Molecular docking of silibinin, oleanolic acid, and ursolic acid against the D13 protein (PDB ID: 6BED). Left Panel: Visualization of the binding position of the ligand within the protein cavity. Central Panel: Ribbon diagram showing the docked complex. Ligands have been shown in red. Right Panel: 2D diagram representing the interactions between the phytochemicals and D13 protein generated by LigPlot. Dashed lines show hydrogen bonds while the red arcs represent hydrophobic interactions. Rifampicin has been used as a control. Docking was performed <italic>via</italic> AutoDock Vina (Version 1.5.7) and the 3D structures were visualized using PyMol tool.</p>
</caption>
<graphic xlink:href="fmicb-13-1073419-g003.tif"/>
</fig>
</sec>
<sec id="sec15">
<label>3.2.2.</label>
<title>Molecular interactions between A26 protein and phytochemicals</title>
<p>The A26 protein is a viral envelope protein known to attach to glycosaminoglycans and extracellular matrix laminin on the host cell surface, thereby promoting endocytosis in the vaccinia virus (<xref ref-type="bibr" rid="ref5">Carter et al., 2005</xref>; <xref ref-type="bibr" rid="ref11">Chiu et al., 2007</xref>). Hence, its direct inhibition would prevent the entry of viruses into the host cell. To exploit A26 as a drug target, its crystal structure was retrieved from PDB, and blind docking was carried out wherein the docking region covered the entire surface of the protein. Results revealed that all the phytochemicals showed exceptionally low binding scores with values below &#x2212;7.0&#x2009;kcal/mol, except rosmarinic acid. Asp267, Thr283, Tyr289, and Asn340 were predicted to be the key amino acid residues involved in mediating interactions with the phytochemicals (<xref rid="fig4" ref-type="fig">Figure 4</xref> and <xref ref-type="supplementary-material" rid="SM2">Supplementary Figure S2</xref>). The amino acid residues on the A26 protein shown to be interacting with the phytochemicals have been summarized in <xref ref-type="supplementary-material" rid="SM10">Supplementary Table S4</xref>. The lead compounds that demonstrated efficient binding were stigmasterol (&#x2212;8.0&#x2009;kcal/mol), oleanolic acid (&#x2212;7.9&#x2009;kcal/mol), ursolic acid (&#x2212;7.8&#x2009;kcal/mol), silibinin (&#x2212;7.8&#x2009;kcal/mol), and withanone (&#x2212;7.8&#x2009;kcal/mol). Stigmasterol was bonded to the A26 protein only <italic>via</italic> hydrophobic bonds, including interactions with the Tyr113, His135, Asp136, Tyr160, and Lys163, Gln164, His167, and Ile 202 residues. On the other hand, ursolic acid formed hydrogen bonds with Arg57, Phe266, and Asp267, and hydrophobic interactions with Lys61, Lys301, Val302, Thr303, Glu304, and Lys306 (<xref rid="fig4" ref-type="fig">Figure 4</xref>). Similarly, silibinin formed five hydrogen bonds with Thr283, Asp310, Arg333, Asn340, and Thr355, along with hydrophobic interactions with Arg46, Asn284, Asn287, Asn334, Asp339, Glu341, and Thr353. Withanone also formed five hydrogen bonds with Thr208, Asp310, Lys311, Asp339, and Glu341, and hydrophobic interactions with Asn287, Tyr289, and Asn340 residues. Apart from these phytochemicals, oleanolic acid showed the best binding pose as it formed one hydrogen bond with Asp267 and multiple hydrophobic bonds with Glu54, Arg57, Asp58, Lys61, Tyr242, Phe266, Val301, Thr303, Glu304, Asp305, and Lys306, residues present in the acid-sensing region at the N-terminal region of A26 protein (<xref rid="fig4" ref-type="fig">Figure 4</xref>). This region plays a quintessential role in the infectivity of the vaccinia virus (<xref ref-type="bibr" rid="ref10">Chang et al., 2019</xref>). Based on the docking position of oleanolic acid, it can be hypothesized that it may act as an inhibitor of the A26 protein due to excellent fitting into the N-terminal domain, thereby lowering viral infectivity.</p>
<fig position="float" id="fig4">
<label>Figure 4</label>
<caption>
<p>Molecular docking of silibinin, oleanolic acid, and ursolic acid against the A26 protein (PDB ID: 6A9S). Left Panel: Visualization of the binding position of the ligand within the protein cavity. Central Panel: Ribbon diagram showing the docked complex. Ligands are shown in red. Right Panel: 2D diagram representing the interactions between phytochemicals and A26 protein generated by LigPlot. Dashed lines depict hydrogen bonds and the red arcs represent hydrophobic interactions. Rifampicin has been used as a control. Docking was performed <italic>via</italic> AutoDock Vina (Version 1.5.7) and the 3D structures were visualized using PyMol tool.</p>
</caption>
<graphic xlink:href="fmicb-13-1073419-g004.tif"/>
</fig>
</sec>
<sec id="sec16">
<label>3.2.3.</label>
<title>Molecular interactions between H3 protein and phytochemicals</title>
<p>The H3 protein is another important envelope protein of poxviruses that binds to heparin sulfate on the host cell surface and assists in the virus attachment (<xref ref-type="bibr" rid="ref43">Singh et al., 2016</xref>). Moreover, it has gained attention among scientists as a potent drug and vaccine target. Hence, we docked bioactive phytochemicals against the H3 viral protein. Our findings revealed that interactions between phytochemicals and the H3 protein were predominantly stabilized by hydrophobic interactions, similar to the case of D13 and A26 proteins (<xref rid="fig5" ref-type="fig">Figure 5</xref> and <xref ref-type="supplementary-material" rid="SM3">Supplementary Figure S3</xref>). All the predicted interactions between the H3 protein and phytochemicals have been mentioned in <xref ref-type="supplementary-material" rid="SM10">Supplementary Table S5</xref>. Thr6, Thr91, Thr94, and Val183 were identified as the crucial amino residues facilitating the molecular interactions. Among the 12 phytochemicals, E-guggulsterone, oleanolic acid, piperidine, quercitin, silibinin, stigmasterol, ursolic acid, and withanone displayed excellent binding energies with values below &#x2212;7.0&#x2009;kcal/mol. However, the top three compounds with a maximum binding affinity toward the H3 protein were ursolic acid (&#x2212;8.0&#x2009;kcal/mol), oleanolic acid (&#x2212;7.9&#x2009;kcal/mol), and withanone (&#x2212;7.2&#x2009;kcal/mol). Ursolic acid formed one hydrogen bond with Tyr181 and multiple hydrophobic bonds with Cys90, Thr94, Asp182, Val183, Gly211, Phe212, Tyr213, and Phe214 (<xref rid="fig5" ref-type="fig">Figure 5</xref>). On similar lines, oleanolic acid formed a lone hydrogen bond with Val183 and numerous hydrophobic bonds with Cys90, Thr94, Ile156, Tyr181, Asp182, Ser209, Ser210, Glu211, Phe212, Tyr213, and Phe214. Withanone also interacted with Thr91 and Thr94 <italic>via</italic> hydrogen bonds and displayed hydrophobic interactions with Cys99, Gly211, Phe212, Tyr213, and Phe214 residues (<xref rid="fig5" ref-type="fig">Figure 5</xref>). Taken together, our results reveal that bioactive phytochemicals can bind effectively to the H3 envelope protein of poxviruses.</p>
<fig position="float" id="fig5">
<label>Figure 5</label>
<caption>
<p>Molecular docking of silibinin, oleanolic acid, and ursolic acid against the H3 protein (PDB ID: 5EJ0). Left Panel: Visualization of the binding position of the ligand within the protein cavity. Central Panel: Ribbon diagram showing the docked complex. Ligands are shown in red. Right Panel: 2D diagram representing the interactions between phytochemicals and H3 protein generated by LigPlot. Dashed lines illustrate hydrogen bonding and the red arcs represent hydrophobic interactions. Rifampicin has been used as a control. Docking was performed <italic>via</italic> AutoDock Vina (Version 1.5.7) and the 3D structures were visualized using PyMol tool.</p>
</caption>
<graphic xlink:href="fmicb-13-1073419-g005.tif"/>
</fig>
</sec>
</sec>
<sec id="sec17">
<label>3.3.</label>
<title>Molecular dynamics simulations</title>
<p>MDsimulations were performed in order to assess the dynamic interactions of the selected phytomolecules with the three proteins under aqueous conditions. Structural deviations during the simulations were described in terms of root mean square deviation (RMSD) of atomic coordinates, root mean square fluctuations (RMSF), and radius of gyration of the protein. The backbone RMSD of the proteins in presence of the three phytomolecules was very close to unliganded protein which suggests that binding of the phytomolecules did not impart any major structural change in the envelope proteins (<xref rid="fig6" ref-type="fig">Figure 6</xref>). Differences were observed in case of ligand RMSD values although all the protein-ligand complexes were stable throughout the 100&#x2009;ns simulation (as depicted in <xref rid="fig7" ref-type="fig">Figure 7</xref>). The radius of gyration values were similar for all three protein in presence of the phytomolecules, except in the case of D13 protein trimer, where ligand binding lead to an increase in the Rg values (<xref rid="fig8" ref-type="fig">Figure 8</xref>). This implies that binding of the phytomolecules reduces the compactness of the trimer. RMSF values in presence of phytomolecules were consistant with the control with variations observed in solvent exposed loop regions (<xref rid="fig9" ref-type="fig">Figures 9</xref>, <xref rid="fig10" ref-type="fig">10</xref>). The binding of oleanolic acid to the H3 envelope protein was very stable throughtout 100&#x2009;ns (<xref rid="fig9" ref-type="fig">Figure 9B</xref>). In case of the A26 protein, the binding of oleanolic acid was less stable as it was displaced from it&#x2019;s docked pose, which was stabilized after 80&#x2009;ns (<xref rid="fig9" ref-type="fig">Figure 9A</xref>). It also showed stable interactions with the hydrophobic cavity between the two monomers of the trimeric D13 (<xref rid="fig10" ref-type="fig">Figure 10</xref>). Although stable hydrogen bonds were not observed, the binding was maintained predominantly through hydrophobic interactions. Similar binding modes were observed for ursolic acid with the three receptors as both are pentacyclic triterpenoids. Ursolic acid showed best binding with the trimeric D13 protein. Silibinin being a flavonolignan, has distinct functional groups compared to oleanoic acid and ursolic acid and thus exhibited more diverse interaction capabilities with the three proteins. It showed stable binding with all three proteins by forming hydrogen bonds, &#x03C0;-stacking interactions and other hydrophobic interactions.</p>
<fig position="float" id="fig6">
<label>Figure 6</label>
<caption>
<p>RMSD profile of protein backbone bound to phytomolecules in case of <bold>(A)</bold> A26 protein <bold>(B)</bold> H3 protein and <bold>(C)</bold> D13 protein.</p>
</caption>
<graphic xlink:href="fmicb-13-1073419-g006.tif"/>
</fig>
<fig position="float" id="fig7">
<label>Figure 7</label>
<caption>
<p>RMSD profile of phytomolecules bound to <bold>(A)</bold> A26 protein <bold>(B)</bold> H3 protein and <bold>(C)</bold> D13 protein.</p>
</caption>
<graphic xlink:href="fmicb-13-1073419-g007.tif"/>
</fig>
<fig position="float" id="fig8">
<label>Figure 8</label>
<caption>
<p>Radius of gyration (Rg) for the proteins bound to phytomolecules <bold>(A)</bold> A26 protein, <bold>(B)</bold> H3 envelop protein, and <bold>(C)</bold> D13 protein.</p>
</caption>
<graphic xlink:href="fmicb-13-1073419-g008.tif"/>
</fig>
<fig position="float" id="fig9">
<label>Figure 9</label>
<caption>
<p>RMSF profile of <bold>(A)</bold> A26 protein and <bold>(B)</bold> H3 protein.</p>
</caption>
<graphic xlink:href="fmicb-13-1073419-g009.tif"/>
</fig>
<fig position="float" id="fig10">
<label>Figure 10</label>
<caption>
<p>RMSF profile of the three monomers of the D13 protein respresented as <bold>(A&#x2013;C)</bold>.</p>
</caption>
<graphic xlink:href="fmicb-13-1073419-g010.tif"/>
</fig>
</sec>
<sec id="sec18">
<label>3.4.</label>
<title>MM/PBSA free energy calculation</title>
<p>To quantify the effect of phytomolecules on the activity of the MPX proteins, the free energy was calculated by MM-PBSA analysis (<xref rid="tab2" ref-type="table">Tables 2</xref>&#x2013;<xref rid="tab4" ref-type="table">4</xref>). The variation of &#x0394;Gbind clearly showed that, among all the test molecules, silibinin demonstrates lowest binding energy values. The difference in &#x0394;Gbind majorly comes from the &#x0394;Evdw (van der Waals) interaction component and &#x0394;Gpsolv (polar part of solvation) free energy components, playing an important role in lower binding free energy. The residue-wise contribution of free energy is depicted in <xref rid="fig11" ref-type="fig">Figure 11</xref>, which highlights the interacting residues of the phytomolecules.</p>
<table-wrap position="float" id="tab2">
<label>Table 2</label>
<caption>
<p>Summary of energy components from MM/PBSA analysis of the lead phytomolecules with the D13 protein.</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="top">Name of phytochemical</th>
<th align="center" valign="top">&#x0394;E<sub>vdw</sub></th>
<th align="center" valign="top">+/&#x2212;</th>
<th align="center" valign="top">&#x0394;E<sub>EEL</sub></th>
<th align="center" valign="top">+/&#x2212;</th>
<th align="center" valign="top">&#x0394;G<sub>psolv</sub></th>
<th align="center" valign="top">+/&#x2212;</th>
<th align="center" valign="top">&#x0394;G<sub>npsolv</sub></th>
<th align="center" valign="top">+/&#x2212;</th>
<th align="center" valign="top">&#x0394;G<sub>bind</sub></th>
<th align="center" valign="top">+/&#x2212;</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="top">Silibinin</td>
<td align="center" valign="top">&#x2212;178.798</td>
<td align="center" valign="top">16.682</td>
<td align="center" valign="top">&#x2212;7.729</td>
<td align="center" valign="top">5.416</td>
<td align="center" valign="top">65.403</td>
<td align="center" valign="top">20.664</td>
<td align="center" valign="top">&#x2212;18.219</td>
<td align="center" valign="top">1.603</td>
<td align="center" valign="top">&#x2212;139.342</td>
<td align="center" valign="top">23.111</td>
</tr>
<tr>
<td align="left" valign="top">Oleanolic acid</td>
<td align="center" valign="top">&#x2212;169.667</td>
<td align="center" valign="top">9.077</td>
<td align="center" valign="top">&#x2212;6.164</td>
<td align="center" valign="top">4.028</td>
<td align="center" valign="top">93.247</td>
<td align="center" valign="top">11.445</td>
<td align="center" valign="top">&#x2212;17.162</td>
<td align="center" valign="top">1.487</td>
<td align="center" valign="top">&#x2212;99.746</td>
<td align="center" valign="top">12.752</td>
</tr>
<tr>
<td align="left" valign="top">Ursolic acid</td>
<td align="center" valign="top">&#x2212;200.753</td>
<td align="center" valign="top">8.962</td>
<td align="center" valign="top">&#x2212;4.753</td>
<td align="center" valign="top">4.489</td>
<td align="center" valign="top">114.744</td>
<td align="center" valign="top">17.493</td>
<td align="center" valign="top">&#x2212;17.970</td>
<td align="center" valign="top">1.078</td>
<td align="center" valign="top">&#x2212;108.732</td>
<td align="center" valign="top">16.563</td>
</tr>
</tbody>
</table>
</table-wrap>
<table-wrap position="float" id="tab3">
<label>Table 3</label>
<caption>
<p>Summary of energy components from MM/PBSA analysis of the lead phytomolecules with the A26 protein.</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="top">Name of phytochemical</th>
<th align="center" valign="top">&#x0394;E<sub>vdw</sub></th>
<th align="center" valign="top">+/&#x2212;</th>
<th align="center" valign="top">&#x0394;E<sub>EEL</sub></th>
<th align="center" valign="top">+/&#x2212;</th>
<th align="center" valign="top">&#x0394;G<sub>psolv</sub></th>
<th align="center" valign="top">+/&#x2212;</th>
<th align="center" valign="top">&#x0394;G<sub>npsolv</sub></th>
<th align="center" valign="top">+/&#x2212;</th>
<th align="center" valign="top">&#x0394;G<sub>bind</sub></th>
<th align="center" valign="top">+/&#x2212;</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="top">Silibinin</td>
<td align="center" valign="top">&#x2212;207.641</td>
<td align="center" valign="top">13.822</td>
<td align="center" valign="top">&#x2212;4.594</td>
<td align="center" valign="top">4.522</td>
<td align="center" valign="top">91.502</td>
<td align="center" valign="top">13.960</td>
<td align="center" valign="top">&#x2212;19.286</td>
<td align="center" valign="top">1.033</td>
<td align="center" valign="top">&#x2212;140.019</td>
<td align="center" valign="top">16.538</td>
</tr>
<tr>
<td align="left" valign="top">Oleanolic acid</td>
<td align="center" valign="top">&#x2212;98.855</td>
<td align="center" valign="top">18.504</td>
<td align="center" valign="top">&#x2212;0.857</td>
<td align="center" valign="top">2.929</td>
<td align="center" valign="top">33.213</td>
<td align="center" valign="top">12.848</td>
<td align="center" valign="top">&#x2212;10.255</td>
<td align="center" valign="top">2.142</td>
<td align="center" valign="top">&#x2212;76.754</td>
<td align="center" valign="top">13.375</td>
</tr>
<tr>
<td align="left" valign="top">Ursolic acid</td>
<td align="center" valign="top">&#x2212;124.155</td>
<td align="center" valign="top">7.587</td>
<td align="center" valign="top">&#x2212;1.594</td>
<td align="center" valign="top">2.501</td>
<td align="center" valign="top">38.955</td>
<td align="center" valign="top">6.111</td>
<td align="center" valign="top">&#x2212;11.191</td>
<td align="center" valign="top">1.060</td>
<td align="center" valign="top">&#x2212;97.985</td>
<td align="center" valign="top">8.263</td>
</tr>
</tbody>
</table>
</table-wrap>
<table-wrap position="float" id="tab4">
<label>Table 4</label>
<caption>
<p>Summary of energy components from MM/PBSA analysis of the lead phytomolecules with the H3 protein.</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="top">Name of phytochemical</th>
<th align="left" valign="top">&#x0394;E<sub>vdw</sub></th>
<th align="center" valign="top">+/&#x2212;</th>
<th align="left" valign="top">&#x0394;E<sub>EEL</sub></th>
<th align="center" valign="top">+/&#x2212;</th>
<th align="center" valign="top">&#x0394;G<sub>psolv</sub></th>
<th align="center" valign="top">+/&#x2212;</th>
<th align="center" valign="top">&#x0394;G<sub>npsolv</sub></th>
<th align="center" valign="top">+/&#x2212;</th>
<th align="center" valign="top">&#x0394;G<sub>bind</sub></th>
<th align="center" valign="top">+/&#x2212;</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="top">Silibinin</td>
<td align="left" valign="top">&#x2212;196.655</td>
<td align="center" valign="top">14.010</td>
<td align="left" valign="top">1.587</td>
<td align="center" valign="top">2.953</td>
<td align="center" valign="top">48.207</td>
<td align="center" valign="top">7.487</td>
<td align="center" valign="top">&#x2212;18.467</td>
<td align="center" valign="top">1.250</td>
<td align="center" valign="top">&#x2212;165.328</td>
<td align="center" valign="top">14.763</td>
</tr>
<tr>
<td align="left" valign="top">Oleanolic acid</td>
<td align="left" valign="top">&#x2212;129.455</td>
<td align="center" valign="top">11.025</td>
<td align="left" valign="top">&#x2212;2.281</td>
<td align="center" valign="top">2.629</td>
<td align="center" valign="top">28.694</td>
<td align="center" valign="top">6.546</td>
<td align="center" valign="top">&#x2212;11.861</td>
<td align="center" valign="top">1.366</td>
<td align="center" valign="top">&#x2212;114.902</td>
<td align="center" valign="top">10.091</td>
</tr>
<tr>
<td align="left" valign="top">Ursolic acid</td>
<td align="left" valign="top">&#x2212;111.655</td>
<td align="center" valign="top">12.511</td>
<td align="left" valign="top">0.013</td>
<td align="center" valign="top">3.350</td>
<td align="center" valign="top">27.926</td>
<td align="center" valign="top">10.039</td>
<td align="center" valign="top">&#x2212;11.217</td>
<td align="center" valign="top">1.418</td>
<td align="center" valign="top">&#x2212;94.933</td>
<td align="center" valign="top">11.705</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig position="float" id="fig11">
<label>Figure 11</label>
<caption>
<p>Free energy contribution of <bold>(A)</bold> H3 residues, <bold>(B)</bold> A26 residues, <bold>(C)</bold> D13 residues (monomer-1), and <bold>(D)</bold> D13 residues (monomer-2) interacting with the three selected phytomolecules. Major residues with negative free energy contribution are highlighted.</p>
</caption>
<graphic xlink:href="fmicb-13-1073419-g011.tif"/>
</fig>
</sec>
<sec id="sec19">
<label>3.5.</label>
<title>Insights from NMA analysis confirms the stability of phytochemical-protein complexes</title>
<p>Based on the results obtained from molecular docking experiments, we shortlisted three lead phytochemicals <italic>viz</italic>. silibinin, oleanolic acid, and ursolic acid, which demonstrated strong binding affinity toward all the three envelope proteins of poxviruses in common. Subsequently, the stability of the docked phytochemical-protein complexes (best binding pose considered) was examined using MD simulations. The main chain deformability for the complexes has been depicted in the first vertical columns of <xref ref-type="supplementary-material" rid="SM4">Supplementary Figure S4</xref> (D13 protein), <xref ref-type="supplementary-material" rid="SM5">Supplementary Figure S5</xref> (A26 protein), and <xref ref-type="supplementary-material" rid="SM6">Supplementary Figure S6</xref> (H3 protein). Deformability represents the capacity of the given molecule to distort at each residue. Hinges within the plot represent high deformability regions in the protein chain. The second vertical columns of <xref ref-type="supplementary-material" rid="SM4">Supplementary Figure S4</xref> (D13 protein), <xref ref-type="supplementary-material" rid="SM5">Supplementary Figure S5</xref> (A26 protein), and <xref ref-type="supplementary-material" rid="SM6">Supplementary Figure S6</xref> (H3 protein) illustrate the eigenvalues of the docked complexes. Silibinin, oleanolic acid, and ursolic acid demonstrated eigenvalues of 1.825063&#x2009;&#x00D7;&#x2009;10<sup>&#x2212;4</sup>, 1.843076&#x2009;&#x00D7;&#x2009;10<sup>&#x2212;4</sup>, and 1.812478&#x2009;&#x00D7;&#x2009;10<sup>&#x2212;4</sup> with D13 protein; 6.057804&#x2009;&#x00D7;&#x2009;10<sup>&#x2212;4</sup>, 6.381300&#x2009;&#x00D7;&#x2009;10<sup>&#x2212;4</sup> and 6.346715&#x2009;&#x00D7;&#x2009;10<sup>&#x2212;4</sup> with A26 protein; 7.358898&#x2009;&#x00D7;&#x2009;10<sup>&#x2212;4</sup>, 7.377467&#x2009;&#x00D7;&#x2009;10<sup>&#x2212;4</sup> and 7.618457&#x2009;&#x00D7;&#x2009;10<sup>&#x2212;4</sup> with H3 protein, respectively. The eigenvalue is a direct measure of the motion stiffness of C&#x03B1; atoms. The low eigenvalues of the silibinin, oleanolic acid, and ursolic acid with the three structural proteins of poxviruses indicated high stability of the docked complexes. Moreover, the co-variance maps of the complexes (first vertical panel: <xref ref-type="supplementary-material" rid="SM7">Supplementary Figures S7</xref>&#x2013;<xref ref-type="supplementary-material" rid="SM9">S9</xref>), which represents the correlated motion between a pair of residues (red color dots), uncorrelated motion (white color dots), and anti-correlated motion (blue color), also confirmed complex stability and stiffness. This was further substantiated by the elastic network models (second vertical panel: <xref ref-type="supplementary-material" rid="SM7">Supplementary Figures S7</xref>&#x2013;<xref ref-type="supplementary-material" rid="SM9">S9</xref>) that compute the normal modes using a linking matrix to depict the pair of atoms that are connected by springs. The intensity of the gray dots represents the stiffness between the interacting regions. In summary, our findings reveal that the phytochemicals associate strongly with the envelope proteins and demonstrate high stability, yielding satisfactory results from the MD simulations.</p>
</sec>
<sec id="sec20">
<label>3.6.</label>
<title>Specificity of the phytochemicals to human protein targets</title>
<p>The LigTMap server screened some potential targets for the tested phytochemicals associated with cellular pathway. LigTMap takes into account the chemical configuration of the ligand and the target, pharmacophore properties drug binding properties. For all the lead phytochemicals, the consensus ligand similarity scores were found to be lower than 0.5 (<xref ref-type="supplementary-material" rid="SM10">Supplementary Tables S6</xref>&#x2013;<xref ref-type="supplementary-material" rid="SM10">S8</xref>), thus indicating their target specificity and lower risk of incurring any off-target effects.</p>
</sec>
</sec>
<sec id="sec21" sec-type="discussions">
<label>4.</label>
<title>Discussion</title>
<p>With the COVID-19 pandemic slowing down, humankind faces another potential threat due to the sudden emergence of MPXV. The multi-country MPX outbreak may not match the scale of the receeding pandemic but demands immediate global attention and surveillance. With the discontinuation of smallpox vaccine (1980) that provided cross-protection against the MPXV, a large fraction of the world population remains susceptible to MPX (<xref ref-type="bibr" rid="ref8">Chadha et al., 2022b</xref>). Moreover, the situation has been aggravated by several cluster outbreaks being reported across non-endemic countries and the identification of 50 genomic substitutions within the 2022 MPXV (<xref ref-type="bibr" rid="ref8">Chadha et al., 2022b</xref>). In view of these events, there is an urgent need to develop therapeutics against the circulating MPXV. In this regard, repurposing bioactive phytochemicals and existing drugs proves to be a promising alternative to combat the wave of MPX infections. Recently, phytochemicals with profound antiviral activity have been successfully repurposed against SARS-CoV-2 (<xref ref-type="bibr" rid="ref1">Balkrishna et al., 2021</xref>; <xref ref-type="bibr" rid="ref30">Lingwan et al., 2021</xref>). Hence, the present study was undertaken to virtually screen bioactive phytochemicals targeting the envelope proteins of poxviruses through comprehensive <italic>in silico</italic> analysis. The envelope proteins of vaccinia virus were used as a template for this study since its structural proteins share high sequence similarity with the MPXV (<xref ref-type="bibr" rid="ref28">Lam et al., 2022</xref>). Although poxviruses possess multiple drug targets, the envelope proteins were chosen specifically in an attempt to prevent virus attachment to the host cell, consequently preventing viral infection at the first instance.</p>
<p>The D13 protein is an essential component of the viral capsid that provides rigidity to the structure of poxviruses, thereby playing an essential role in viral morphogenesis. It is a scaffold protein that forms a honeycomb-like lattice on a viral membrane and aids in forming pleomorphic immature virions in conjunction with A14 and A17 proteins (<xref ref-type="bibr" rid="ref47">Van Vliet et al., 2009</xref>). Interestingly, the downregulation of the D13 protein results in irregular single bilayer membranes instead of immature virions (<xref ref-type="bibr" rid="ref51">Zhang and Moss, 1992</xref>). It has been experimentally demonstrated that the D13 protein associates with the viral membrane assembly proteins to form membrane precursors (<xref ref-type="bibr" rid="ref33">Maruri-Avidal et al., 2013</xref>). Vaccinia virus adopts two modes of entry into the host cell, which include (i) direct fusion with the plasma membrane and (ii) endocytosis. The A26 protein is an envelope protein known to interact with the components of the extracellular matrix of the host cell, thereby promoting acid-dependent endocytosis in vaccinia virus (<xref ref-type="bibr" rid="ref3">Bengali et al., 2009</xref>). Studies have demonstrated that the loss of the A26 protein blocks the endocytic pathway, thereby triggering the alternate plasma membrane fusion-based route for viral entry into the host cells (<xref ref-type="bibr" rid="ref10">Chang et al., 2019</xref>). Also, experimental findings by Chang et al. revealed that the loss of the first 75 amino acids from the N-terminal region of the A26 protein in vaccinia virus significantly reduced infectivity as compared to the wild-type strain (<xref ref-type="bibr" rid="ref10">Chang et al., 2019</xref>). Another important envelope protein of poxviruses is the H3 protein that binds to heparin sulfate residues on the host cell surface, directly facilitating virus attachment (<xref ref-type="bibr" rid="ref43">Singh et al., 2016</xref>). The crystal structure of the H3 protein shows that it harbors a glycosyltransferase fold and thus binds to UDP glucose (<xref ref-type="bibr" rid="ref43">Singh et al., 2016</xref>). It is highly conserved and is an immunodominant antigen that is a major target of neutralizing antibodies. Considering the pivotal role of these envelope proteins in the pathogenesis of poxviruses, we selected them as suitable drug targets for this study. Rifampicin was used as a control ligand since it is known to arrest viral morphogenesis before membrane assembly (<xref ref-type="bibr" rid="ref17">Grimley et al., 1970</xref>). Additionally, rifampicin has been shown to occupy the phenylalanine-rich pocket of the D13 trimer, thereby inhibiting its binding to the A17 protein and affecting the formation of the immature virions (<xref ref-type="bibr" rid="ref16">Garriga et al., 2018</xref>). However, this binding is reversible, and the assembly of membrane proteins resumes upon withdrawal of the drug. Rifampicin has widely been used as an antibacterial drug as it targets DNA-dependent RNA polymerases. However, its use as an antiviral drug is skeptical in clinical settings because of its low potency and the emergence of resistant viruses.</p>
<p>In the present study, phytochemicals (ligands) were reported to bind effectively to the envelope proteins of MPXV with significantly lower binding energies. The phytochemicals could easily occupy and fit into the active site (pocket) of vaccinia virus A26, D13, and H3 envelope proteins, thus acting as a roadblock for the virus attachment to the host cell. Additionally, the test phytochemicals and control (rifampicin) shared similar molecular interactions with the amino acid residues on the viral proteins. Moreover, with respect to a single protein, certain amino acid residues displayed common interactions with few phytochemicals. Among the 12 compounds screened, 11 displayed binding energies lower than-6.0&#x2009;kcal/mol against all the three target proteins, except rosmarinic acid, which showed binding energy of &#x2212;6.0&#x2009;kcal/mol against the H3 protein. Among the lead phytochemicals, oleanolic acid and ursolic acid belong to the category of pentacyclic triterpenoids. Numerous studies have demonstrated the antiviral, antitumor, and antimicrobial prospects of triterpenoids (<xref ref-type="bibr" rid="ref50">Yu et al., 2014</xref>; <xref ref-type="bibr" rid="ref49">Xiao et al., 2018</xref>). Oleanolic acid, an essential component of olive oil, exhibits immunomodulatory effects on the human immune system and has been extensively used in cancer prevention by abrogating various signaling pathways (<xref ref-type="bibr" rid="ref39">R&#x00ED;os, 2010</xref>). Numerous studies have also illustrated the application of oleanolic acid and ursolic acid as antiviral agents against different viruses (<xref ref-type="bibr" rid="ref21">Khwaza et al., 2018</xref>; <xref ref-type="bibr" rid="ref42">Shan et al., 2021</xref>). Moreover, silibinin is a flavonolignan known to exhibit anti-neoplastic and antioxidative properties. It also inhibits prostate cancer by targeting multiple signaling pathways, inducing apoptosis, and inhibiting metastasis and angiogenesis (<xref ref-type="bibr" rid="ref14">Deep and Agarwal, 2010</xref>). The existing literature and upcoming investigations collectively increase the foothold for repurposing phytochemicals as alternate therapy against emerging pathogens in clinical settings.</p>
<p>In a recently published study, FDA-approved drugs have also been evaluated for repurposing them against the viral capsid and replication proteins of the vaccinia virus through a comprehensive <italic>in silico</italic> approach (<xref ref-type="bibr" rid="ref28">Lam et al., 2022</xref>). The role of hydrophobic interactions in stabilizing the protein-ligand complex has been well documented across various investigations (<xref ref-type="bibr" rid="ref37">Patil et al., 2010</xref>; <xref ref-type="bibr" rid="ref15">Ferreira De Freitas and Schapira, 2017</xref>). Moreover, hydrogen bondings at the drug-ligand interface are known to play a vital role in stabilizing intermolecular interactions within the complex, thereby promoting complex stiffness (<xref ref-type="bibr" rid="ref37">Patil et al., 2010</xref>). In this study, a relatively greater number of hydrogen and hydrophobic bonds were formed by the test phytochemicals as compared to the rifampicin. Multiple hydrogen bonds between the phytochemicals and target protein(s) can possibly enhance the stability of the interacting complex. The findings from the molecular docking and counter verification by MD simulation and NMA studies studies collectively provided evidence that bioactive phytochemicals can effectively associate with viral envelope proteins and stimulate conformational changes that may interfere with binding to their cognate ligands on the host cell surface.</p>
<p>Altogether, our data sheds light on the possible molecular interaction(s) between bioactive phytochemicals and the envelope proteins of MPXV, thereby interfering with the functional mechanisms of viral entry into the host cell. The investigation generates preliminary findings that indicate the ability of oleanolic acid, ursolic acid, and silibinin to effectively target the D13, A26, and H3 envelope proteins of MPXV. Nevertheless, these <italic>in silico</italic> findings warrant further experimental validation for their application as alternative medicine against MPX in laboratory settings.</p>
</sec>
<sec id="sec22" sec-type="conclusions">
<label>5.</label>
<title>Conclusion</title>
<p>In summary, of all the phytochemicals tested, oleanolic acid, ursolic acid, and silibinin can be envisioned as the lead compounds for repurposing phytochemicals against poxviruses, including the MPXV. Apart from displaying efficacious binding through a plethora of hydrogen and hydrophobic bonds, these 3 phytochemicals were well-fitted in the pockets of all the three envelope proteins taken into consideration. Thus, these compounds dwell into the criteria of a multi-targeted drug delivery approach which utilizes distinct drug candidates that act synergistically against a particular target through different molecular mechanisms. Owing to their natural existence, these bioactive phytochemicals can be used as alternatives to existing chemical drugs that inflict adverse drug reactions and off-target effects. Hence, this study points toward the possible use of phytochemicals to combat the ongoing multi-country outbreak of MPX. However, this demands further validation through <italic>in vitro</italic> and <italic>in vivo</italic> testing.</p>
</sec>
<sec id="sec23" sec-type="data-availability">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/<xref rid="sec26" ref-type="sec">Supplementary material</xref>, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="sec24">
<title>Author contributions</title>
<p>PG: idea conceptualization, perform molecular docking and MD simulation experiments, data analysis and curation, and manuscript writing. JC: idea conceptualization, experimental design, data curation, manuscript writing and editing, and critical analysis. KH: continuous motivation and formal supervision and analysis. SS: continuous motivation, supervision of the study, and manuscript editing. All authors read and approved the final version of this manuscript.</p>
</sec>
<sec id="conf1" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="sec100" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
</body>
<back>
<ack>
<p>JC would like to thank the Indian Council of Medical Research (ICMR), New Delhi for providing financial assistance in the form of fellowship (JRF &#x0026; SRF). PG and JC would also like to acknowledge the relentless and generous efforts of Tushar Kushwaha and Govinda Rao Dabburu for their cooperation and assistance in molecular docking and simulation studies.</p>
</ack>
<sec id="sec26" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary material for this article can be found online at: <ext-link xlink:href="https://www.frontiersin.org/articles/10.3389/fmicb.2022.1073419/full#supplementary-material" ext-link-type="uri">https://www.frontiersin.org/articles/10.3389/fmicb.2022.1073419/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Data_Sheet_1.PDF" id="SM1" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"><label>SUPPLEMENTARY FIGURE S1</label><caption><p>Visualization of molecular docking and interactions between other test phytochemicals docked against the D13 protein (PDB ID: 6BED).</p></caption></supplementary-material>
<supplementary-material xlink:href="Data_Sheet_1.PDF" id="SM2" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"><label>SUPPLEMENTARY FIGURE S2</label><caption><p>Visualization of molecular docking and interactions between other test phytochemicals docked against the A26 protein (PDB ID: 6A9S).</p></caption></supplementary-material>
<supplementary-material xlink:href="Data_Sheet_1.PDF" id="SM3" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"><label>SUPPLEMENTARY FIGURE S3</label><caption><p>Visualization of molecular docking and interactions between other test phytochemicals docked against the H3 protein (PDB ID:5EJ0).</p></caption></supplementary-material>
<supplementary-material xlink:href="Data_Sheet_1.PDF" id="SM4" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"><label>SUPPLEMENTARY FIGURE S4</label><caption><p>Molecular dynamics simulation analysis of silibinin, oleanolic acid, and ursolic acid docked against the D13 protein (PDB ID: 6BED). The left and right panels represent the deformability plot and eigenvalues (plots), respectively.</p></caption></supplementary-material>
<supplementary-material xlink:href="Data_Sheet_1.PDF" id="SM5" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"><label>SUPPLEMENTARY FIGURE S5</label><caption><p>Molecular dynamics simulation analysis of silibinin, oleanolic acid, and ursolic acid docked against the A26 protein (PDB ID: 6A9S). The left and right panels represent the deformability plot and eigenvalues (plots), respectively.</p></caption></supplementary-material>
<supplementary-material xlink:href="Data_Sheet_1.PDF" id="SM6" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"><label>SUPPLEMENTARY FIGURE S6</label><caption><p>Molecular dynamics simulation analysis of silibinin, oleanolic acid, and ursolic acid docked against the H3 protein (PDB ID: 5EJ0). The left and right panels represent the deformability plot and eigenvalues (plots), respectively.</p></caption></supplementary-material>
<supplementary-material xlink:href="Data_Sheet_1.PDF" id="SM7" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"><label>SUPPLEMENTARY FIGURE S7</label><caption><p>Molecular dynamics simulation analysis of silibinin, oleanolic acid, and ursolic acid docked against the D13 protein (PDB ID: 6BED). Left vertical panel shows the covariance matrix. Correlated region has been depicted in red, uncorrelated region in white, and anticorrelated region in blue. Right vertical panel shows elastic network of the docked complexes (dark gray regions indicates more stiffer regions).</p></caption></supplementary-material>
<supplementary-material xlink:href="Data_Sheet_1.PDF" id="SM8" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"><label>SUPPLEMENTARY FIGURE S8</label><caption><p>Molecular dynamics simulation analysis of silibinin, oleanolic acid, and ursolic acid docked against the A26 protein (PDB ID: 6A9S). Left vertical panel shows the covariance matrix. Correlated region has been depicted in red, uncorrelated region in white, and anticorrelated region in blue. Right vertical panel shows elastic network of the docked complexes (dark gray regions indicates more stiffer regions).</p></caption></supplementary-material>
<supplementary-material xlink:href="Data_Sheet_1.PDF" id="SM9" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"><label>SUPPLEMENTARY FIGURE S9</label><caption><p>Molecular dynamics simulation analysis of silibinin, oleanolic acid, and ursolic acid docked against the H3 protein (PDB ID: 5EJ0). Left vertical panel shows the covariance matrix. Correlated region has been depicted in red, uncorrelated region in white, and anticorrelated region in blue. Right vertical panel shows elastic network of the docked complexes (dark gray regions indicates more stiffer regions).</p></caption></supplementary-material>
<supplementary-material xlink:href="Table_1.DOC" id="SM10" mimetype="application/msword" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="ref1"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Balkrishna</surname> <given-names>A.</given-names></name> <name><surname>Pokhrel</surname> <given-names>S.</given-names></name> <name><surname>Singh</surname> <given-names>H.</given-names></name> <name><surname>Joshi</surname> <given-names>M.</given-names></name> <name><surname>Mulay</surname> <given-names>V. P.</given-names></name> <name><surname>Haldar</surname> <given-names>S.</given-names></name> <etal/></person-group>. (<year>2021</year>). <article-title>Withanone from <italic>Withania somnifera</italic> attenuates SARS-CoV-2 RBD and host ACE2 interactions to rescue spike protein induced pathologies in humanized zebrafish model</article-title>. <source>Drug Des. Devel. Ther.</source> <volume>15</volume>, <fpage>1111</fpage>&#x2013;<lpage>1133</lpage>. doi: <pub-id pub-id-type="doi">10.2147/DDDT.S292805</pub-id>, PMID: <pub-id pub-id-type="pmid">33737804</pub-id></citation></ref>
<ref id="ref2"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Benet</surname> <given-names>L. Z.</given-names></name> <name><surname>Hosey</surname> <given-names>C. M.</given-names></name> <name><surname>Ursu</surname> <given-names>O.</given-names></name> <name><surname>Oprea</surname> <given-names>T. I.</given-names></name></person-group> (<year>2016</year>). <article-title>BDDCS, the rule of 5 and drugability</article-title>. <source>Adv. Drug Deliv. Rev.</source> <volume>101</volume>, <fpage>89</fpage>&#x2013;<lpage>98</lpage>. doi: <pub-id pub-id-type="doi">10.1016/J.ADDR.2016.05.007</pub-id>, PMID: <pub-id pub-id-type="pmid">27182629</pub-id></citation></ref>
<ref id="ref3"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bengali</surname> <given-names>Z.</given-names></name> <name><surname>Townsley</surname> <given-names>A. C.</given-names></name> <name><surname>Moss</surname> <given-names>B.</given-names></name></person-group> (<year>2009</year>). <article-title>Vaccinia virus strain differences in cell attachment and entry</article-title>. <source>Virology</source> <volume>389</volume>, <fpage>132</fpage>&#x2013;<lpage>140</lpage>. doi: <pub-id pub-id-type="doi">10.1016/J.VIROL.2009.04.012</pub-id>, PMID: <pub-id pub-id-type="pmid">19428041</pub-id></citation></ref>
<ref id="ref4"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bunge</surname> <given-names>E. M.</given-names></name> <name><surname>Hoet</surname> <given-names>B.</given-names></name> <name><surname>Chen</surname> <given-names>L.</given-names></name> <name><surname>Lienert</surname> <given-names>F.</given-names></name> <name><surname>Weidenthaler</surname> <given-names>H.</given-names></name> <name><surname>Baer</surname> <given-names>L. R.</given-names></name> <etal/></person-group>. (<year>2022</year>). <article-title>The changing epidemiology of human monkeypox&#x2014;a potential threat? A systematic review</article-title>. <source>PLoS Negl. Trop. Dis.</source> <volume>16</volume>:<fpage>e0010141</fpage>. doi: <pub-id pub-id-type="doi">10.1371/JOURNAL.PNTD.0010141</pub-id>, PMID: <pub-id pub-id-type="pmid">35148313</pub-id></citation></ref>
<ref id="ref5"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Carter</surname> <given-names>G. C.</given-names></name> <name><surname>Law</surname> <given-names>M.</given-names></name> <name><surname>Hollinshead</surname> <given-names>M.</given-names></name> <name><surname>Smith</surname> <given-names>G. L.</given-names></name></person-group> (<year>2005</year>). <article-title>Entry of the vaccinia virus intracellular mature virion and its interactions with glycosaminoglycans</article-title>. <source>J. Gen. Virol.</source> <volume>86</volume>, <fpage>1279</fpage>&#x2013;<lpage>1290</lpage>. doi: <pub-id pub-id-type="doi">10.1099/VIR.0.80831-0</pub-id>, PMID: <pub-id pub-id-type="pmid">15831938</pub-id></citation></ref>
<ref id="ref6"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chadha</surname> <given-names>J.</given-names></name> <name><surname>Gupta</surname> <given-names>M.</given-names></name> <name><surname>Nagpal</surname> <given-names>N.</given-names></name> <name><surname>Sharma</surname> <given-names>M.</given-names></name> <name><surname>Adarsh</surname> <given-names>T.</given-names></name> <name><surname>Joshi</surname> <given-names>V.</given-names></name> <etal/></person-group>. (<year>2021</year>). <article-title>Antibacterial potential of indigenous plant extracts against multidrug-resistant bacterial strains isolated from New Delhi region</article-title>. <source>GSC Biol. Pharm. Sci.</source> <volume>14</volume>, <fpage>185</fpage>&#x2013;<lpage>196</lpage>. doi: <pub-id pub-id-type="doi">10.30574/GSCBPS.2021.14.2.0053</pub-id></citation></ref>
<ref id="ref7"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chadha</surname> <given-names>J.</given-names></name> <name><surname>Harjai</surname> <given-names>K.</given-names></name> <name><surname>Chhibber</surname> <given-names>S.</given-names></name></person-group> (<year>2022a</year>). <article-title>Repurposing phytochemicals as anti-virulent agents to attenuate quorum sensing-regulated virulence factors and biofilm formation in <italic>Pseudomonas aeruginosa</italic></article-title>. <source>Microb. Biotechnol.</source> <volume>15</volume>, <fpage>1695</fpage>&#x2013;<lpage>1718</lpage>. doi: <pub-id pub-id-type="doi">10.1111/1751-7915.13981</pub-id>, PMID: <pub-id pub-id-type="pmid">34843159</pub-id></citation></ref>
<ref id="ref8"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chadha</surname> <given-names>J.</given-names></name> <name><surname>Khullar</surname> <given-names>L.</given-names></name> <name><surname>Gulati</surname> <given-names>P.</given-names></name> <name><surname>Chhibber</surname> <given-names>S.</given-names></name> <name><surname>Harjai</surname> <given-names>K.</given-names></name></person-group> (<year>2022b</year>). <article-title>Insights into the monkeypox virus: making of another pandemic within the pandemic?</article-title> <source>Environ. Microbiol.</source> <volume>24</volume>, <fpage>4547</fpage>&#x2013;<lpage>4560</lpage>. doi: <pub-id pub-id-type="doi">10.1111/1462-2920.16174</pub-id>, PMID: <pub-id pub-id-type="pmid">35974453</pub-id></citation></ref>
<ref id="ref9"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chadha</surname> <given-names>J.</given-names></name> <name><surname>Khullar</surname> <given-names>L.</given-names></name> <name><surname>Mittal</surname> <given-names>N.</given-names></name></person-group> (<year>2022c</year>). <article-title>Facing the wrath of enigmatic mutations: a review on the emergence of severe acute respiratory syndrome coronavirus 2 variants amid coronavirus disease-19 pandemic</article-title>. <source>Environ. Microbiol.</source> <volume>24</volume>, <fpage>2615</fpage>&#x2013;<lpage>2629</lpage>. doi: <pub-id pub-id-type="doi">10.1111/1462-2920.15687</pub-id>, PMID: <pub-id pub-id-type="pmid">34320263</pub-id></citation></ref>
<ref id="ref10"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chang</surname> <given-names>H. W.</given-names></name> <name><surname>Yang</surname> <given-names>C. H.</given-names></name> <name><surname>Luo</surname> <given-names>Y. C.</given-names></name> <name><surname>Su</surname> <given-names>B. G.</given-names></name> <name><surname>Cheng</surname> <given-names>H. Y.</given-names></name> <name><surname>Tung</surname> <given-names>S. Y.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>Vaccinia viral A26 protein is a fusion suppressor of mature virus and triggers membrane fusion through conformational change at low pH</article-title>. <source>PLoS Pathog.</source> <volume>15</volume>:<fpage>e1007826</fpage>. doi: <pub-id pub-id-type="doi">10.1371/JOURNAL.PPAT.1007826</pub-id>, PMID: <pub-id pub-id-type="pmid">31220181</pub-id></citation></ref>
<ref id="ref11"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chiu</surname> <given-names>W.-L.</given-names></name> <name><surname>Lin</surname> <given-names>C.-L.</given-names></name> <name><surname>Yang</surname> <given-names>M.-H.</given-names></name> <name><surname>Tzou</surname> <given-names>D.-L. M.</given-names></name> <name><surname>Chang</surname> <given-names>W.</given-names></name></person-group> (<year>2007</year>). <article-title>Vaccinia virus 4c (A26L) protein on intracellular mature virus binds to the extracellular cellular matrix laminin</article-title>. <source>J. Virol.</source> <volume>81</volume>, <fpage>2149</fpage>&#x2013;<lpage>2157</lpage>. doi: <pub-id pub-id-type="doi">10.1128/JVI.02302-06</pub-id>, PMID: <pub-id pub-id-type="pmid">17166913</pub-id></citation></ref>
<ref id="ref12"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Choudhury</surname> <given-names>A.</given-names></name> <name><surname>Das</surname> <given-names>N. C.</given-names></name> <name><surname>Patra</surname> <given-names>R.</given-names></name> <name><surname>Bhattacharya</surname> <given-names>M.</given-names></name> <name><surname>Ghosh</surname> <given-names>P.</given-names></name> <name><surname>Patra</surname> <given-names>B. C.</given-names></name> <etal/></person-group>. (<year>2021</year>). <article-title>Exploring the binding efficacy of ivermectin against the key proteins of SARS-CoV-2 pathogenesis: an approach</article-title>. <source>Future Virol.</source> <volume>16</volume>, <fpage>277</fpage>&#x2013;<lpage>291</lpage>. doi: <pub-id pub-id-type="doi">10.2217/FVL-2020-0342</pub-id></citation></ref>
<ref id="ref13"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Daina</surname> <given-names>A.</given-names></name> <name><surname>Michielin</surname> <given-names>O.</given-names></name> <name><surname>Zoete</surname> <given-names>V.</given-names></name></person-group> (<year>2017</year>). <article-title>SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules</article-title>. <source>Sci. Rep.</source> <volume>7</volume>, <fpage>1</fpage>&#x2013;<lpage>13</lpage>. doi: <pub-id pub-id-type="doi">10.1038/srep42717</pub-id>, PMID: <pub-id pub-id-type="pmid">28256516</pub-id></citation></ref>
<ref id="ref14"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Deep</surname> <given-names>G.</given-names></name> <name><surname>Agarwal</surname> <given-names>R.</given-names></name></person-group> (<year>2010</year>). <article-title>Antimetastatic efficacy of silibinin: molecular mechanisms and therapeutic potential against cancer</article-title>. <source>Cancer Metastasis Rev.</source> <volume>2010</volume> 293, <fpage>447</fpage>&#x2013;<lpage>463</lpage>. doi: <pub-id pub-id-type="doi">10.1007/S10555-010-9237-0</pub-id></citation></ref>
<ref id="ref15"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ferreira De Freitas</surname> <given-names>R.</given-names></name> <name><surname>Schapira</surname> <given-names>M.</given-names></name></person-group> (<year>2017</year>). <article-title>A systematic analysis of atomic protein&#x2013;ligand interactions in the PDB</article-title>. <source>Medchemcomm</source> <volume>8</volume>, <fpage>1970</fpage>&#x2013;<lpage>1981</lpage>. doi: <pub-id pub-id-type="doi">10.1039/C7MD00381A</pub-id>, PMID: <pub-id pub-id-type="pmid">29308120</pub-id></citation></ref>
<ref id="ref16"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Garriga</surname> <given-names>D.</given-names></name> <name><surname>Headey</surname> <given-names>S.</given-names></name> <name><surname>Accurso</surname> <given-names>C.</given-names></name> <name><surname>Gunzburg</surname> <given-names>M.</given-names></name> <name><surname>Scanlon</surname> <given-names>M.</given-names></name> <name><surname>Coulibaly</surname> <given-names>F.</given-names></name></person-group> (<year>2018</year>). <article-title>Structural basis for the inhibition of poxvirus assembly by the antibiotic rifampicin</article-title>. <source>Proc. Natl. Acad. Sci. U. S. A.</source> <volume>115</volume>, <fpage>8424</fpage>&#x2013;<lpage>8429</lpage>. doi: <pub-id pub-id-type="doi">10.1073/PNAS.1810398115</pub-id>, PMID: <pub-id pub-id-type="pmid">30068608</pub-id></citation></ref>
<ref id="ref17"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Grimley</surname> <given-names>P. M.</given-names></name> <name><surname>Rosenblum</surname> <given-names>E. N.</given-names></name> <name><surname>Mims</surname> <given-names>S. J.</given-names></name> <name><surname>Moss</surname> <given-names>B.</given-names></name></person-group> (<year>1970</year>). <article-title>Interruption by rifampin of an early stage in vaccinia virus morphogenesis: accumulation of membranes which are precursors of virus envelopes</article-title>. <source>J. Virol.</source> <volume>6</volume>, <fpage>519</fpage>&#x2013;<lpage>533</lpage>. doi: <pub-id pub-id-type="doi">10.1128/JVI.6.4.519-533.1970</pub-id>, PMID: <pub-id pub-id-type="pmid">5497899</pub-id></citation></ref>
<ref id="ref18"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gulati</surname> <given-names>P.</given-names></name> <name><surname>Yadav</surname> <given-names>A.</given-names></name> <name><surname>Chadha</surname> <given-names>J.</given-names></name> <name><surname>Singh</surname> <given-names>S.</given-names></name></person-group> (<year>2021</year>). <article-title>Virtual screening of phytochemicals targeting the Main protease and spike protein of SARS-CoV-2: an in silico approach</article-title>. <source>J. Biomed. Res. Environ. Sci.</source> <volume>2</volume>, <fpage>1121</fpage>&#x2013;<lpage>1131</lpage>. doi: <pub-id pub-id-type="doi">10.37871/JBRES1357</pub-id></citation></ref>
<ref id="ref19"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Isidro</surname> <given-names>J.</given-names></name> <name><surname>Borges</surname> <given-names>V.</given-names></name> <name><surname>Pinto</surname> <given-names>M.</given-names></name> <name><surname>Sobral</surname> <given-names>D.</given-names></name> <name><surname>Santos</surname> <given-names>J. D.</given-names></name> <name><surname>Nunes</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2022</year>). <article-title>Phylogenomic characterization and signs of microevolution in the 2022 multi-country outbreak of monkeypox virus</article-title>. <source>Nat. Med.</source> <volume>28</volume>, <fpage>1569</fpage>&#x2013;<lpage>1572</lpage>. doi: <pub-id pub-id-type="doi">10.1038/s41591-022-01907-y</pub-id>, PMID: <pub-id pub-id-type="pmid">35750157</pub-id></citation></ref>
<ref id="ref20"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jos</surname> <given-names>J.</given-names></name> <name><surname>Ramon</surname> <given-names>J.</given-names></name> <name><surname>Ramon</surname> <given-names>R.</given-names></name> <name><surname>Opez-Blanco</surname> <given-names>L.</given-names></name> <name><surname>Aliaga</surname> <given-names>J. I.</given-names></name> <name><surname>Quintana-Ort&#x00ED;</surname> <given-names>E. S.</given-names></name> <etal/></person-group>. (<year>2014</year>). <article-title>iMODS: internal coordinates normal mode analysis server</article-title>. <source>Nucleic Acids Res.</source> <volume>42</volume>, <fpage>W271</fpage>&#x2013;<lpage>W276</lpage>. doi: <pub-id pub-id-type="doi">10.1093/NAR/GKU339</pub-id>, PMID: <pub-id pub-id-type="pmid">24771341</pub-id></citation></ref>
<ref id="ref21"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Khwaza</surname> <given-names>V.</given-names></name> <name><surname>Oyedeji</surname> <given-names>O. O.</given-names></name> <name><surname>Aderibigbe</surname> <given-names>B. A.</given-names></name> <name><surname>Ciminiello</surname> <given-names>P.</given-names></name> <name><surname>Mangoni</surname> <given-names>A.</given-names></name> <name><surname>Menna</surname> <given-names>M.</given-names></name> <etal/></person-group>. (<year>2018</year>). <article-title>Antiviral activities of oleanolic acid and its analogues</article-title>. <source>Molecules</source> <volume>23</volume>:<fpage>2300</fpage>. doi: <pub-id pub-id-type="doi">10.3390/MOLECULES23092300</pub-id>, PMID: <pub-id pub-id-type="pmid">30205592</pub-id></citation></ref>
<ref id="ref22"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kugelman</surname> <given-names>J. R.</given-names></name> <name><surname>Johnston</surname> <given-names>S. C.</given-names></name> <name><surname>Mulembakani</surname> <given-names>P. M.</given-names></name> <name><surname>Kisalu</surname> <given-names>N.</given-names></name> <name><surname>Lee</surname> <given-names>M. S.</given-names></name> <name><surname>Koroleva</surname> <given-names>G.</given-names></name> <etal/></person-group>. (<year>2014</year>). <article-title>Genomic variability of monkeypox virus among humans, Democratic Republic of the Congo</article-title>. <source>Emerg. Infect. Dis.</source> <volume>20</volume>, <fpage>232</fpage>&#x2013;<lpage>239</lpage>. doi: <pub-id pub-id-type="doi">10.3201/EID2002.130118</pub-id>, PMID: <pub-id pub-id-type="pmid">24457084</pub-id></citation></ref>
<ref id="ref23"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kumar</surname> <given-names>V.</given-names></name> <name><surname>Dhanjal</surname> <given-names>J. K.</given-names></name> <name><surname>Bhargava</surname> <given-names>P.</given-names></name> <name><surname>Kaul</surname> <given-names>A.</given-names></name> <name><surname>Wang</surname> <given-names>J.</given-names></name> <name><surname>Zhang</surname> <given-names>H.</given-names></name> <etal/></person-group>. (<year>2020</year>). <article-title>Withanone and Withaferin-a are predicted to interact with transmembrane protease serine 2 (TMPRSS2) and block entry of SARS-CoV-2 into cells</article-title>. <source>J. Biomol. Struct. Dyn.</source> <volume>40</volume>, <fpage>1</fpage>&#x2013;<lpage>13</lpage>. doi: <pub-id pub-id-type="doi">10.1080/07391102.2020.1775704</pub-id>, PMID: <pub-id pub-id-type="pmid">32469279</pub-id></citation></ref>
<ref id="ref24"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kumar</surname> <given-names>S.</given-names></name> <name><surname>Garg</surname> <given-names>C.</given-names></name> <name><surname>Kaushik</surname> <given-names>S.</given-names></name> <name><surname>Buttar</surname> <given-names>H.</given-names></name> <name><surname>Garg</surname> <given-names>M.</given-names></name></person-group> (<year>2021</year>). <article-title>Demystifying therapeutic potential of medicinal plants against chikungunya virus</article-title>. <source>Indian J. Pharmacol.</source> <volume>53</volume>, <fpage>403</fpage>&#x2013;<lpage>411</lpage>. doi: <pub-id pub-id-type="doi">10.4103/IJP.IJP_81_20</pub-id>, PMID: <pub-id pub-id-type="pmid">34854411</pub-id></citation></ref>
<ref id="ref25"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kumar</surname> <given-names>A.</given-names></name> <name><surname>Mishra</surname> <given-names>D. C.</given-names></name> <name><surname>Angadi</surname> <given-names>U. B.</given-names></name> <name><surname>Yadav</surname> <given-names>R.</given-names></name> <name><surname>Rai</surname> <given-names>A.</given-names></name> <name><surname>Kumar</surname> <given-names>D.</given-names></name></person-group> (<year>2021</year>). <article-title>Inhibition potencies of phytochemicals derived from sesame against SARS-CoV-2 main protease: a molecular docking and simulation study</article-title>. <source>Front. Chem.</source> <volume>9</volume>:<fpage>773</fpage>. doi: <pub-id pub-id-type="doi">10.3389/FCHEM.2021.744376</pub-id>, PMID: <pub-id pub-id-type="pmid">34692642</pub-id></citation></ref>
<ref id="ref26"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kumari</surname> <given-names>R.</given-names></name> <name><surname>Kumar</surname> <given-names>R</given-names></name> <collab id="coll1">Open Source Drug Discovery Consortium</collab> <name><surname>Lynn</surname> <given-names>A</given-names></name></person-group>. (<year>2014</year>). <article-title>g_mmpbsa--a GROMACS tool for high-throughput MM-PBSA calculations</article-title>. <source>J. Chem. Inf. Model.</source> <volume>54</volume>:<fpage>1951</fpage>&#x2013;<lpage>1962</lpage>. doi: <pub-id pub-id-type="doi">10.1021/ci500020m</pub-id></citation></ref>
<ref id="ref27"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kwofie</surname> <given-names>S. K.</given-names></name> <name><surname>Broni</surname> <given-names>E.</given-names></name> <name><surname>Teye</surname> <given-names>J.</given-names></name> <name><surname>Quansah</surname> <given-names>E.</given-names></name> <name><surname>Issah</surname> <given-names>I.</given-names></name> <name><surname>Wilson</surname> <given-names>M. D.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>Pharmacoinformatics-based identification of potential bioactive compounds against Ebola virus protein VP24</article-title>. <source>Comput. Biol. Med.</source> <volume>113</volume>:<fpage>103414</fpage>. doi: <pub-id pub-id-type="doi">10.1016/J.COMPBIOMED.2019.103414</pub-id>, PMID: <pub-id pub-id-type="pmid">31536833</pub-id></citation></ref>
<ref id="ref28"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lam</surname> <given-names>H. Y. I.</given-names></name> <name><surname>Guan</surname> <given-names>J. S.</given-names></name> <name><surname>Mu</surname> <given-names>Y.</given-names></name></person-group> (<year>2022</year>). <article-title>In silico repurposed drugs against monkeypox virus</article-title>. <source>Molecules</source> <volume>27</volume>, <fpage>1</fpage>&#x2013;<lpage>20</lpage>. doi: <pub-id pub-id-type="doi">10.3390/MOLECULES27165277/S1</pub-id>, PMID: <pub-id pub-id-type="pmid">36014515</pub-id></citation></ref>
<ref id="ref29"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ligon</surname> <given-names>B. L.</given-names></name></person-group> (<year>2004</year>). <article-title>Monkeypox: a review of the history and emergence in the Western hemisphere</article-title>. <source>Semin. Pediatr. Infect. Dis.</source> <volume>15</volume>, <fpage>280</fpage>&#x2013;<lpage>287</lpage>. doi: <pub-id pub-id-type="doi">10.1053/J.SPID.2004.09.001</pub-id>, PMID: <pub-id pub-id-type="pmid">15494953</pub-id></citation></ref>
<ref id="ref30"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lingwan</surname> <given-names>M.</given-names></name> <name><surname>Shagun</surname> <given-names>S.</given-names></name> <name><surname>Pahwa</surname> <given-names>F.</given-names></name> <name><surname>Kumar</surname> <given-names>A.</given-names></name> <name><surname>Verma</surname> <given-names>D. K.</given-names></name> <name><surname>Pant</surname> <given-names>Y.</given-names></name> <etal/></person-group>. (<year>2021</year>). <article-title>Phytochemical rich Himalayan Rhododendron arboreum petals inhibit SARS-CoV-2 infection in vitro</article-title>. <source>J. Biomol. Struct. Dyn.</source>, <fpage>1</fpage>&#x2013;<lpage>11</lpage>. doi: <pub-id pub-id-type="doi">10.1080/07391102.2021.2021287</pub-id>, PMID: <pub-id pub-id-type="pmid">34961411</pub-id></citation></ref>
<ref id="ref31"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lipinski</surname> <given-names>C. A.</given-names></name></person-group> (<year>2004</year>). <article-title>Lead-and drug-like compounds: the rule-of-five revolution</article-title>. <source>Drug Discov. Today Technol.</source> <volume>1</volume>, <fpage>337</fpage>&#x2013;<lpage>341</lpage>. doi: <pub-id pub-id-type="doi">10.1016/J.DDTEC.2004.11.007</pub-id>, PMID: <pub-id pub-id-type="pmid">24981612</pub-id></citation></ref>
<ref id="ref32"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Marennikova</surname> <given-names>S. S.</given-names></name> <name><surname>Seluhina</surname> <given-names>E. M.</given-names></name> <name><surname>Malceva</surname> <given-names>N. N.</given-names></name> <name><surname>Cimiskjan</surname> <given-names>K. L.</given-names></name> <name><surname>Macevic</surname> <given-names>G. R.</given-names></name></person-group> (<year>1972</year>). <article-title>Isolation and properties of the causal agent of a new variola-like disease (monkeypox) in man</article-title>. <source>Bull. World Health Organ.</source> <volume>46</volume>, <fpage>599</fpage>&#x2013;<lpage>611</lpage>. PMID: <pub-id pub-id-type="pmid">4340219</pub-id></citation></ref>
<ref id="ref33"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Maruri-Avidal</surname> <given-names>L.</given-names></name> <name><surname>Weisberg</surname> <given-names>A. S.</given-names></name> <name><surname>Moss</surname> <given-names>B.</given-names></name></person-group> (<year>2013</year>). <article-title>Direct formation of vaccinia virus membranes from the endoplasmic reticulum in the absence of the newly characterized L2-interacting protein A30.5</article-title>. <source>J. Virol.</source> <volume>87</volume>, <fpage>12313</fpage>&#x2013;<lpage>12326</lpage>. doi: <pub-id pub-id-type="doi">10.1128/JVI.02137-13</pub-id>, PMID: <pub-id pub-id-type="pmid">24027302</pub-id></citation></ref>
<ref id="ref34"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mourya</surname> <given-names>D.</given-names></name> <name><surname>Yadav</surname> <given-names>P.</given-names></name> <name><surname>Ullas</surname> <given-names>P.</given-names></name> <name><surname>Bhardwaj</surname> <given-names>S.</given-names></name> <name><surname>Sahay</surname> <given-names>R.</given-names></name> <name><surname>Chadha</surname> <given-names>M.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>Emerging/re-emerging viral diseases &#x0026; new viruses on the Indian horizon</article-title>. <source>Indian J. Med. Res.</source> <volume>149</volume>, <fpage>447</fpage>&#x2013;<lpage>467</lpage>. doi: <pub-id pub-id-type="doi">10.4103/IJMR.IJMR_1239_18</pub-id>, PMID: <pub-id pub-id-type="pmid">31411169</pub-id></citation></ref>
<ref id="ref35"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ortiz</surname> <given-names>C. L. D.</given-names></name> <name><surname>Completo</surname> <given-names>G. C.</given-names></name> <name><surname>Nacario</surname> <given-names>R. C.</given-names></name> <name><surname>Nellas</surname> <given-names>R. B.</given-names></name></person-group> (<year>2019</year>). <article-title>Potential inhibitors of Galactofuranosyltransferase 2 (GlfT2): molecular docking, 3D-QSAR, and in silico ADMETox studies</article-title>. <source>Sci. Rep.</source> <volume>9</volume>, <fpage>17096</fpage>&#x2013;<lpage>17028</lpage>. doi: <pub-id pub-id-type="doi">10.1038/s41598-019-52764-8</pub-id>, PMID: <pub-id pub-id-type="pmid">31745103</pub-id></citation></ref>
<ref id="ref36"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pathania</surname> <given-names>S.</given-names></name> <name><surname>Randhawa</surname> <given-names>V.</given-names></name> <name><surname>Kumar</surname> <given-names>M.</given-names></name></person-group> (<year>2019</year>). <article-title>Identifying potential entry inhibitors for emerging Nipah virus by molecular docking and chemical-protein interaction network</article-title>. <source>J. Biomol. Struct. Dyn.</source> <volume>38</volume>, <fpage>5108</fpage>&#x2013;<lpage>5125</lpage>. doi: <pub-id pub-id-type="doi">10.1080/07391102.2019.1696705</pub-id>, PMID: <pub-id pub-id-type="pmid">31771426</pub-id></citation></ref>
<ref id="ref37"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Patil</surname> <given-names>R.</given-names></name> <name><surname>Das</surname> <given-names>S.</given-names></name> <name><surname>Stanley</surname> <given-names>A.</given-names></name> <name><surname>Yadav</surname> <given-names>L.</given-names></name> <name><surname>Sudhakar</surname> <given-names>A.</given-names></name> <name><surname>Varma</surname> <given-names>A. K.</given-names></name></person-group> (<year>2010</year>). <article-title>Optimized hydrophobic interactions and hydrogen bonding at the target-ligand Interface leads the pathways of drug-designing</article-title>. <source>PLoS One</source> <volume>5</volume>:<fpage>e12029</fpage>. doi: <pub-id pub-id-type="doi">10.1371/JOURNAL.PONE.0012029</pub-id>, PMID: <pub-id pub-id-type="pmid">20808434</pub-id></citation></ref>
<ref id="ref38"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Reynolds</surname> <given-names>M. G.</given-names></name> <name><surname>Carroll</surname> <given-names>D. S.</given-names></name> <name><surname>Karem</surname> <given-names>K. L.</given-names></name></person-group> (<year>2012</year>). <article-title>Factors affecting the likelihood of monkeypox&#x2019;s emergence and spread in the post-smallpox era</article-title>. <source>Curr. Opin. Virol.</source> <volume>2</volume>, <fpage>335</fpage>&#x2013;<lpage>343</lpage>. doi: <pub-id pub-id-type="doi">10.1016/J.COVIRO.2012.02.004</pub-id>, PMID: <pub-id pub-id-type="pmid">22709519</pub-id></citation></ref>
<ref id="ref39"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>R&#x00ED;os</surname> <given-names>J. L.</given-names></name></person-group> (<year>2010</year>). <article-title>Effects of triterpenes on the immune system</article-title>. <source>J. Ethnopharmacol.</source> <volume>128</volume>, <fpage>1</fpage>&#x2013;<lpage>14</lpage>. doi: <pub-id pub-id-type="doi">10.1016/J.JEP.2009.12.045</pub-id></citation></ref>
<ref id="ref40"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Schmid</surname> <given-names>N.</given-names></name> <name><surname>Eichenberger</surname> <given-names>A. P.</given-names></name> <name><surname>Choutko</surname> <given-names>A.</given-names></name> <name><surname>Riniker</surname> <given-names>S.</given-names></name> <name><surname>Winger</surname> <given-names>M.</given-names></name> <name><surname>Mark</surname> <given-names>A. E.</given-names></name> <etal/></person-group>. (<year>2011</year>). <article-title>Definition and testing of the GROMOS force-field versions 54A7 and 54B7</article-title>. <source>Eur. Biophys. J.</source> <volume>40</volume>, <fpage>843</fpage>&#x2013;<lpage>856</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s00249-011-0700-9</pub-id>, PMID: <pub-id pub-id-type="pmid">21533652</pub-id></citation></ref>
<ref id="ref41"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sch&#x00FC;ttelkopf</surname> <given-names>A. W.</given-names></name> <name><surname>van Aalten</surname> <given-names>D. M.</given-names></name></person-group> (<year>2004</year>). <article-title>PRODRG: a tool for high-throughput crystallography of protein-ligand complexes</article-title>. <source>Acta Crystallogr. D Biol. Crystallogr.</source> <volume>60</volume>, <fpage>1355</fpage>&#x2013;<lpage>1363</lpage>. doi: <pub-id pub-id-type="doi">10.1107/S0907444904011679</pub-id>, PMID: <pub-id pub-id-type="pmid">15272157</pub-id></citation></ref>
<ref id="ref42"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shan</surname> <given-names>T.</given-names></name> <name><surname>Ye</surname> <given-names>J.</given-names></name> <name><surname>Jia</surname> <given-names>J.</given-names></name> <name><surname>Wang</surname> <given-names>Z.</given-names></name> <name><surname>Jiang</surname> <given-names>Y.</given-names></name> <name><surname>Wang</surname> <given-names>Y.</given-names></name> <etal/></person-group>. (<year>2021</year>). <article-title>Viral UL8 is involved in the antiviral activity of oleanolic acid against HSV-1 infection</article-title>. <source>Front. Microbiol.</source> <volume>12</volume>:<fpage>2046</fpage>. doi: <pub-id pub-id-type="doi">10.3389/FMICB.2021.689607</pub-id>, PMID: <pub-id pub-id-type="pmid">34354687</pub-id></citation></ref>
<ref id="ref43"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Singh</surname> <given-names>K.</given-names></name> <name><surname>Gittis</surname> <given-names>A. G.</given-names></name> <name><surname>Gitti</surname> <given-names>R. K.</given-names></name> <name><surname>Ostazeski</surname> <given-names>S. A.</given-names></name> <name><surname>Su</surname> <given-names>H.-P.</given-names></name> <name><surname>Garboczi</surname> <given-names>D. N.</given-names></name></person-group> (<year>2016</year>). <article-title>The vaccinia virus H3 envelope protein, a major target of neutralizing antibodies, exhibits a glycosyltransferase fold and binds UDP-glucose</article-title>. <source>J. Virol.</source> <volume>90</volume>, <fpage>5020</fpage>&#x2013;<lpage>5030</lpage>. doi: <pub-id pub-id-type="doi">10.1128/JVI.02933-15</pub-id>, PMID: <pub-id pub-id-type="pmid">26937025</pub-id></citation></ref>
<ref id="ref44"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sklenovsk&#x00E1;</surname> <given-names>N.</given-names></name> <name><surname>Van Ranst</surname> <given-names>M.</given-names></name></person-group> (<year>2018</year>). <article-title>Emergence of Monkeypox as the most important Orthopoxvirus infection in humans</article-title>. <source>Front. Public Heal.</source> <volume>6</volume>:<fpage>241</fpage>. doi: <pub-id pub-id-type="doi">10.3389/FPUBH.2018.0024</pub-id>, PMID: <pub-id pub-id-type="pmid">30234087</pub-id></citation></ref>
<ref id="ref45"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Titanji</surname> <given-names>B. K.</given-names></name> <name><surname>Tegomoh</surname> <given-names>B.</given-names></name> <name><surname>Nematollahi</surname> <given-names>S.</given-names></name> <name><surname>Konomos</surname> <given-names>M.</given-names></name> <name><surname>Kulkarni</surname> <given-names>P. A.</given-names></name></person-group> (<year>2022</year>). <article-title>Monkeypox: a contemporary review for healthcare professionals</article-title>. <source>Open Forum Infect. Dis.</source> <volume>9</volume>:<fpage>ofac310</fpage>. doi: <pub-id pub-id-type="doi">10.1093/OFID/OFAC310</pub-id>, PMID: <pub-id pub-id-type="pmid">35891689</pub-id></citation></ref>
<ref id="ref46"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Van Der Spoel</surname> <given-names>D.</given-names></name> <name><surname>Lindahl</surname> <given-names>E.</given-names></name> <name><surname>Hess</surname> <given-names>B.</given-names></name> <name><surname>Groenhof</surname> <given-names>G.</given-names></name> <name><surname>Mark</surname> <given-names>A. E.</given-names></name> <name><surname>Berendsen</surname> <given-names>H. J.</given-names></name></person-group> (<year>2005</year>). <article-title>GROMACS: fast, flexible, and free</article-title>. <source>J. Comput. Chem.</source> <volume>26</volume>, <fpage>1701</fpage>&#x2013;<lpage>1718</lpage>. doi: <pub-id pub-id-type="doi">10.1002/jcc.20291</pub-id>, PMID: <pub-id pub-id-type="pmid">16211538</pub-id></citation></ref>
<ref id="ref47"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Van Vliet</surname> <given-names>K.</given-names></name> <name><surname>Mohamed</surname> <given-names>M. R.</given-names></name> <name><surname>Zhang</surname> <given-names>L.</given-names></name> <name><surname>Villa</surname> <given-names>N. Y.</given-names></name> <name><surname>Werden</surname> <given-names>S. J.</given-names></name> <name><surname>Liu</surname> <given-names>J.</given-names></name> <etal/></person-group>. (<year>2009</year>). <article-title>Poxvirus proteomics and virus-host protein interactions</article-title>. <source>Microbiol. Mol. Biol. Rev.</source> <volume>73</volume>, <fpage>730</fpage>&#x2013;<lpage>749</lpage>. doi: <pub-id pub-id-type="doi">10.1128/MMBR.00026-09</pub-id>, PMID: <pub-id pub-id-type="pmid">19946139</pub-id></citation></ref>
<ref id="ref48"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wallace</surname> <given-names>A. C.</given-names></name> <name><surname>Laskowski</surname> <given-names>R. A.</given-names></name> <name><surname>Thornton</surname> <given-names>J. M.</given-names></name></person-group> (<year>1995</year>). <article-title>LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions</article-title>. <source>Protein Eng. Des. Sel.</source> <volume>8</volume>, <fpage>127</fpage>&#x2013;<lpage>134</lpage>. doi: <pub-id pub-id-type="doi">10.1093/PROTEIN/8.2.127</pub-id></citation></ref>
<ref id="ref90"><citation citation-type="other"><person-group person-group-type="author"><collab id="coll22">WHO</collab></person-group> (<year>2022</year>). <article-title>Monkeypox. World Health Organization</article-title>. <ext-link xlink:href="https://www.who.int/news-room/fact-sheets/detail/monkeypox" ext-link-type="uri">https://www.who.int/news-room/fact-sheets/detail/monkeypox</ext-link> [Accessed 10th October 2022].</citation></ref>
<ref id="ref49"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Xiao</surname> <given-names>S.</given-names></name> <name><surname>Tian</surname> <given-names>Z.</given-names></name> <name><surname>Wang</surname> <given-names>Y.</given-names></name> <name><surname>Si</surname> <given-names>L.</given-names></name> <name><surname>Zhang</surname> <given-names>L.</given-names></name> <name><surname>Zhou</surname> <given-names>D.</given-names></name></person-group> (<year>2018</year>). <article-title>Recent progress in the antiviral activity and mechanism study of pentacyclic triterpenoids and their derivatives</article-title>. <source>Med. Res. Rev.</source> <volume>38</volume>, <fpage>951</fpage>&#x2013;<lpage>976</lpage>. doi: <pub-id pub-id-type="doi">10.1002/MED.21484</pub-id>, PMID: <pub-id pub-id-type="pmid">29350407</pub-id></citation></ref>
<ref id="ref50"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yu</surname> <given-names>M.</given-names></name> <name><surname>Si</surname> <given-names>L.</given-names></name> <name><surname>Wang</surname> <given-names>Y.</given-names></name> <name><surname>Wu</surname> <given-names>Y.</given-names></name> <name><surname>Yu</surname> <given-names>F.</given-names></name> <name><surname>Jiao</surname> <given-names>P.</given-names></name> <etal/></person-group>. (<year>2014</year>). <article-title>Discovery of pentacyclic triterpenoids as potential entry inhibitors of influenza viruses</article-title>. <source>J. Med. Chem.</source> <volume>57</volume>, <fpage>10058</fpage>&#x2013;<lpage>10071</lpage>. doi: <pub-id pub-id-type="doi">10.1021/JM5014067</pub-id>, PMID: <pub-id pub-id-type="pmid">25383779</pub-id></citation></ref>
<ref id="ref51"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>Y.</given-names></name> <name><surname>Moss</surname> <given-names>B.</given-names></name></person-group> (<year>1992</year>). <article-title>Immature viral envelope formation is interrupted at the same stage by lac operator-mediated repression of the vaccinia virus D13L gene and by the drug rifampicin</article-title>. <source>Virology</source> <volume>187</volume>, <fpage>643</fpage>&#x2013;<lpage>653</lpage>. doi: <pub-id pub-id-type="doi">10.1016/0042-6822(92)90467-4</pub-id>, PMID: <pub-id pub-id-type="pmid">1546459</pub-id></citation></ref>
</ref-list>
<fn-group>
<fn id="fn0004"><p><sup>1</sup><ext-link xlink:href="http://www.swissadme.ch" ext-link-type="uri">http://www.swissadme.ch</ext-link></p></fn>
<fn id="fn0005"><p><sup>2</sup><ext-link xlink:href="https://biosig.lab.uq.edu.au/pkcsm/prediction" ext-link-type="uri">https://biosig.lab.uq.edu.au/pkcsm/prediction</ext-link></p></fn>
<fn id="fn0006"><p><sup>3</sup><ext-link xlink:href="https://www.rcsb.org/" ext-link-type="uri">https://www.rcsb.org/</ext-link></p></fn>
<fn id="fn0007"><p><sup>4</sup><ext-link xlink:href="https://pubchem.ncbi.nlm.nih.gov/" ext-link-type="uri">https://pubchem.ncbi.nlm.nih.gov/</ext-link></p></fn>
<fn id="fn0008"><p><sup>5</sup><ext-link xlink:href="http://www.pymol.org/pymol" ext-link-type="uri">http://www.pymol.org/pymol</ext-link></p></fn>
<fn id="fn0009"><p><sup>6</sup><ext-link xlink:href="https://imods.iqfr.csic.es/" ext-link-type="uri">https://imods.iqfr.csic.es/</ext-link></p></fn>
<fn id="fn0010"><p><sup>7</sup><ext-link xlink:href="https://cbbio.online/LigTMap/" ext-link-type="uri">https://cbbio.online/LigTMap/</ext-link></p></fn>
</fn-group>
</back>
</article>