Next Article in Journal
Pro-Inflammatory and Pro-Apoptotic Effects of the Non-Protein Amino Acid L-Azetidine-2-Carboxylic Acid in BV2 Microglial Cells
Next Article in Special Issue
Protective Effects of One 2,4-Dihydro-3H-Pyrazol-3-one Derivative against Posterior Capsular Opacification by Regulation of TGF-β2/SMADs and Non-SMAD Signaling, Collagen I, and Fibronectin Proteins
Previous Article in Journal
Dynamic Expression of Palmitoylation Regulators across Human Organ Development and Cancers Based on Bioinformatics
Previous Article in Special Issue
Evaluation of the Regenerative Potential of Platelet-Lysate and Platelet-Poor Plasma Derived from the Cord Blood Units in Corneal Wound Healing Applications: An In Vitro Comparative Study on Corneal Epithelial Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Proposed Mechanism for Emodin as Agent for Methicillin Resistant Staphylococcus Aureus: In Vitro Testing and In Silico Study

by
Mohammed M. Ghoneim
1,2,†,
Rasha Hamed Al-Serwi
3,†,
Mohamed El-Sherbiny
4,
El-sayed M. El-ghaly
1,
Amal E. Hamad
5,
Mohamed A. Abdelgawad
6,
Ehab A. Ragab
1,
Sarah I. Bukhari
7,
Khulud Bukhari
8,
Khaled Elokely
9,* and
Manal A. Nael
9,10,*
1
Pharmacognosy and Medicinal Plants Department, Faculty of Pharmacy, Al-Azhar University, Cairo 11884, Egypt
2
Department of Pharmacy Practice, College of Pharmacy, AlMaarefa University, Riyadh 11597, Saudi Arabia
3
Department of Basic Dental Sciences, College of Dentistry, Princess Nourah bint Abdulrahman University, Riyadh 11671, Saudi Arabia
4
Department of Basic Medical Sciences, College of Medicine, AlMaarefa University, Riyadh 11597, Saudi Arabia
5
Department of Pharmaceutical Analytical Chemistry, Faculty of Pharmacy, Menoufia University, Shebin El-Koum 32511, Egypt
6
Department of Pharmaceutical Chemistry, College of Pharmacy, Jouf University, Sakaka 72341, Saudi Arabia
7
Department of Pharmaceutics, College of Pharmacy, King Saud University, Riyadh 11451, Saudi Arabia
8
Department of Microbiology and Parasitology, College of Veterinary Medicine, King Faisal University, Hofuf, Al-Ahsa 36362, Saudi Arabia
9
Department of Chemistry, Institute for Computational Molecular Science, Temple University, Philadelphia, PA 19122, USA
10
Department of Pharmaceutical Chemistry, Faculty of Pharmacy, Tanta University, Tanta 31527, Egypt
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Curr. Issues Mol. Biol. 2022, 44(10), 4490-4499; https://doi.org/10.3390/cimb44100307
Submission received: 5 September 2022 / Revised: 19 September 2022 / Accepted: 24 September 2022 / Published: 27 September 2022

Abstract

:
In the search for a new anti-MRSA lead compound, emodin was identified as a good lead against methicillin-resistant Staphylococcus aureus (MRSA). Emodin serves as a new scaffold to design novel and effective anti-MRSA agents. Because rational drug discovery is limited by the knowledge of the drug target, α-hemolysin of Staphylococcus aureus was used in this study because it has an essential role in Staphylococcus infections and because emodin shares structural features with compounds that target this enzyme. In order to explore emodin’s interactions with α-hemolysin, all possible ligand binding pockets were identified and investigated. Two ligand pockets were detected based on bound ligands and other reports. The third pocket was identified as a cryptic site after molecular dynamics (MD) simulations. MD simulations were conducted for emodin in each pocket to identify the most plausible ligand site and to aid in the design of potent anti-MRSA agents. Binding of emodin to site 1 was most stable (RMSD changes within 1 Å), while in site 2, the binding pose of emodin fluctuated, and it left after 20 ns. In site 3, it was stable during the first 50 ns, and then it started to move out of the binding site. Site 1 is a possible ligand binding pocket, and this study sheds more light on interaction types, binding mode, and key amino acids involved in ligand binding essential for better lead design. Emodin showed an IC50 value of 6.3 μg/mL, while 1, 6, and 8 triacetyl emodin showed no activity against MRSA. A molecular modeling study was pursued to better understand effective binding requirements for a lead.

1. Introduction

Humanity is facing many challenges that threaten its own existence. At the time of this study, the world was facing one of its worst nightmares—a global pandemic that disrupted human activities and increased hospitalization and intensive care admission to such a level that it surpassed maximum capacity many times [1,2,3]. This increased hospitalization brings forward the looming risks of nosocomial infections that could complicate cases, prolong hospital stays and even increase fatalities. Amongst the most notorious nosocomial infections, MRSA (methicillin-resistant Staphylococcus aureus) is a top assailant [4]. Preventive measures have led to a continuous decline in healthcare-related MRSA invasive infections in developed countries; nevertheless, this decline is now slowing down. Moreover, community-based infections, which constitute a reservoir, are not decreasing, and many developing countries still have a stubbornly high prevalence [5]. The high virulence of MRSA and the increased difficulty of treating it leads to a high fatality rate [6]. Another worrying feature is the ease with which it acquires resistance to new antibiotics making its treatment more difficult by the day and increasing its mortality [7].
Amongst its vast and versatile arsenal, cytolytic toxins are especially important for its virulence. α-toxin is an important one that helps MRSA kill leukocytes [8]. This toxin, also referred to as α-hemolysin [9], is a host-specific protein that binds to the cell walls of leukocytes and perforates them, causing cell lysis, thus representing an important therapeutic target in the fight against MRSA. The structure and mechanism of action of α-hemolysin, therefore, have been widely studied [10,11,12], and it was selected for our molecular modeling study because it is the target of emodin-like compounds [13]. Other anthraquinones have been reported for their anti-MRSA activity [14,15].
Emodin was identified as the active ingredient from Rhizoma Polygoni Cuspidati, which is responsible for significant anti-MRSA activity in the in vitro antibacterial experiments. The study suggested that emodin could damage cell wall integrity [16]. Emodin did not show significant cytotoxicity against human cells. The synthesis and antibacterial evaluations were reported for emodin and several analogs, indicating the importance of the hydroxyl and methyl groups in the emodin scaffold for the anti-MRSA activity [17]. The mechanism of action of emodin is not clear. Proteomic studies were used to explore the molecular mode of action of emodin, showing significant changes in several MRSA proteins, which indicates that emodin could exert its effects by multiple mechanisms [18].
Molecular modeling is used in this work to investigate the binding of emodin and related compounds to three proposed ligand sites in α-hemolysin. Molecular Docking and Molecular Dynamics (MD) simulations were conducted to gain a better understanding of what makes a good fit for α-hemolysin active sites. This is an essential step if we are to successfully develop effective anti α-hemolysin therapeutics that increases activity against MRSA.

2. Materials and Methods

Antimicrobial Assay

Emodin (Compound 1) and 1, 6, 8 triacetyl emodin (Compound 2) were tested for antimicrobial activity. All organisms were obtained from the American Type Culture Collection (Manassas, VA, USA). Fungal strains: Candida albicans ATCC 90028, Cryptococcus neoformans ATCC 90113, Aspergillus fumigatus ATCC 90906, methicillin-resistant Staphylococcus aureus ATCC 43300 (MRSA), Escherichia coli ATCC 35218, Pseudomonas aeruginosa ATCC 27853 and Mycobacterium intracellulare ATCC 23068. Ciprofloxacin and amphotericin B were used as positive controls for bacteria and fungi, respectively [19].
The following protocol was followed: The American Type Culture Collection was used to obtain all of the organisms (Manassas, VA, USA). Optical density was employed to track growth for all species, and a modified version of the Clinical and Laboratory Standards Institute (formerly NCCLS) procedures was utilized for susceptibility testing. For growth detection of M. intracellulare and A. fumigatus, 5% Alamar Blue (Biosource International, Camarillo, CA, USA) was added to media (Biosource International, Camarillo, CA, USA). Samples were duplicated and transferred to 96-well flat bottom microplates after being serially diluted in 20% DMSO/saline. By adjusting the OD630 of the cell/spore suspension in the RPMI at pH 4.5 incubation broth for Candida albicans, microbial inocula were created. For Cryptococcus neoformans, Sabouraud dextrose was used, followed by cation-adjusted Muller–Hinton at a pH 7.3 for MRSA, 5% Alamar Blue (Biosource International, Camarillo, CA, USA) in Middlebrook 7H9 broth with OADC enrichment, pH = 7.3 for M. intracellulare, and 5% Alamar Blue in RPMI at pH 7.3 for A. fumigatus. Each assay contained drug controls, ciprofloxacin (ICN Biomedicals, Solon, OH, USA) for bacteria and amphotericin B (ICN Biomedicals, OH) for fungi. Prior to and following incubation, all organisms were read at 630 nm using a Polarstar Galaxy plate reader (BMG Lab Technologies, Ortenberg, Germany) or 544 ex/590 cm (M. intracellulare, A. fumigatus) using a Biotek Powerwave XS plate reader (Bio-Tek Instruments, Winooski, VT, USA). By taking 5 L from each clear well, putting it to agar, and letting it sit there until growth appeared, it was possible to calculate the minimum fungicidal or bactericidal concentrations. The MFC/MBC is defined as the test concentration at which the organism dies at the lowest level (allows no growth on agar). Samples are first examined in duplicate at a concentration of 50 µg/mL, and percent inhibitions (% inh.) are estimated in relation to the negative and positive controls. Extracts that exhibit inhibition of 50% move on to the secondary assay. Samples dissolved to 20 mg/mL are tested in the Secondary Assay at 200, 40, 8 µg/mL, while samples dissolved to 2 mg/mL are tested at 20, 4, 0.8 g/mL. The Tertiary Assay is used for pure substances with an IC50 of less than 7 µg/mL in the Secondary Assay. The Tertiary Assay involves testing pure chemicals against microorganisms at concentrations of 20, 10, 5.0, …, 0.02 µg/mL and calculating IC50s. Using the XLFit fit curve fitting program, all IC50s are determined. The antifungal drug control is amphotericin B, and the antibacterial drug control is ciprofloxacin.

3. Molecular Modeling

Protein Preparation. The crystal structure of α-hemolysin of Staphylococcus aureus (PDB accession code: 6U3T) was obtained from the RSCB protein data bank repository [20]. The Protein Preparation Wizard (PrepWizard) module of Schrodinger suite [21,22] was used to prepare the protein structure. Missing hydrogen atoms were added, and atom bond orders were assigned. Prime was used to fill in missing residues and loops. Protonation and tautomerization states of amino acid residues were assigned for pH 7.4. Water molecules that are not forming at least three hydrogen bonds to non-waters were deleted. Hydrogen bond networks were optimized for pH 7.4. Then, the protein structure was minimized with the OPLS3 forcefield [23].
Molecular Docking. Three possible ligand binding sites were used to prepare the receptor grids for docking [24]. Site 1 was defined using bound ligand coordinates. The center of which was determined by the following amino acids: Asn176, Gln177, Trp179, Tyr182, Gln194, Leu195, Met197, Lys198, Thr199 and Agr200. Site 2 was defined according to published reports by the center of the following amino acids: Tyr102, Ar104, Asn105, Ile107, Asp108, Thr109, Glu111, Tyr112, Met113, Ser114, Gly126, Asp127, Asp128, and Ile132. Site 3, in close proximity to site 2, was determined from the molecular dynamics runs, and it was defined as the center of: Ala62, Gly63, Thr60, Gln64, Arg66, Tyr28, Asp29, Lys30, Gly33, Met34, Ala35, Ser221, Ser222 and Gly223. The docking grids were built by Glide [24,25,26,27]. LigPrep [28] was used to prepare ligand structures for docking. Standard precision (SP) docking algorithm in Glide was used, and softening of the receptor potential was allowed to add more receptor flexibility during the docking step. The per-atom van der Waals radii and charges were scaled to 0.85. The five best poses were kept for further calculations. Prime [29,30,31] was used to refine the docking poses, and more receptor flexibility was allowed for residues within 10 Å of docked ligand. Prime MM-GBSA was then used to calculate the ligand binding free energy (∆G). Another run of induced fit docking (IFD) was performed using the same amino acids as the center for each site [32,33,34,35].
Molecular Dynamics Simulations. MD simulations were carried out in triplicates. Three MD systems were built (Table 1). Selected docking poses were solvated in the TIP4P water box [36] and then neutralized with Na+ ions. DESMOND [37] with the OPLS3 force field was used to run the simulations. NPT ensemble at a constant pressure of 1 bar and temperature of 300 K was set for the calculations by using the Nosé–Hoover chain and Martyna–Tobias–Klein coupling schemes [38,39,40]. Long-range electrostatic energy and forces were calculated using particle-mesh-based Ewald method [41], and a tolerance was set to 1e−9. A variant of the SHAKE algorithm was used to allow the time step to be increased. A time-scale splitting (RESPA-based) integrator was utilized for the numerical integration [37]. The short range and bonded interaction were updated every 2ps. The long-range and non-bonded interactions were revised every 6ps. A cutoff of 9.0 Å was defined for the short-range Coulomb interactions.

4. Results

α-Hemolysin was identified as a target for MRSA because it plays an essential role in S. aureus infections. The X-ray crystal structure of the α-hemolysin protein was resolved as a monomer, which assembles to form a heptameric pore (Figure 1). The assembly is a mushroom-shaped structure. The secondary structure of each monomer is made up of twelve β-sheets and one α-helix.
Compound 1 (Figure 2): emodin (Sigma Aldrich) was evaluated for its antimicrobial activity against MRSA, and it exhibited an inhibitory activity with an IC50 value of 6.3 μg/mL (Table 2). Emodin was subjected to acetylation using pyridine and acetic anhydride to give 1, 6, 8 triacetyl emodin (compound 2), which led to loss of the activity against MRSA.
To explore a good understanding of the role of hydroxyl groups; a computational study has been performed for both emodin and 1, 6, 8 triacetyl emodin.
The ligand binding pockets in α-hemolysin were explored and defined based on the atomic data of protein complexed with ligands and previously studied sites. Site 1 was determined using ligand coordinates and the center of amino acids that surround the ligand. Site 2 was determined from previous reports [20]. Emodin docked with high affinity into site 1 and site 2. Because emodin is a rigid structure and the identified binding cavities are not wide, an induced fit docking protocol was attempted to provide an insight into the most reliable binding mode in each site. In the case of site 1 (Figure 3), emodin was docked with an IFD score of −661.05 kcal/mol. It exhibited hydrogen bonds with Gln194, Met197, and Arg200. It showed π–π stacking with Trp179, and hydrophobic contacts with several amino acids. Eomdin was docked with an IFD score of −660.45 kcal/mol at site 2. It displayed hydrogen bonds with Asn105, Met113, Ser114 and Asp127, and hydrophobic contacts with several other amino acids. Emodin, in site 3, docked with an IFD score of −524.37 kcal/mol, and it exhibited hydrogen bonds with Gln64, Ser221 and Ser222.
Three systems of molecular dynamics (MD) simulations were constructed. The first MD system was built for the best binding mode of emodin in site 1, and the second system was for site 2. Molecular dynamics (MD) simulations were carried out for 100 ns to study the stability of emodin’s binding poses. The root mean square deviations (RMSD) were computed according to the following equation.
R M S D x = 1 N i = 1 N ( r i ( t x ) r i ( t r e f ) ) 2
N: the system’s number of atoms, tref: the reference time, r’i: the position of protein backbone and ligand atoms (i) in frame x, tx: the record time of frame x.
The average change in displacement of emodin’s atoms throughout the MD simulation in site 1 with respect to the original starting point indicated a stable pose in the binding pocket over the course of the MD time. It showed an RMSD value of ~2.0 Å compared to the 1.5 Å of the protein peptide chain (Figure 4A). While in site 2, the binding pose of emodin fluctuated with an average RMSD value of ~40.0 Å with respect to that of the protein backbone of 2.4 Å (Figure 4B). Emodin left site 2 after 20 ns. Analysis of the MD simulations showed that the loop, which is composed of amino acids 118 through 143, fluctuated the most during the simulation. Emodin left site 2 and occupied another cryptic cavity that may accommodate ligands. Therefore, emodin was docked into site 3 and the best pose was used to construct the third MD system. Emodin fluctuated throughout the MD simulations with an average RMSD of 3.0 Å with respect to the protein backbone of 2.0 Å. It was stable during the first 50 ns, and then it started to move out of the binding site.
Emodin is a planar molecule that contains several polar atoms. Interactions of emodin with protein were monitored throughout the simulation. Emodin interacts with the amino acids in site 1 through hydrogen bonds, π–π contact and water-mediated interactions (Figure 5). On average, emodin has three contacts (Figure 6). The most stable contacts were indicated with Met 197 and Arg200, followed by Asn178, Trp179, and Gly180 (Figure 6).
In site 2, emodin did not show stable contact with the amino acids in this site. Instead, it displayed interactions with amino acids in site 2 and site 3 as well (Figure S7). During the first 20 ns, only Met113 showed a stable contact, which was lost during the second 80 ns. Emodin started to form more contacts with Tyr65, Val67, Asp212, and Pro213 during the last 80 ns (Figure S8). Emodin was more stable in site 3 in comparison to site 2 (Figure S9). Ser222, Ala62, Gly63, and Tyr28 demonstrated the most stable interactions (Figure S10).
Compound 2 did not show any good binding mode in either site. This could be attributed to extra molecular volume added by the acetyl groups and the block of the interactions brought by the hydroxyl groups of emodin. The planar features of emodin do not make the compound a perfect lead for modification. Understanding the interaction types, binding mode, and recognizing the key amino acids involved in ligand binding would lead to a better lead design. Modifying the scaffold and keeping the important functional groups of emodin would increase ligand affinity and would provide good candidates to manage resistant Staphylococcus aureus.

5. Conclusions

In this study, emodin was identified as a potential antimicrobial agent against MRSA (IC50 of 6.3 μg/mL). However, being a planar compound, it is not a good candidate for a lead modification. Instead, studying the differences between its binding to α-hemolysin and that of the inactive 1, 6, 8 triacetyl emodin serve to better understand binding modes, sites and interaction types. This is an important step in the quest to develop successful anti-MRSA leads. Three possible binding sites were explored for α-hemolysin using MD simulations. Only site 1 is expected to be the possible ligand binding pocket. Finding out the ligand pocket is essential in the stage of potent lead design and modification to develop potent anti-MRSA agents.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cimb44100307/s1, NMR data and MS data of compounds (1 and 2) are included as Supplementary Material.

Author Contributions

Conceptualization, M.M.G., K.E. and M.A.N.; methodology, All authors; software, K.E. and M.A.N.; validation, All authors; formal analysis, All authors; investigation, All authors; resources, All authors; data curation, All authors; writing—original draft preparation, M.M.G., A.E.H., K.E. and M.A.N.; writing—review and editing, All authors; visualization, All authors. All authors have read and agreed to the published version of the manuscript.

Funding

The work was supported by Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2022R199), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia. This publication was supported by AlMaarefa University researchers supporting program (grant number: MA-006), AlMaarefa University, Riyadh, Saudi Arabia. This project includes calculations carried out on HPC resources supported in part by the National Science Foundation through major research instrumentation grant number 1625061 and by the US Army Research Laboratory under contract number W911NF-16-2-0189.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Marks, K.J.; Whitaker, M.; Anglin, O.; Milucky, J.; Patel, K.; Pham, H.; Chai, S.J.; Kirley, P.D.; Armistead, I.; McLafferty, S. Hospitalizations of children and adolescents with laboratory-confirmed COVID-19—COVID-NET, 14 States, July 2021–January 2022. Morb. Mortal. Wkly. Rep. 2022, 71, 271. [Google Scholar] [CrossRef] [PubMed]
  2. Xie, Y.; Xu, E.; Bowe, B.; Al-Aly, Z. Long-term cardiovascular outcomes of COVID-19. Nat. Med. 2022, 28, 583–590. [Google Scholar] [CrossRef] [PubMed]
  3. Murray, C.J. COVID-19 will continue but the end of the pandemic is near. Lancet 2022, 399, 417–419. [Google Scholar] [CrossRef]
  4. Magill, S.S.; O’Leary, E.; Janelle, S.J.; Thompson, D.L.; Dumyati, G.; Nadle, J.; Wilson, L.E.; Kainer, M.A.; Lynfield, R.; Greissman, S. Changes in prevalence of health care–associated infections in US hospitals. N. Engl. J. Med. 2018, 379, 1732–1744. [Google Scholar] [CrossRef]
  5. Kourtis, A.P.; Hatfield, K.; Baggs, J.; Mu, Y.; See, I.; Epson, E.; Nadle, J.; Kainer, M.A.; Dumyati, G.; Petit, S. Vital signs: Epidemiology and recent trends in methicillin-resistant and in methicillin-susceptible Staphylococcus aureus bloodstream infections—United States. Morb. Mortal. Wkly. Rep. 2019, 68, 214. [Google Scholar] [CrossRef]
  6. Purrello, S.; Garau, J.; Giamarellos, E.; Mazzei, T.; Pea, F.; Soriano, A.; Stefani, S. Methicillin-resistant Staphylococcus aureus infections: A review of the currently available treatment options. J. Glob. Antimicrob. Resist. 2016, 7, 178–186. [Google Scholar] [CrossRef]
  7. Chambers, H.F.; DeLeo, F.R. Waves of resistance: Staphylococcus aureus in the antibiotic era. Nat. Rev. Microbiol. 2009, 7, 629–641. [Google Scholar] [CrossRef]
  8. Lowy, F.D. Staphylococcus aureus infections. N. Engl. J. Med. 1998, 339, 520–532. [Google Scholar] [CrossRef]
  9. Bhakdi, S.; Tranum-Jensen, J. Alpha-toxin of Staphylococcus aureus. Microbiol. Mol. Biol. Rev. 1991, 55, 733–751. [Google Scholar] [CrossRef]
  10. Meesters, C.; Brack, A.; Hellmann, N.; Decker, H. Structural characterization of the α-hemolysin monomer from Staphylococcus aureus. Protein Struct. Funct. Bioinform. 2009, 75, 118–126. [Google Scholar] [CrossRef]
  11. Gouaux, E. α-Hemolysin fromStaphylococcus aureus: An archetype of β-barrel, channel-forming toxins. J. Struct. Biol. 1998, 121, 110–122. [Google Scholar] [CrossRef] [PubMed]
  12. Hua, L.; Hilliard, J.J.; Shi, Y.; Tkaczyk, C.; Cheng, L.I.; Yu, X.; Datta, V.; Ren, S.; Feng, H.; Zinsou, R.; et al. Assessment of an anti-alpha-toxin monoclonal antibody for prevention and treatment of Staphylococcus aureus-induced pneumonia. Antimicrob. Agents Chemother. 2014, 58, 1108–1117. [Google Scholar] [CrossRef]
  13. Tang, F.; Li, L.; Meng, X.-M.; Li, B.; Wang, C.-Q.; Wang, S.-Q.; Wang, T.-L.; Tian, Y.-M. Inhibition of alpha-hemolysin expression by resveratrol attenuates Staphylococcus aureus virulence. Microb. Pathog. 2019, 127, 85–90. [Google Scholar] [CrossRef] [PubMed]
  14. Ghoneim, M.M.; Elokely, K.M.; Atef, A.; Mohammad, A.-E.I.; Jacob, M.; Radwan, M.M.; Doerksen, R.J.; Cutler, S.J.; Ross, S.A. Asphodosides AE, anti-MRSA metabolites from Asphodelus microcarpus. Phytochemistry 2014, 105, 79–84. [Google Scholar] [CrossRef] [PubMed]
  15. Ghoneim, M.M.; Elokely, K.M.; Atef, A.; Abd Elsalam, I.M.; Jacob, M.; Cutler, S.J.; Doerksen, R.J.; Ross, S.A. Isolation and characterization of new secondary metabolites from Asphodelus microcarpus. Med. Chem. Res. 2014, 23, 3510–3515. [Google Scholar] [CrossRef] [PubMed]
  16. Cao, F.; Peng, W.; Li, X.; Liu, M.; Li, B.; Qin, R.; Jiang, W.; Cen, Y.; Pan, X.; Yan, Z. Emodin is identified as the active component of ether extracts from Rhizoma Polygoni Cuspidati, for anti-MRSA activity. Can. J. Physiol. Pharmacol. 2015, 93, 485–493. [Google Scholar] [CrossRef]
  17. Chalothorn, T.; Rukachaisirikul, V.; Phongpaichit, S.; Pannara, S.; Tansakul, C. Synthesis and antibacterial activity of emodin and its derivatives against methicillin-resistant Staphylococcus aureus. Tetrahedron Lett. 2019, 60, 151004. [Google Scholar] [CrossRef]
  18. Ji, X.; Liu, X.; Peng, Y.; Zhan, R.; Xu, H.; Ge, X. Comparative analysis of methicillin-sensitive and resistant Staphylococcus aureus exposed to emodin based on proteomic profiling. Biochem. Biophys. Res. Commun. 2017, 494, 318–324. [Google Scholar] [CrossRef]
  19. Bharate, S.B.; Khan, S.I.; Yunus, N.A.; Chauthe, S.K.; Jacob, M.R.; Tekwani, B.L.; Khan, I.A.; Singh, I.P. Antiprotozoal and antimicrobial activities of O-alkylated and formylated acylphloroglucinols. Bioorg. Med. Chem. 2007, 15, 87–96. [Google Scholar] [CrossRef]
  20. Liu, J.; Kozhaya, L.; Torres, V.J.; Unutmaz, D.; Lu, M. Structure-based discovery of a small-molecule inhibitor of methicillin-resistant Staphylococcus aureus virulence. J. Biol. Chem. 2020, 295, 5944–5959. [Google Scholar] [CrossRef] [Green Version]
  21. Sastry, G.M.; Adzhigirey, M.; Day, T.; Annabhimoju, R.; Sherman, W. Protein and ligand preparation: Parameters, protocols, and influence on virtual screening enrichments. J. Comput.-Aided Mol. Des. 2013, 27, 221–234. [Google Scholar] [CrossRef] [PubMed]
  22. Schrödinger Release 2018-4: Protein Preparation Wizard; Epik, Schrödinger, LLC: New York, NY, USA, 2016; Impact, Schrödinger, LLC: New York, NY, USA, 2016; Prime, Schrödinger, LLC: New York, NY, USA, 2018.
  23. Harder, E.; Damm, W.; Maple, J.; Wu, C.; Reboul, M.; Xiang, J.Y.; Wang, L.; Lupyan, D.; Dahlgren, M.K.; Knight, J.L. OPLS3: A force field providing broad coverage of drug-like small molecules and proteins. J. Chem. Theory Comput. 2016, 12, 281–296. [Google Scholar] [CrossRef]
  24. Schrödinger Release 2018-4: Glide; Schrödinger, LLC: New York, NY, USA, 2018.
  25. Friesner, R.A.; Banks, J.L.; Murphy, R.B.; Halgren, T.A.; Klicic, J.J.; Mainz, D.T.; Repasky, M.P.; Knoll, E.H.; Shelley, M.; Perry, J.K. Glide: A new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. J. Med. Chem. 2004, 47, 1739–1749. [Google Scholar] [CrossRef] [PubMed]
  26. Friesner, R.A.; Murphy, R.B.; Repasky, M.P.; Frye, L.L.; Greenwood, J.R.; Halgren, T.A.; Sanschagrin, P.C.; Mainz, D.T. Extra precision glide: Docking and scoring incorporating a model of hydrophobic enclosure for protein− ligand complexes. J. Med. Chem. 2006, 49, 6177–6196. [Google Scholar] [CrossRef] [PubMed]
  27. Halgren, T.A.; Murphy, R.B.; Friesner, R.A.; Beard, H.S.; Frye, L.L.; Pollard, W.T.; Banks, J.L. Glide: A new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening. J. Med. Chem. 2004, 47, 1750–1759. [Google Scholar] [CrossRef]
  28. Schrödinger Release 2018-4: LigPrep; Schrödinger, LLC: New York, NY, USA, 2018.
  29. Jacobson, M.P.; Pincus, D.L.; Rapp, C.S.; Day, T.J.; Honig, B.; Shaw, D.E.; Friesner, R.A. A hierarchical approach to all-atom protein loop prediction. Proteins Struct. Funct. Bioinform. 2004, 55, 351–367. [Google Scholar] [CrossRef] [PubMed]
  30. Jacobson, M.P.; Friesner, R.A.; Xiang, Z.; Honig, B. On the role of the crystal environment in determining protein side-chain conformations. J. Mol. Biol. 2002, 320, 597–608. [Google Scholar] [CrossRef]
  31. Schrödinger Release 2018-4: Prime; Schrödinger, LLC: New York, NY, USA, 2018.
  32. Farid, R.; Day, T.; Friesner, R.A.; Pearlstein, R.A. New insights about HERG blockade obtained from protein modeling, potential energy mapping, and docking studies. Bioorg. Med. Chem. 2006, 14, 3160–3173. [Google Scholar] [CrossRef]
  33. Sherman, W.; Day, T.; Jacobson, M.P.; Friesner, R.A.; Farid, R. Novel procedure for modeling ligand/receptor induced fit effects. J. Med. Chem. 2006, 49, 534–553. [Google Scholar] [CrossRef]
  34. Sherman, W.; Beard, H.S.; Farid, R. Use of an induced fit receptor structure in virtual screening. Chem. Biol. Drug Des. 2006, 67, 83–84. [Google Scholar] [CrossRef]
  35. Schrödinger Release 2018-4: Induced Fit Docking Protocol; Glide, Schrödinger, LLC: New York, NY, USA, 2018; Prime, Schrödinger, LLC: New York, NY, USA, 2018.
  36. Horn, H.W.; Swope, W.C.; Pitera, J.W.; Madura, J.D.; Dick, T.J.; Hura, G.L.; Head-Gordon, T. Development of an improved four-site water model for biomolecular simulations: TIP4P-Ew. J. Chem. Phys. 2004, 120, 9665–9678. [Google Scholar] [CrossRef] [PubMed]
  37. Bowers, K.J.; Chow, D.E.; Xu, H.; Dror, R.O.; Eastwood, M.P.; Gregersen, B.A.; Klepeis, J.L.; Kolossvary, I.; Moraes, M.A.; Sacerdoti, F.D. Scalable algorithms for molecular dynamics simulations on commodity clusters. In Proceedings of SC’06: Proceedings of the 2006 ACM/IEEE Conference on Supercomputing, Tampa, FL, USA, 11–17 November 2006; p. 43. [Google Scholar]
  38. Melchionna, S.; Ciccotti, G.; Lee Holian, B. Hoover NPT dynamics for systems varying in shape and size. Mol. Phys. 1993, 78, 533–544. [Google Scholar] [CrossRef]
  39. Martyna, G.J.; Tobias, D.J.; Klein, M.L. Constant pressure molecular dynamics algorithms. J. Chem. Phys. 1994, 101, 4177–4189. [Google Scholar] [CrossRef]
  40. Martyna, G.J.; Klein, M.L.; Tuckerman, M. Nosé–Hoover chains: The canonical ensemble via continuous dynamics. J. Chem. Phys. 1992, 97, 2635–2643. [Google Scholar] [CrossRef]
  41. Toukmaji, A.Y.; Board Jr, J.A. Ewald summation techniques in perspective: A survey. Comput. Phys. Commun. 1996, 95, 73–92. [Google Scholar] [CrossRef]
Figure 1. α-Hemolysin structure. (A). A side view. (B). A top view. (C). The possible ligand binding sites occupied by emodin in a single monomer. Protein is shown as a cartoon and ligand as spheres.
Figure 1. α-Hemolysin structure. (A). A side view. (B). A top view. (C). The possible ligand binding sites occupied by emodin in a single monomer. Protein is shown as a cartoon and ligand as spheres.
Cimb 44 00307 g001
Figure 2. Structure of Compounds 1 and 2.
Figure 2. Structure of Compounds 1 and 2.
Cimb 44 00307 g002
Figure 3. The interaction of emodin with the amino acids in site 1 (A), site 2 (B) and site 3 (C). Hydrogen bonds are shown as arrows and π–π stacking as green lines.
Figure 3. The interaction of emodin with the amino acids in site 1 (A), site 2 (B) and site 3 (C). Hydrogen bonds are shown as arrows and π–π stacking as green lines.
Cimb 44 00307 g003
Figure 4. RMSD of emodin pose in site 1 (A), site 2 (B), and site 3 (C). The left Y axis is for protein RMSD, and the right scale is for ligand RMSD. How stable a ligand is in relation to a protein and its binding pocket is shown by the ligand RMSD (right Y-axis). Red for the protein backbone and blue for the ligand RMSD. MD simulations time in ns (X-axis).
Figure 4. RMSD of emodin pose in site 1 (A), site 2 (B), and site 3 (C). The left Y axis is for protein RMSD, and the right scale is for ligand RMSD. How stable a ligand is in relation to a protein and its binding pocket is shown by the ligand RMSD (right Y-axis). Red for the protein backbone and blue for the ligand RMSD. MD simulations time in ns (X-axis).
Cimb 44 00307 g004
Figure 5. Emodin contacts with amino acids in site 1.
Figure 5. Emodin contacts with amino acids in site 1.
Cimb 44 00307 g005
Figure 6. Average contacts of emodin in site 1.
Figure 6. Average contacts of emodin in site 1.
Cimb 44 00307 g006
Table 1. System details for MD simulation.
Table 1. System details for MD simulation.
Binding SiteSimulation Time (ns)No of AtomsNo of Water AtomsCounter Ion Conc. mM
Site 110046,68314,0223.890
Site 210035,65810,3475.272
Site 310040,92612,1034.507
Table 2. Antibacterial activity of compounds 1 and 2.
Table 2. Antibacterial activity of compounds 1 and 2.
Compound 1
IC50 (μg/mL)
Compound 2
IC50 (μg/mL)
Ciprofloxacin *
IC50 (μg/mL)
Staphylococcus aureus4.2NA0.1
methicillin-resistant S. aureus6.3NA0.1
Test concentrations (µg/mL)20, 4, 0.820, 4, 0.820, 4, 0.8
* Ciprofloxacin was used as positive control; NA = not active up to the maximum dose tested 20 µg/mL. Compounds 1 and 2 revealed no activity against Candida albicans, Cryptococcus neoformans, Aspergillus fumigatus, Escherichia coli, Pseudomonas aeruginosa, and Mycobacterium intracellulare.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ghoneim, M.M.; Al-Serwi, R.H.; El-Sherbiny, M.; El-ghaly, E.-s.M.; Hamad, A.E.; Abdelgawad, M.A.; Ragab, E.A.; Bukhari, S.I.; Bukhari, K.; Elokely, K.; et al. Proposed Mechanism for Emodin as Agent for Methicillin Resistant Staphylococcus Aureus: In Vitro Testing and In Silico Study. Curr. Issues Mol. Biol. 2022, 44, 4490-4499. https://doi.org/10.3390/cimb44100307

AMA Style

Ghoneim MM, Al-Serwi RH, El-Sherbiny M, El-ghaly E-sM, Hamad AE, Abdelgawad MA, Ragab EA, Bukhari SI, Bukhari K, Elokely K, et al. Proposed Mechanism for Emodin as Agent for Methicillin Resistant Staphylococcus Aureus: In Vitro Testing and In Silico Study. Current Issues in Molecular Biology. 2022; 44(10):4490-4499. https://doi.org/10.3390/cimb44100307

Chicago/Turabian Style

Ghoneim, Mohammed M., Rasha Hamed Al-Serwi, Mohamed El-Sherbiny, El-sayed M. El-ghaly, Amal E. Hamad, Mohamed A. Abdelgawad, Ehab A. Ragab, Sarah I. Bukhari, Khulud Bukhari, Khaled Elokely, and et al. 2022. "Proposed Mechanism for Emodin as Agent for Methicillin Resistant Staphylococcus Aureus: In Vitro Testing and In Silico Study" Current Issues in Molecular Biology 44, no. 10: 4490-4499. https://doi.org/10.3390/cimb44100307

APA Style

Ghoneim, M. M., Al-Serwi, R. H., El-Sherbiny, M., El-ghaly, E. -s. M., Hamad, A. E., Abdelgawad, M. A., Ragab, E. A., Bukhari, S. I., Bukhari, K., Elokely, K., & Nael, M. A. (2022). Proposed Mechanism for Emodin as Agent for Methicillin Resistant Staphylococcus Aureus: In Vitro Testing and In Silico Study. Current Issues in Molecular Biology, 44(10), 4490-4499. https://doi.org/10.3390/cimb44100307

Article Metrics

Back to TopTop