Investigating the Antiviral Properties of Nyctanthes arbor-tristis Linn against the Ebola, SARS-CoV-2, Nipah, and Chikungunya Viruses: A Computational Simulation Study

Background: The hunt for naturally occurring antiviral compounds to combat viral infection was expedited when COVID-19 and Ebola spread rapidly. Phytochemicals from Nyctanthes arbor-tristis Linn were evaluated as significant inhibitors of these viruses. Methods: Computational tools and techniques were used to assess the binding pattern of phytochemicals from Nyctanthes arbor-tristis Linn to Ebola virus VP35, SARS-CoV-2 protease, Nipah virus glycoprotein, and chikungunya virus. Results: Virtual screening and AutoDock analysis revealed that arborside-C, beta amyrin, and beta-sitosterol exhibited a substantial binding affinity for specific viral targets. The arborside-C and beta-sitosterol molecules were shown to have binding energies of −8.65 and −9.11 kcal/mol, respectively, when interacting with the major protease. Simultaneously, the medication remdesivir exhibited a control value of −6.18 kcal/mol. The measured affinity of phytochemicals for the other investigated targets was −7.52 for beta-amyrin against Ebola and −6.33 kcal/mol for nicotiflorin against Nipah virus targets. Additional molecular dynamics simulation (MDS) conducted on the molecules with significant antiviral potential, specifically the beta-amyrin-VP35 complex showing a stable RMSD pattern, yielded encouraging outcomes. Conclusions: Arborside-C, beta-sitosterol, beta-amyrin, and nicotiflorin could be established as excellent natural antiviral compounds derived from Nyctanthes arbor-tristis Linn. The virus-suppressing phytochemicals in this plant make it a compelling target for both in vitro and in vivo research in the future.


Introduction
Nyctanthes arbor-tristis (N.arbortristis) is a plant belonging to the Oleaceae family and the Nyctanthes genus.It is traditionally used for its therapeutic properties in tropical and subtropical regions of South and Southeast Asia.It grows in the outer Himalayas of India.It can be found in areas such as Jammu and Kashmir, Assam's east coast, Bengal, and Tripura, which stretch across the country's central region and down to the Godavari River in the south, as well as in Nepal.It is typically a tiny shrub or tree with beautiful, incredibly ease.The Ebola virus is a highly deadly infectious agent capable of infecting people.The Ebola virus protein VP35 hinders the synthesis of host interferon (IFN)-α/β by disrupting the host's immunological response to viral invasion.As a result, it is seen as a potential therapeutic target.The viral protein VP35 inhibits the activity of human-type I IFNs to suppress the response of dendritic cells (DCs) to the viral infection.This mechanism helps the virus avoid the buildup of signal transducers, specifically activators of transcription 1 (STAT1), in the infected cells.Moreover, VP35 inhibits dsRNA-dependent protein kinase receptor (PKR) activation, which is responsible for the production of IFN [10].
The glycoprotein NiV-G attached to the exterior of the infectious agent is crucial for identifying the ephrin-B2 (EFNB2) and ephrin-B3 (EFNB3) sensors on the surface of host cells.This interaction facilitates the binding of the pathogen to the cells, resulting in the fusion of their membranes and viral entry into the host.Therefore, the molecular suppression of NiV-G is a key approach for the control of NiV infection [11].The nsP2 cysteine protease has been identified as a vital target for drugs in the fight against illnesses associated with Alphaviruses, such as the chikungunya virus, since it plays a key role in viral replication.
Thus, in the search for antiviral molecules derived from Nyctanthes arbor-tristis, phytochemicals were analyzed against targets from different viruses through computational screening, interaction analysis, and molecular dynamics simulation (MDS) studies.

Results
We used computational techniques like molecular docking, ADMET analysis, and MDS for the study.The results are discussed in the following sections.

Docking Data Interpretation
After performing the AutoDock tool's docking analysis, we selected the top two natural compounds based on binding energies for further ADMET and interaction modeling analyses.The results showed that the natural compounds had interaction energies ranging between −9.11 and −3.73 kcal/mol compared with remdesivir (−6.8 kcal/mol) (see Table 1 and Supporting Information (Table S1)).
Table 1.Molecular interaction analysis details the screened tops of two natural compounds and the main protease of a SARS-CoV-2 crystal structure (PDB: 6LU7).In the "number of hydrogen bonds" column, UNK1 = ligand molecules selected for the interaction.The binding energy values for selected natural compounds arborside-C and betasitosterol with the main protease were found to be −8.85 and −9.11 kcal/mol, respectively.At the same time, the control drug, remdesivir, had a value of −6.18 kcal/mol.
Table 2. Molecular interaction analysis details the screened top natural compound with control drug remdesivir and Ebola virus VP35 bound to a small molecule (PDB: 4IBB).In the "number of hydrogen bonds" column, UNK = ligand molecules selected for the interaction.Beta-amyrin demonstrated a more favorable binding energy of −7.52 kcal/mol compared to remdesivir.It formed two hydrogen bonds.Beta-amyrin had a lower inhibition constant of 3.08 µM, suggesting potentially higher efficacy in inhibiting Ebola virus VP35.Amino acid residues Cys247, Lys248, Ile295, and Pro304 formed alkyl interactions, while amino acid residues Ala221, Gln241, Phe328, Asp302, Ile303, and Ile297 contributed to the binding via hydrophobic interactions (Figure 5A,B).Further, Table 3 details the molecular interaction between the control drug remdesivir and the natural compound nicotiflorin, which interacts with the Nipah G attachment glycoprotein (PDB ID: 3D11).
Table 3.Molecular interaction analysis details the screened top natural compound with the control drug remdesivir and the crystal structure of the Nipah G attachment glycoprotein (PDB ID: 3D11).In the "number of hydrogen bonds" column, UNK = ligand molecules selected for the interaction.Remdesivir showed a binding energy of −6.27 kcal/mol, indicating its effectiveness in binding with the protein.It formed five hydrogen bonds.The inhibition constant for remdesivir was 25.29 µM, suggesting moderate binding affinity.Hydrophobic interactions were formed by amino acid residues.Pro353, Gly352, Cys282, Tyr351, Phe458, Trp504, Leu305, Arg242, Asp219, Pro220, Val507, Thr218, Ser241, and Leu234 were involved in hydrophobic interactions.Other interactions also formed, like Pi-Alkyl with Tyr280 and Pro441, Pi-Sigma with His281, Pi-Cation with Lys560, Pi-Anion with Asp302, and attractive charge with Asp302 (Table 3).The natural compound nicotiflorin exhibited a slightly higher binding energy of −6.39 kcal/mol, which could imply a stronger binding affinity.It formed eight hydrogen bonds during the interaction.Nicotiflorin had a lower inhibition constant of 20.63 µM, further supporting the notion of it having a stronger binding affinity.Amino acid residues Asn586, Asp302, Gly301, Tyr280, Leu305, Thr218, and Pro220 were involved in hydrophobic interactions.Alkyl/Pi-Alkyl contact was formed by Ile588 and Tyr581 (Figure 6A,B).While both compounds exhibited significant interactions with the Nipah G attachment glycoprotein, nicotiflorin may offer enhanced binding characteristics, as evidenced by its higher binding energy and lower inhibition constant.Further investigation into nicotiflorin's antiviral potential is warranted.
The molecular interaction analysis summarized in Table 4 provides insights into the interactions between the natural compound beta-amyrin, the control drug remdesivir, and the chikungunya virus nsP2 protease (PDB ID: 3TRK).
Table 4. Molecular interaction analysis detailing the screened top natural compound with the control drug remdesivir and the structure of the chikungunya virus nsP2 protease (PDB ID: 3TRK).In the "number of hydrogen bonds" column, UNK = ligand molecules selected for the interaction.Remdesivir exhibited a binding energy of −5.48 kcal/mol and formed three hydrogen bonds.Its inhibition constant was 95.82 µM, suggesting a moderate binding affinity.The amino acid residues His1236, Ala1237, Lys1045, Tyr1047, Tyr1079, Leu1205, Trp1084, Pro1049, Lys1091, Glu1050, Arg1271, Thr1268, Leu1243, and Glu1204 were involved in hydrophobic interactions.Lys1239, Ala1046, and Val1051 were involved in Alkyl/Pi-alkyl interactions (Table 4).The natural compound beta-amyrin exhibited a significantly lower binding energy of −8.54 kcal/mol, indicating a potentially stronger binding affinity.It formed one hydrogen bond with a length of 2.11674 Angstrom.Beta-amyrin had a lower inhibition constant of 551.06 nM, suggesting that it might be more effective in inhibiting the chikungunya virus nsP2 protease.Amino acid residues Val1051, Gln1241, Leu1205, Gly1206, Glu1204, Leu1243, and Leu1203 created hydrophobic interactions, while Tyr1079, Trp1084, Ala1046, Tyr1047, Ala1040, Lys1239, and Lys1045 were involved in Alkyl/Pi-Alkyl interactions (Figure 7A,B).The data suggest that beta-amyrin could potentially be more effective than remdesivir in inhibiting the chikungunya virus nsP2 protease due to its lower Ki and higher binding energy.Both compounds exhibited hydrophobic interactions but differed in the residues involved, which could play a role in their respective efficacies and specificities.Further experimental validation is required to confirm these findings and explore potential therapeutic applications.

Inhibition
Overall, each compound interacted differently and efficiently with the viral proteins, and their binding energies and inhibition constants varied.All the selected compounds showed promising results, justifying further investigation for potential antiviral applications.

Drug-Likeness and ADMET Analysis
After examining several parameters, including blood-brain barrier (BBB) permeability, gastrointestinal absorption, inhibition of cytochrome, substrate interaction with P glycoprotein, and log Kp value for skin permeation, arborside-C showed significant results and nearly identical data to those of the control drug remdesivir based on the ADME data that were obtained for beta-sitosterol, arborside-C, nicotiflorin, and beta-amyrin (Figure 8A-D) from the SwissADME server (see Supporting Information (Table S2)).According to a drug-likeness analysis, the compound arborside-C showed two violations of Lipinski's rule of five requirements, similar to the control drug [15].Furthermore, the synthetic accessibility value and the bioavailability score fell within the predetermined cutoff range (see Supporting Information (Table S3)).
Regarding the Lipinski rule of five, arborside-C violated two criteria: it had a molecular weight greater than 500 g/mol and more than nine rotatable bonds.However, it did not infringe on any of the other criteria.Therefore, arborside-C was considered a "drug-like" molecule.
Arborside-C had a bioavailability score of 0.11, which indicated that it was poorly absorbed into the bloodstream.This is likely due to its high lipophilicity.Arborside-C also had a synthetic accessibility score of 6.14, which indicated that it was relatively easy to synthesize.
Beta-sitosterol did not violate any criteria of the Lipinski rule of five.Therefore, it was considered a "drug-like" molecule.
Beta-sitosterol had a bioavailability score of 0.55, indicating that it was well absorbed into the bloodstream.This is likely due to its low lipophilicity.Beta-sitosterol also had a synthetic accessibility score of 6.3, suggesting that it was relatively easy to synthesize.
Overall, remdesivir, arborside-C, and beta-sitosterol are all small-molecule drugs shown to have potential for treating COVID-19.Remdesivir and arborside-C are rela-tively easy to synthesize but have low bioavailability.Beta-sitosterol is more bioavailable but more challenging to synthesize.Further research is needed to determine which of these compounds is most effective and safe for the treatment of COVID-19.Additionally, the toxicity observation of this study was completed using the server pkCSM (http://biosig.unimelb.edu.au/pkcsm/theory,accessed on 1 March 2024).Arborside-C is a non-toxic compound.It met the requirements established by factors like AMES toxicity, liver toxicity, sensitization of skin, toxicity of T. pyriformis, and minnow toxicity (see Supporting Information (Table S4)).
Arborside-C is not mutagenic, as indicated by its negative AMES test result.It had an acceptable dose of −0.397 mg/kg/day in humans.Arborside-C is not a hERG I or hERG II inhibitor.It had an oral acute toxicity (LD50) of 2.701 mg/kg and oral chronic toxicity (LOAEL) of 4.723 mg/kg in rats.Arborside-C does not cause hepatotoxicity or skin sensitization.It had a T. pyriformis toxicity of 0.285 log µg/L and a minnow toxicity of 3.731 log LC50.
Beta-sitosterol has been shown to have several health benefits, including reducing cholesterol levels, improving heart health, and relieving the symptoms of benign prostatic hyperplasia (BPH).Beta-sitosterol is mutagenic above 0.559 mg/kg/day in humans, as reflected by its positive AMES test.Beta-sitosterol is not an hERG I or hERG II inhibitor.It has an LD50 of 2.482 mg/kg and a LOAEL of −0.515 mg/kg.Beta-sitosterol does not cause hepatotoxicity or skin sensitization.It has a T. pyriformis toxicity of 0.285 log µg/L and a minnow toxicity of 4.614 log LC50.
Overall, remdesivir and arborside-C are both relatively safe compounds.They are not mutagenic and do not cause hepatotoxicity or skin sensitization.However, they are both moderately toxic, with oral acute toxicity (LD50) values of 2.329 mg/kg and 2.701 mg/kg, respectively, in rats.Beta-sitosterol is less harmful than remdesivir or arborside-C, with an oral acute toxicity (LD50) value of 2.482 mg/kg in rats.Additional investigation is required to ascertain the extended-term safety of arborside-C and beta-sitosterol.
We also predicted the ADMET properties of beta-amyrin and nicotiflorin.Both betaamyrin and nicotiflorin were predicted to have low gastrointestinal absorption.They were predicted to be non-permeable to the BBB.Neither was expected to be a substrate for P-glycoprotein, a transporter protein that can limit drug absorption.Neither compound was predicted to inhibit any of the major enzymes of the cytochrome P450 (CYP) system involved in drug metabolism.The prediction of BBB non-permeability suggests that these compounds may not pass the BBB to reach the central nervous system.The value was −2.41, indicating low skin permeation.Negative log Kp values correspond to less skin permeation of the molecule.The value was significantly lower at −9.91, indicating even less skin permeation than beta-amyrin.This could be advantageous for some therapeutic applications, as it would reduce the risk of side effects on the central nervous system.
Beta-amyrin and nicotiflorin showed a greater binding affinity for other viral receptors and were also screened for drug-likeness using the SwissADME server.Beta-amyrin had a molecular weight (MW) of 426.72 g/mol and a total polar surface area (TPSA) of 20.23 A 2 .Nicotiflorin had an MW of 594.52 g/mol and a TPSA of 249.20 A 2 .Both compounds violated the Muegge size principle.Beta-amyrin had zero rotatable bonds and satisfied all other filtering criteria from the SwissADME server.Nicotiflorin had six rotatable bonds and violated filters for the Ghose, Veber, and Egan rules.The bioavailability score was 0.55 for beta-amyrin and 0.17 for nicotiflorin.
Both compounds exhibited distinct drug-likeness profiles.Beta-amyrin leaned toward lipophilicity and moderate oral absorption potential, while nicotiflorin was more polar, water-soluble, and challenging to synthesize.These insights can guide researchers in selecting promising drug candidates.Beta-amyrin and nicotiflorin were observed to be safe, as predicted by the results from the pkCSM server.A comparison of both compounds revealed no AMES mutagenicity.
The maximum tolerated dose (human) of beta-amyrin was −0.312, and that of nicotiflorin was 0.221.Nicotiflorin had a slightly higher tolerated dose.Neither compound inhibits hERG I.Both compounds inhibit hERG II.
The oral acute toxicity (LD50) parameter value for beta-amyrin was 2.296 in rats, suggesting a relatively low acute toxicity, while that of nicotiflorin was 2.518.Both compounds had comparable acute toxicity.Beta-amyrin's LOAEL value was 0.922, indicating chronic exposure tolerance, while nicotiflorin's LOAEL value was significantly higher at 5.488, suggesting better tolerance in chronic exposure scenarios.Nicotiflorin showed better chronic tolerance.Neither compound exhibited hepatotoxicity, and both compounds were non-sensitizers.The T. pyriformis toxicity value for beta-amyrin was 0.388, which is low.Nicotiflorin had a value of 0.285, slightly lower than that of beta-amyrin.Thus, both compounds had low toxicity in T. pyriformis.
The minnow toxicity value for beta-amyrin was −1.747, indicating low toxicity to minnows.Nicotiflorin's value was significantly higher at 7.767, suggesting significant tolerance to minnows.Thus, nicotiflorin showed better tolerance to minnows.
Both compounds exhibited low acute toxicity, were not hepatotoxic, and did not cause skin sensitization.Nicotiflorin showed higher tolerance in chronic exposure scenarios and significantly lower minnow toxicity.

MDS Analysis
After a successful run, the XMGRACE tool was used to generate GROMACS trajectory files [16] (Turner, 2005).The developed 2D plots were analyzed for the radius of gyration (Rg), root mean square deviation (RMSD), root mean square fluctuation (RMSF), and hydrogen bond during a 100-ns simulation.The beta-sitosterol-protease complex showed higher fluctuation than other complexes.
The deviation of the protease simulation in water and protease simulation in the presence of arborside-C, beta-sitosterol, and remdesivir was between 0.15 and 0.25 nm (Figure 9A).The Arborside-C_Protease complex showed a similar value (0.15) to that of the protease simulation in water and a lower RMSD value than the remdesivir-protease simulation value (more than 0.2 nm).
RMSF fluctuation values for complexes were between 0.1 and 0.35 nm (Figure 9B).All the chosen compounds had average RMSF values of around 0.1 nm at the 50-55, 120, 150-160, 188-190, 225-230, and 275-280 amino acid residue regions.A protein's tertiary structure's compactness and stability are indicated by its Rg values.When testing the stability of proteins in the presence of certain chemicals, this computation is essential, which was observed between 2.2 and 2.25 nm, as shown in the plot.At the same time, protease simulation in water and arborside-C and remdesivir presence in the system displayed numbers that were quite close (≥2.2nm), while beta-sitosterol showed values near 2.25 nm (Figure 9C).
Hydrogen bonds are crucial in the interaction between ligands and proteins [17].We generated a hydrogen bond plot.The Remdesivir_Protease complex formed eight hydrogen bonds compared to the Arborside-C_Protease complex, which included a maximum of seven bonds.In contrast, beta-sitosterol formed one hydrogen bond during the simulation period (Figure 9D).
Both beta-amyrin_VP35 and remdesivir_VP35 exhibited fluctuations in RMSD over time (Figure 10A).The RMSD values remained below 0.15 nm throughout the 100-ns simulation.Both compounds maintained structural stability during the simulation.Remdesivir_VP35 showed more significant fluctuations in RMSF than beta-amyrin_VP35, with an average value between 0.05 and 0.1 nm (Figure 10B).Remdesivir and beta-amyrin experienced flexibility or variability in specific amino acid residues and showed a similar pattern throughout the simulation.Further, both compounds exhibited identical fluctuations in the Rg over time.Rg values remained between 1.46 nm and 1.48 nm (Figure 10C).The overall compactness of the protein structures was maintained during the simulation.Remdesivir_VP35 consistently exceeded the formation of hydrogen bonds 1-9, while beta-amyrin_VP35 formed 1-2 hydrogen bonds (Figure 10D).Remdesivir formed more hydrogen bonds with the protein, potentially influencing its stability and interactions.
We also analyzed similar MDS parameters to those of the MERS-CoV-2 protease receptor for Nicotiflorin_Glycoprotein and Remdesivir_Glycoprotein complexes of the Nipah glycoprotein.Both compounds exhibited similar behavior regarding RMSD over time.The RMSD values remained below 0.2 nm throughout the 100-ns simulation (Figure 11A).Both Nicotiflorin_Glycoprotein, shown at approximately 0.15 nm, and Remdesivir_Glycoprotein at 0.1 nm maintained structural stability during the simulation.As shown in Figure 11B, the Nicotiflorin_Glycoprotein complex showed slightly higher fluctuations in RMSF across various amino acid residues compared to the Remde-sivir_Glycoprotein.However, both complexes showed similar average values of 0.1 nm until the end of the simulation.Nicotiflorin_Glycoprotein exhibited greater flexibility or variability in specific protein structure regions (470-580 amino acids).
Both compounds exhibited similar fluctuations in Rg over time.Rg values remained between 2.10 nm and 2.13 nm (Figure 11C).The overall compactness of the protein structures was maintained during the simulation.Both Nicotiflorin_Glycoprotein and Remde-sivir_Glycoprotein formed a significant number of hydrogen bonds.A total of 1-9 hydrogen bonds formed during the 100 ns simulation (Figure 11D).The fluctuations in hydrogen bonding patterns suggested dynamic interactions.Both compounds engaged in essential hydrogen bonding with the surrounding environment.
Overall, both compounds exhibited stable structures (low RMSD) and maintained compactness (consistent Rg).Nicotiflorin_Glycoprotein showed greater flexibility (higher RMSF) but formed a comparable number of hydrogen bonds.Remdesivir_Glycoprotein maintained stability and formed a similar number of hydrogen bonds.The differences in RMSF and hydrogen bonding may influence the binding affinities and interactions with the protein.In summary, Nicotiflorin_Glycoprotein and Remdesivir_Glycoprotein interacted differently with the protein, and their distinct features may impact their antiviral potential.Further studies are needed to understand their mechanisms and fully guide drug development.
Furthermore, we analyzed the time evolution of various properties of beta-amyrin's and remdesivir's interactions with chikungunya protease during MDS.
The RMSD of Beta-Amyrin_Protease and Remdesivir_Protease relative to their initial structures over time (ns) were assessed.The RMSD of both complexes increased slowly over the course of the simulation, with remdesivir showing a slightly lower RMSD than beta-amyrin, with values ranging from 0.2 to 0.4 nm.This suggests that beta-amyrin undergoes slightly larger structural fluctuations than remdesivir during the simulation (Figure 12A).The RMSF of both molecules fluctuated around 0.1-0.2nm throughout the simulation.This suggests that beta-amyrin and remdesivir have relatively stable backbone structures with minor fluctuations (Figure 12B).The Rg of Beta-Amyrin_Protease and Remdesivir_Protease over time (100 ns) was assessed.The Rg of both molecules fluctuated slightly over the course of the simulation, with remdesivir having a marginally higher average Rg compared to beta-amyrin at the end of the simulation, with a value between 2.1 and 2.2 nm.This suggests that remdesivir and beta-amyrin may have adopted a more extended conformation during the simulation (Figure 12C).Hydrogen bonds formed in the interaction between the protein and beta-amyrin and protease and remdesivir protease over time (ns).The number of hydrogen bonds between remdesivir and the protein fluctuated around 1-6 throughout the simulation, while the number of hydrogen bonds between beta-amyrin and the protein fluctuated around 2-3 (Figure 12D).Overall, the results suggest that both beta-amyrin and remdesivir interact with the protein target.However, beta-amyrin appears to form a more stable interaction with the protein, likely due to a combination of factors such as smaller structural fluctuations and a higher number of hydrogen bonds.
MDS trajectory files were employed in the MMPBSA analysis process.Tables 5 and 6 contain the results of the Poisson-Boltzmann complex energy and ligand-receptor energy component calculations, respectively.The total free binding energy for the given complexes is based on the provided table.Each component contributes to the overall binding energy and understanding these values sheds light on the stability and interactions within the complexes, as shown in Table 5.
Here, ∆Vdwaals = Van der Waals energy; ∆EPB = Polar contribution to the solvation energy; ∆EEL = Electrostatic molecular energy; ∆ENPOLAR = Nonpolar contribution of repulsive solute-solvent interactions to the solvation energy; ∆G Gas = Total gas phase molecular energy; ∆EDISPER = Nonpolar contribution of attractive solute-solvent interactions to the solvation energy; ∆G Solv = Total solvation energy, and ∆G Total = Total binding energy.
The total binding energy (∆GTotal) was negative for all complexes, indicating that the binding was favorable.However, the binding energy's efficacy varied between the complexes.

Discussion
Beta-amyrin showed binding energy, suggesting interactions with the chikungunya protease.It can be concluded that the presence of arborside-C and remdesivir in the target system stabilizes the protein, as evidenced by the lower RMSD and Rg values.In contrast, beta-sitosterol destabilizes the protein, as evidenced by the higher RMSD and Rg values.The hydrogen bond analysis also supports this conclusion.The Remdesivir_Protease complex formed eight hydrogen bonds, while the Arborside-C_Protease complex formed a maximum of seven bonds.In contrast, the Beta-Sitosterol_Protease complex formed only one hydrogen bond.These results suggest that arborside-C and remdesivir may be more effective protease inhibitors than beta-sitosterol.Further research is needed to confirm these findings.
More and more plant-based medications are used because of their therapeutic benefits and essential roles in achieving positive health outcomes [18].Many plant compounds have been analyzed through different approaches to explore their interacting potential with viral proteases [5].Numerous medicinal plants are rich in flavonoids and polyphenolic chemicals that have been shown to have anti-protease efficacy against multiple infections like hepatitis, flu, HIV, cytomegalovirus, herpes, smallpox, and Ebola [19][20][21].Their antiviral activity could involve inhibiting reverse transcriptase enzyme proteases and DNA polymerase.
Compared to the native ligand, Mpro and other natural compounds studied, such as 2,2-di(3-indolyl)-3-indolone, terrequinone A, and betulinic acid, have demonstrated a greater affinity for binding.An anthraquinone from the observed list, alterporriol-Q, was discovered to be the most effective inhibitor of SARS-CoV-2 Mpro.Additionally, MD modeling experiments for the Mpro complex with alterporriol-Q indicate that alterporriol-Q does not significantly change the structure of Mpro [22][23][24].
In diverse Nyctanthes arbortristis tissues, steroids, flavonoids, terpenes, aliphatic chemicals, and alkaloids were found-the two primary types of chemicals made by this medicinal plant are glycosides and alkaloids.Several studies have focused on plant extracts from different parts of the plant for their therapeutic effects, such as anti-inflammatory, anti-allergic [25], antibacterial, antioxidant [26], anti-diabetic, anti-cancer, antiviral, and antimalarial properties [27].The antiviral potential of arbortristoside A and C molecules against encephalomyocarditis virus and Semliki forest virus has been explored [28].The studied bioactive molecules were observed to be more potent than the antiviral medication remdesivir.Arbrorside-C and beta-sitosterol have very good interaction energies, and these molecules significantly decreased protease-mediated infectivity.Recently, dietary bioactive compounds from three ayurvedic medicinal herbs, vasicoline (∆G = −7.34kcal/mol), solanocapsine (∆G = −8.59kcal/mol), and 18 β-glycyrrhetinic acid (∆G = −8.86kcal/mol) were reported to have significant binding interactions to Mpro enzyme compared to the native N3 inhibitor (∆G = −5.41kcal/mol) remdesivir.This supported our observed results, as the binding energies were much higher than those of the reported flavonoids.The binding of the molecules could be at the receptor binding domain sites at the His41 and Cys145.Quercetin and protease interactions were reported to be influenced by the presence of the 7-OH group and the acetylation of the sugar moiety.Moreover, the binding was based on their computer study of catechin gallate and quercetin 3-O-malonylglucoside, which are reported to inhibit the main protease by interacting with the active sites of the target protein (His41, Met49, Cys145, Met165, and Thr190) [29].Researchers have recently explored the antiviral potential of beta-amyrin.This study has also shown that beta-amyrin interacted with viral protease, albeit with less potency than VP35, 3tkr chikungunya protease.Nicotiflorin interacts strongly with the Nipah G attachment glycoprotein.This suggests that nicotiflorin has medical uses, such as cardiac activity, which could be important in protecting the life of a heart patient who is infected with a virus [30,31].
Drug-likeness, molecular dynamics, and ADMET analysis of selected compounds represented important data more significantly than the remdesivir.The docking ADME results indicate the therapeutic significance of the studied antiviral molecules, which need further investigation.

Ligand Preparation
The 2D conformational information, including SMILES IDs, of the most significant 26 natural compounds of the plant N. arbor-tristis and control drug remdesivir was downloaded from the PubChem database [32] and utilized by converting it into 3D structural files for further molecular docking and simulation purposes.Furthermore, generated ligand files were submitted to the Discovery Studio Visualizer version 24.1.0.23298 for the energy minimization process using the CHARMm forcefield [33], which signifies the functional empirical energy to treat macromolecular systems as a model [34][35][36][37].

Viral Receptor Preparation
The target COVID-19 main protease (PDB ID: 6LU7) [38], Ebola virus VP35 (PDB ID: 4IBB) [39], and Nipah G attachment glycoprotein (PDB ID: 3D11) [40] and the structure of the chikungunya virus nsP2 protease (PDB ID: 3TRK) [41] were downloaded from the Protein Data Bank (PDB) database [42].We removed the published structures' HETATM and water molecules and applied the CHARMm forcefield [33].The native 3D structural conformation of selected receptors was also examined, and amino acid residues of active sites involved in binding with pre-bounded ligands were discovered.These amino acid residues were then considered for docking studies.With the aid of Discovery Studio Visualizer version 24.1.0.23298,PDB files were rendered in 3D [34,35].

Molecular Docking Studies
To estimate the binding significance between the receptor structures of the different viruses and natural compounds, the MGL tool AutoDock 4.2, the Lamarckian genetic algorithm (LGA), and energy ∆G were calculated, as described by Morris [43,44]: Here, the dispersion of two Gaussian functions in this instance has an attractive term of ∆G gauss.Repulsion: if the distance is less than a certain threshold, the square of the distance Ramp function for ∆Ghbond: also applied to interactions between metal ions ∆Ghydrophobic: ramp characteristic, ∆Gtors: represents the bond numbers that can be rotated.
The target molecules were dehydrated, and Kollman united, Gasteiger charges, and hydrogen were added to perform a docking with a point spacing of 0.375 by default in the grid box specification of 60 × 60 × 60 • for X, Y, and Z with a grid center: for COVID-19 protease (PDB ID: 6LU7), −13.721, 16.337, and 69.721; for Ebola virus VP35 (PDB ID: 4IBB), 1.637, 21.039, and 1.327; for the crystal structure of the Nipah G attachment glycoprotein (PDB ID: 3D11), 35.436, 11.123, and 102.638; and for the structure of the chikungunya virus nsP2 protease (PDB ID: 3TRK), 12.76, 21.994, and 23.723, respectively.The parameters like the number of individuals in the population for Genetic Algorithm (150), the maximum number of energy evaluations (2,500,000), the maximum number of generations (27,000), the number of top individuals to survive to the next generation (1), the rate of gene mutation (0.02), the rate of crossover (0.8), the iterations of Solis and Wets local search (300), the consecutive successes before changing rho (0.1), the consecutive failures before changing rho (4), the probability of performing a local search on an individual (0.6), and the total 10 run were finally set for execution.The last obtained 3D conformations of the complexes were visualized in Discovery Studio Visualizer version 24.1.0.23298 [34,45].

Drug-Likeness and ADMET
The SwissADME and pkCSM online servers were used (http://www.swissadme.com,accessed on 1 March 2024), and drug-likeness and ADMET analyses of selected natural compounds in the computational mode were completed through the online resource provided by the Swiss Institute of Bioinformatics (SIB), Lausanne, Switzerland [45][46][47].

MDS
Natural compound-receptor complexes' docking interactions were further assessed through MDS analysis; consequently, the MDS environment was configured to run simulations for these complexes for 100 ns using GROMACS [48].
The pdb2gmx module was used to create receptor topology files, and then CHARMM27 all-atom force fields were chosen.In the following step, the SwissParam server-generated topology files for the selected natural compounds and control drug remdesivir were used [49].A water-filled triclinic box model of the unit cell was constructed for the solvation process.The network was stabilized after including Na + and Cl − ions, leading to a decrease in potential.After reaching harmony, the system (chosen complexes) was subjected to the NVT and NPT two-step arrays (constant particles, temperature, and pressure).Both ensembles' pressure and temperature coupling can be controlled, allowing for more accurate simulations and excellent system stability [50].The RMSD calculation tool integrated into the GROMACS, that is, the gmx rms package, was utilized to evaluate complex trajectory datasets [51], gmxrmsf for the calculation of RMSF, gmxgyrate for the Rg [51], and gmxhbond to determine how fast hydrogen bonds form throughout the exchange.The Xmgrace tool was used to make 2D plots [16].
Furthermore, the Molecular Mechanics-Poisson-Boltzmann Surface Area (MM-PBSA) [52] method-based free energy calculation was executed after utilizing the MMPBSA.pyPython script developed by Miller et al. [53].The free energy calculation was finally run by the gmx_MMPBSA program [54].

Conclusions
Effective antiviral therapies are urgently needed for the management of viral diseases due to viruses' high transmission rates, mutagenicity, and genomic reorganization.The major Mpro enzymes, helicases, RNA-dependent RNA polymerase, and glycoproteins are examples of viral molecules that can be targeted with antiviral drugs to stop the replication and spread of viral infections.Viruses are extremely mutagenic, although any genetic change to the highly conserved protease enzyme frequently has adverse effects on them.Thus, using protease enzyme inhibitors may lower the probability of mutation-mediated medication resistance and offer an efficient antiviral defense.Several currently available antiviral medications and dietary bioactive agents are used against viruses.We used computational tools and methods to determine how arborside-C, beta-sitosterol, and betaamyrin from Nyctanthes arbortristis bound to targets on Ebola virus VP35, SARS-CoV-2 protease, Nipah virus glycoprotein, and chikungunya virus.According to ADMET data, the characteristics of arborside-C and the reference medication remdesivir are similar.For the complexes of phytochemicals with targets, additional MDSs produced encouraging outcomes.RMSD, as observed in nm, signified the stability of the complexes.Based on our data, the phytochemicals of Nyctanthes arbor tristis could be recognized as excellent broad-spectrum natural antiviral chemicals.Our findings suggest that these substances could be used as drugs to treat viral diseases of various origins.These chemicals may be candidates for treatment regimens following in vitro and in vivo clinical trials.

Figure 1 .
Figure 1.(A) Interacting site of arborside-C (central green ball and stick pattern) with the SARS-CoV-2 main protease (surrounding grey flat ribbon pattern).The blue stick pattern shows the interacting amino acid residues of the binding pocket, and green dotted lines show hydrogen bond formation.(B) The 2D graphic of amino acid residues involved in the interactions, shown by spheres.

Figure 2 .
Figure 2. (A) Interacting site of the binding site of beta-sitosterol (central ball and stick pattern) with the SARS-CoV-2 main protease (surrounding grey flat ribbon pattern).The blue stick pattern shows the interacting amino acid residues of the binding pocket, and green dotted lines show hydrogen bond formation.(B) The 2D graphic of amino acid residues involved in the interactions, shown by spheres.

Figure 3 .
Figure 3. (A) The interacting site of remdesivir (central light blue ball and stick pattern) with the SARS-CoV-2 main protease (surrounding grey flat ribbon pattern).The blue stick pattern shows the interacting amino acid residues of the binding pocket, and green dotted lines show hydrogen bond formation.(B) The 2D graphic of amino interactions, shown by spheres.

Figure 4 .
Figure 4. 3D visualization of the SARS-CoV-2 main protease (solid grey surface) binding pocket containing arborside-C (green), beta-sitosterol (red), and the control drug remdesivir (blue), shown as ball and stick patterns in the center.

Figure 5 .
Figure 5. (A) The interacting site of beta-amyrin (central purple ball and stick pattern) with Ebola virus VP35 bound to a small molecule (PDB ID: 4IBB, surrounding grey flat ribbon pattern).The blue stick pattern shows the interacting amino acid residues of the binding pocket, and green dotted lines show hydrogen bond formation.(B) The 2D graphic of amino interactions, shown by spheres.

Figure 6 .
Figure 6.(A) The interacting site of nicotiflorin (central olive ball and stick pattern) with Nipah G attachment glycoprotein (PDB ID: 3D11, surrounding grey flat ribbon pattern).The blue stick pattern shows the interacting amino acid residues of the binding pocket, and green dotted lines show hydrogen bond formation.(B) The 2D graphic of amino interactions, shown by spheres.

Figure 7 .
Figure 7. (A) The interacting site of beta-amyrin (central purple ball and stick pattern) with the chikungunya virus nsP2 protease (PDB ID: 3TRK, surrounding grey flat ribbon pattern).The blue stick pattern shows the interacting amino acid residues of the binding pocket, and green dotted lines show hydrogen bond formation.(B) The 2D graphic of amino acid interactions, shown by spheres.

Figure 9 .
Figure 9.The plots generated by trajectory files.(A) RMSD plot of Arborside-C_Protease (green), Beta-Sitosterol_Protease (red), and Remdesivir_Protease (blue) complexes.(B) Root mean square fluctuation (RMSF) plot illustrating the volatility of each amino acid residue.(C) The radius of gyration (Rg) plot explaining the protease's compactness, tightening in the presence of arborside-C, beta-sitosterol, and remdesivir.(D) Plot representing the interacting hydrogen bonds during a 100-ns simulation of selected complexes.

Figure 10 .
Figure 10.Plots generated by trajectory files.(A) RMSD plot of Beta-Amyrin_VP35 (green) and Remdesivir_VP35 (blue) complexes.(B) RMSF explaining each amino acid residue fluctuation.(C) Rg plot representing Ebola VP35 compactness, tightening in the presence of beta-amyrin and remdesivir.(D) Plot showing hydrogen bond interactions during a 100-ns simulation of selected complexes.

Figure 11 .
Figure 11.Plots generated by trajectory files.(A) RMSD graph of Nicotiflorin_Glycoprotein (pink) and Remdesivir_Glycoprotein (blue) complexes.(B) RMSF graph representing the fluctuation per amino acid residue.(C) Rg plot showing the Nipah glycoprotein tightening and compactness in the presence of nicotiflorin and remdesivir.(D) The number of hydrogen bond interactions during a 100-ns simulation of selected complexes.

Figure 12 .
Figure 12.Plots generated by trajectory files.(A) RMSD plot of Beta-Amyrin_Protease (green) and Remdesivir_Protease (blue) complexes.(B) RMSF plot showing the per amino acid residue fluctuation.(C) Rg graph showing the chikungunya protease tightening and compactness in the presence of beta-amyrin and remdesivir.(D) Interacting hydrogen bonds during a 100-ns simulation of selected complexes.