Skip to main content

ORIGINAL RESEARCH article

Front. Cell. Infect. Microbiol., 24 May 2023
Sec. Virus and Host
This article is part of the Research Topic Evolution, Pathogenesis, Host Interactions and Therapeutic Strategies Against Monkeypox Virus View all 10 articles

Mechanistic inhibition of Monkeypox and Marburg virus infection by O-rhamnosides and Kaempferol-o-rhamnosides derivatives: a new-fangled computational approach

  • 1Biophysics and Biomedicine Research Lab, Department of Electrical & Electronic Engineering, Islamic University, Kushtia, Bangladesh
  • 2Laboratory of Computational Research for Drug Design and Material Science, Department of Chemistry, European University of Bangladesh, Dhaka, Bangladesh
  • 3Department of Microbiology, West Bengal State University, West Bengal, Kolkata, India
  • 4Department of Health Sciences, Novel Global Community Educational Foundation, Habersham, NSW, Australia
  • 5Department of Poultry Science, Faculty of Animal Science & Veterinary Medicine, Sher-e-Bangla Agricultural University, Dhaka, Bangladesh
  • 6Department of Microbiology, Adamas University, West Bengal, Kolkata, India
  • 7School of Electronic Science and Engineering, Southeast University, Nanjing, China
  • 8Department of Life Sciences, Presidency University, Kolkata, West Bengal, India
  • 9Department of Pharmacy, Daffodil International University, Sukrabad, Dhaka, Bangladesh
  • 10Department of Science and Engineering, Novel Global Community Educational Foundation, Habersham, NSW, Australia
  • 11Department of Neuroscience, AFNP Med, Wien, Austria
  • 12Pharmaceutical Biotechnology Laboratory, Department of Pharmaceutical Chemistry, College of Pharmacy, King Saud University, Riyadh, Saudi Arabia
  • 13Microbiology Division, Department of Botany, Gauhati University, Assam, India
  • 14Department of Occupational Medicine, Tainan Municipal Hospital (managed by Show Chwan Medical Care Corporation), Tainan, Taiwan
  • 15Department of Public Health, College of Medicine, National Cheng Kung University, Tainan, Taiwan
  • 16Department of Rasa Shastra and Bhaishajya Kalpana, Faculty of Ayurveda, Institute of Medical Sciences, Banaras Hindu University, Varanasi, Uttar Pradesh, India

The increasing incidence of Monkeypox virus (Mpox) and Marburg virus (MARV) infections worldwide presents a significant challenge to global health, as limited treatment options are currently available. This study investigates the potential of several O-rhamnosides and Kaempferol-O-rhamnosides as Mpox and MARV inhibitors using molecular modeling methods, including ADMET, molecular docking, and molecular dynamics/MD simulation. The effectiveness of these compounds against the viruses was assessed using the Prediction of Activity Spectra for Substances (PASS) prediction. The study’s primary focus is molecular docking prediction, which demonstrated that ligands (L07, L08, and L09) bind to Mpox (PDB ID: 4QWO) and MARV (PDB ID: 4OR8) with binding affinities ranging from -8.00 kcal/mol to -9.5 kcal/mol. HOMO-LUMO based quantum calculations were employed to determine the HOMO-LUMO gap of frontier molecular orbitals (FMOs) and to estimate chemical potential, electronegativity, hardness, and softness. Drug similarity and ADMET prediction assessments of pharmacokinetic properties revealed that the compounds were likely non-carcinogenic, non-hepatotoxic, and rapidly soluble. Molecular dynamic (MD) modeling was used to identify the most favorable docked complexes involving bioactive chemicals. MD simulations indicate that varying types of kaempferol-O-rhamnoside are necessary for successful docking validation and maintaining the stability of the docked complex. These findings could facilitate the discovery of novel therapeutic agents for treating illnesses caused by the Mpox and MARV viruses.

Introduction

In the middle of the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection that emerged as a novel virus in 2019, and has since caused the ongoing coronavirus disease 2019 (COVID-19). Two additional zoonotic viruses, Monkeypox virus (Mpox) virus and Marburg virus (MARV), have recently re-emerged in the current years (Islam et al., 2023; Zaeck et al., 2023). With over 60,799 confirmed cases of Mpox reported from 99 countries and locations, along with 20 deaths as of September 16, 2022, the re-emergence of Mpox, a member of the family Poxviridae and genus Orthopoxvirus that also contains smallpox virus, causing monkeypox (Mpox) in humans, has posed a public health emergency of international concern (Mohapatra et al., 2022; Khattak et al., 2023). MARV, a member of Filoviridae family that also contains deadly Ebola virus, which causes Marburg virus disease (MVD), and has a fatality rate of up to 90%, re-emerged in Ghana, an African country, in July 2022, as well as it was also reported from Guinea in 2021; there have been a total of 15 MVD outbreaks, the vast majority of which have occurred in Africa (Abir et al., 2022; Aborode et al., 2022; Zhao et al., 2022). MARV was initially documented to infect people in 1967, while Mpox was first identified in 1958, with the first zoonotic human case confirmed in 1970. These two viruses have been rare and forgotten human diseases for almost 50 years. Among these, Mpox has recently become a global health emergency (Farahat and Memish, 2022; Anwar et al., 2023).

The viruses are obligate intracellular parasites that can act as living, non-loving objects (Claverie, 2006). They are very tiny in size and have a simple structure. They have the lack any metabolic activity and consist of either DNA or RNA as a nucleic acid, protein, and membrane lipids which surround them. Generally, the proteins provide the structural shape of the virus, which include capsid proteins, matrix proteins, and membrane glycoproteins (Whittaker and Helenius, 1998; Claverie, 2006). Due to obligate parasites in nature, they can act as non-living objects outside of host cells, and when they reach the physiological system or any living host cell, they are able to cross cell membranes even several other barriers in order to transport their genetics into the cytosol or nucleus of the host cell (shown in Figure 1), which ultimately results in pathogenicity or infectious disease to the host cells (Marie, 1957). The viral infection spreads quickly, and host-cell propagation occurs in three primary stages: viral particle formation in infected host cells; virus escape into the extracellular space; and finally, virus entrance into a newly exposed cell membrane (Boulant et al., 2015). Diseases caused by viruses often start in the periphery of the body, and may migrate to the internal organs, liver, kidney, heart, gastrointestinal system, neurological system of mammals, where they can directly impact the vital organs, central nervous system (CNS) and the peripheral nervous system (PNS) (Becker, 2021).

FIGURE 1
www.frontiersin.org

Figure 1 General pathogenesis and life cycle of virus (Created with BioRender.com).

According to the recent outbreaks and findings, viruses are still the cause of a vast variety of infectious illnesses in humans and animals, as well as a substantial proportion of newly developing and reemerging infections that are considered to be of significant health concern (Fauci, 2001). Some emerging/re-emerging notable pathogenic viral infections include SARS CoV-2 (Islam et al., 2022; Shanmugaraj et al., 2022), influenza viruses including bird flu virus (H5N1), Crimean-Congo hemorrhagic fever virus (CCHFV), Zika (ZIKAV), Nipah (NiV), Hendra virus (HeV), Ebola (Coltart et al., 2017), Monkeypox (Beer and Rao, 2019), Langya Henipavirus (Akash et al., 2022) and Marburg virus (Akash et al., 2023). Due to the rarity of Mpox, it has infected a huge number of people during the current COVID-19 pandemic still has a lot of unknowns (León-Figueroa et al., 2022; Taseen et al., 2022). Infections can spread from animals to humans, making this a zoonotic virus. Rapid human-to-human transmission has been primarily noticed with the recent re-emergence of Mpox in multiple nations, causing a worldwide health emergency. This virus is the most commonly found in Western and Central Africa (Chakraborty C. et al., 2022; Velavan and Meyer, 2022). The clinically significant illness produced by Mpox is comparable to those seen in smallpox but is not as devastating (Harapan et al., 2020; Khattak et al., 2022). This pathogenic virus was first reported in 1970 in an infant who had previously been affected by smallpox (Ladnyj et al., 1972). Despite the fact that this disease was discovered more than 60 years ago, it has lately been reported to infect SARS CoV-2 patients who have been hospitalized since the emergence of the COVID-19 pandemic (Shanmugaraj et al., 2022). Since January 2022, the World Health Organization has received reports of more than 47000 patients diagnosed with Mpox (Sharma et al., 2022). Researchers and healthcare providers are getting more concerned about the possibility of a newer pandemic of Mpox (Adalja and Inglesby, 2022; Ahmed et al., 2022). Secondly, MARV as another zoonotic pathogen has recently attracted attention of the scientists owing to its recent year’s outbreaks in Ghana and Guinea (Nyakarahuka et al., 2019). In Uganda, there were reported for three separate MVD epidemics in 2007, 2012, and 2014, and thereafter from other African countries and beyond Africa, 6-8 MVD epidemics were also identified (Adjemian et al., 2011; Knust et al., 2015; Nyakarahuka et al., 2017; Zhao et al., 2021). The MARV produces severe or even life-threatening effects owing to its very high case fatality rates (Balter, 2000; Zhao et al., 2021).

There are very few approved drugs, licensed vaccine, or therapeutic without vast ranging of potential activity yet available against both the Mpox and MVD, though efforts have been made, as well currently being advanced for finding out effective prophylaxis and treatment regimens for both of the diseases (Dulin et al., 2021; Chakraborty S. et al., 2022; Hickman et al., 2022; Sah et al., 2022). However, European medical agency has approved tecovirimat drug against Mpox (Ali et al., 2023). Many researchers are trying to develop potential medication from the natural sources against mentioned disease. Currently, a number of different drugs have been shown to be effective against the Mpox when tested in vitro, such as Adamantane derivatives, Resveratol, Tecovirimat, Mitoxantrone, and Ribavirin. Nevertheless, their usefulness and safety in the context of Mpox have not been shown, and further basic and clinical research is needed (Ali et al., 2023; Khan et al., 2023; Khani et al., 2023).

In addition to novel drugs, chemical ligands, broad-spectrum antivirals, and broadly neutralizing antibodies, it is worth noting that the prophylactic and therapeutic aspects of medicinal plants and herbs, their extracts and metabolites, phytochemicals, immunity-enhancing foods, and dietary nutrition have shown potentials for countering several important infectious pathogens and deadly viruses, including Mpox and MARV (Khan et al., 2023; Nguyen, 2023; Pourhajibagher and Bahador, 2023). O-rhamnosides, including Kaempferol-O-rhamnosides, are a class of polyphenolic compounds with a variety of beneficial biological effects, such as antimicrobial, antioxidant, and anti-inflammatory properties (León-Figueroa et al., 2022; Velavan and Meyer, 2022). Their propensity to bind with p-hydroxybenzoic acid in DNA may also explain why they suppress polymerase chain reaction. Virus replication is stifled as a result of this interaction, which blocks the transcription of viral DNA into RNA (Zhao et al., 2021). Therefore, Kaempferol-o-rhamnoside and its eight derivatives were chosen as the primary molecules in this investigation, and their antiviral capability was investigated against Mpox and MARV infections for their utility in treating patients with Mpox and MVD. Although these two pathogenic viruses have the feasible potential to cause another pandemic, it is unfortunate that no effective or potentially effective medication has been developed to date (Kortepeter et al., 2020). So, this investigation has been performed to identify a potential medication which may inhibit Mpox and Mpox with lower adverse effects by using Computer-Aided Drug Design (CADD) techniques and quantum chemistry of chemical descriptors. In comparison to traditional procedures, these strategies have been shown to be productive across the whole drug development process, resulting in a reduction in the amount of money and time required to produce a medicine (Rahman et al., 2012).

Materials and methods

Optimization and ligand preparation

O-rhamnoside is a type of carbohydrate and contains the active methyl group where one, two and three groups are attached to form the L02, L03, and L04 compounds. Kaempferol, on the other hand, is a phytochemical derived from the Talipariti elatum tree and was employed as the mother compound in this investigation. Next, the Kaempferol is attached with L01, L02, L03, and L04 to form L06, L07, L08 and L09 as derivatives of Kaempferol-o-rhamnoside.

DFT functional methods were employed towards complete molecular optimization using DMol3 code from Material Studio 08 package (Kumer and Khan, 2021). To ensure accuracy, the DMol3 code was put up using the functional of B3LYP and Gaussian double zeta plus polarization function basis set (DNP). Following geometric optimization, the HOMO and LUMO molecular border orbital diagrams were split, and calculated the electron affinity, electron negativity, energy gap, chemical potential, hardness, softness, and electrophilicity index using the following equations (a-h), respectively. Finally, the optimized molecule was stored as a PDB file for future in-silico investigation, such as molecular docking, molecular dynamics, and ADMET analysis.

Egap =( ELUMOEHOMO).(a);I=EHOMO..(b)A=ELUMO…………(c);(x)=I+A2………….(d);(ω)=μ22η………….(e)(μ)=I+A2…………(f);(η)=I+A2……….(g);(s)=1η…………(h)

termination of the data of ADMET and lipinski rule

It is critical to forecast ADMET traits during the discovery or development of a new pharmaceutical in order to prevent treatment failure throughout clinical investigations (Cheng et al., 2013). As a consequence, in-silico pharmacokinetic prediction plays a significant role. The pkCSM was utilized to create this prediction so that there was less risk of the new drug molecules failing during the trial phases and greater possibility of them making it to the final step as prospective drug candidates (Adalja and Inglesby, 2022). This prediction is based on key factors, such as molecule uptake in the human intestine, percolation capabilities of the blood-brain barrier, central nervous system, metabolism capabilities, total clearance, bioavailability, and toxicity levels.

PASS prediction

The PASS prediction has been made to determine the capability of the ant-viral drug. The online web program PASS (http://www.pharmaexpert.ru/passonline/) was used to predict the anti-bacterial, anti-fungal, and anti-parasitic spectrum (Lagunin et al., 2000). The configurations of Kaempferol-o-rhamnoside derivatives were illustrated first, and then they were transformed into smile forms by the addition of the free online programs provided by SwissADME (http://www.swissadme.ch) (Daina et al., 2017). These programs are well-known for their ability to determine PASS spectrum by making use of the PASS web tool. PASS findings are represented by the probabilities Pa (probability for an active molecule) and Pi (probability for an inactive molecule). Pa and Pi grades may vary from 0.00 to 1.00, and Pa and Pi must be less than 1, since potentialities can be anticipated in whatever way the researcher chooses. The drug candidate should be potential and bioactive if the score of Pa > Pi (Poroikov and Filimonov, 2005).

Preparation of protein and method for molecular docking

Molecular docking is one of the most effective ways of virtual screening, particularly when there is a computer-based drug design and three-dimensional structure of the receptor and ligand accessible. Firstly, two targeted receptor proteins, such as Mpox (PDB ID: 4QWO), and MARV (PDB ID: 4OR8) were download from Protein Data Bank (PDB) (https://www.rcsb.org/) which were isolated by X-ray diffraction method with high stable configuration. After that, PyMol version 2021 was applied for removing all the excess heteroatom and minimized the energy by Swiss PDB viewer (Joseph and Namboothiri, 2015; Mooers and Brown, 2021). Finally, docking analysis was completed with the help of AutoDock tool and during this progres (Dallakyan and Olson, 2015), the grid box parameter was used as center X =12.4697, Y =15.9818, Z=16.0634, and Dimension (Å) X = 35.144, Y= 37.645Å, and Z= 36.966 Å for Mpox (PDB ID 4QWO). In case of MARV, center is in X = 14.267 Å, Y = 4.929, Z= 14.557, and Dimension (Å) X = 53.382 Å, Y= 48.667 Å, and Z= 66.040 Å. Three-dimensional structures of Mpox and MARV proteins are displayed in Figure 2.

FIGURE 2
www.frontiersin.org

Figure 2 Three-dimensional protein structure of Monkeypox and Marburg viruses (Zhang et al., 2014; Minasov et al., 2022).

Molecular dynamics simulation

On a high-configuration computer, molecular dynamics simulation was used in live view according to the capabilities afforded by NAMD (Kumer et al., 2022b). The docking results for the most powerful drugs were validated by molecular dynamics simulations up to 100 nanoseconds for the hollow tube (drug-protein). In addition, NAMD is particularly well suited to the increasingly popular Beowulf-class PC clusters, which are quite similar to the workstation clusters for which it was originally designed. MD simulation determines the stability of molecules and verifies the docking fitting. In this study, the MD simulation was performed at 100ns and applying AMBER14 force field (Nath et al., 2021). Before that, the whole system was in the presence of a water solvent equilibrated with 0.9% NaCl at 298 K temperature in presence of liquid system with NVT ensemble having the simulation box size as boundary 1.5 and X=5, Y=5 and Z=5 which is standard for these proteins for solvation, as well as, its vector coordinate of box is at -25.06, 18.79 and -7.06, respectively. VMD evaluated the performance of this result using RMSD and RMSF. Next, the Ramachandran Plot and B-factor were evaluated to explain the validation and stability of docked complex.

Results and analysis

Chemistry of O-rhamnosides and Kaempferol-O-rhamnosides derivatives

In case of O-rhamnoside (L01), it is revealed that the L02, L03 and L04 are formed by attaching one, two and three O-rhamnoside rings with the original compound of O-rhamnoside, as shown in Figure 3. Moreover, the compound L01, L02, and L03 is attached with Kaempferol to form L07, L08 and L09 compounds. Finally, the computational and in silico analysis as well as SAR and comparative studies are performed on them.

FIGURE 3
www.frontiersin.org

Figure 3 Chemical structure of derivatives of O-rhamnosides and Kaempferol-O-rhamnosides.

Optimized structure of derivative of O-rhamnosides and Kaempferol-o-rhamnosides

In computational chemistry, the use of quantum mechanical methods for the calculation of thermodynamic, molecular orbital, and molecule electrostatic characteristics is prevalent. The program known as Gaussian 09 has been run in order to improve the geometry and further enhance each generated derivative. The density functional theory (DFT), functional of B3LYP and Gaussian double zeta plus polarization function basis set (DNP), was utilized in order to improve and predict the molecular orbital and thermal properties of the molecules. All compounds had their dipole moment, enthalpy, free energy, and electronic energy determined. The optimized molecular structures are shown in Figure 4.

FIGURE 4
www.frontiersin.org

Figure 4 Optimized structure of derivative of O-rhamnosides and Kaempferol-o-rhamnosides.

PASS prediction

The probable biological spectrum for Kaempferol-o-rhamnoside derivatives has been predicted by applying the web server PASS. The PASS data is summarized as Pa and Pi, which are shown in Table 1. According to the presupposition in Table 1, Kaempferol-o-rhamnoside derivatives 02–09 demonstrated 0.379<Pa< 0.712 for antiviral (Influenza), 0.423<Pa<0.655 for antibacterial, 0.632< Pa< 0.755 for antifungal, and 0.214<Pa<0.391 for anti-parasitic, while Ligand 07 demonstrated Pa 0.634<pa. Consequently, the PASS assessment showed that 0.379< Pa< 0.712 for antiviral (Influenza), which indicated that the Kaempferol-O-rhamnoside derivatives have a greater potential as antiviral agent compared to other antibacterial and anti-parasitic characteristics. Although the antifungal Pa score is higher, it has been largely ignored in favor of the antiviral Pa score since researchers are more interested in discovering new antiviral drugs.

TABLE 1
www.frontiersin.org

Table 1 Data of PASS prediction.

Lipinski’s rule, pharmacokinetics and drug likeness

Due to a variety of factors, the most medication candidates never become commercially available or cannot pass the clinical or preclinical stages. It is the utmost significance to create trustworthy computational approaches for the estimation of the drug-likeness of new drug candidates that increase the percentage of successful drug discovery and development attempts during the trial phase. The drug-likeness prediction model was developed by Christopher A. Lipinski in 1997, which included molecular descriptors connected to the numbers of factors that can differentiate between possible drugs and non-drugs. These factors include topological polar surface area, number of rotatable bonds, hydrogen bond acceptor, hydrogen bond donor, and molecular weight. The molecular weight of the reported Kaempferol-o-rhamnoside derivatives (L01-L09) was at 164.16-754.69, the number of rotatable bonds was 0.0-10, the hydrogen bond acceptor was 05-19, the hydrogen bond donor was 04-13, and the topological polar surface area was 90.15Ų -318.37Ų. So, only Ligand of L01, L02, L06, and L07 satisfied the overall Lipinski rule (as shown in Table 2). However, the other molecules do not satisfy because they have a larger topological polar surface area, have a hydrogen bond acceptor, or are a hydrogen bond donor. Finally, the bioavailability score for ligands L01, L02, L04, and L05 is greater than 0.55, and the remaining molecules’ bioavailability scores are reported at 0.17. So, the topological polar surface area or hydrogen bond acceptor and hydrogen bond donor have been abandoned.

TABLE 2
www.frontiersin.org

Table 2 Data of lipinski rule, pharmacokinetics and drug likeness.

Quantum calculation and HOMO-LUMO and frontier molecular orbital (FMO)

Frontier molecular orbitals (FMO) are the most important orbitals in a molecule. They are used by scientists to study chemical reactivity and kinetic stability. People often talk about the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) when talking about these molecular orbitals on the edge of the atom (LUMO). When you say “electronic absorption,” you mean a change from the ground state to the first excited state. This is usually thought of as one electron moving from the HOMO state to the LUMO state (Perepichka and Bryce, 2005; Kumer et al., 2021a). A molecule’s high level of chemical stability and low level of chemical reactivity are directly linked to how long its energy difference lasts (Kumer et al., 2022a). Because it makes it easier for electrons to move from one electron to another, a low energy gap is linked to low chemical stability and high chemical reactivity. Because of this, it takes a lot more energy for an electron to move from its ground state, called HOMO, to its excited state, called LUMO. Table 3 highlights the HOMO and LUMO energies, the HOMO-LUMO gap (Δ), and the hardness (η) and softness (S) indexes of all derivatives of Kaempferol-o-rhamnoside, while Figure 5 highlights, for better understanding, the frontier orbital diagrams of both HOMO and LUMO, as shown using a variety of color schemes to facilitate better comprehension. In the context of HOMO, the shade of deep radish signifies the positive nodes of orbitals, while the shade of yellow corresponds to the negative node of orbitals.

TABLE 3
www.frontiersin.org

Table 3 Data of chemical descriptors.

FIGURE 5
www.frontiersin.org

Figure 5 Frontier Molecular orbital for derivative of O-rhamnosides and Kaempferol-o-rhamnosides.

As indicated in Table 3 and Figure 5, the energy gap value of Kaempferol-O-rhamnoside derivative 01-04 is consistently larger (9.524, 9.311, 9.095 and 9.303 eV) compared to other considered molecules. The range of chemical hardness lies between 4.762 and 3.477, while the softness values of the Kaempferol-O-rhamnoside derivatives extend from 0.2100 to 0.2876. It's noticeable that hardness and softness share a reciprocal relationship; a decrease in hardness consequently corresponds to an increase in softness. In the case of ligands, the least hardness is observed alongside the highest softness.

From the Figure 5, it is observed that the LUMO is scattered throughout in whole molecules for L01 and L02, but it is different for the other ligands (L03-L09). The LUMO is stayed in the O-rhamnoside for L01 to L09 but the HOMO is found in the kaempferol ring.

Molecular electrostatic potential charge distribution mapping

The Molecular Electrostatic Potential, also known as MEP, is utilized as a global reactivity map. This map illustrates the amount of an organic molecule that is most equipped for the electrophilic and nucleophilic onslaught of charged point-like reagents (Rahman et al., 2022). It is helpful in understanding the biological recognition system and the coupling of hydrogen bonds. This MEP counter map offers a straightforward method for forecasting how various types of geometry could interact with one another. The significance of molecular electrostatic potential (MEP) arises from the fact that it concurrently demonstrates a molecular shape and size and also positive, negative, and neutral electrostatic potential provinces in terms of visual assessment. In addition to this, MEP is very effective in the investigation of macromolecules with physical and chemical features. On the basis of the optimized structure of the Kaempferol-O-rhamnoside derivatives (L01-L09), MEP was utilized to make a prediction about the active sites on the surface that would be susceptible to electrophilic and nucleophilic interactions. The different electrostatic potential has been shown by the red and blue hue. The red color indicates the maximum negative area, which is a significant site for electrophilic interactions; the blue color represents the largest positive region, which is an acceptable site for nucleophilic interactions; and the green color represents the neutral potential area. (See Figure 6).

FIGURE 6
www.frontiersin.org

Figure 6 MEP mapping for the ligands after DFT optimization.

Molecular docking

The forces between molecules that occur between the ligand and the receptor during the formation of a protein-ligand complex are called binding affinity. This method has grown in popularity as a result of the research and investigation of innovative drug design (Meng et al., 2011). A binding energy of -6.00 kcal/mole has been speculated as a potential drug candidate (Kumer et al., 2021b). In this experiment, the Auto-Dock Vina program was used to investigate a set of Kaempferol-o-rhamnoside derivatives by the in-silico method in order to emphasize their putative binding energy and interaction configurations with the targeted protein of Mpox (PDB: 4QWO) and MARV (PDB: 4OR8). In light of the findings that were acquired from the docking inspection, the three derivatives with the highest binding energies were found. The results of computed analysis showed that three derivatives L07, L08, and L09 had potential docking scores of -9.2 (> -6.0), -9.0(>-6.0) and -9.4(>-0.6) against protein Mpox (PDB: 4QWO). Also, these three derivatives L07, L08, and L09 had potential docking scores -8.1 (> -6.0), -8.4(>-0.6) and -9.0(>-0.6) against Crystal Structure of the MARV (PDB: 4OR8). Where the traditional drug Cidofovir had the docking score -6.0(=-6.0) and -5.2(<-6.0) against protein Mpox (PDB: 4QWO) and Crystal Structure of the MARV (PDB: 4OR8) respectively. It demonstrates that L07, L08, and L09 substances can firmly bind to the protein receptors’ pockets than Cidofovir (as shown in Table 4). Besides, it is seen that by increasing the side chain, the binding energy is equally increased. It has been said that compared with standard antiviral drugs like Cidofovir, Ligand 03–09 shows better effectiveness.

TABLE 4
www.frontiersin.org

Table 4 Binding affinity (kcal/mol) against Monkeypox and Marburg virus.

Protein-ligand interaction

The Kaempferol-o-rhamnoside derivatives generated a wide variety of different types of non-covalent interactions with the active sites of the microbial proteins. These observed bonds included pi-alkyl, pi-anion, pi-sigma, amide pi-stacked, and pi-donor hydrogen bonds. Based on these results, it is abundantly obvious that aromatic substituent may quickly enhance the binding ability of uridine esters as well as their potential to inhibit the growth of reported viral pathogens due to the high electron density of aromatic substituent. The influence of H-bonds had a significant impact on the selectivity of ligand binding to the receptor, drug design in chemical and biological processes, and molecular recognition and biological activity. Therefore, the H-bond surface of all derivative ligands that interact with both proteins is presented in Figure 7. The 2D structure of protein ligand interaction also illustrates how to determine how many bioactive sides are present.

FIGURE 7
www.frontiersin.org

Figure 7 Lead picture of ligand and standard.

Molecular dynamics simulation

Molecular dynamics (MD) simulations are effective tools for comprehending the physiological underpinnings that encompass the structure and function of bioactive molecules (Barroso da Silva et al., 2020). The initial conception of proteins as relatively inflexible arrangements has been superseded by a dynamic model according to which the internal movements of macromolecules and the conformational changes that follow from those movements play a significant key role in the development of new drug molecules. This MD simulation is a helpful tool for predicting the stability of drug-like molecules (Hickman et al., 2022). The outcome of MD simulation is described by RMSD and RMSF, which mean root-mean-square deviation and root-mean-square fluctuation. Stable molecules are biological molecules with a RMSD and RMSF score of less than 2.0 A° (Kawsar et al., 2022). So, based on maximum docking outcomes, Ligands of L07, and L09 have been conducted in MD simulation as well as standard Cidofovir (D1). According to the finding of MD simulation against Mpox (PDB ID 4QWO), the RMSD is about 0.9 Å –1.2 Å and similar outcomes were found in terms of RMSD in terms of backbone with amino acid residues. But, in RMSF, little impact has been seen in the addition of backbone with amino acid residues and there is fluctuation 0.7 Å – 0.8 Å (as shown in Figures 8A, B).

FIGURE 8
www.frontiersin.org

Figure 8 Analysis of (A) RMSD; (B) RMSF for L07, L09 and Cidofovir drugs in complex with Monkeypox virus (C) RMSD and (D) RMSF for L08, L09 and Cidofovir drugs in complex with Marburg virus.

Secondly, the MD simulation has been completed for ligands 08 and 09 towards the MARV (PDB 4OR8). In this case, the RMSD and RMSF have been obtained with equivalent outcomes, and the RMSD is about 0.9 Å – 1.2 Å both time vs. backbone and backbone with amino acid residues, while the RMSF score at 0.7 Å – 0.8 Å (as shown in Figures 8C, D). In summary, the RMSD and RMSF scores are accepted and less than 2.00 Å which means these drugs could be stable after reaching biological systems.

Ramachandran plot from molecular dynamic study

The Ramachandran plot is one of the most important parameters for checking chemical stability of the protein or enzyme structure satisfying both stereochemistry and steriochemical structure. In addition, it is also used for validation of docking and stability of docked complex of protein ligand. Also, this analysis illustrates the energy minimization of Botherombin and most favorable region. From the Figure 9, the Ramachandran plot is showed for the top three drugs (L07, L08 and L09), and compare with standard drugs (D1). There are 225 amino acids rasidues where 216 amino acid residues (96.0%) are out of favorable region, and 9 amino acids (4.00%) in unfavorable region for L07 which is almost similar trend for the L08, L09 and D1, means that there are no change after entering the ligand in protein. So that, it can be concluded that the molecular docking is valid and highly stable docked complexes were obtained after docking.

FIGURE 9
www.frontiersin.org

Figure 9 Ramachandran Plot from Molecular Dynamic Study for L07, L08, L09 and D1 for Monkeypox virus.

B-factor in case of normal mode analysis (NMA) from molecular dynamic study

The expression B-factor, occasionally entitled the Debye-Waller factor, temperature factor, or atomic displacement parameter is to show by curve. In gerenal, the most significance of it is to express the higher flexibility results in larger displacements or lower electron density, indicating the lower stability. So after molecular docking, the B-factor of docked protein complex can be calculated that helps to say its deformation or stability at each of its residues by Normal Mode Analysis (NMA) of protein. From the Figure 10, the NMA of docked complex has been shown where the value of NMA is almost 0.8 or less than 1.0 for all residues for L07, L08, L09, and D1, that ensure the stability and validation of docking procedure.

FIGURE 10
www.frontiersin.org

Figure 10 B-factor from Molecular Dynamic Study.

ADME and toxicity data

The present system of advanced manufacturing drug development includes the use of in silico techniques as an essential element. These tools are used to explore the absorption, distribution, metabolism, and excretion (ADME-PK) features of novel chemical substances. So, this ADMET analysis was performed to examine both the pharmacokinetic (PK) features and the toxicity of the compound. Table 5 provides an overview of the compiled data for ADMET prediction. The Water solubility Log S standard score is considered to be from -4.0 to -6.0 and -2.0 to -4.0, respectively, for minimum and maximum solubility substances (Rout et al., 2022). Our finding log S is predicted to be -0.152 to -3.543, which falls in the range of maximum aqueous solubility. Substances (L01, L04, L06, L07, and L08) have a high rate of absorption, according to the results of the human intestinal absorption test (HIA) (HIA of more than 30% is considered a high absorption rate) (Pires et al., 2015). Additionally, the Caco-2 cell permeability measure was used to assess the absorption of the substances that were investigated. According to the findings, only three of the molecules exhibit Caco-2 cell permeability in a favorable manner. The VDSs (human) have been obtained from -0.095 to 0.519, which indicates better distribution. Consequently, no single compound can cross the BBB. Furthermore, the inhibitory or substrate response of the cytochrome P450 enzymes was used to confirm the metabolic pathways of the substances that were under investigation (CYPs). This enzyme is crucial to the oxidation process and helps ease the removal of foreign organic molecules like medicines. It also serves a crucial function in the production of energy in the cell. Ligands 04, 06 and 07 found to be actively inhibited the CYP450 1A2 Inhibitor, while no compounds were found to be CYP450 2C9 Substrate. In excretion prediction, the rate of total clearance is about 0.066 ml/min/kg to -6.618 ml/min/kg, which indicates a better clearance rate than clearance rate or negative score. The results of the toxicity tests, which included the AMES toxicity test and the hepatotoxicity test, suggest that the anticipated chemicals are safe, excluding 04, 06, and 07. On the other hand, the maximum tolerated dose in humans was in the range of 0.295–2.193 mg/kg/day.

TABLE 5
www.frontiersin.org

Table 5 ADME and Toxicity prediction.

Conclusion

Kaempferol-O-rhamnoside derivatives (L01-L09) have been studied for their potential as agonists to suppress Mpox and MARV infections; these compounds have been modeled using computational tools using DFT method, as well as molecular docking, molecular dynamics simulation, and in silico study.

In addition, investigations are conducted on pharmacokinetics, ADMET, drug-likeness, HOMO/LUMO gap, molecular electrostatic potential, and PASS prediction in an effort to offer them with an effective antiviral drug against Mpox and MARV infections. With beginning, the PASS prediction investigation outcome obtained at 0.379< Pa<0.712 for anti-viral (Influenza), 0.423< Pa<0.655 for anti-bacterial, 0.632< Pa<0.755 for antifungal, and 0.214< Pa<0.391 for anti-parasitic, so it can be said that they have a high possibility against viral pathogens. Next, molecular docking has been described as the most important part of this investigation, and it has achieved an acceptable docking score, which was even higher than the FDA-approved anti-viral Cidofovir drug. In accumulation, there are formed the hydrogen bonds after docking that conveys the more stability of protein as inhibitor. All of designed drugs showed the negative value of BBB permeability, CYP4501A2 inhibitor, CYP4502C9 substrate and Renal OCT2 substrate even they have no AMES toxicity. In case of quantum calculation, the energy gap is about 9.00 for 1st four ligands (L01-L04) which is more reactive, and it is slightly lower for next class that is about 7.0 or below. However, the softness of L01 –L04 is about 0.21 where the softness for other is about 0.25 to 0.27. However, the 2nd series of ligand (L05-L09) after joining is more active in case of docking but chemical more stable where L07, L08 and L09 show higher binding score. So, it can be said that with increasing the chain length, the binding capacity increases.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding authors.

Author contributions

AK, NM- Conceptualization, study design, writing, data acquisition. MA- Doing the simulation, DFT calculation and data analysis. AC, SM, AG, UC, SA, AK, AAla- Writing and data acquisition. AD, RS, AG, AAle and K-TC – Review, editing, and supervision. All authors contributed to the article and approved the submitted version.

Acknowledgments

This work was funded by the Researchers Supporting Project Number (RSP2023R339) at King Saud University, Riyadh 11451, Saudi Arabia. This study was supported by a grant (no. RA112001) from the Tainan Municipal Hospital (managed by Show Chwan Medical Care Corporation), Tainan, Taiwan. The authors extend their appreciation to their respective institutes for providing facilities to conduct this research.

Conflict of interest

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.

Publisher’s note

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.

References

Abir, M. H., Rahman, T., Das, A., Etu, S. N., Nafiz, I. H., Rakib, A., et al. (2022). Pathogenicity and virulence of marburg virus. Virulence 13, 609–633. doi: 10.1080/21505594.2022.2054760

PubMed Abstract | CrossRef Full Text | Google Scholar

Aborode, A. T., Wireko, A. A., Bel-Nono, K. N., Quarshie, L. S., Allison, M., Bello, M. A. (2022). Marburg virus amidst COVID-19 pandemic in Guinea: fighting within the looming cases. Int. J. Health Plann. Manage. 37, 553. doi: 10.1002/hpm.3332

PubMed Abstract | CrossRef Full Text | Google Scholar

Adalja, A., Inglesby, T. (2022). A novel international monkeypox outbreak. Am. Coll. Physicians 175, 1175–1176. doi: 10.7326/M22-1581

CrossRef Full Text | Google Scholar

Adjemian, J., Farnon, E. C., Tschioko, F., Wamala, J. F., Byaruhanga, E., Bwire, G. S., et al. (2011). Outbreak of marburg hemorrhagic fever among miners in kamwenge and ibanda districts, Uganda, 2007. J. Infect. Dis. 204, S796–S799. doi: 10.1093/infdis/jir312

PubMed Abstract | CrossRef Full Text | Google Scholar

Ahmed, M., Naseer, H., Arshad, M., Ahmad, A. (2022). Monkeypox in 2022: a new threat in developing. Ann. Med. Surg. 78, 103975. doi: 10.1016/j.amsu.2022.103975

CrossRef Full Text | Google Scholar

Akash, S., Emran, T. B., Chopra, H., Dhama, K. (2023). Re-emerging of marburg virus: warning about its virulence and potential impact on world’s health. Int. J. Surg. 109, 165–166. doi: 10.1097/JS9.0000000000000162

PubMed Abstract | CrossRef Full Text | Google Scholar

Akash, S., Rahman, M. M., Islam, M. R., Sharma, R. (2022). Emerging global concern of langya henipavirus: pathogenicity, virulence, genomic features, and future perspectives. J. Med. Virol. 95 (1), e28127. doi: 10.1002/jmv.28127

PubMed Abstract | CrossRef Full Text | Google Scholar

Ali, Y., Imtiaz, H., Tahir, M. M., Gul, F., Saddozai, U. A. K., ur Rehman, A., et al. (2023). "Fragment-based approaches identified tecovirimat-competitive novel drug candidate for targeting the F13 protein of the monkeypox virus,". Viruses 15, 570. doi: 10.3390/v15020570

PubMed Abstract | CrossRef Full Text | Google Scholar

Anwar, F., Haider, F., Khan, S., Ahmad, I., Ahmed, N., Imran, M., et al. (2023). Clinical manifestation, transmission, pathogenesis, and diagnosis of monkeypox virus: a comprehensive review. Life 13, 522. doi: 10.3390/life13020522

PubMed Abstract | CrossRef Full Text | Google Scholar

Balter, M. (2000). Emerging diseases. On the trail of Ebola and Marburg viruses. Science 290 (5493), 923–925. doi: 10.1126/science.290.5493.923

CrossRef Full Text | Google Scholar

Barroso da Silva, F. L., Carloni, P., Cheung, D., Cottone, G., Donnini, S., Foegeding, E. A., et al. (2020). Understanding and controlling food protein structure and function in foods: perspectives from experiments and computer simulations. Annu. Rev. Food Sci. Technol. 11, 365–387. doi: 10.1146/annurev-food-032519-051640

PubMed Abstract | CrossRef Full Text | Google Scholar

Becker, S. C. (2021). Untersuchung zur transmission, virusverbreitung und pathogenität der zecken-übertragenen flaviviren tick-borne encephalitis virus und langat virus in mäusen. doi: 10.25673/39254

CrossRef Full Text | Google Scholar

Beer, E. M., Rao, V. B. (2019). A systematic review of the epidemiology of human monkeypox outbreaks and implications for outbreak strategy. PloS Negl. Trop. Dis. 13, e0007791. doi: 10.1371/journal.pntd.0007791

PubMed Abstract | CrossRef Full Text | Google Scholar

Boulant, S., Stanifer, M., Lozach, P.-Y. (2015). Dynamics of virus-receptor interactions in virus binding, signaling, and endocytosis. Viruses 7, 2794–2815. doi: 10.3390/v7062747

PubMed Abstract | CrossRef Full Text | Google Scholar

Chakraborty, C., Bhattacharya, M., Nandi, S. S., Mohapatra, R. K., Dhama, K., Agoramoorthy, G. (2022). Appearance and re-appearance of zoonotic disease during the pandemic period: long-term monitoring and analysis of zoonosis is crucial to confirm the animal origin of SARS-CoV-2 and monkeypox virus. Veterinary Q. 42, 119–124. doi: 10.1080/01652176.2022.2086718

CrossRef Full Text | Google Scholar

Chakraborty, S., Chandran, D., Mohapatra, R. K., Alagawany, M., Nahed, A., Sharma, A. K., et al. (2022). Clinical management, antiviral drugs and immunotherapeutics for treating monkeypox. an update on current knowledge and futuristic prospects. Int. J. Surg. 105, 106847. doi: 10.1016/j.ijsu.2022.106847

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, F., Li, W., Liu, G., Tang, Y. (2013). In silico ADMET prediction: recent advances, current challenges and future trends. Curr. topics medicinal Chem. 13, 1273–1289. doi: 10.2174/15680266113139990033

CrossRef Full Text | Google Scholar

Claverie, J.-M. (2006). Viruses take center stage in cellular evolution. Genome Biol. 7, 1–5. doi: 10.1186/gb-2006-7-6-110

CrossRef Full Text | Google Scholar

Coltart, C. E., Lindsey, B., Ghinai, I., Johnson, A. M., Heymann, D. L. (2017). The Ebola outbreak, 2013–2016: old lessons for new epidemics. Philos. Trans. R. Soc. B: Biol. Sci. 372, 20160297. doi: 10.1098/rstb.2016.0297

CrossRef Full Text | Google Scholar

Daina, A., Michielin, O., Zoete, V. (2017). SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci. Rep. 7, 42717. doi: 10.1038/srep42717

PubMed Abstract | CrossRef Full Text | Google Scholar

Dallakyan, S., Olson, A. J. (2015). Small-molecule library screening by docking with PyRx. Chem. biology: Methods Protoc., 243–250. doi: 10.1007/978-1-4939-2269-7_19

CrossRef Full Text | Google Scholar

Dulin, N., Spanier, A., Merino, K., Hutter, J. N., Waterman, P. E., Lee, C., et al. (2021). Systematic review of marburg virus vaccine nonhuman primate studies and human clinical trials. Vaccine 39, 202–208. doi: 10.1016/j.vaccine.2020.11.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Farahat, R. A., Memish, Z. A. (2022). Re-emergence of marburg virus and monkeypox in the shadow of COVID-19 pandemic: current situation and implications–correspondence. Int. J. Surg. 106, 106923. doi: 10.1016/j.ijsu.2022.106923

PubMed Abstract | CrossRef Full Text | Google Scholar

Fauci, A. S. (2001). Infectious diseases: considerations for the 21st century. Clin. Infect. Dis. 32, 675–685. doi: 10.1086/319235

PubMed Abstract | CrossRef Full Text | Google Scholar

Harapan, H., Setiawan, A. M., Yufika, A., Anwar, S., Wahyuni, S., Asrizal, F. W., et al. (2020). Knowledge of human monkeypox viral infection among general practitioners: a cross-sectional study in Indonesia. Pathog. Global Health 114, 68–75. doi: 10.1080/20477724.2020.1743037

CrossRef Full Text | Google Scholar

Hickman, M. R., Saunders, D. L., Bigger, C. A., Kane, C. D., Iversen, P. L. (2022). The development of broad-spectrum antiviral medical countermeasures to treat viral hemorrhagic fevers caused by natural or weaponized virus infections. PloS Negl. Trop. Dis. 16, e0010220. doi: 10.1371/journal.pntd.0010220

PubMed Abstract | CrossRef Full Text | Google Scholar

Islam, M. R., Akash, S., Rahman, M. M., Sharma, R. (2023). Epidemiology, pathophysiology, transmission, genomic structure, treatment, and future perspectives of the novel marburg virus outbreak. Int. J. Surg. 109, 36–38. doi: 10.1097/JS9.0000000000000096

PubMed Abstract | CrossRef Full Text | Google Scholar

Islam, M. R., Rahman, M. M., Ahasan, M. T., Sarkar, N., Akash, S., Islam, M., et al. (2022). The impact of mucormycosis (black fungus) on SARS-CoV-2-infected patients: at a glance. Environ. Sci. pollut. Res. 29, 69341–69366. doi: 10.1007/s11356-022-22204-8

CrossRef Full Text | Google Scholar

Joseph, D., Namboothiri, R. P. (2015). Molecular docking analysis of natural ligands against infectious proteins of human immunodeficiency virus and hepatitis b virus. AQUINAS J. OF Multidiscip. Res. 1, 60.

Google Scholar

Kawsar, S. M., Kumer, A., Munia, N. S., Hosen, M. A., Chakma, U., Akash, S. (2022). Chemical descriptors, PASS, molecular docking, molecular dynamics and ADMET predictions of glucopyranoside derivatives as inhibitors to bacteria and fungi growth. Organic Commun. 15, 203. doi: 10.25135/acg.oc.122.2203.2397

CrossRef Full Text | Google Scholar

Khan, A., Adil, S., Qudsia, H. A., Waheed, Y., Alshabrmi, F. M., Wei, D.-Q. (2023). Structure-based design of promising natural products to inhibit thymidylate kinase from monkeypox virus and validation using free energy calculations. Comput. Biol. Med., 106797. doi: 10.1016/j.compbiomed.2023.106797

PubMed Abstract | CrossRef Full Text | Google Scholar

Khani, E., Afsharirad, B., Entezari-Maleki, T. (2023). Monkeypox treatment: current evidence and future perspectives. J. Med. Virol. 95, e28229. doi: 10.1002/jmv.28229

PubMed Abstract | CrossRef Full Text | Google Scholar

Khattak, S., Qaisar, M., Zaman, S., Khan, T. A., Ali, Y., Wu, D.-D., et al. (2022). Monkeypox virus preparation in Pakistan-next viral zoonotic disease outbreak after COVID-19? Biomed. Lett. 8, 196–201. doi: 10.47262/BL/8.2.20220610

CrossRef Full Text | Google Scholar

Khattak, S., Rauf, M. A., Ali, Y., Yousaf, M. T., Liu, Z., Wu, D.-D., et al. (2023). The monkeypox diagnosis, treatments and prevention: a review. Front. Cell. Infection Microbiol. 12, 2005. doi: 10.3389/fcimb.2022.1088471

CrossRef Full Text | Google Scholar

Knust, B., Schafer, I. J., Wamala, J., Nyakarahuka, L., Okot, C., Shoemaker, T., et al. (2015). Multidistrict outbreak of marburg virus disease–Uganda, 2012. J. Infect. Dis. 212, S119–S128. doi: 10.1093/infdis/jiv351

PubMed Abstract | CrossRef Full Text | Google Scholar

Kortepeter, M. G., Dierberg, K., Shenoy, E. S., Cieslak, T. J. (2020). Marburg virus disease: a summary for clinicians. Int. J. Infect. Dis. 99, 233–242. doi: 10.1016/j.ijid.2020.07.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumer, A., Chakma, U., Islam, M., Howlader, D., Hossain, T. (2021a). The computational investigation of sixteen antiviral drugs against main protease (mpro) and spike protease (spro) of sars-cov-2. J. Chilean Chem. Soc. 66, 5339–5351. doi: 10.4067/s0717-97072021000405339

CrossRef Full Text | Google Scholar

Kumer, A., Chakma, U., Matin, M. M. (2022a). Bilastine based drugs as SARS-CoV-2 protease inhibitors: molecular docking, dynamics, and ADMET related studies. Orbital: Electronic J. Chem., 15–23. doi: 10.17807/orbital.v14i1.1642

CrossRef Full Text | Google Scholar

Kumer, A., Chakma, U., Matin, M. M., Akash, S., Chando, A., Howlader, D. (2021b). The computational screening of inhibitor for black fungus and white fungus by d-glucofuranose derivatives using in silico and SAR study. Organic Commun. 14, 305–322. doi: 10.25135/acg.oc.116.2108.2188

CrossRef Full Text | Google Scholar

Kumer, A., Chakma, U., Rana, M. M., Chandro, A., Akash, S., Elseehy, M. M., et al. (2022b). Investigation of the new inhibitors by sulfadiazine and modified derivatives of α-d-glucopyranoside for white spot syndrome virus disease of shrimp by in silico: quantum calculations, molecular docking, ADMET and molecular dynamics study. Molecules 27, 3694. doi: 10.3390/molecules27123694

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumer, A., Khan, M. W. (2021). The effect of alkyl chain and electronegative atoms in anion on biological activity of anilinium carboxylate bioactive ionic liquids and computational approaches by DFT functional and molecular docking. Heliyon 7, e07509. doi: 10.1016/j.heliyon.2021.e07509

PubMed Abstract | CrossRef Full Text | Google Scholar

Ladnyj, I., Ziegler, P., Kima, E. (1972). A human infection caused by monkeypox virus in basankusu territory, democratic republic of the Congo. Bull. World Health Organ. 46, 593.

PubMed Abstract | Google Scholar

Lagunin, A., Stepanchikova, A., Filimonov, D., Poroikov, V. (2000). PASS: prediction of activity spectra for biologically active substances. Bioinformatics 16, 747–748. doi: 10.1093/bioinformatics/16.8.747

PubMed Abstract | CrossRef Full Text | Google Scholar

León-Figueroa, D. A., Bonilla-Aldana, D. K., Pachar, M., Romaní, L., Saldaña-Cumpa, H. M., Anchay-Zuloeta, C., et al. (2022). The never-ending global emergence of viral zoonoses after COVID-19? the rising concern of monkeypox in Europe, north America and beyond. Travel Med. Infect. Dis. 49, 102362. doi: 10.1016/j.tmaid.2022.102362

PubMed Abstract | CrossRef Full Text | Google Scholar

Marie, S. A. (1957). "Viruses: are they alive? New Scholasticism 31, 297–316. doi: 10.5840/newscholas195731343

CrossRef Full Text | Google Scholar

Meng, X.-Y., Zhang, H.-X., Mezei, M., Cui, M. (2011). Molecular docking: a powerful approach for structure-based drug discovery. Curr. computer-aided Drug design 7, 146–157. doi: 10.2174/157340911795677602

CrossRef Full Text | Google Scholar

Minasov, G., Inniss, N. L., Shuvalova, L., Anderson, W. F., Satchell, K. J. (2022). Structure of the monkeypox virus profilin-like protein A42R reveals potential functional differences from cellular profilins. Acta Crystallographica Section F: Struct. Biol. Commun. 78, 371–377. doi: 10.1107/S2053230X22009128

CrossRef Full Text | Google Scholar

Mohapatra, R. K., Tuli, H. S., Sarangi, A. K., Chakraborty, S., Chandran, D., Chakraborty, C., et al. (2022). Unexpected sudden rise of human monkeypox cases in multiple non-endemic countries amid COVID-19 pandemic and salient counteracting strategies: another potential global threat? Int. J. Surg. (London England) 103, 106705. doi: 10.1016/j.ijsu.2022.106705

CrossRef Full Text | Google Scholar

Mooers, B. H., Brown, M. E. (2021). Templates for writing PyMOL scripts. Protein Sci. 30, 262–269. doi: 10.1002/pro.3997

PubMed Abstract | CrossRef Full Text | Google Scholar

Nath, A., Kumer, A., Zaben, F., Khan, M. (2021). Investigating the binding affinity, molecular dynamics, and ADMET properties of 2, 3-dihydrobenzofuran derivatives as an inhibitor of fungi, bacteria, and virus protein. Beni-Suef Univ. J. Basic Appl. Sci. 10, 1–13. doi: 10.1186/s43088-021-00117-8

CrossRef Full Text | Google Scholar

Nguyen, H. D. (2023). In silico identification of potential effects of natural compounds on monkeypox. Phytotherapy Res. doi: 10.1002/ptr.7733

CrossRef Full Text | Google Scholar

Nyakarahuka, L., Ojwang, J., Tumusiime, A., Balinandi, S., Whitmer, S., Kyazze, S., et al. (2017). Isolated case of marburg virus disease, Kampala, Uganda, 2014. Emerging Infect. Dis. 23, 1001. doi: 10.3201/eid2306.170047

CrossRef Full Text | Google Scholar

Nyakarahuka, L., Shoemaker, T. R., Balinandi, S., Chemos, G., Kwesiga, B., Mulei, S., et al. (2019). Marburg virus disease outbreak in kween district Uganda, 2017: epidemiological and laboratory findings. PloS Negl. Trop. Dis. 13, e0007257. doi: 10.1371/journal.pntd.0007257

PubMed Abstract | CrossRef Full Text | Google Scholar

Perepichka, D. F., Bryce, M. R. (2005). Molecules with exceptionally small HOMO–LUMO gaps. Angewandte Chemie Int. Edition 44, 5370–5373. doi: 10.1002/anie.200500413

CrossRef Full Text | Google Scholar

Pires, D. E., Blundell, T. L., Ascher, D. B. (2015). pkCSM: predicting small-molecule pharmacokinetic and toxicity properties using graph-based signatures. J. medicinal Chem. 58, 4066–4072. doi: 10.1021/acs.jmedchem.5b00104

CrossRef Full Text | Google Scholar

Poroikov, V., Filimonov, D. (2005). “PASS: prediction of biological activity spectra for substances,” in Predictive toxicology (CRC Press), 471–490.

Google Scholar

Pourhajibagher, M., Bahador, A. (2023). Virtual screening and computational simulation analysis of antimicrobial photodynamic therapy using propolis-benzofuran a to control of monkeypox. Photodiagnosis Photodyn. Ther. 41, 103208. doi: 10.1016/j.pdpdt.2022.103208

PubMed Abstract | CrossRef Full Text | Google Scholar

Rahman, M. A., Chakma, U., Kumer, A., Rahman, M. R., Matin, M. M. (2022). Uridine-derived 4-aminophenyl 1-thioglucosides: DFT optimized FMO, ADME, and antiviral activities study. Biointerface Res. Appl. Chem. 13, 1–15. doi: 10.33263/BRIAC131.052

CrossRef Full Text | Google Scholar

Rahman, M. M., Karim, M. R., Ahsan, M. Q., Khalipha, A. B. R., Chowdhury, M. R., Saifuzzaman, M. (2012). Use of computer in drug design and drug discovery: a review. Int. J. Pharm. Life Sci. 1. doi: 10.3329/ijpls.v1i2.12955

CrossRef Full Text | Google Scholar

Rout, J., Swain, B. C., Tripathy, U. (2022). In silico investigation of spice molecules as potent inhibitor of SARS-CoV-2. J. Biomolecular Structure Dynamics 40, 860–874. doi: 10.1080/07391102.2020.1819879

CrossRef Full Text | Google Scholar

Sah, R., Abdelaal, A., Asija, A., Basnyat, S., Sedhai, Y. R., Ghimire, S., et al. (2022). Monkeypox virus containment: the application of ring vaccination and possible challenges. J. Travel Med. 29, taac085. doi: 10.1093/jtm/taac085

PubMed Abstract | CrossRef Full Text | Google Scholar

Shanmugaraj, B., Khorattanakulchai, N., Phoolcharoen, W. (2022). Emergence of monkeypox: another concern amidst COVID-19 crisis. Asian Pacific J. Trop. Med. 15, 193. doi: 10.4103/1995-7645.346081

CrossRef Full Text | Google Scholar

Sharma, A., Fahrni, M. L., Choudhary, O. P. (2022). Monkeypox outbreak: new zoonotic alert after the COVID-19 pandemic. Int. J. Surg. (London England) 104, 106812. doi: 10.1016/j.ijsu.2022.106812

CrossRef Full Text | Google Scholar

Taseen, S., Nasir, F., Abbas, M., Altaf, M., Asghar, M. S., Tahir, M. J. (2022). Post-pandemic world at the mercy of monkeypox virus outbreak: time to worry or not? J. Med. Virol. 95 (1), e27948. doi: 10.1002/jmv.27948

PubMed Abstract | CrossRef Full Text | Google Scholar

Velavan, T. P., Meyer, C. G. (2022). Monkeypox 2022 outbreak: an update. Trop. Med. Int. Health 27, 604–605. doi: 10.1111/tmi.13785

PubMed Abstract | CrossRef Full Text | Google Scholar

Whittaker, G. R., Helenius, A. (1998). Nuclear import and export of viruses and virus genomes. Virology 246, 1–23. doi: 10.1006/viro.1998.9165

PubMed Abstract | CrossRef Full Text | Google Scholar

Zaeck, L. M., Lamers, M. M., Verstrepen, B. E., Bestebroer, T. M., van Royen, M. E., Götz, H., et al. (2023). Low levels of monkeypox virus-neutralizing antibodies after MVA-BN vaccination in healthy individuals. Nat. Med. 29, 270–278. doi: 10.1038/s41591-022-02090-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, A. P., Bornholdt, Z. A., Abelson, D. M., Saphire, E. O. (2014). Crystal structure of marburg virus VP24. J. Virol. 88, 5859–5863. doi: 10.1128/JVI.03565-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, F., He, Y., Lu, H. (2022). Marburg virus disease: a deadly rare virus is coming. BioScience Trends 16, 312–316. doi: 10.5582/bst.2022.01333

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, C., Wang, F., Tang, B., Han, J., Li, X., Lian, G., et al. (2021). Anti-inflammatory effects of kaempferol-3-O-rhamnoside on HSV-1 encephalitis in vivo and in vitro. Neurosci. Lett. 765, 136172. doi: 10.1016/j.neulet.2021.136172

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Monkeypox virus, Marburg virus, drug development, O-rhamnosides, Kaempferol-o-rhamnosides, admet, molecular docking, molecular dynamics simulations

Citation: Al Mashud MA, Kumer A, Mukerjee N, Chandro A, Maitra S, Chakma U, Dey A, Akash S, Alexiou A, Khan AA, Alanazi AM, Ghosh A, Chen K-T and Sharma R (2023) Mechanistic inhibition of Monkeypox and Marburg virus infection by O-rhamnosides and Kaempferol-o-rhamnosides derivatives: a new-fangled computational approach. Front. Cell. Infect. Microbiol. 13:1188763. doi: 10.3389/fcimb.2023.1188763

Received: 17 March 2023; Accepted: 26 April 2023;
Published: 24 May 2023.

Edited by:

Saadullah Khattak, Henan University, China

Reviewed by:

Niraj Kumar Jha, Sharda University, India
Daniel Uti, Federal University of Health Sciences Otukpo Benue State, Nigeria
Andrew Awuah Wireko, Sumy State University, Ukraine

Copyright © 2023 Al Mashud, Kumer, Mukerjee, Chandro, Maitra, Chakma, Dey, Akash, Alexiou, Khan, Alanazi, Ghosh, Chen and Sharma. 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.

*Correspondence: Rohit Sharma, rohitsharma@bhu.ac.in; Kow-Tong Chen, ktchen@mail.ncku.edu.tw

These authors have contributed equally to this work and share first authorship

Disclaimer: 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.