Ginger (Zingiber officinale) phytochemicals-gingerenone-A and shogaol inhibit SaHPPK: molecular docking, molecular dynamics simulations and in vitro approaches.

BACKGROUND: 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. METHODS: 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. RESULTS: 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. CONCLUSIONS: 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.

Comparative Structural Dynamics of tRNA(Phe) with Respect to Hinge Region Methylated Guanosine: A Computational Approach.

Transfer RNAs (tRNAs) contain various uniquely modified nucleosides thought to be useful for maintaining the structural stability of tRNAs. However, their significance for upholding the tRNA structure has not been investigated in detail at the atomic level. In this study, molecular dynamic simulations have been performed to assess the effects of methylated nucleic acid bases, N (2)-methylguanosine (m(2)G) and N (2)-N (2)-dimethylguanosine (m 2 (2) G) at position 26, i.e., the hinge region of E. coli tRNA(Phe) on its structure and dynamics. The results revealed that tRNA(Phe) having unmodified guanosine in the hinge region (G26) shows structural rearrangement in the core of the molecule, resulting in lack of base stacking interactions, U-turn feature of the anticodon loop, and TΨC loop. We show that in the presence of the unmodified guanosine, the overall fold of tRNA(Phe) is essentially not the same as that of m(2)G26 and m 2 (2) G26 containing tRNA(Phe). This structural rearrangement arises due to intrinsic factors associated with the weak hydrogen-bonding patterns observed in the base triples of the tRNA(Phe) molecule. The m(2)G26 and m 2 (2) G26 containing tRNA(Phe) retain proper three-dimensional fold through tertiary interactions. Single-point energy and molecular electrostatics potential calculation studies confirmed the structural significance of tRNAs containing m(2)G26 and m 2 (2) G26 compared to tRNA with normal G26, showing that the mono-methylated (m(2)G26) and dimethylated (m 2 (2) G26) modifications are required to provide structural stability not only in the hinge region but also in the other parts of tRNA(Phe). Thus, the present study allows us to better understand the effects of modified nucleosides and ionic environment on tRNA folding.

Conformational Preferences of Modified Nucleoside 5-Taurinomethyluridine, τm(5)U Occur at 'wobble' 34th Position in the Anticodon Loop of tRNA.

Conformational preferences of hypermodified nucleoside 5-taurinomethyluridine 5'-monophoshate 'p-τm(5)U' (-CH2-NH2(+)-CH2-CH2-SO3(-)) have been investigated using semi-empirical RM1 method. Automated geometry optimization using ab initio molecular orbital HF-SCF (6-31G**) and DFT (B3LYP/6-31G**) calculations have also been made to compare the salient features. The RM1 preferred most stable conformation of 'p-τm(5)U' has been stabilized by hydrogen bonding interactions between O(11a)…HN(8), O1P(34)…HN(8), and O1P(34)…HC(10). Another conformational study of 5-taurinomethyluridine side chain has also been performed in context of anticodon loop bases of E. coli tRNA(Leu). The atom O(11a) of τm(5)U(34) side chain interacts with adenosine (A35) as well as ribose-phosphate backbone which might provide structural stability to the anticodon loop. The glycosyl torsion angle of τm(5)U retains 'anti'-conformation. The solvent accessible surface area calculations revealed the role of τm(5)U in tRNA(Leu) anticodon loop. MD simulation results are found in agreement with RM1 preferred stable structure. The MEPs calculations of τm(5)U(34):G3 model show unique potential tunnels between the hydrogen bond donor and acceptor atoms as compared to τm(5)U(34):A3 model. Thus, these results could pave the way to understand the role of τm(5)U(34) to recognize UUG/UUA codons at atomic level in the mitochondrial disease, MELAS.

MD simulation studies to investigate iso-energetic conformational behaviour of modified nucleosides m(2)G and m(2) 2G present in tRNA.

Modified nucleic acid bases are most commonly found in tRNA. These may contain modifications from simple methylation to addition of bulky groups. Methylation of the four canonical nucleotide bases at a wide variety of positions is particularly prominent among the known modification. Methylation of N2 group of guanine is a relatively common modification in tRNA and rRNA. N2-methylguanosine (m(2)G) is the second most often encountered nucleoside in E. coli tRNAs. N2, N2- dimethylguanosine (m(2) 2G) is found in the majority of eukaryotic tRNAs and involved in forming base pair interactions with adjacent bases. Hence, in order to understand the structural significance of these methylated nucleic acid bases we have carried out molecular dynamics simulation to see the salvation effect. The results obtained shows iso-energetic conformational behaviors for m(2)G and m(2) 2G. The simulation trajectory of m(2)G shows regular periodical fluctuations suggesting that m(2)G is equally stable as either s-cis or s-trans rotamers. The two rotamers of m(2)G may interact canonically or non-canonically with opposite base as s-trans m(2)G26:C/A/U44 and s-cis m(2)G26:A/U44. The free rotations around the C-N bond could be the possible reason for these iso-energetic conformations. Dimethylation of G has almost no influence on base pairing with either A or U. Thus, these results reveal that modified nucleosides m(2)G and m(2) 2G may play an important role to prevent tRNA from adopting the unusual mitochondrial like conformation.

Conformational preferences of modified nucleoside N(2)-methylguanosine (m(2)G) and its derivative N(2), N(2)-dimethylguanosine (m(2)(2)G) occur at 26th position (hinge region) in tRNA.

Conformational preferences of the modified nucleosides N(2)-methylguanosine (m(2)G) and N(2), N(2)-dimethylguanosine (m(2)(2)G) have been studied theoretically by using quantum chemical perturbative configuration interaction with localized orbitals (PCILO) method. Automated complete geometry optimization using semiempirical quantum chemical RM1, along with ab initio molecular orbital Hartree-Fock (HF-SCF), and density functional theory (DFT) calculations has also been made to compare the salient features. Single-point energy calculation studies have been made on various models of m(2)G26:C/A/U44 and m(2)(2)G26:C/A/U44. The glycosyl torsion angle prefers "syn" (χ = 286°) conformation for m(2)G and m(2)(2)G molecules. These conformations are stabilized by N(3)-HC2' and N(3)-HC3' by replacing weak interaction between O5'-HC(8). The N(2)-methyl substituent of (m(2)G26) prefers "proximal" or s-trans conformation. It may also prefer "distal" or s-cis conformation that allows base pairing with A/U44 instead of C at the hinge region. Thus, N(2)-methyl group of m(2)G may have energetically two stable s-trans m(2)G:C/A/U or s-cis m(2)G:A/U rotamers. This could be because of free rotations around C-N bond. Similarly, N(2), N(2)-dimethyl substituent of (m(2)(2)G) prefers "distal" conformation that may allow base pairing with A/U instead of C at 44th position. Such orientations of m(2)G and m(2)(2)G could play an important role in base-stacking interactions at the hinge region of tRNA during protein biosynthesis process.