- Open Access
Ginger (Zingiber officinale) phytochemicals—gingerenone-A and shogaol inhibit SaHPPK: molecular docking, molecular dynamics simulations and in vitro approaches
Annals of Clinical Microbiology and Antimicrobials volume 17, Article number: 16 (2018)
Antibiotic resistance is a defense mechanism, harbored by pathogens to survive under unfavorable conditions. Among several antibiotic resistant microbial consortium, Staphylococcus aureus is one of the most havoc microorganisms. Staphylococcus aureus encodes a unique enzyme 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (SaHPPK), against which, none of existing antibiotics have been reported.
Computational approaches have been instrumental in designing and discovering new drugs for several diseases. The present study highlights the impact of ginger phytochemicals on Staphylococcus aureus SaHPPK. Herein, we have retrieved eight ginger phytochemicals from published literature and investigated their inhibitory interactions with SaHPPK. To authenticate our work, the investigation proceeds considering the known antibiotics alongside the phytochemicals. Molecular docking was performed employing GOLD and CDOCKER. The compounds with the highest dock score from both the docking programmes were tested for their inhibitory capability in vitro. The binding conformations that were seated within the binding pocket showing strong interactions with the active sites residues rendered by highest dock score were forwarded towards the molecular dynamic (MD) simulation analysis.
Based on molecular dock scores, molecular interaction with catalytic active residues and MD simulations studies, two ginger phytochemicals, gingerenone-A and shogaol have been proposed as candidate inhibitors against Staphylococcus aureus. They have demonstrated higher dock scores than the known antibiotics and have represented interactions with the key residues within the active site. Furthermore, these compounds have rendered considerable inhibitory activity when tested in vitro. Additionally, their superiority was corroborated by stable MD results conducted for 100 ns employing GROMACS package.
Finally, we suggest that gingerenone-A and shogaol may either be potential SaHPPK inhibitors or can be used as fundamental platforms for novel SaHPPK inhibitor development.
Staphylococcus aureus has evolved as one of the most devastating pathogens, demonstrating a wide range of antibiotic resistance . Staphylococcus aureus is a gram positive, non-motile bacterium. This facultative anaerobe is a gram positive, non-motile bacterium hailing from Staphylococcaceae family, powered to infect every known mammalian species causing food poisoning [2, 3]. This is an ectopic commensal and is niched on mucosal membranes and skin of humans . It is transmitted to foods via air, dust, and the lids covering the food containers [5, 6] and the food handlers carry the bacteria on their heads and noses, hence, has an ability to colonize on the normal humans and transmit through direct contact with the bacteria-colonized person. Staphylococcus intoxication occurs due to toxin-contaminated food consumption. Such condition is symptomized very quickly (2–8 h) and is associated with vomiting, abdominal cramps, nausea and/or diarrhea [7, 8]. Even though, staphylococcus intoxication subsides within 48 h, nevertheless, it becomes severe in children and elders  and causes several life threatening infections like, impetigo, ritter disease, osteomyelitis, septic arthritis, endocarditis, toxic shock syndrome, pneumonia, thrombophlebitis and deep skin abscess and infection [10, 11]. Several antibiotics have been used to combat the bacteria [12,13,14,15,16,17,18,19,20] such as, penicillin, methicillin, oxacillin, various vancomycins and glycopeptides, daptomycin, tetracyclines, aminoglycosides, linezolid, chloramphenicol, florfenciol, macrolides, and streptogramins . However, Staphylococcus aureus exerts resistance by several mechanisms that could be broadly categorized into mutations that occur at the chromosomal genes and by horizontally acquired resistance . Specifically, gaining resistance through mutations can happen when the inhibitor is unable to bind to the accurate drug target, derepressing the drug resistance efflux pumps and by mutations that can amend the structure and composition of the drug targets . On the other hand, the horizontally acquired resistance may occur by alteration and inactivation of enzymatic drug, change in the drug binding site, dislocating the drug from its appropriate position and by drug efflux . Adapting either of the mechanisms, the organism endeavours to survive avoiding the encounter with the drug/antibiotic or neutralizes them . Besides these, antibiotic abuse can also add to the raise the resistance . Consequently, the effective treatment is hampered and promotes the infection and enhances the economic burden [23, 24].
Nevertheless, concerns for this bacterium rise due to its resistance against methicillin, often called the Methicillin resistance S. aureus (MRSA) that is prevalent currently by exhibiting diverse phenotypes . This ‘superbug’ was responsible for 19,000 deaths in USA in 1 year  and can be classified as hospital MRSA (haMRSA), referring to those originating from the hospitals and the community MRSA (caMRSA), indicating to those prevalent in the community . Besides these there is another MRSA called as livestock-associated methicillin resistant Staphylococcus aureus . However, caMRSAs are increasingly virulent because of the presence of elevated levels of alpha toxins and phenol-soluble modulin . Not only to methicillin, this bacterium has also been shown resistance against several antibiotics including wonder drug penicillin [1, 22, 30, 31]. Pathogenicity and antibiotic resistance potential of S. aureus drives our interest to develop new drugs that can effectively challenge this bacterium.
Choosing an appropriate target for the discovery of novel antimicrobial drugs is a very important aspect . It will be ideal to identify a target that is confined to pathogen and is not present in host such that the drugs can effectively render its efficacy on the target alone, doing no harm to the host . Accordingly, 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK, EC 18.104.22.168) has been chosen as the drug target for the present investigation  as it is not present in humans . This monomeric protein comprises of 158 residues, which catalyzes the transfer of pyrophosphoryl group from ATP to 6-hydroxymethyl-7, 8-dihydropterin (HMDP), the substrate. It has a molecular weight of 18 kDa pictured by three-layered α-β-α fold. Typically, the HPPK has three loops comprising of loop 1 with residues from 12 to 14, loop 2 consisting of residues from 45 to 51 and loop 3 consisting of residues from 82 to 94, which demonstrate a significant change in the conformation as compared with the protein structure during catalysis. However, loop 3 is known to display major conformational changes . The rigidity of the ternary structure is attributed when the pterin substrate binds and thus loop 3 closes over the binding site. SaHPPK is an ideal target for novel drug designing due to its expression in pathogen only, and is not targeted by existing antibiotics. Alternatively, this enzyme has also been found in other bacteria such as Yersinia pestis, Haemophilus influenzae, Streptococcus pneumonia, Francisella tularensis, Saccharomyces cerevisiae and Escherichia coli, demonstrating conserved active sites [36, 37]. Consequently, it can be suggested that the discovery of drugs against SaHPPK might have inhibitory effects on its homologues on other pathogens .
An overwhelming significance of phytochemicals to be employed as drug molecules is their capability to induce nutritional values besides acting as a medicine and thus, redirecting towards the formulation of nutraceuticals. Nutraceuticals are the food that have both the nutritional and the pharmaceutical value  and refers to as a food or part of the food that induces health benefits while offering medicinal values and thus participates in improving the health of an individual. Additionally, they are cost effective in their formulation and produce no side effects and serious toxicities upon prolonged usage . Due to these properties, the phytochemicals have high advantage to be labelled as drugs [39, 40]. Ginger (Zingiber officinale) is an aromatic, pungent and spicy herb, enriched with the natural phytochemicals. For a time, this herbal species has been used as a flavoring agent and has been placed on the top list of folk medicines against common cold , sore throat  etc. It is therefore understood that ginger is regarded to be safe , however, little is known regarding its mechanism of action and hence, careful assessment is required before considering ginger phytochemicals for any therapy . Moreover, ginger extracts were known to showcase its inhibitory effect on Staphylococcus aureus [44, 45]. All these intriguing factors focus our research to investigate ginger phytochemicals against SaHPPK and elaborate their mechanism of interaction.
In order to accomplish this goal, we have retrieved eight selected ginger phytochemicals from published literature based upon their therapeutic ability and antimicrobial activity [44, 46,47,48,49,50,51,52]. The predictive inhibitory effects of these selected phytochemicals have been evaluated against SaHPPK by molecular docking. To further infer on the mode of binding and interaction with catalytic active residues, successful candidates have been subjected to MD simulations. For the accomplishment of this objective, docking of seven antibiotics, such as streptomycin, ampicillin, amoxicillin, methicillin, penicillin, trimethoprim, and sulfamethoxazole were also considered for comparative analysis. Finally, gingerenone-A and shogaol propounded as potential ginger-phytochemicals against SaHPPK.
Selection of the protein and its preparation
SaHPPK structure used for current investigation, has been obtained from RCSB Protein Data Bank (PDB) with PDB code 3QBC , which is in complex with 2-amino-8-sulfanyl-1,9-dihydro-6H-purin-6-one (8MG). The protein was prepared employing Discovery Studio v4.5 (hereinafter D.S v4.5) by removing all heteroatoms and the addition of hydrogen atoms . The structure was energy minimized, until the convergence gradient satisfied was obtained. The active site was evaluated 10.0 Å around 8MG and the key residues were identified as Ala44, Thr43, Val46, and Asn56 . Interrogating the active site revealed the presence of two Phe residues; Phe54 and Phe123, that are located on either side of the 8MG . Moreover, the histidine residues of protein were oriented in accordance with crystal structure to ND1H protonation state.
Selection and preparation of the ligands
Ginger has several active compounds  and to the best of our knowledge, they have not been tested against SaHPPK as no reports were retrieved upon performing a systematic search. This triggered our interest to understand how these phytochemicals effect the SaHPPK and therefore, for the current study, phytochemicals namely 6-dehydrogingerdion, gingerenone-A, gingerol, paradol, shogaol, zingerone, trans-1,8-cineole-3,6-dihydroxy-3-O-β-d-glucopyranoside and trans-3-hydroxy-1,8-cineole-O-β-d-glucopyranoside were selected that have not been assessed against SaHPPK [54, 55]. More specifically the selection of these phytochemicals was done based upon their therapeutic ability as reported earlier [46,47,48,49,50]. The 2D structures of the selected phytochemicals were represented in (Fig. 1). The corresponding 2D structures of the eight phytochemicals were sketched on ChemSketch (http://www.acdlabs.com/resources/freeware/ChemSketch/) and were subsequently, imported onto D.S v4.5 to generate their 3D structures.
Molecular docking mechanism
Molecular docking is a promising strategy to mimic intermolecular binding modes and interactions. Particularly, molecular docking relies on binding site topology, intermolecular affinity and interaction of key residues with the ligand. For the current study, Genetic Optimization for Ligand Docking (GOLD) v5.2.2 was employed to perform the docking studies. Goldscore was recruited to compute the binding affinities between the protein and ligands, whereas the Chemscore was used for the rescoring purpose. Goldscore comprises of external H-bond, external vdW, internal vdW and internal torsion. Moreover, to obtain an appropriate binding pattern of ligands, 30 docking poses were allowed to generate. Additionally, for identifying the best pose, the Goldscore, interactions between the protein’s active site residues and ligand and the binding modes were considered.
To further validate the obtained results, a second docking programme, CDOCKER implemented on D.S v4.5 has been employed. The results were evaluated based upon CDOCKER interaction energy; higher CDOCKER interaction energy implies greater favourable binding . This is a grid based docking operates by employing CHARMm and facilities the generation of random ligand conformations retrieved form the initial structure.
In-vitro antimicrobial analysis of phytochemicals
In-vitro evaluation of the phytochemicals was performed to infer the results obtained from the computational approach. The phytochemicals and the Staphylococcus aeureus were the generous gifts from the Osmania University, Department of Botany, Hyderabad. To maintain the prepared medium with no contamination, it was autoclaved for 15 min at 15 lb pressure along with petri dishes, spreader, 4–25 ml conical flasks, forceps, inoculation loops and cotton balls. The agar media was then transferred into the petri dishes and was allowed to solidify.
Preparation of the culture media
The nutrient agar media was prepared by suspending 28 g of nutrient agar in 1000 ml distilled water according to Mueller and Hinton , and was maintained at pH 7.0 and at room temperature.
Preparation of the inoculum
20 ml of the above-prepared media was transferred onto the petri-dishes and was allowed to solidify. A loopful of bacteria was transferred to 10 ml of distilled water in a test tube and the addition is continued until the turbidity is equal to standard 0.5 McFarland. Employing the cotton swabs the inoculum was gently swabbed on the surface of the media and were then allowed to dry.
Preparation of the disks
The disks (Whatman 1 filter paper) were prepared with the help of the punch machine of 6 mm in diameter. For the current experiment, the two phytochemicals that have produced the highest dock score along with the antibiotic amoxicillin were considered for current in vitro test. The samples were prepared in the concentration of 1 mg/1 ml.
The bacterial strain, overnight culture was grown in broth was adjusted to an inoculum size of 106 CFU/ml for inoculation of the agar plates. The sterilized disks were carefully impregnated with the sample, allowed to dry for 1 min, and then transferred onto the petri dishes containing 20 ml on nutrient agar. The plates were allowed to incubate for 24 h at 35 ± 2 °C and was followed by measuring the zone of inhibition expressed in mm and was performed in triplicates .
Minimum inhibitory concentration
To quantify the minimum inhibitory concentration (MIC), different concentration of the phytochemicals were tested against the Staphylococcus aureus (MTCCB 737) ranging between 0.05 and 2 mg/ml. The MIC is defined as the lowest concentration of the phytochemical (highest dilution) that can inhibit the growth of the bacteria.
In vitro results were analyzed employing the GraphPad Prism v7.02 and were expressed as mean ± standard deviation recruiting the Turkey’s method and the correlation analysis was executed by two-way ANOVA, P value less than 0.005 was considered as significant.
Molecular dynamics simulations to assess the binding modes of the hits against the reference
To gain further insight into the protein–ligand interactions, the selected ginger phytochemicals were subjected to MD simulations along with reference compound (8MG) and amoxicillin. Parameters for protein’s topology and coordinates were developed by CHARMm27 ff [59,60,61,62] in GROMACS 5.0.7 . The topologies of the ligands and the cofactor were extracted from the SwissParam . The parameters for topology and coordinates of protein and for corresponding ligand were merged, and ten independent systems (one for each phytochemical, one for amoxicillin, and one for 8MG) were designed. Each system was solvated in a dodecahedron box, using TIP3P water model, and neutralized with counter ions. Each solvated system was energy minimized by employing steepest descent algorithm for 10,000 steps and an upper limit of force being lower than 1000 kJ/mol was employed to remove any bad contacts and steric clashes of protein–ligand complexes. Every minimized system was grouped into protein–ligand and solvent-ions to escape collapse, and subsequently, subjected to equilibration. The equilibration of each system was comprised of two components. First, equilibration was conducted at constant volume (NVT) for 1 ns at constant temperature of 300 K using Berendsen thermostat algorithm . Following this, second equilibration was executed for 1 ns at constant pressure (NPT) of 1 bar maintained by Parrinello–Rahman barostat  and LINCS  was employed to constrain all bonds. Particle Mesh Ewald (PME)  was used to calculate long-range electrostatic interactions with a cut-off of 1.2 nm. All short-range non-bonded interactions were calculated within a cut-off of 1.2 nm. A cut-off distance of 12 Å was attributed for Coulombic and van der Waals interactions. All simulations were executed using the NPT ensemble for 100 ns, and coordinates were saved after each 2 fs intervals. The results were examined recruiting visual molecular dynamics  and D.S v4.5.
Molecular docking results showed that the eight ginger phytochemicals have strong interactions with SaHPPK protein, however, gingerenone-A and shogaol displayed the highest Goldscore of 63.62 and 55.48 respectively (Additional file 1: Table S1). On the other hand, it was noted that the antibiotics showed lower Goldscore than the phytochemicals. Among them, amoxicillin has generated the highest dock score of 41.98 (Additional file 1: Table S2), and therefore, this antibiotic was considered for further studies. Hereinafter, the co-crystal, 8MG in SaHPPK crystal structure, was designated as the reference compound. Characteristically, this is a co-crystal located at the active site of the protein and imparts knowledge on the location where the chosen ligands should be anchored at the proteins binding groove. Additionally, this guides the key residues that are involved in the inhibition and marked at 10.0 Å around its location. Moreover, an interaction with these residues labels prospective drug candidates as effective. Furthermore, the reference molecule should logically determine an appropriate binding mode of the candidate molecules. Therefore, the 8MG has been represented as the reference compound. To further evaluate the best conformation, the dock scores rendered by GOLD and CDOCKER, catalytic active residue(s) interactions and the binding modes were opted as the determinant factors. The study was proceeded in the presence of Mg2+-AMPCC.
In-vitro antimicrobial analysis of phytochemicals
The ginger phytochemicals have rendered remarkable results and were comparable with the standard reference antibiotic and a control into which no inoculum was added. The mean zone of inhibition of the phytochemical shogaol has recorded to be between 6 and 12 mm and gingerenone-A was found to be between 2 and 8 mm, respectively. These reading are in harmony with that of the reference antibiotic and were observed to be 4–16 mm. The minimum zone of inhibition was observed at 25 µg/ml for the phytochemicals and the reference antibiotic. The in vitro results further state that the phytochemicals could induce the inhibitory effect in par with amoxicillin (Fig. 2).
Molecular dynamics simulations
To gain insight into the interaction mechanism of ligands (selected phytochemicals, amoxicillin, and reference compound), MD simulations were carried out [70,71,72] for 100 ns and the results were analyzed. Throughout the MD run, dynamic behavior and conformational changes of all the ligands were monitored. The root mean square deviation (RMSD) of backbone atoms was evaluated to estimate the stability of each protein–ligand complex. The reference molecule displayed an average RMSD of ~ 0.18 nm, while gingerenone-A and shogaol projected ~ 0.15 and ~ 0.16 nm respectively (Fig. 3), while amoxicillin has demonstrated a RMSD of ~ 0.17 nm. Additionally, all RMSD plots were noticed to be in 0.1–0.27 nm range (Additional file 1: Figure S1) suggesting their stability during the simulation . Furthermore, the potential energy analysis notified that all the compounds have ranged between − 390,000 and − 395,000 kJ/mol, (Fig. 4). The root mean square fluctuation (RMSF) profiles of backbone atoms of the corresponding systems shed light on their fluctuations. Accordingly, for all the systems, backbone fluctuation was observed within 0.49 nm and was found to be in similar manner. However, moderate deviations were noticed for amoxicillin between 0.1 and 0.4 nm as depicted in blue box with reference to the residues that lie between 90 and 110 (Fig. 5a) and their corresponding atoms (Fig. 5b). This deviation might be because of the non-bonded water molecule that is present in the active site.
The binding mode assessment was executed utilizing the last 20 ns trajectories. Upon superimposition of the representative structure of each system, it was observed that all the ligands have occupied the binding pocket in a similar manner as was observed for reference compound (Fig. 6) and displayed similar binding pattern. The binding pocket is present towards the loop 2 and the cofactor site is located near to loop 3. Since the target protein is devoid of cofactor AMPCC and two Mg2+ ions, we have imported their coordinates from SaHPPK protein with the PDB code 5ETR. Additionally, the key residues that shaped the active site were found to interact with the reference compound as well as for the selected phytochemicals. Inspecting the molecular interactions revealed that the reference compound has generated four hydrogen bonds with key residues of SaHPPK (Fig. 7a). The carboxylic oxygen (hereinafter O) of Ala44 has interacted with N3 of reference compound, whereas, O of Val46 has H-bond interaction with N5 of reference compound. One H-bond was detected between the OG1 of Thr43 residue and N3 of reference compound, while another H-bond was formed between ND2 of Asn56 residue and O1 of reference compound. OD1 atom of Asn56 additionally participated in H-bond interaction with the H13 atom of the ligand. Furthermore, it was observed, that all the H-bonds displayed a distance of ~ 2.8 Å.
Amoxicillin formed one hydrogen bond between N atom of Tyr48 has formed the H-bond with the O2 atom of the ligand with a distance of 2.9 Å, while another H-bond was observed for NE22 of Gln51 and O2 of amoxicillin with a bond distance of 2.9 Å (Fig. 7b).
Hydrogen bond interactions with key residues were noticed with all the phytochemicals, however their number varied significantly (Fig. 7). To further authenticate our results, we have done a comparison with the co-crystal inhibitor and the known antibiotic. The active site and the key residues (Ala44, Thr43, Val46, and Asn56) were defined 10.0 Å around the inhibitor. The prospective drug molecules were examined critically for their interactions with these residues. Amongst them, phytochemicals, shogaol, gingerol, trans-1,8-cineole-3,6-dihydroxy-3-O-β-d-glucopyranoside and trans-1,8-cineole-3,6-dihydroxy-3-O-β-d-glucopyranoside have displayed 2 H-bonds with protein, 6-dehydrogingerdione and zingerone have formed one H-bond each while gingerenone-A displayed 4 H-bonds (Fig. 7).
Further delineating, it was noted that one H-bond was observed between the N of Val46 and O2 of gingerenone-A with a bond distance of 2.3 Å. The NE2 of Gln51 has interacted with the O4 of the ligand with a distance of 2.9 Å. Another H-bond was noticed between the H44 of the gingerenone-A and OD1 of Asn56 with a distance of 2.0 Å, while the other H-bond was detected between NH2 of Arg121 and O3 of gingerenone-A with bond distance of 1.3 Å, (Fig. 7c). Gingerenone-A also formed water-mediated interaction with Asp95 of SaHPPK. On the other hand, gingerol formed two hydrogen bonds with Val46 and Gln51 each. The O atom of the Val46 has interacted with H37 of gingerol with a distance of 1.9 Å. Gln51 has participated in the H-bond formation with its NE2 atom and the O2 of the ligand with a distance of 2.9 Å (Fig. 7d). Moreover, shogaol formed two H-bond with Val46 and Tyr48, respectively. The H41 atom of the ligand has formed a bond with the O of Val46 with a distance of 2.9 Å, while the N atom of Tyr48 has interacted with the O3 of the ligand with a distance of 2.9 Å (Fig. 7e). Additionally, a water-mediated bond with Asp95 stabilized shogaol. Zingerone has also generated an H-bond between N of Val46 and O3 of ligand with a bond distance of 2.9 Å, (Fig. 7f). When 6-dehydrogingerdion was evaluated for H-bond analysis, it was observed that O3 of the inhibitor and N of Val46 formed a single H-bond with a bond length of 2.5 Å, (Fig. 7g). Phytochemical paradol did not render any hydrogen bond; however, it had shown van der Waals interactions Table 1. Phytochemical trans-1,8-cineole-3,6-dihydroxy-3-O-β-d-glucopyranoside demonstrated two H-bonds with Val46 and Gln51. The O atom of Val46 has involved in the hydrogen bond with H48 of the ligand with a distance of ~ 1.8 Å. The second bond was formed between the HE22 of Gln51 and O20 of the ligand represented by a bond distance of 2.1 Å, (Fig. 7h). The compound trans-3-hydroxy-1,8-cineole-O-β-d-glucopyranoside has rendered 2 H-bonds anchored by Val46 and Gln51, respectively. The O atom of Val46 and H47 of the ligand joined by a H-bond displaying a bond length of 1.8 Å. The second H-bond was formed between the HE22 atom of Gln51 and O19 of the ligand with a distance of 2.1 Å, (Fig. 7i).
Focusing on the importance of water molecule in SaHPPK stability and augmenting its reactivity, it was speculated that water molecule at active site plays a crucial rule (Fig. 7). The water mediated bond with Asp95 was noticed in the presence of all the phytochemicals as was observed with the reference molecule, however this interaction was absent with the antibiotic amoxicillin, (Fig. 7b). Further details of the interactions are recorded in Table 1. The resultant docked poses were validated by MD simulation analysis and it was confirmed that the binding stability of the selected poses remained unaltered during the simulation.
The interactions of the cofactor and the Mg2+ with the protein were additionally evaluated that the benzene ring of the adenine group has interacted with three hydrogen bonds formed by Ile98 and Ser112 residues. Leu71 additionally holds the adenine group by the hydrophobic interactions. The ribose moiety has interacted with Lys110 demonstrated by a hydrogen bond. Arg121 has interacted with O1G of the cofactor by electrostatic bond. Furthermore, the electrostatic bonds hold the Mg2+ ions represented by Glu78 and Asp97. Additionally, the exposed O atoms of the cofactor firmly hold the Mg2+ ions (Additional file 1: Figure S2).
Despite tremendous progress in medical sciences, the effective treatment against infectious microorganisms remains a major challenge. The primary reason behind this failure is the ability of the microorganisms to gain resistance against antibiotics, which is conferred by a variety of mechanisms. Consequently, it is essential to develop new drugs that can effectively combat the microorganisms and to overcome their pathogenicity. Staphylococcus aureus is one of the widely known pathogenic microorganisms that has gained resistance against several antibiotics. In this study, ginger phytochemicals were employed to evaluate their inhibitory effects when challenged against microbial pathogenicity. Since, HPPK plays a key role in microbial folate pathway and hence, SaHPPK might be an ideal target for novel inhibitors. Additionally, an ideal target should possess the following attributes, such as having no homolog in humans, should exist in large range of bacteria performing characteristics role, should be specifically druggable and should possess a low cross-resistance potential (https://www.ncbi.nlm.nih.gov/books/NBK200811/#sec_17). Because SaHPPK represents all these characteristic features, we have relied on SaHPPK for our current investigation.
Besides, the selection of the phytochemicals have been performed based upon the literature search taking into consideration that the phytochemicals portray and are embedded with therapeutic activities [46,47,48,49,50]. We aimed at understanding how these compounds that have a similar structure act when challenged against SaHPPK and further which phytochemical is potential against the targeted protein. Such a study was conducted earlier considering different phytochemicals against different targets and diseases [74,75,76,77]. The study was conducted in the presence of Mg2+-AMPCC to determine the effect of the chosen phytochemicals.
Our investigation of eight ginger phytochemicals against SaHPPK demonstrated that two candidate phytochemicals gingerenone-A and shogaol have higher inhibitory effects. Furthermore, highest dock score, stable orientation of selected phytochemicals in SaHPPK active site and stable interactions with key residues, support our investigation. Furthermore, to validate the dock results, we have performed the in silico investigation employing second docking tool, CDOCKER. These docking results reaffirm the superiority of gingerenone-A and shogaol rendered by highest CDOCKER interaction energy (Additional file 1: Table S1) and the interaction with key residues located at the active site of the protein. Additionally, different interactions also affirmed that gingerenone-A and shogaol might be potential scaffolds to be developed has novel SaHPPK inhibitors.
In the current study, all the phytochemical have displayed a score greater than the reference, however, amongst them, since, gingerenone-A and shogaol projected remarkable scores by both the molecular docking programmes and further displayed stable MD results, we therefore have considered only the top two dock scored compounds for further evaluation and subsequently for the in vitro experiments. Delineating on in vitro results, it can be observed that the lead candidates have exhibited remarkable inhibition at all the concentrations including at 25 µg/ml. However shogaol has demonstrated an overall greater inhibition at 50 µg/ml while gingerenone-A rendered a marginally higher degree of inhibition at 75 µg/ml. Nevertheless, both the phytochemicals have conferred with the inhibitory activities in par with the antibiotic thus; affirm the inhibitory potential of ginger phytochemicals against SaHPPK.
Since, it is widely accepted that S. aureus is resistant to existing antibiotics; we speculate that none of these antibiotics have potential to inhibit SaHPPK. One of the possible reasons for this failure is the missing of H-bonding of amoxicillin with Val46. Valine at position 46 of SaHPPK has been reported as key catalytic residue and forms stable H-bond with reference compound . Our findings showed that the phytochemicals as well as the reference compound formed stable H-bond with Val46 after 100 ns MD simulation, while such a binding pattern was not observed for amoxicillin. Delineation on Val46 and its significance, we scrupulously monitored the interaction with their respective atoms of the ligands throughout the MD run. Val46 was seen to anchor with the ligands within the acceptable hydrogen bond length of < 3 Å. Conversely, amoxicillin failed to represent the substantial interaction. This comparison and MD examination led us to the conclusion that Val46 is crucial in any SaHPPK targeted small molecule therapy. Apart from that, water molecule plays very important role in increasing the accuracy and feasibility of chemical reactions [78, 79]. The co-crystal structure of SaHPPK also proclaimed a crucial role of the water molecule . SaHPPK crystal structure revealed that the reference compound is stabilized in enzyme’s active site by a water-mediated bond present at Asp95 of the enzyme. Our results also showed the same pattern of water-mediated stability for all phytochemicals as well as the reference compound (Fig. 7). Conversely, amoxicillin could not gain water-mediated stability in SaHPPK active site. Based on our results we emphasize that the interaction of Val46 residue and water-mediated interaction of Asp95 are the pre-requisite for SaHPPK functional exploration and/or inhibition.
The crystal structure of protein SaHPPK PBD code: 3QBC, has three loops loop 1, loop 2, loop 3; comprising of 12–14, 45–51 and 82–94 residues, respectively . As in the crystal structure, the co-crystal substrate 8MG has been sandwiched between aromatic residues Phe54 and Phe123. Our results also confirmed similar π–π stacked interactions between the ligand molecules and the Phe54 and Phe123 of the crystal structure.
We further investigated the structural topology and the behaviour of the three loops upon the interaction with the phytochemicals and the evaluation was based upon the findings of Kaifu et al. . Loop 1 showed the semi-open conformation, while loop 2 and loop 3 demonstrated the closed conformations as represented in (Fig. 8 and Table 2). From the results, it was evident that loop 2 predominantly displayed closed conformation and loop 3 on the other hand showed closed conformation except for amoxicillin. These results further impart information that the presence of AMPCC in the cofactor site has induced the closed conformation of loop 3. We further speculated that Asp95 that lies in close proximity with loop 3 might also have played a role in providing closed conformation while it shows semi-closed conformation with amoxicillin, in which case the Asp95 binding was absent (Fig. 8).
We further assessed to comprehend the binding ability of the phytochemicals across different homologues. It is well reported that the HPPK structures are highly conserved with six-strands in α-β-α fold with Mg2+ ions that recognize the ATP and HMDP substrates. Additionally, among the present HPPK proteins, the S. aureus homologue shares the identities with E. coli and Y. pestis by 39%, H. influenza and S. cerevisiae by 37 and 34% with S. pneumonia [36, 81]. It can be speculated that, since the organisms share a conserved active site and high structural similarity of the ternary complexes, the identified phytochemicals could also exert their antimicrobial effect on its homologues as was reported earlier .
Based upon the results obtained from molecular docking, MD simulation and in vitro studies it can be indicated that gingerenone-A and shogaol might be more effective against SaHPPK. Consequently, we suggest these two ginger-phytochemicals as foundation scaffolds for the development of novel SaHPPK inhibitors.
The best way to increase the antibiotic activity especially for the multidrug resistant bacteria is the inclusion of natural products into the drug formulation, as they offer a plethora of advantages. In the present investigation, ginger phytochemicals were evaluated for the prospective drugs. Out of the eight phytochemicals, shogaol, and gingerenone-A were potential drug candidates demonstrating highest dock scores and strong active site residue interactions. Furthermore, the MD simulation results confirmed their stable orientation and strong interactions with catalytic active residues of SaHPPK catalytic pocket. We therefore speculate that the two phytochemicals can be effective against SaHPPK.
Chambers HF, Deleo FR. Waves of resistance: Staphylococcus aureus in the antibiotic era. Nat Rev Microbiol. 2009;7:629–41.
Le Loir Y, Baron F, Gautier M. Staphylococcus aureus and food poisoning. Genet Mol Res. 2003;2:63–76.
Johler S, Giannini P, Jermini M, Hummerjohann J, Baumgartner A, Stephan R. Further evidence for staphylococcal food poisoning outbreaks caused by egc-encoded enterotoxins. Toxins (Basel). 2015;7:997–1004.
Otto M. Staphylococcus epidermidis—the “accidental” pathogen. Nat Rev Microbiol. 2009;7:555–67.
Weese JS. Methicillin-resistant Staphylococcus aureus in animals. ILAR J England. 2010;51:233–44.
Lambrechts AA, Human IS, Doughari JH, Lues JFR. Bacterial contamination of the hands of food handlers as indicator of hand washing efficacy in some convenient food industries. Pak J Med Sci. 2014;30:755.
Argudín MÁ, Mendoza MC, Rodicio MR. Food poisoning and Staphylococcus aureus enterotoxins. Toxins (Basel). 2010;2:1751–73.
Thakkar S, Agrawal R. A case of Staphylococcus aureus enterocolitis: a rare entity. Gastroenterol Hepatol. 2010;6:115–7.
Kadariya J, Smith TC, Thapaliya D. Staphylococcus aureus and staphylococcal food-borne disease: an ongoing challenge in public health. Biomed Res Int. 2014;2014:827965.
Hartman-Adams H, Banvard C, Juckett G. Impetigo: diagnosis and treatment. Am Fam Physician. 2014;90:229–35.
McCaig LF, McDonald LC, Mandal S, Jernigan DB. Staphylococcus aureus-associated skin and soft tissue infections in ambulatory care. Emerg Infect Dis. 2006;12:1715–23.
Rayner C, Rayner C, Munckhof WJ, Munckhof WJ. Antibiotics currently used in the treatment of infections caused by Staphylococcus aureus. Intern Med J. 2005;35(Suppl 2):S3–16.
Sandberg A, Jensen KS, Baudoux P, Van Bambeke F, Tulkens PM, Frimodt-Møller N. Intra- and extracellular activities of dicloxacillin against Staphylococcus aureus in vivo and in vitro. Antimicrob Agents Chemother. 2010;54:2391–400.
Garzoni C, Uçkay I, Belaieff W, Breilh D, Suvà D, Huggler E, et al. In vivo interactions of continuous flucloxacillin infusion and high-dose oral rifampicin in the serum of 15 patients with bone and soft tissue infections due to Staphylococcus aureus—a methodological and pilot study. Springerplus. 2014;3:287.
Rubin JE, Ball KR, Chirino-Trejo M. Antimicrobial susceptibility of Staphylococcus aureus and Staphylococcus pseudintermedius isolated from various animals. Can Vet J. 2011;52:153–7.
Soares GMS, Figueiredo LC, Faveri M, Cortelli SC, Duarte PM, Feres M. Mechanisms of action of systemic antibiotics used in periodontal treatment and mechanisms of bacterial resistance to these drugs. J Appl Oral Sci. 2012;20:295–304.
Sakoulas G, Olson J, Yim J, Singh NB, Kumaraswamy M, Quach DT, et al. Cefazolin and ertapenem, a synergistic combination used to clear persistent Staphylococcus aureus bacteremia. Antimicrob Agents Chemother. 2016;60:6609–18.
Hu Q, Peng H, Rao X. Molecular events for promotion of vancomycin resistance in vancomycin intermediate Staphylococcus aureus. Front Microbiol. 2016;7:1601.
Matsumoto K, Watanabe E, Kanazawa N, Fukamizu T, Shigemi A, Yokoyama Y, et al. Pharmacokinetic/pharmacodynamic analysis of teicoplanin in patients with MRSA infections. Clin Pharmacol. 2016;8:15–8.
Delgado A, Zaman S, Muthaiyan A, Nagarajan V, Elasri MO, Wilkinson BJ, et al. The fusidic acid stimulon of Staphylococcus aureus. J Antimicrob Chemother. 2008;62:1207–14.
Foster TJ. Antibiotic resistance in Staphylococcus aureus. Current status and future prospects. FEMS Microbiol Rev. 2017;41:430–49.
Lowy FD. Antimicrobial resistance: the example of Staphylococcus aureus. J Clin Investig. 2003;111:1265–73.
McAdam AJ, Hooper DC, DeMaria A, Limbago BM, O’Brien TF, McCaughey B. Antibiotic resistance: how serious is the problem, and what can be done? Clin Chem. 2012;58:1182–6.
Ventola CL. The antibiotic resistance crisis: part 1: causes and threats. P T A Peer-Reviewed J Formul Manag. 2015;40:277–83.
Deplano A, Vandendriessche S, Nonhoff C, Denis O. Genetic diversity among methicillin-resistant Staphylococcus aureus isolates carrying the mecC gene in Belgium. J Antimicrob Chemother. 2014;69:1457–60.
Chhabra S, Dolezal O, Collins BM, Newman J, Simpson JS, Macreadie IG, et al. Structure of S. aureus HPPK and the discovery of a new substrate site inhibitor. PLoS ONE. 2012;7:e29444.
Huang H, Flynn NM, King JH, Monchaud C, Morita M, Cohen SH. Comparisons of community-associated methicillin-resistant Staphylococcus aureus (MRSA) and hospital-associated MSRA infections in Sacramento, California. J Clin Microbiol. 2006;44:2423–7.
Sharma M, Nunez-Garcia J, Kearns AM, Doumith M, Butaye PR, Angeles Argudín M, et al. Livestock-associated methicillin resistant Staphylococcus aureus (LA-MRSA) clonal complex (CC) 398 isolated from UK animals belong to European lineages. Front Microbiol. 2016;7:1741.
David MZ, Daum RS. Community-associated methicillin-resistant Staphylococcus aureus: epidemiology and clinical consequences of an emerging epidemic. Clin Microbiol Rev. 2010;23:616–87.
Pantosti A, Sanchini A, Monaco M. Mechanisms of antibiotic resistance in Staphylococcus aureus. Future Microbiol. 2007;2:323–34.
Hiramatsu K, Katayama Y, Matsuo M, Sasaki T, Morimoto Y, Sekiguchi A, et al. Multi-drug-resistant Staphylococcus aureus and future chemotherapy. J Infect Chemother. 2014;20:593–601.
Pasdaran A, Hamedi A, Mamedov N. Antibacterial and insecticidal activity of volatile compounds of three algae species of Oman Sea. Int J Sec Metab. 2016;3:66–73.
Haag NL, Velk KK, Wu C. Potential antibacterial targets in bacterial central metabolism. Int J Adv Life Sci. 2012;4:21–32.
Shi G, Shaw G, Liang YH, Subburaman P, Li Y, Wu Y, et al. Bisubstrate analogue inhibitors of 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase: new design with improved properties. Bioorg Med Chem. 2012;20:47–57.
Xiao B, Shi G, Gao J, Blaszczyk J, Liu Q, Ji X, et al. Unusual conformational changes in 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase as revealed by x-ray crystallography and NMR. J Biol Chem. 2001;276:40274–81.
Chhabra S, Newman J, Peat TS, Fernley RT, Caine J, Simpson JS, et al. Crystallization and preliminary X-ray analysis of 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase from Staphylococcus aureus. Acta Crystallogr Sect F Struct Biol Cryst Commun. 2010;66:575–8.
Dennis ML, Pitcher NP, Lee MD, DeBono AJ, Wang Z-C, Harjani JR, et al. Structural basis for the selective binding of inhibitors to 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase from Staphylococcus aureus and Escherichia coli. J Med Chem. 2016;59:5248–63.
Nasri H, Baradaran A, Shirzad H, Kopaei MR. New concepts in nutraceuticals as alternative for pharmaceuticals. Int J Prev Med. 2014;5:1487–99.
Kaul PN, Joshi BS. Alternative medicine: herbal drugs and their critical appraisal–part II. Prog Drug Res Switzerland. 2001;57:1–75.
Molinari G. Natural products in drug discovery: present status and perspectives. Adv Exp Med Biol. 2009;655:13–27.
Raal A, Volmer D, Sõukand R, Hratkevitš S, Kalle R. Complementary treatment of the common cold and flu with medicinal plants—results from two samples of pharmacy customers in Estonia. PLoS ONE. 2013;8:e58642.
Khayat S, Kheirkhah M, Behboodi Moghadam Z, Fanaei H, Kasaeian A, Javadimehr M. Effect of treatment with ginger on the severity of premenstrual syndrome symptoms. ISRN Obstet Gynecol. 2014;2014:792708.
Wilkinson JM. Effect of ginger tea on the fetal development of Sprague-Dawley rats. Reprod Toxicol. 2000;14:507–12.
Gull I, Saeed M, Shaukat H, Aslam SM, Samra Z, Athar AM. Inhibitory effect of Allium sativum and Zingiber officinale extracts on clinically important drug resistant pathogenic bacteria. Ann Clin Microbiol Antimicrob. 2012;11:8.
Sharma PK, Singh V, Ali M. Chemical composition and antimicrobial activity of fresh rhizome essential oil of Zingiber officinale Roscoe. Pharmacogn J. 2016;8:185–90.
Huang S-H, Lee C-H, Wang H-M, Chang Y-W, Lin C-Y, Chen C-Y, et al. 6-Dehydrogingerdione restrains lipopolysaccharide-induced inflammatory responses in RAW 264.7 macrophages. J Agric Food Chem. 2014;62:9171–9.
Antony P, Vijayan R. Identification of novel aldose reductase inhibitors from spices: a molecular docking and simulation study. PLoS ONE. 2015;10:e0138186.
Semwal RB, Semwal DK, Combrinck S, Viljoen AM. Gingerols and shogaols: important nutraceutical principles from ginger. Phytochemistry. 2015;117:554–68.
Wei C-K, Tsai Y-H, Korinek M, Hung P-H, El-Shazly M, Cheng Y-B, et al. 6-Paradol and 6-shogaol, the pungent compounds of ginger, promote glucose utilization in adipocytes and myotubes, and 6-paradol reduces blood glucose in high-fat diet-fed mice. Int J Mol Sci. 2017;18:168.
Ahmad B, Rehman MU, Amin I, Arif A, Rasool S, Bhat SA, et al. A review on pharmacological properties of zingerone (4-(4-hydroxy-3-methoxyphenyl)-2-butanone). Sci World J. 2015;2015:1–6.
Guo T, Tan S-B, Wang Y, Chang J. Two new monoterpenoid glycosides from the fresh rhizome of Tongling White Ginger (Zingiber officinale). Nat Prod Res. 2018;32:71–6.
Chakotiya AS, Tanwar A, Narula A, Sharma RK. Zingiber officinale: its antibacterial activity on Pseudomonas aeruginosa and mode of action evaluated by flow cytometry. Microb Pathog. 2017;107:254–60.
Fu W, Chen L, Wang Z, Zhao C, Chen G, Liu X, et al. Determination of the binding mode for anti-inflammatory natural product xanthohumol with myeloid differentiation protein 2. Drug Des Dev Ther. 2016;10:455–63.
Rahmani AH, Shabrmi FM, Aly SM. Active ingredients of ginger as potential candidates in the prevention and treatment of diseases via modulation of biological activities. Int J Physiol Pathophysiol Pharmacol. 2014;6:125–36.
Ghasemzadeh A, Jaafar HZE, Rahmat A. Changes in antioxidant and antibacterial activities as well as phytochemical constituents associated with ginger storage and polyphenol oxidase activity. BMC Complement Altern Med. 2016;16:382.
Rampogu S, Rampogu Lemuel M. Network based approach in the establishment of the relationship between type 2 diabetes mellitus and its complications at the molecular level coupled with molecular docking mechanism. Biomed Res Int. 2016;2016:6068437.
Mueller JH, Hinton J. A protein-free medium for primary isolation of the gonococcus and meningococcus. Exp Biol Med. 1941;48:330–3.
Meriga B, Mopuri R, MuraliKrishna T. Insecticidal, antimicrobial and antioxidant activities of bulb extracts of Allium sativum. Asian Pac J Trop Med. 2012;5:391–5.
Mackerell AD. Empirical force fields for biological macromolecules: overview and issues. J Comput Chem. 2004;25:1584–604.
MacKerell AD, Bashford D, Bellott M, Dunbrack RL, Evanseck JD, Field MJ, et al. All-atom empirical potential for molecular modeling and dynamics studies of proteins. J Phys Chem B. 1998;102:3586–616.
Zhu X, Lopes PEM, Mackerell AD. Recent developments and applications of the CHARMM force fields. Wiley Interdiscip Rev Comput Mol Sci. 2012;2:167–85.
Mallajosyula SS, Jo S, Im W, MacKerell AD. Molecular dynamics simulations of glycoproteins using CHARMM. Methods Mol Biol. 2015;1273:407–29.
Van Der Spoel D, Lindahl E, Hess B, Groenhof G, Mark AE, Berendsen HJC. GROMACS: fast, flexible, and free. J Comput Chem. 2005;26:1701–18.
Zoete V, Cuendet MA, Grosdidier A, Michielin O. SwissParam: a fast force field generation tool for small organic molecules. J Comput Chem. 2011;32:2359–68.
Berendsen HJC, Postma JPM, van Gunsteren WF, DiNola A, Haak JR. Molecular dynamics with coupling to an external bath. J Chem Phys. 1984;81:3684–90.
Parrinello M. Polymorphic transitions in single crystals: a new molecular dynamics method. J Appl Phys. 1981;52:7182.
Hess B, Bekker H, Berendsen HJC, Fraaije JGEM. LINCS: a linear constraint solver for molecular simulations. J Comput Chem. 1997;18:1463–72.
Darden T, York D, Pedersen L. Particle mesh Ewald: an N·log(N) method for Ewald sums in large systems. J Chem Phys. 1993;98:10089.
Humphrey W, Dalke A, Schulten KVMD. Visual molecular dynamics. J Mol Graph. 1996;14:33–8.
Rampogu S, Baek A, Son M, Zeb A, Park C, Kumar R, et al. Computational exploration for lead compounds that can reverse the nuclear morphology in Progeria. Biomed Res Int. 2017;2017:1–15.
Rampogu S, Son M, Park C, Kim H-H, Suh J-K, Lee K. Sulfonanilide derivatives in identifying novel aromatase inhibitors by applying docking, virtual screening, and MD simulations studies. Biomed Res Int. 2017;2017:1–17.
Kumar R, Bavi R, Jo MG, Arulalapperumal V, Baek A, Rampogu S, et al. New compounds identified through in silico approaches reduce the α-synuclein expression by inhibiting prolyl oligopeptidase in vitro. Sci Rep. 2017;7:10827.
Verma S, Grover S, Tyagi C, Goyal S, Jamal S, Singh A, et al. Hydrophobic interactions are a key to MDM2 inhibition by polyphenols as revealed by molecular dynamics simulations and MM/PBSA free energy calculations. PLoS ONE. 2016;11:e0149014.
Seyedi SS, Shukri M, Hassandarvish P, Oo A, Shankar EM, Abubakar S, et al. Computational approach towards exploring potential anti-Chikungunya activity of selected flavonoids. Sci Rep. 2016;6:24027.
Kong D, Zhang Y, Yamori T, Duan H, Jin M. Inhibitory activity of flavonoids against class I phosphatidylinositol 3-kinase isoforms. Molecules. 2011;16:5159–67.
Antony P, Vijayan R. Acetogenins from Annona muricata as potential inhibitors of antiapoptotic proteins: a molecular modeling study. Drug Des Dev Ther. 2016;10:1399–410.
Middleton E, Kandaswami C, Theoharides TC. The effects of plant flavonoids on mammalian cells: implications for inflammation, heart disease, and cancer. Pharmacol Rev. 2000;52:673–751.
Roberts BC, Mancera RL. Ligand-protein docking with water molecules. J Chem Inf Model. 2008;48:397–408.
Thilagavathi R, Mancera RL. Ligand-protein cross-docking with water molecules. J Chem Inf Model. 2010;50:415–21.
Gao K, Jia Y, Yang M. A network of conformational transitions revealed by molecular dynamics simulations of the binary complex of Escherichia coli 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase with MgATP. Biochemistry. 2016;55:6931–9.
Derrick JP. The structure and mechanism of 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase. Vitam Horm. 2008;79:411–33.
SR and KWL conceived the idea of the project. SR, AB have conducted the computational works. SR, AB and AZ wrote the manuscript. RGG performed the in vitro experiments, RSB and RK performed the MD simulations of phytochemicals. SR, YSK, YJK and KWL approved the analysis of the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
All the data is available with the manuscript and online.
Consent for publication
Ethics approval and consent to participate
This research was supported by Pioneer Research Center Program through the National Research Foundation (NRF) funded by the Ministry of Science, ICT and Future Planning (NRF-2015M3C1A3023028). Next-Generation BioGreen 21 Program (PJ01106202) from Rural Development Administration (RDA) of Republic of Korea also supported this work. This material is based upon work supported by the Ministry of Trade, Industry & Energy (MOTIE, Korea) under Industrial Technology Innovation Program (No. 10038744), ‘Establishment of Drug Development System Using the Drug Repositioning Technology’.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.