Drug Design and Bioinformatics Unit, Medical Biotechnology Department, Biotechnology Research Center, Pasteur Institute of Iran, Tehran, 13169-43551, Iran; Department of Chemistry, Faculty of Science, Tarbiat Modares University, Tehran, 14115-175, Iran.
Department of Chemistry, Faculty of Science, Tarbiat Modares University, Tehran, 14115-175, Iran.
Solid State Nucl Magn Reson. 2024 Oct;133:101960. doi: 10.1016/j.ssnmr.2024.101960. Epub 2024 Aug 22.
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 MraY-capuramycin, MraY-carbacaprazamycin, MraY-3'-hydroxymureidomycin A, and MraY-muraymycin D2 complexes by QTAIM and NBO analyses from DFT QM/MM MD calculations, on the O chemical shielding (CS) and electric field gradient (EFG) tensors of carboxylate (Oδ), carbonyl (C═O), and hydroxyl (O-H) oxygens in these models. The O 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 O CS and EFG tensors in these models, it was found that the nuclear shielding, σ, for carboxylate or carbonyl oxygen increases (shielding effect) as the H-bond length decreases and the percentage p-character of n/n lone pair partner in the CT interaction enhances. In contrast, the σ (O-H) decreases (deshielding effect) with a reduction in the H-bond length as well as with an enhancement in percentage s-character of the n lone pair/σ* antibond. By reducing the H-bond length or by increasing p-character of the n/n lone pair, the Oδ/O═C quadrupole coupling constant smoothly decreases, while the Oδ/O═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 O NMR spectroscopic and H-bonding characteristics in the protein-ligand complexes.
本研究建立在前人研究的基础上,旨在通过 QTAIM 和 NBO 分析,从 DFT QM/MM MD 计算中,研究和阐明氢键(H 键)和电荷转移(CT)相互作用的各种特征,这些特征存在于 MraY-capuramycin、MraY-carbacaprazamycin、MraY-3'-hydroxymureidomycin A 和 MraY-muraymycin D2 复合物的抑制剂结合口袋(标记为 QM 模型 I-IV)中,对这些模型中羧酸盐(Oδ)、羰基(C═O)和羟基(O-H)氧的 O 化学位移(CS)和电场梯度(EFG)张量的影响。在 M06-2X/6-31G**水平上,通过包括使用极化连续体模型的溶剂效应,计算了 QM 模型 I-IV 中这三种类型氧的 O CS 和 EFG 张量。从这些模型中计算出的 O CS 和 EFG 张量中发现,随着 H 键长度的减小和 CT 相互作用中 n/n 孤对的 p-特征百分比的增加,羧酸盐或羰基氧的核屏蔽,σ,增加(屏蔽效应)。相比之下,随着 H 键长度的减小以及 n 孤对的 s-特征百分比的增加,O-H 的σ(O-H)减小(去屏蔽效应)。通过减小 H 键长度或增加 n/n 孤对的 p-特征,Oδ/O═C 四极耦合常数平滑减小,而 Oδ/O═C 不对称参数平滑增加。此外,这些计算出的参数与实验值吻合较好。这里获得的信息对于进一步理解蛋白质-配体复合物中 O NMR 光谱和氢键特征之间的经验相关性特别有价值。