Pharmaco-phytochemistry of esters isolated from leaf extracts of Spondias mombin as potential antiviral agents

The present work reports on the isolation and characterisation of two novel antiviral ester compounds from dichloromethane leaf extracts of Spondias mombin ( SM ). The characterisation and structural elucidation were established from spectroscopic evidence of nuclear magnetic resonance, Fourier transform infrared and mass spectroscopy. The compounds identified were 6-methylheptyl pentadecanoate and 6-methylheptyl-15-(1,2,3,4,4a,8a-hexahydronaphthalen-1-yl)pentadecanoate. The novel isolated ester compounds were reported to have anti-rhino virus activity in silico against a known biological target (HsNMT1) that plays a key role in developing therapeutics against the common cold. Molecular docking analysis revealed the binding affinity across all targets within the range of − 4.6 to − 8.2 kcal/mol, whilst molecular dynamic simulation showed that systems attained good stability due to the maintenance of mean root-mean-square deviation values within the acceptable range of 1.5 − 2.5 Å. It can be concluded that the novel compounds are potential inhibitory candidates against rhinovirus protein target HsNMT1. However, in vitro and in vivo experiments are further required to validate the possible inhibitory candidates against rhinovirus disease (common cold).


Introduction
Novel phytochemical compounds have diverse phytochemical and pharmacological properties.These phytochemicals are abundant in natural products 1 , with some used as new drugs, leading to drug discovery 2 .4][5] Natural products, such as Spondias mombin (Anacardiaceae; SM) leaf extracts, have been used to treat several infectious diseases. 6The pharmacological activities of leaf extracts of S. mombin have been attributed to some bioactive compounds isolated from the medicinal plant.Some of these isolated compounds from leaf extracts of S. mombin include esters such as 3β-olean-12-en-3-yl (9Z)-hexadec-9-enoate 7 , chlorogenic acid butyl ester 8 and caffeoyl ester 9,10 .The leaves are reported to be part of the medicinal plant, and are mostly used for treating viral respiratory infections such as rhinovirus in traditional African healing systems. 11man rhinoviruses affect the upper and lower respiratory tract and cause common colds associated with pneumonia, wheezing and asthma. 12olated compounds from natural products are identified by several chromatographic and spectroscopic methods, such as thin layer chromatography (TLC), column chromatography (CC), Fourier transform infrared (FTIR) spectroscopy, gas chromatography and mass spectroscopy (GC-MS) and nuclear magnetic resonance (NMR), among others.These methods for the separation, identification and structural determination of phytochemicals are becoming increasingly powerful. 13TLC might be the simplest of all chromatographic methods, but it provides critical information for identifying compounds separated by other methods during the phytochemistry analysis of natural products. 14IR is known to identify only the types of functional groups in a compound, most commonly CH 2 , CH 3 , =CH, ≡CH, O-H, C=O, C-O, C=C, C≡C, C-O-C and C-C-O. 15FTIR and NMR analysis, coupled with mass spectroscopy (MS), are helpful tools in the structural elucidation of an isolated compound.
In order to assess the therapeutic potential of the selected compounds via in silico methods, it was necessary to identify a peculiar biological target.The prediction of the biological activities of the compounds was determined by utilising the PASSonline software. 16This software is used for the prediction of different physiological activities for multiple compounds, both natural and synthetic, based on their chemical formula.Additionally, PASSonline predicts pharmacological effects, mechanisms of action, adverse effects, interaction with metabolic enzymes and transporters, and influence on gene expression.It uses the 2D molecular fragments known as multilevel neighbours of atom descriptors, which postulates that a compound's molecular structure may determine whether it is biological. 17With this software, the evaluated activity of a compound is estimated as probable activity (Pa) and probable inactivity (Pi). 18 aimed to isolate, identify, characterise and predict the antiviral properties of novel esters using computer-aided drug design methods on novel inhibitors of human rhinovirus protease through molecular docking and molecular dynamic simulation.The results reveal, for the first time, two ester compounds from S. mombin leaf extracts that possess antiviral properties in an in silico molecular target prediction.

Material processing and extraction
Fresh leaves of S. mombin Linn were collected from Cape Coast, Ghana, and authenticated.Leaves dried at room temperature were pulverised by using a hammer mill.The leaf powder of mass 100 g was initially defatted with 1 L of hexane and then extracted with 2 L DCM by cold maceration for 72 h until the solvent was clear.The extracts were filtered with filter paper and concentrated using a rotary evaporator under reduced pressure at 40 °C.The concentrate was completely dried with a weight of 5.19 g (5.17%) and denoted as SM-DCM.

General analytical information
1 H and 13 C NMR spectra were recorded on a Bruker AV 400 MHz instrument at 400 MHz ( 1 H NMR) and 100 MHz ( 13 C NMR).All 1 H NMR spectra were measured in parts per million (ppm) downfield or relative to the residual proton signals of d1-chloroform (CDCl 3 ,7.26ppm).All 13 C NMR spectra were reported in ppm relative to the residual carbon signals of CDCl 3 (77.16ppm).Coupling constants (J) are reported in Hertz (Hz).Multiplicity is indicated as follows: s (singlet), d (doublet), t (triplet), q (quartet), p (pentet) and m (multiplet). 15Thin layer chromatography (TLC) was performed on precoated Merck Silica gel 60 F254 plates using different polarities of hexane-ethyl acetate solvent systems, and compounds were visualised with UV light at 254 nm. 19he retention factor (R f ) values of the different spots that were observed were calculated. 20e R f values were calculated using Equation 1 FTIR spectroscopy was performed using a PerkinElmer Spectrum 100 spectrophotometer at room temperature, whilst the MS of isolates was determined using a CombiFlash Purlon Mass Spectrometer (2000 Da Polarity Auto Switching).

Preliminary phytochemical screening of SM-DCM extract
Preliminary phytochemical screening was performed as a qualitative process to investigate the presence of different classes of phytochemicals according to standard procedures as reported by other authors. 21,22Briefly, crude DCM leaf extract of S. mombin was used for the qualitative analysis to determine the presence of alkaloids, steroid flavonoids, saponins, terpenoid tannins, anthraquinone derivatives and cardiac glycosides.

Column fractionation of SM-DCM mixture
The S. mombin DCM extract was loaded onto a glass column packed with silica gel.It was then eluted with mixtures of ethyl acetate and hexane of increasing gradient polarity, starting with 100% hexane to 100% ethyl acetate.One hundred and seventy-one fractions were collected in 50 mL aliquots.Based on their TLC analysis, aliquots 55-100 with a similar R f value (retention factor), were bulked together (denoted C) for further separation of the two compounds.It was further separated using column chromatography with silica gel using a solvent mixture of gradient, ethyl acetate and hexane.Seventy fractions were collected in 10 mL aliquots, and based on the TLC, aliquots 1-23 were bulked into CS1 and 24 to 70 into CS2.

Characterisation of isolated compounds
Chemical shifts are reported for DSS-trimethyl singlet resonance at 0.0000 ppm and multiplicity.

Biological activity prediction via PASSonline
PASSonline software 16 is used to predict physiological activities, pharmacological effects, mechanisms of action, toxic and adverse effects, interaction with metabolic enzymes and transporters and influence on gene expression for multiple compounds, both natural and synthetic, based on their chemical formula.
The evaluated activity of a compound is estimated as Pa and Pi. 18ompounds presenting Pa higher than Pi relative to a particular activity are considered feasible for that specific medical activity, and those with Pi higher than Pa are, therefore, eliminated.To this end, the selected compounds were assessed for their biological activities on PASSonline.

Molecular docking
The X-ray crystal structures of some selected rhinovirus antiviral targets (PDB ID: 5FX6, 5MU6, 4C2X and 1CQQ) 23 were retrieved from the Protein Data Bank 24 .These structures were co-crystallised with native inhibitors that defined their respective binding sites.The structures, 5FX6, 5MU6, 4C2X and 1CQQ, were then prepared by using UCSF Chimera version 1.13.1 25 to remove all non-standard residues, and Modeller 9.25 version 26 was employed to fix missing residues.The binding site residues were obtained by zoning the native inhibitors and selecting residues that lie within 5 Å for each target protein.Subsequently, the selected compounds were optimised using Avogadro 2.0 software and saved.Molecular docking was carried out for the three selected compounds against each of the rhinovirus target proteins using Prix software.The target that showed the best docking properties against all compounds was selected for molecular dynamics simulation.

Molecular dynamic simulation
Molecular dynamic (MD) simulations were performed using the AMBER18 GPU package for the best-docked ligand, CS1 and CS2, and IMP-1088 for the target (HsNMT1).The ligand and receptor were both defined and optimised using the AMBER force fields using the Antechamber and LEAP modules, respectively.Solvation and neutralisation were carried out on the receptor prior to its combination with the ligand.Partial minimisation of the receptor in the system was conducted for 2500 steps with a restraint potential of 500 kcal/mol Å2, followed by complete minimisation of 10 000.The system underwent heating at 300 K using a Langevin thermostat in a canonical ensemble (NVT).Equilibration of the system was carried out to ensure that AMBER rechecks the system, and it was at 300 K. MD simulation was run for 12 h at 100 ns, and results were obtained in the form of trajectories and analysed using statistics.The trajectories generated allow for the measurement of the binding energies of the association of the ligand with the receptor.The visualisation of the interactions was produced by Snapshots and Discovery Studio.

Binding free energy analysis via MM/GBSA method
The molecular mechanics/generalised-born surface area (MM/ GBSA) 27,28 method was employed in estimating the binding free energy for each of the inhibitor-bound systems.The binding free energy (ΔG bind ) was calculated from the following equations: ΔG bind = E gas + ΔG sol − TS, Equation 3where ΔG bind is the summation of the gas phase and solvation energy terms less the entropy (TS) term.
E gas is the sum of the AMBER force field internal energy terms E int (bond, angle and torsion), the covalent van der Waals (E vdw ) and the non-bonded electrostatic energy component (E lec ).The solvation energy is calculated from the following equation: where ΔG bind is taken to be the sum of the gas phase and solvation energy terms less the entropy (TΔS) term and G Complex represents the energy of the receptor-ligand complex, whilst G receptor and G ligand represent the energies of receptor and ligand, respectively.E gas denotes gasphase energy; E int signifies internal energy; and E ele and E vdw indicate the electrostatic and van der Waals contributions, respectively.E gas is the gas phase, elevated directly from the FF14SB force terms.G asol denotes solvation-free energy and can be decomposed into polar and non-polar contribution states.The polar solvation contribution, G GB , is determined by solving the GB equation, whereas G SA , the non-polar solvation contribution, is estimated from the solvent-accessible surface area (SASA) determined using a water probe radius of 1.4 Å. T and S correspond to temperature and total solute entropy, respectively.γ is a constant. 29Per-residue decomposition analyses were also carried out to estimate the individual energy contribution of residues in the substrate pocket towards the affinity and stabilisation of each target.

Preliminary phytochemical screening and TLC test of compounds CS1 and CS2
Table 1 indicates a phytochemical test to examine the qualitative chemical constituents contained in leaf extracts of S. mombin.The phytochemical test revealed the presence of anthraquinone derivatives, steroids, tannins and cardiac glucosides.The results of the preliminary phytochemical screening are in line with reports by other authors. 30,31e percentage yield in this study was calculated using the weight of the extracted sample divided by the total sample used and found to be 11.14% (Table 2).
During the TLC test, the retention factor values obtained in this experiment (Table 3) did not give many clues as to the type of compounds in the extract, but they did suggest the polarity of the compounds as reported by Talukdar et al. 14 The authors indicated that a high R f value in a less polar solvent system possesses low polarity. 14

Structural elucidation of 6-methylheptyl pentadecanoate (CS1)
The FTIR spectrum was used to identify the functional groups of the active components present in the extract based on the peak values in the region of IR radiation.When the extract was analysed by FTIR, the functional groups of the components were separated based on their peak ratio.
Electron ionisation mass spectroscopy (EIMS) of the isolated compound CS1 showed a mass ion peak at m/z 355 [M + H] (Supplementary figure 3), from which a molecular formula of C 23 H 46 O 2 was assigned.Typically, molecule CS1, at a retention time (Rt) of 0.714 min (Table 3), produces a precursor ion at m/z 355 [M + H], and the fragmentation of this molecule (Figure 2) generated product ions at m/z 298.These were derived from the loss of the isobutyl side chain (-57 Da) after a possible 1,3 methyl rearrangement of the isopropyl derivative of methylheptyl pentadecanoate to a more stable butyl pentadecanoate derivative (Figure 2).Product ions at m/z 284, due to the neutral loss of methyl (−14 Da), and at m/z 266 (loss of propyl molecule) were also observed.Based on these data, CS1 was identified as 6-methylheptyl pentadecanoate.
We, therefore, propose the structure and IUPAC name for compound CS1 based on the information obtained as 6-methylheptyl pentadecanoate.
The fragmentation pattern of 6-methylheptyl pentadecanoate was based on the analysis of mass spectroscopy in Supplementary figure 3. The fragmentation pattern is indicated below (Figure 2).
Proton on carbon C8 confirms the ester nature of CS2.  3).
The 2D NMR of the proposed compound, CS2, demonstrates varying coupling dimensions.A typical example is the two protons on carbon 10 close to the carbonyl moiety coupled with the adjacent proton on carbon 9 (Supplementary figure 13) as confirmed by HSQC Supplementary figure 14.Similarly, protons on carbons 12 and 13 can also be confirmed by the HSQC spectrum, which also coupled with each other to give triplets at δ 2.29 and 1.65.The methylene protons  The mass spectrometric analysis of CS2, showed a mass ion peak at m/z 489 (M + H), from which a molecular formula of C 33 H 60 O 2 was assigned (Supplementary figure 4).A retention time of (Rt) of 0.464 min (Table 3) produced a precursor ion at m/z 489 [M + H], and fragmentation of this molecule (Figure 4) generated product ions at m/z 414, derived from the loss of the isopentyl side chain (-75 Da), and m/z 359 due to the loss of propanol (-58 Da), m/z 300 also due to loss of the second propanol (-59 Da).Based on these data, in addition to the NMR and FTIR data, molecule CS2 was identified as 6-methylheptyl-15-(1,2,3,4,4a,8a-hexahydronaphthalen-1-yl)pentadecanoate.The fragmentation pattern of CS2 was based on the analysis of mass spectroscopy in Supplementary figure 4. The fragmentation pattern is indicated in Figure 4.

Biological activity prediction
A biological activity spectrum for a substance is a list of biological activity types for which the probability to be revealed (Pa) and the probability not to be revealed (Pi) are calculated.Pa and Pi values are independent, and their values vary from 0 to 1. Biological activity spectra were predicted for the two isolated structures of 6-methylheptyl pentadecanoate and 6-methylheptyl-15-(1,2,3,4,4a,8a-hexahydronaphthalen-1-yl)pentadecanoate via PASSonline 2005 version. 56Generally, in predicting the desired biological activity, Pa > Pi is considered feasible as there is a high chance of the compound revealing that activity.If Pa > 0.7, the compound is likely to reveal its activity in experiments, but in this case, the chance of being the analogue of the known pharmaceutical agent is high.If 0.5 < Pa < 0.7, the compound is likely to reveal this activity in experiments, but this is less, and the compound is not so similar to the known pharmaceutical agent.If Pa < 0.5, the compound is unlikely to reveal this activity in experiments, but if the presence of this activity is confirmed in the experiment, the compound might be a new chemical entity.
The biological activities predicted for each of the compounds herein include antieczematic, antiphobic, and antipruritic for CS1, as shown in Table 4, and antieczematic and antiulcerative for CS2, as shown in Table 4. Findings from the biological activity prediction show that both compounds had diverse activities towards different biological processes.However, the selected compounds were predicted to have a common antiviral property, particularly against rhinovirus.
In this study, special attention was given to certain reported activities of S. mombin that actively have antiviral properties. 57,58Hence, the selection of a suitable biological activity related to the antiviral activity of its isolated compounds, CS1 and CS2, was feasible.Additionally, the desired novelty of a chemical compound is important as well.The predicted Pa values for CS1 (0.655) and CS2 (0.643) both correlated to antiviral activity (rhinovirus), which falls within the 0.5 < Pa < 0.7 threshold correlating to a novel compound that has no known similarity to a known pharmaceutical agent.Subsequently, various antiviral macromolecules were selected to test the efficiency of CS1 and CS2 via in silico molecular docking.

Molecular docking
Molecular docking of selected rhinovirus targets used was human rhinovirus HRV (5FX6), human N-myristoyltransferase-HsNMT1 (5MU6), human N-myristoyltransferases-HsNMT2 (4C2X), and rhinovirus 3C protease (1CQQ).The compounds showed good binding towards the selected targets, as evidenced by obtaining an overall binding affinity in the range of −4.6 to −8.2 kcal/mol across all targets, as shown in Table 5.However, CS1 and CS2 proved to have the best binding affinity when docked to HsNMT1 (5MU6), suggesting they may have potential activity towards the HsMNT1 micromolecule, which is a potential drug target in developing therapeutics against the common cold.

Analysis of molecular dynamic simulation
Molecular dynamic simulations were conducted to assess the conformation dynamics as well as the spatial distribution of atoms in the backbone structure of HsMNT1 upon binding of the compounds.MD simulations were also employed to further validate findings from molecular docking by showing the most stable conformations of the complexed structures across time.Post-MD analysis protocols, including root-mean-square deviation (RMSD), and root-mean-square fluctuation (RMSF), radius of gyration (RoG) and solvent-accessible surface area (SASA), were employed to provide insights on the structural impact of the phytochemical compounds on HsNMT1.An error assessment was also established in analysing all MD trajectories to consider technical and biological variability.Eliminating these systematic errors lowers experimental variability and makes it possible to determine the underlying dynamics of protein motion in cellular signalling with greater accuracy.

Structural stability of HsNMT1
A 150-ns MD simulation trajectory was established to analyse the conformational dynamics of the c-α atoms in the backbone structure of HsNMT1 in all the simulated systems.The root means square deviation gives an estimation of the protein convergence and stability of the simulated system.Furthermore, the RMSD value estimates the average variation in atomic displacement over a given period of time compared to a reference time. 59The acceptable threshold for an average change in RMSD of a protein-ligand complex is between 1 and 3 Å.If the RMSD average is more significant than this threshold, it implies there is an extensive conformational alteration in the structure of the protein.
Findings show that systems converged early during the simulation and maintained steady atomic motions till the 150-ns simulation run, as shown in Figure 5A.The mean RMSD estimated for all the simulated systems were 1.88, 2.15, 1.54 and 1.83 Å for the unbound HsNMT1, 6-methylheptyl pentadecanoate, 6-methylheptyl-15-( within the acceptable range of 1.5−2.5 Å during the simulation.Also, good stability highlights the reliability of the simulated systems for further conformational analysis.

Structural flexibility of HsNMT1
The

Radius of gyration
The spatial arrangement of atoms in a protein-ligand complex system around its axis is known as the radius of gyration (RoG). 56,60stimating RoG is one of the most crucial indicators for predicting a macromolecule's structural activity, and it provides insights into variations in the compactness of the protein complex.Therefore, the stability of the unbound HsNMT1, CS1, CS2 and IMP-

Solvent-accessible surface area
Solvent-accessible surface area (SASA) impacts the structure and activity of biological macromolecules.SASA analysis provides important insights into residual exposure to surrounding solvent molecules during the simulation.Furthermore, due to the location of active site residues at the surface of the protein, greater insight into residue accessibility to solvent would be important in understanding the solvent-like behaviour (hydrophilic or hydrophobic) of a molecule as well as the proteinligand complex. 61,62SASA analysis can also be used to describe protein folding and unfolding. 61As such, the SASA for the simulated systems was computed, as shown in Figure 5D.The averages estimated for the simulated systems were 18570.40,17877.74,17707.02and 18000.32Å 2 for the Apo (iron-free), CS1, CS2 and Imp-1088, respectively.The SASA values of the complexed systems were slightly lower than those of the unbound HsNMT1 system, indicating a lower surface area exposed to solvent.The binding of the inhibitors induces rigidity in the amino acids in the structure of HsNMT1 upon binding.Findings further highlight the similarity in the structural impact of the compounds and the reference inhibitor of HsNMT1.

Binding free energy
The mechanics/generalised-born surface area (MM/GBSA) method was employed to estimate the binding free energetics of the complexed systems of CS1 and CS2, including the reference IMP-1088 compound.It is well recognised that the MM/GBSA method for predicting binding energy is more accurate than most molecular docking scoring functions and computationally less complex than alchemical free energy techniques. 63-66The computed binding free energies for the complexed systems of HsNMT1 were estimated to be −35.20 kcal/mol for CS1, −44.55 kcal/mol for CS2 and −47.06 kcal/mol for IMP-1088.Findings show that CS2 had the strongest binding free energy among the two compounds; however, both compounds demonstrated overall stronger energies than the reference compound used in the study.The results indicate that these compounds can be considered potential inhibitors of HsNMT1.Table 6 indicates the energy terms that contribute to the binding free energy, the most favourable components being the Δ  , Δ  and Δ gas , whilst ΔG sol is unfavourable.The MM/GBSA method is a well-known technique that demonstrates computational effectiveness using implicit solvent and also offers a transparent environment for determining the physical causes of observed effects in protein-ligand interactions. 28,66Taken together, the energies presented by these compounds suggest the spontaneity, permeation and a measure of the reaction kinetics that characterise their complexing with the target protein.

Binding interactions
The types of interactions a molecule has in a target protein's binding pocket emphasise how therapeutically effective it is for the protein. 67he binding interactions of CS1, CS2 and the reference IMP-1088 compound bound to HsNMT1 were assessed.The CS1 and CS2 as potential inhibitors were observed to engage in a variety of interactions involving conventional and carbon-hydrogen bonds, van der Waals and pi-alkyl, alkyl interaction as depicted in Figure 5.The variation in interaction types between the potential inhibitors and the binding site residues was attributed to the different molecular features.Assessing the interaction profile of the reference (IMP-1088) compounds showed similar interaction types, as observed in Figure 6.The interactions observed herein include conventional and carbon-hydrogen bonds, van der Waals and pi-alkyl, alkyl interaction, pi-pi stacked, and pi-pi T-shaped.Findings revealed similar interactions with binding site residue between the compounds and the reference compound, suggesting CS1  The outcome of this investigation highlighted several possible biological activities.However, the selection of suitable biological activity was considered based on a higher Pa value than a Pi value.Special attention was given to the reported biological activity associated with the S. mombin leaf extracts known to have antiviral activity.Thus, the suitable biological activity predicted for the two isolated novel esters was antiviral activity, particularly towards rhinovirus.Amongst the selected targets, CS1 and CS2 showed a higher binding potential for HsMNT1, an essential enzyme in treating the common cold.The MD simulation employed to test the effect of the compounds against HsNMT1 enzymes revealed that the compounds exhibited good stability, flexibility, structural rigidity, and reduced surface area exposed to solvents.These structural effects of the compounds towards HsNMT1 were similar to the structural effects of the reference inhibitor, suggesting the potential inhibitory effects of the compounds towards HsNMT1.
In silico molecular recognition protocols were employed to assess the pharmacological effects of the compounds CS1 and CS2 on the S. mombin leaf.The predicted biological activity for the two isolated novel esters was anti-rhinovirus activity.
Molecular docking analysis indicated that CS1 and CS2 showed a higher binding potential towards HsMNT1.The MM/GBSA method revealed stronger binding free energy in CS1 and CS2 than the reference compound.Assessment of binding interactions also shows similarity in interactions CS1, CS2 and the reference IMP-1088 inhibitor, indicating the potential to elicit similar therapeutic effects against HsNMT.

Conclusion
The current study of the phytochemical analysis of DCM leaf extracts of S. mombin led to the identification of two esters that had previously not been reported in the plant.These compounds, 6-methylheptyl pentadecanoate and 6-methylheptyl-15-(1,2,3,4,4a,8a-hexahydronaphthalen-1-yl)pentadecanoate, possess anti-rhino virus (HsNMT1) properties as indicated through an in silico molecular targeting prediction.
Further in vitro validation is required to optimise a potential drug candidate.
: R f = Distance traveled by the solute ________________________ Distance traveled by the solvent Equation 1
root means square fluctuations were assessed to determine the relative flexibility of the c-α atoms in the backbone structure of HsNMT1 upon binding of the inhibitors.As such, the RMSF values of the unbound HsNMT1, CS1, CS2 and IMP-1088 complexed with HsNMT1 were estimated to observe the change in protein structural flexibility during the simulation run.As shown in Figure 5B, all the selected compounds, including the reference IMP-1088 compound, show a peak area of the protein at Glu130, Leu175, Lys240, Ser315 and Thr395 residual positions that fluctuate the most during the simulation.It was observed that the amino acid residues where the reference IMP-1088 bound have similar structural behaviour as that of the phytochemically bound systems of HsMNT1.The mean RMSF values estimated were 0.98 ± 0.03, 1.01 ± 0.04, 0.87 ± 0.02 and 0.95 ± 0.03 Å for unbound HsNMT1, CS1, CS2 and IMP-1088, respectively, showing that the values are very close to each other.However, in comparing the relative flexibilities of the simulated systems, the complexed HsNMT1 systems show lower fluctuations in contrast to the native unbound system of HsNMT1, indicating that the bound inhibitors enact rigidity on the protein structure.
1088 complexes was estimated by measuring RoG over the 150-ns simulation as shown in Figure 5C.The respective RoG averages computed were 21.85, 21.75, 21.77 and 21.78 Å for the Apo (HsNMT1), CS1, CS2 and IMP-1088 systems.The similarity in mean values of the native unbound state (apo) of HsNMT1 and the bound complexes indicates that the selected compounds do not induce major conformational changes to the active site upon binding.

Figure 5 :
Figure 5: Comparative C-α RMSD, RMSF, RoG and SASA plots showing conformational alterations upon binding of the compounds and reference compound to HsNMT1 over the 150-ns MD simulation time.(A) The RMSD plots, which indicate the compounds induced relative stability on the HsNMT1 enzyme upon binding.(B) The RMSF plots indicating peak regions of residual fluctuations.(C) Relative compactness of all simulated systems of complex structures and the unbound (Apo) system.(D) The surface area exposed to solvent between the simulated systems.

Figure 6 :
Figure 6: 2D molecular interactions of inhibitors (A) CS1, (B) CS2 and (C) IMP-1088 within the binding site of HsNMT1 showing similar interactions with the binding site residues, suggesting the compounds have the potential to elicit similar therapeutic effects as reference IMP-1088.

table 2 :
Physical appearance, weight of dry crude extract per solvent of extraction and percentage (%) yield from 60 g SM dry powder

table 4 :
-yl)pentadecanoate and IMP-1088 complex systems, respectively.As observed from the findings, all systems attained good stability due to the maintenance of mean RMSD values Predicted biological activity via PASSonline 0.655 0.004 Antiviral (rhinovirus) 0.643 0.013 Antiviral (rhinovirus) Key: Pa, probability to be revealed; Pi, probability not revealed