Theoretical investigation of the effects of diverse hydrogen-bonding characteristics on the 17O chemical shielding and electric field gradient tensors within the active sites of MraYAA bound to nucleoside antibiotics capuramycin, carbacaprazamycin, 3'-Hydroxymureidomycin A, and muraymycin D2.

This study builds upon our prior researches and seeks to investigate and clarify the influences of various characteristics of hydrogen bonds (H-bonds) and charge transfer (CT) interactions, which were detected within the inhibitor binding pockets (labeled as the QM models I-IV) of MraYAA-capuramycin, MraYAA-carbacaprazamycin, MraYAA-3'-hydroxymureidomycin A, and MraYAA-muraymycin D2 complexes by QTAIM and NBO analyses from DFT QM/MM MD calculations, on the 17O chemical shielding (CS) and electric field gradient (EFG) tensors of carboxylate (Oδ), carbonyl (C═O), and hydroxyl (O-H) oxygens in these models. The 17O CS and EFG tensors of these three types of oxygens in QM models I-IV were calculated at the M06-2X/6-31G** level by including the solvent effects using the polarizable continuum model. From the computed 17O CS and EFG tensors in these models, it was found that the nuclear shielding, σiso, for carboxylate or carbonyl oxygen increases (shielding effect) as the H-bond length decreases and the percentage p-character of nOδ/nC═O lone pair partner in the CT interaction enhances. In contrast, the σiso (17O-H) decreases (deshielding effect) with a reduction in the H-bond length as well as with an enhancement in percentage s-character of the nOH lone pair/σ*O-H antibond. By reducing the H-bond length or by increasing p-character of the nOδ/nC═O lone pair, the 17Oδ/17O═C quadrupole coupling constant smoothly decreases, while the 17Oδ/17O═C asymmetry parameter smoothly increases. Moreover, these calculated parameters are in a good agreement with the experimental values. The information garnered here is valuable particularly for further understanding of empirical correlations between 17O NMR spectroscopic and H-bonding characteristics in the protein-ligand complexes.

Intermolecular Interactions of Nucleoside Antibiotic Tunicamycin with On-Target MraYCB-TUN and Off-Target DPAGT1-TUN in the Active Sites Delineated by Quantum Mechanics/Molecular Mechanics Calculations.

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.

ß-Hydroxylation of α-amino-ß-hydroxylbutanoyl-glycyluridine catalyzed by a nonheme hydroxylase ensures the maturation of caprazamycin.

Caprazamycin is a nucleoside antibiotic that inhibits phospho-N-acetylmuramyl-pentapeptide translocase (MraY). The biosynthesis of nucleoside antibiotics has been studied but is still far from completion. The present study characterized enzymes Cpz10, Cpz15, Cpz27, Mur17, Mur23 out of caprazamycin/muraymycin biosynthetic gene cluster, particularly the nonheme αKG-dependent enzyme Cpz10. Cpz15 is a ß-hydroxylase converting uridine mono-phosphate to uridine 5' aldehyde, then incorporating with threonine by Mur17 (Cpz14) to form 5'-C-glycyluridine. Cpz10 hydroxylates synthetic 11 to 12 in vitro. Major product 13 derived from mutant Δcpz10 is phosphorylated by Cpz27. ß-Hydroxylation of 11 by Cpz10 permits the maturation of caprazamycin, but decarboxylation of 11 by Mur23 oriented to muraymycin formation. Cpz10 recruits two iron atoms to activate dioxygen with regio-/stereo-specificity and commit electron/charge transfer, respectively. The chemo-physical interrogations should greatly advance our understanding of caprazamycin biosynthesis, which is conducive to pathway/protein engineering for developing more effective nucleoside antibiotics.

Theoretical Study of Intermolecular Interactions between Critical Residues of Membrane Protein MraYAA and Promising Antibiotic Muraymycin D2.

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.

Determination and evaluation of secondary structure content derived from calcium-induced conformational changes in wild-type and mutant mnemiopsin 2 by synchrotron-based Fourier-transform infrared spectroscopy.

Mnemiopsin 2 from a luminous ctenophore with two functional EF-hand motifs is a calcium-regulated photoprotein that is responsible for emitting a bright blue bioluminescence upon reacting with coelenterazine and calcium ions in Mnemiopsis leidyi. Synchrotron radiation-based Fourier-transform infrared (SR-FTIR) spectroscopy was applied to analyze the distribution of secondary structures, the conformational changes resulting from calcium binding and the structural stabilities in wild-type mnemiopsin 2, as well as its mutant type that possesses three EF-hand motifs. The distribution of secondary structures of these proteins indicates that mutant apo-mnemiopsin 2 has a more stable secondary structure than the wild-type. Analyses of the SR-FTIR spectra revealed that the conformational changes at the secondary structures of both mnemiopsin 2 depend on the calcium concentrations, such that the most noticeable changes in structures of wild-type and mutant mnemiopsin 2 occur at optimum concentrations 6 and 2 mM of calcium chloride, respectively. The addition of calcium to both proteins increases the proportions of their secondary structures in the amide I and II regions. The major amide I bands in the IR spectra of both mnemiopsin­calcium complexes shift towards smaller wavenumbers, whereas their main amide II bands are identified at larger wavenumbers.

Characterization of Dimeric Interactions within Protrusion-Domain Interfaces of Parallel and X-Shaped Conformations of Macrobrachium rosenbergii Nodavirus: A Theoretical Study Using the DFT Method along with QTAIM and NBO Analyses.

The protrusion domain (P-domain; MrNVPd) of Macrobrachium rosenbergii nodavirus (MrNV) exists in two conformations, parallel and X-shaped. We have performed a theoretical study to gain insight into the nature of the dimeric interactions involving the dimeric interfaces within parallel and X-shaped conformations of MrNVPd by applying the quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) analyses in the framework of the density functional theory (DFT) approach. The results reveal that the dimer-dimer interfaces of MrNVPd have hydrogen bonds of common types. Leu255-Lys287, Tyr257-Lys287, Lys287-Ser253, Met294-Cys328, Asp295-Lys327, Ser298-Ser324, Ile326-Asp295, and Cys328-Met294 are the key residue pairs of the dimer-dimer interfaces to maintain the dimer-dimer structures of MrNVPd through charge-charge, charge-dipole, dipole-dipole, hydrophobic, and hydrogen bonding interactions. The strengths of these intermolecular dimer-dimer interactions in the parallel conformation are much greater than those in the X-shaped conformation. The parallel trimeric interface is held basically by electrostatic and hydrophobic interactions. The electrostatic interactions accompanying a strong hydrogen bond of Oγ1-Hγ1···Oγ1 in the Thr276 A-Thr276 D pair maintain the intermolecular interface of two X-shaped MrNVPd dimers.

DFT, QTAIM, and NBO studies on the trimeric interactions in the protrusion domain of a piscine betanodavirus.

Crystal structure of the protrusion domain (P-domain) of the grouper nervous necrosis virus (GNNV) shows the presence of three-fold trimeric protrusions with two asymmetrical calcium cations along the non-crystallographic three-fold axis. The trimeric interaction natures of the interacting residues and the calcium cations with the neighboring residues within the trimeric interface have been studied by the quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) analyses in the framework of the density-functional theory (DFT) approach. The results revealed that residues Leu259, Val274, Trp280, and Gln322 of subunit A, Arg261, Asp275, Ala277, and Gln322 of subunit B, Leu259, Asp260, Arg261, Ala277, Val278, and Leu324 of subunit C are the main residues involved in the trimeric interactions. Charge-dipole, dipole-dipole, and hydrogen bonding interactions make the significant contributions to these trimeric interactions. Among different interacting residues within trimeric interface, residue pair Arg261 B-Leu259C forms the strongest hydrogen bond inside the interface between subunits B and C. It was also found that calcium cations interact with residues Asp273, Val274, and Asp275 of subunits A, B, and C through charge-charge and charge transfer interactions.

A theoretical study on the characteristics of the intermolecular interactions in the active site of human androsterone sulphotransferase: DFT calculations of NQR and NMR parameters and QTAIM analysis.

A theoretical study at the level of density functional theory (DFT) was performed to characterize noncovalent intermolecular interactions, especially hydrogen bond interactions, in the active site of enzyme human androsterone sulphotransferase (SULT2A1/ADT). Geometry optimization, interaction energy, (2)H, (14)N, and (17)O electric field gradient (EFG) tensors, (1)H, (13)C, (17)O, and (15)N chemical shielding (CS) tensors, Natural Bonding Orbital (NBO) analysis, and quantum theory of atoms in molecules (QTAIM) analysis of this active site were investigated. It was found that androsterone (ADT) is able to form hydrogen bonds with residues Ser80, Ile82, and His99 of the active site. The interaction energy calculations and NBO analysis revealed that the ADT molecule forms the strongest hydrogen bond with Ser80. Results revealed that ADT interacts with the other residues through electrostatic and Van der Waals interactions. Results showed that these hydrogen bonds influence on the calculated (2)H, (14)N, and (17)O quadrupole coupling constants (QCCs), as well as (1)H, (13)C, (17)O, and (15)N CS tensors. The magnitude of the QCC and CS changes at each nucleus depends directly on its amount of contribution to the hydrogen bond interaction.