ABSTRACT
Tunicamycin (TUN) is a nucleoside antibiotic with a complex structure comprising uracil, tunicamine sugar, N-acetylglucosamine (GlcNAc), and fatty acyl tail moieties. TUN, known as a canonical inhibitor, blocks vital functions of certain transmembrane protein families, for example, the insect enzyme dolichyl phosphate α-N-acetylglucosaminylphosphotransferase (DPAGT1) of Spodoptera frugiperda and the bacterial enzyme phospho-N-acetylmuramoylpentapeptide translocase (MraYCB) of Clostridium bolteae. Accurate description of protein-drug interactions has an immense impact on structure-based drug design, while the main challenge is to create proper topology and parameter entries for TUN in modeling protein-TUN interactions given the structural complexity. Starting from DPAGT1-TUN and MraYCB-TUN crystal structures, we first sketched these structural complexes on the basis of the CHARMM36 force field and optimized each of them using quantum mechanics/molecular mechanics (QM/MM) calculations. By continuing calculations on the active site (QM region) of each optimized structure, we specified the characteristics of intermolecular interactions contributing to the binding of TUN to each active site by quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) analyses at the M06-2X/6-31G** level. The results outlined that TUN insertion into each active site requires multiple weak, moderate, and strong hydrogen bonds accompanying charge-dipole, dipole-dipole, and hydrophobic interactions among different TUN moieties and adjacent residues. The water-mediated interactions also play central roles in situating the uracil and tunicamine moieties of TUN within the DPAGT1 active site as well as in preserving the uracil-binding pocket in the MraYCB active site. The TUN binds more strongly to DPAGT1 than to MraYCB. The information garnered here is valuable particularly for better understanding mode of action at the molecular level, as it is conducive to developing next generations of nucleoside antibiotics.
ABSTRACT
Phospho-N-acetylmuramoyl-pentapeptide translocase (MraYAA) from Aquifex aeolicus is the binding target for the nucleotide antibiotic muraymycin D2 (MD2). MraYAA in the presence of the MD2 ligand has been crystallized and released, while the interactions between the ligand and active-site residues remain less quantitatively and qualitatively defined. We characterized theoretically the key residues involved in noncovalent interactions with MD2 in the MraYAA active site. We applied the quantum theory of atoms in molecules and natural bond orbital analyses based on the density functional theory method on the solved crystal structure of MraY with the MD2 to quantitatively estimate the intermolecular interactions. The obtained results revealed the presence of multiple hydrogen bonds in the investigated active site with strength ranging from van der Waals to covalent limits. Lys70, Asp193, Gly194, Asp196, Gly264, Ala321, Gln305, and His325 are key active-site residues interacting with MD2. Conventional and unconventional hydrogen bonds in addition with charge-dipole and dipole-dipole interactions contribute significantly to stabilize the MD2 binding to the MraYAA active site. It was also found that water molecules inside the active site have substantial effects on its structure stability through hydrogen-bonding interactions with MD2 and the interacting residues.