洞察结构、光谱和氢键相互作用的模式nicotinamide-oxalic酸盐(我)用实验和理论方法
- 1勒克瑙大学物理系,印度勒克瑙
- 2部门的材料和环境化学、斯德哥尔摩大学、瑞典斯德哥尔摩
在目前的工作,nicotinamide-oxalic酸(NIC-OXA,形成I)被缓慢蒸发结晶盐的水溶液中。了解分子结构和光谱性质的网卡与OXA co-crystallization后,实验红外(IR)、拉曼光谱签名、x射线粉末衍射(XRPD)和差示扫描量热法(DSC)技术被用来描述和验证的盐。采用密度泛函理论(DFT)方法执行所有理论计算采用B3LYP / 6 - 311 + + G (d, p)功能/基底集合。实验几何参数匹配良好的相关性比单体和二聚体的理论参数,由于覆盖最近的氢键相互作用出现在盐的晶体结构。二聚物的红外光谱和拉曼光谱显示,红色(向下)转移和扩大的乐队(N15-H16), (N38-H39)和(C13 = O14)债券NIC (C26 = O24)、(C3 = O1), OXA (C26 = O25)组,因此参与NIC-OXA盐的形成。原子分子(AIM)分析显示,(N8-H9···O24)是最强的(传统)分子间氢键相互作用的二聚体模型盐相互作用能的最大值−12.1千卡摩尔−1。此外,福克的自然键轨道(NBO)分析矩阵表明,在二聚体模型中,(N8-H9···O24)债券负责盐的稳定的能源价值13.44千卡摩尔−1。前线分子轨道(fmo)分析表明,NIC-OXA(我)盐比NIC活性和更不稳定,能源缺口的NIC-OXA(我)盐小于网卡。全球和本地反应描述符参数计算单体和二聚体模型的盐。亲电、亲核和中性活性的网卡,OXA,单体和二聚体模型的盐被绘制可视化分子静电势(MESP)表面。这项研究提供了有价值的见解结合实验和理论结果可以定义分子的物理化学性质。
1介绍
Pepinsky提出“晶体工程”这一术语表示的过程构建晶体材料使用非共价分子或离子合成的目的,现在,评估或操纵组织性能的关系(Pepinsky 1955;Desiraju 1989;巴塔查里亚et al ., 2018)。这种方法也适用于医药科学,作为原料药(api)拥有不同的氢键可能性(通常负责生物活性)(Desiraju 1995;Desiraju 2010;Duggirala et al ., 2015)。api的物化性能,如力学性能和热学和化学稳定性,可以增强或改变与盐的形成,多晶型物、溶剂化物和cocrystals (库马尔,2018;Kumar et al ., 2018)。制药盐被定义为一种得药物,可以形成完整的API之间的质子转移和coformer (Bhogala et al ., 2005;蔡尔兹et al ., 2007)。盐和cocrystals都可以由多种非共价相互作用,如离子相互作用,π-π叠加,氢键和范德华力(马丁et al ., 2013;骏马,2013;王et al ., 2015;陈et al ., 2016)。盐形成晶体工程是用来修改配方中的药物与化学改性产品,发展,生物制药和治疗属性没有改变化学结构和药理学。晶体工程的概念是理解基本的氢键网络讨论分子和超分子结构之间的关系默尔顿和Zaworotko, 2001;弗尼和Pilati, 2002;山和Zaworotko, 2008)。据报道,超过60%的药物临床试验的失败由于溶解度低,由于分子量增加,亲油性(Sareen et al ., 2012)。因此,药物的发展和有效性一直是一个热门话题的兴趣不同的科学家和研究人员研究药物的超分子结构和增强属性通过晶体工程策略(Vioglio et al ., 2017)。它也证明了co-crystallization过程的药物可以显著提高api的不同物理化学性质(对乙酰氨基酚(斯利瓦斯塔瓦et al ., 2019),加巴喷丁(2008年温格和伯恩斯坦),enoxacin (刘et al ., 2020)、卡维地洛(Hiendrawan et al ., 2017),等等)。
烟酰胺(水溶性维生素酮;NIC)吸引了很多研究人员的关注,这是一个代数余子式和前体的烟酰胺腺嘌呤二核苷酸(NAD),负责细胞功能,生长,生存和新陈代谢(Soma Lalam, 2022)。网卡是NAD的口服活性形式,可以获得来自植物和动物来源如花椰菜、牛肉,黄瓜,和牛奶,可以用来治疗糙皮病,糖尿病,肥胖和饮食,aging-induced (米尔斯et al ., 2016)。网卡已被证明是人类潜在的安全容许(老大et al ., 2020;廖et al ., 2021)。草酸盐(OXA)是用于制药的准备,它已经显示出良好的耐受性和低毒性(贝尼托·et al ., 2022 a)。它的主要来源是植物、蔬菜、水果、坚果和种子(Chukwuebuk Chinenye, 2015)。OXA包含两个羰基(C = O)和羟基(OH)组、中(C = O)组与质子给予体可以形成氢键,而哦组与质子受体可以形成氢键,所以很容易形成盐或cocrystals (贝尼托·et al ., 2022 b)。《盐的形成发生在一个质子从OXA转移到网卡的吡啶环上的氮原子。
在延续先前发表的工作nicotinamide-oxalic酸(II)盐(时et al ., 2022)和cocrystal nicotinamide-citric酸(时et al ., 2019),在目前的研究中,特征,光谱分析、热分析、氢键图案,反应分析NIC-OXA(我)盐进行。为此,我们采用了一种结合实验和密度泛函理论(DFT)方法来描述单体和二聚体的结构和物理化学性质的模型NIC-OXA(我)盐。分子结构和光谱获得的B3LYP功能采用6 - 311 + + G (d, p)基组与实验结果比较,得到一个洞察API (NIC)之前和之后的模式与coformer的盐的形成(OXA)。x射线粉末衍射(XRPD) (Spiliopoulou et al ., 2020)模式得到了讨论《盐的相组成。差示扫描量热法(DSC) (Amarachi 2020)技术被用来研究热行为和评估physiosorbed水分子的痕迹。目的和NBO分析证明了氢键在盐的性质和强度,分别。此外,全球和当地的反应描述符计算检查盐的化学性质。这项研究揭示了几何和物理化学性质NIC-OXA(我)盐。
2实验细节
NIC-OXA(我)盐制备混合NIC的克分子数相等的比例和OXA微孔水。最初,过度饱和溶液在室温下制备不饱和,加水。解决方案是使用一种高级绘画纸滤纸过滤,和水晶产品被隔离在室温下缓慢的蒸发。NIC-OXA(1:1或形式I)由缓慢的蒸发盐结晶从水(Athimoolam Natarajan, 2007)。分析了秋石用x射线粉末衍射(XRPD)和差示扫描量热法(DSC)来验证产品的材料。它是在单斜结晶P21/ c空间群与草酸NIC阳离子和阴离子(单deprotonated)。
2.1红外(IR)光谱
红外测量进行使用力量Alpha-P ATR-spectrometer配有钻石晶体(德国卡尔斯鲁厄)衰减全反射配置。数据evathluation与作品进行了程序(力量、Ettlingen、德国)。单梁背景没有样品和驱动的单光束光谱样本范围400 - 4000厘米−1通过平均20扫描的光学分辨率4厘米吗−1。
2.2拉曼光谱学
拉曼光谱收集使用LabRAM HR 800拉曼光谱仪(Horiba、法国)配备Nd: YAG激光操作在532海里。样本被放置在显微镜玻片和暴露在激光功率为0.8兆瓦通过物镜与×50放大。通常,50扫描是积累1 - 2 s曝光时间从100年到4000厘米−1光谱分辨率为2厘米−1获取光谱。
2.3 x射线粉末衍射(XRPD)
XRPD模式样本收集使用一个苍天x射线衍射仪(PANalytical、荷兰)配备PIXcel检测器和一个单色Cu-Kα1辐射x射线管(λ= 1.54059)。管电压和电流安培值被设定在45 kV和40 mA,分别。样本以硅wafer-based zero-background持有人。白色的产品在2θ固体样本扫描范围有些°,增加步长为0.02。
2.4差示扫描量热法(DSC)
DSC分析样品进行TA Q1000 DSC的仪器。样品(1 - 2毫克)卷曲在non-hermetic铝锅和扫描的升温速率下10°C /分钟连续干燥氮气气氛(流量50毫升/分钟)使用一个类似的空锅作为参考。
3量子化学计算的细节
网卡的电子结构和能量优化,OXA, NIC-OXA盐(单体和二聚体模型)计算DFT / B3LYP / 6 - 311 + + G (d, p)(理论/功能/基础集)方法(承德et al ., 1988;帕尔和阳,1989年;帕尔和阳,1995年使用高斯09年套件()弗里希et al ., 2009)。振动作业和波数的所有正常模式定义基于Pulay的建议(Fogarasi et al ., 1992)使用Gar2Ped程序(马丁和Alsenoy, 1995年)。获得的波数是由波数比例的线性扩展(WLS)过程包括非简谐振动的效果(吉田et al ., 2002)。分子结构的图形表示和计算数据的几何图形绘制与GaussView (弗里希et al ., 2000)。几何和拓扑参数的相互作用的原子的债券在债券临界点(bcp)研究了巴德的原子在分子的量子理论(巴德,1990)使用AIMALL程序包(巴德和Cheeseman, 2000)。全球反应描述符(硬度(η),柔软(S),化学势(μ),亲电性指数(ω),和电负性(χ))计算最高占据分子轨道的能量(人类)和最低未占据分子轨道(LUMO)所有分子(皮尔森,1963;皮尔森,1989;陆et al ., 2022)。当地的反应描述符(fk,年代k,ωk)可以预测原子适合亲核,亲电和激进的攻击(Flores-Holguin et al ., 2019)。
4结果与讨论
4.1几何优化和能量
网卡的结构,OXA NIC-OXA(我)盐的晶体状态(Athimoolam Natarajan, 2007;Jarzembska et al ., 2014;巴塔查里亚,2020)所示补充数据S1-S3分别和网卡的基态结构优化,OXA,单体和二聚体的模型NIC-OXA(我)盐所示补充数据S4, S5;图1,2,分别。网卡的基态几何参数优化,OXA,单体和二聚体的模型NIC-OXA(我)盐及其实验值来源于结晶数据给出了补充表S1。NIC的基态能量优化,OXA,盐和单体和二聚体模型(我)−261740.6278−237480.5176−499192.8929−998362.5581千卡摩尔−1,分别。
比较了网卡的优化参数和单体之间NIC-OXA模型(我)盐,表明计算不超过0.0030/0.0054不同,3.6/2.4°,1.6/1.0°实验/理论键长,键角,和二面角角度,分别,除了连接与氢原子的债券。然而,0.0044/0.0101在实验/理论的主要变化值(C10-N8)债券中注意到NIC-OXA(我)盐,导演氢键的形成:(N8-H9···O2 = C3)和(N8-H9···O6-C4),分别在NIC (API)和OXA (coformer)(如所示图1)。的键角和二面角角度的变化相应地这些债券被发现。
同样,在比较实验/理论NIC-OXA OXA和单体模型的几何参数(我)盐,0.0406/0.0292的重大变化,8.2/6.7°和2.6/2.0°的键长,键角,和二面角角度,分别被发现,除了参与的债券持有氢原子和分子间氢键。明显不同的0.0862/0.1183的观察实验/优化债券的长度参数(O1群= C3)债券从OXA单体NIC-OXA(我)盐。这是由于这样的事实:个人OXA氢原子H8连着O5一个羧基,在的NIC-monomer OXA(我)盐、氢被转移到网卡的吡啶环上的氮原子,导致(N8-H9)债券。同样,相应的键角变化观察到和二面角角度,导致分子间氢键的存在(N8-H9···O2 = C3)在盐单体模型(我)(如所示图1)。
单体和二聚体的优化参数模型NIC-OXA(我)盐(Athimoolam Natarajan, 2007)也与实验值相比(表补充表S1)。所示补充表S1从实验结果,计算值复制键长,键角,和二面角0.0021/0.0135,1.1°/ 1.0°,和4.5°的单体、二聚体/ 0.7°(C22-N8),(同样的= O5),和(C18-C12-C13 = O14) / (C12-C18-C20-C22)债券,分别观察,除了债券参与氢键的形成。然而,0.0014的差异(C3 = O2)和0.0076 (C4-O6)债券被认为由于分子间氢键的存在(N8-H9···O6-C4)和(C10-H11···O2 = C3),分别(所示图1)。理论的几何参数是在良好的协议与实验结果在盐二聚体模型(我)由于公司最近的氢键相互作用。
可以看出C-N-C键的键角的变化从实验/单体和二聚体模型被注意到123.04°/ 123.29°,123.60°,分别,这证实了电荷转移现象NIC-OXA(我)盐的报道艾伦(2002)。
4.2振动作业和波数
网卡,OXA NIC-OXA(单体/二聚体)盐包含15日8日,23/46的原子,分别;因此,他们拥有39个,18日,63/132正常的振动模式,活跃在红外和拉曼光谱。结合实验和理论振动作业和波数的两个成分NIC, OXA, NIC-OXA(单体/二聚体)盐了补充表S2-S4,分别。拉曼强度,直接正比于喇曼散射截面,可以导出了喇曼散射振幅DFT-based计算直接不给这些强度(龙格,总值1984;Guirgis et al ., 2003)。分子的红外光谱(相关振动)通常被用来识别官能团和化学键分子系统。然而,定义这些振动红外波段实验是一项艰巨的任务,需要理论计算。温度和非谐效应通常不被认为是在计算,因此实验和计算谱之间的差异被发现(Buelna-Garcia et al ., 2020;Buelna-Garcia et al ., 2021)。已经提出,实验红外光谱通常观察到室温(293 K)和理论上模拟在0 K。观察到的差异,计算红外光谱的波数NIC-OXA(我)盐也可能是由于这样的温度效应。
对比实验和计算红外光谱和拉曼光谱网卡/ OXA所示补充数据S6-S9。对比观察的红外和拉曼光谱NIC-OXA(我)盐单体和二聚体的模拟光谱模型描述了图3,4,分别。
表1代表债券长度和拉伸波数的显著变化由于NIC-OXA的形成(我)盐的两个成分:网卡,OXA。图1,2显示,在单体/二聚体,(C3 = O2和C4-O6) /四(C26 = O24), (C3 = O1), (C26 = O25),和(C3 = O2)组OXA形成氢键(N8-H9、C10-H11和N8-H9) / (N8-H9), (N38-H40) / (N15-H16) / (C41-H42), (C22-H23)和(C41-H42)的网卡,分别。
4.2.1戒指准备NIC的振动
芳香胺的氮伸展振动被分配在该地区的1382−1266厘米−1(西尔弗斯坦et al ., 1981)。个人NIC的伸缩振动(C6-N3)估计1272厘米−1观察到的峰值对应的,1253/1256厘米−1中所描绘的一样,在红外/拉曼光谱补充表S2。同样,在单体和二聚体模型中,相同的拉伸振动(C22-N8)计算1116/1116厘米−1,红外/拉曼峰1126/1129厘米−1分别给出了补充表S4。很明显,碳氮从结合网卡是免费的,但是在单体和二聚体模型中,质子化作用发生在从OXA吡啶氮氢的转移。因此,一个大的区别是NIC中观察到的波数和单体和二聚体。见表1,减少计算拉伸波数(氮)债券的156/156厘米1键长增量的0.007/0.0095从网卡单体/二聚体盐证实了其参与的分子间氢键NIC-OXA(我)盐。在NIC (C6-H12)伸展振动计算为3010厘米−1在3013/3010厘米,相应地分配−1在红外/拉曼光谱。(碳氢键)平面,平面外弯曲变形计算为1341和980 cm−1,相应的观察到1340/1346和979/972厘米−1,分别。单体和二聚体,(C10-H11)伸展振动计算3038/3056厘米−1对应的拉曼峰3023厘米−1。不对称的扭力,环估计418厘米−1并观察到411厘米−1在红外和415厘米−1拉曼光谱地区(中给出补充表S4)。
4.2.2酰胺(C = O-NH2)集团NIC的振动
NIC由酰胺基团的吡啶环碳原子。个人NIC的羰基伸缩振动(C4 = O1)组织估计1712厘米−1。此外,同样的羰基(C13 = O14)单体和二聚体拉伸模式估计1726/1,668厘米−1并在1689/1688厘米观察实验−1在红外/拉曼光谱。很明显从上述波数这羰基是免费的NIC和单体的结合;因此,有一个大的区别在模拟和观察到的波数的值没有网卡和单体的波数的减少。二聚体模型中,有一个减少44厘米−1的增量计算债券长度值0.0072。因此,这一组的二聚体参与分子间氢键分子与周边NIC(所示图2)。拉伸振动(h)债券的酰胺基影响氢键的存在(曹et al ., 2017)。转让NIC,伸展的振动(N2-H10) / (N2-H11)计算3523/3410厘米−1,分别。单体,拉伸(N15-H16)和(N15-H17)的振动模式计算为3520和3411 cm−1,分别。的拉伸振动(N38-H40)和(N15-H16) / (N38-H39)计算在3469和3392/3141厘米−1并观察到3335厘米−1红外光谱和3139厘米−1分别在拉曼光谱。很明显,NH2集团在孤立的NIC和单体自由结合。同时,二聚体,分别为54和269厘米的减少−1,增加债券长度为0.0016和0.0225,分别在(N38-H40)和(N16-H17) / (N38-H39)债券,表明分子间氢键的存在(N38-H40···O1 = C3), (N15-H16···O1 = C3)和(N38-H39···O14 = C13)。因此,二聚体值与实验值。
4.2.3羧基OXA的振动
的拉伸振动(C = O)集团1780 - 1680厘米的范围内−1(施立夫et al ., 1951;Koczon et al ., 1999)。OXA分子,伸缩振动(C2 = O3)集团的估计1788厘米−1。在单体和二聚体模型中,相同的拉伸振动(C3 = O2) / (C26 = O24)计算1722/1,652厘米−1,观察到的峰值对应1647/1646厘米−1在红外/拉曼光谱,这表明,在这两种盐单体/二聚体模型,在计算波数有减少66/136厘米−1从OXA以及债券长度的0.0693/0.0568的增量,分别,因此揭示分子间氢键的存在(N8-H9···O2 = C3)和(N8-H9···O24 = C26),分别指导盐形成(如中给出补充表S4)。氢气从OXA从O1群转移到N8(吡啶环)和O1和N31 O24 N8,分别在单体和二聚体。这个(C2-O1)债券在个别OXA计算1127厘米−1和拉曼峰观察到1110厘米−1。
4.3 x射线粉末衍射(XRPD)
XRPD分析(Klitou et al ., 2020)进行标识隔离材料的相组成。XRPD数据收集初始组件(烟酰胺和草酸)并与实验数据(图5 c)和模拟XRPD (图5 d)获得单晶x射线衍射。图5说明所得NIC-OXA(我)盐发现符合模拟XRPD模式,并没有额外的山峰从最初的组件。因此,孤立的材料是一种单相所需的材料和不包含未反应的初始组件。
4.4差示扫描量热法(DSC)
评估存在的痕迹physiosorbed水分子没有发现XRPD数据、差示扫描量热法(DSC)进行了分析。图6显示了DSC曲线(升温速率:10°C /分钟下氮净化气体)。很明显,准备盐有一个吸热峰在207°C。没有其他自由水或相变吸热峰。尽管从水中材料准备,没有一丝physiosorbed水从DSC热谱。图6说明了明显不同NIC-OXA盐的熔点比初始组件。DSC和XRPD数据分析表明,准备NIC-OXA(我)盐是一个水晶阶段和自由从初始组件和水分子的痕迹(Shimpi et al ., 2017)。
4.5在债券关键路径拓扑和能量参数(bcp)
巴德的原子在分子的量子理论(QTAIM)表明,一双是原子间键的路径显示了分子间氢键相互作用(Jabłoński 2015;Jabłoński 2019)。DFT或波函数的方法可以用来执行拓扑分析理论只有通过知识的电子密度(Tognetti Joubert, 2014)。QTAIM,成键分子中两个原子之间的交互验证的位置(3−1)债券临界点(bcp)和相关的债券路径(原子相互作用线)之间的两个或两个以上相互作用的原子。几何和拓扑参数,如电子密度(ρb)和电子密度(∇拉普拉斯算子2ρb在bcp),用于提供的强度和性质的存在氢键。科赫和Popelier表明氢键的存在(质子供体···受体)依赖于ρb的范围内(0.002 - -0.040 a u。)和(∇2ρb0.024 - -0.139 a u。) (科赫和Popelier, 1995年)。此外,Rozas等人的标准来描述分子的氢键相互作用如下:(i) (∇2ρb)< 0和Hb< 0强劲和共价性质;(2)(∇2ρb)> 0和Hb< 0为媒介和部分共价性质;和(3)(∇2ρb)> 0和Hb> 0为弱和静电性质(Alkorta et al ., 1998)。此外,强和弱相互作用的概念与传统和非传统氢键(Desiraju 2002)。
QTAIM分子图给出了单体和二聚体模型的优化几何图形盐(我)给出补充图S10;图7,分别。氢键的几何和拓扑参数(分子内和分子间)为单体和二聚体模型的盐(我)给出补充表S5;表2,分别。
表2。几何和拓扑参数(键长)的债券的分子间氢键相互作用的原子NIC-OXA(我)盐:电子密度(ρb),电子密度的拉普拉斯算子(∇2ρb),电子的动能密度(Gb潜在的能量密度(V),电子b),总电子能量密度(Hb)债券临界点(bcp)和估计的相互作用能(Eint)。
氢键的标准如下:债券距离的两个相互作用的原子必须小于这些原子的范德华半径之和为单体和二聚体模型的盐(我),列出补充表S6、S7,分别。在盐单体模型(我),最强的分子间氢键(传统)(N8-H9···O2)相互作用能−8.0006千卡摩尔−1拥有中、部分共价性质基于科赫和Popelier标准。分子间氢键的相互作用(C10-H11···O2)和能源价值−5.8986千卡摩尔−1被发现的上述和非常规性质相同。另一个氢键(N8-H9···O6)计算相互作用能的价值−2.9493千卡摩尔−1还拥有一个媒介(部分共价)(补充表S5)。
在盐二聚体模型(我),传统的基于能量的分子间氢键的顺序值(N8-H9···O24) > (N38-H40···O1) > (N15-H16···O1) > (N38-H39···O14),这显示了媒体和部分共价性质,而(C22-H23···O25) > (C41-H42···O2) > (C41-H42···O1)是非常规的形成分子间氢键和直接NIC-OXA(我)盐(给出表2)。
4.6自然键轨道(NBO)分析
NBO分析是一个强大的工具来检查共轭或hyper-conjugative电荷转移相互作用现象,提供了一个详细的研究intramolecular-and分子间氢键的分子系统(由漫画家et al ., 2020)。(即电子基的趋势。,the intensity of charge transfer between electron donor and acceptor) is generally dependent on the extent of stabilization energy E (2) and decides the strength of the bond in a particular molecule. The electron density delocalization between an unoccupied (Rydberg or anti-bonded) and an occupied (lone pair or bonded) orbital corresponds to a stabilized donor–acceptor interaction (Anugrah et al ., 2022)。亲水相互作用的二阶micro-disturbance理论的NBO依据单体和二聚体模型NIC-OXA(我)盐了补充表S8 S9,分别。
在盐单体模型,最强的分子内氢键相互作用在单位1 (OXA)单位2 (NIC) [LP (3) O2→π* (O1群= C3)]和[LP (1) N15→π* (C13 = O14)]从最大E(2)值82.84和63.45千卡摩尔−1分别分别稳定分子和环R1。的电荷转移[单元1单元2]和[单元2单元1];[LP (2) O2→σ* (N8-H9)]和[σ(C10-H11)→σ* (O2 = C3)]与E(2)4.50和0.25千卡的摩尔−1分别证实了很强的分子间氢键的存在(N8-H9···O2)和(C10-H11···分别为O2 = C3)(中描述补充表S8)。
同样,在盐的二聚体模型、电荷转移相互作用在单位1 (OXA), 2 (NIC), 3 (OXA)和4 (NIC) [LP (3) O1→π* (O2 = C3)], [LP (1) C12→π* (N8-C10)], [LP (3) O24→π* (O25 = C26)],和[LP (1) C35→π* (N31-C33)]与E(2)值为88.62,194.72,101.58,和297.50千卡摩尔−1分别确认的分子内电荷转移相互作用与稳定NIC-OXA(我)盐。而其他分子间电荷转移相互作用,如[LP (1) O1→σ* (N15-H16)], [LP (1) O1→σ* (N38-H40)],[σ(N15-H16)→σ* (O1群= C3)],[σ(C22-H23)→σ* (O25 = C26)], [LP (2) O14→σ* (N38-H39)], [LP (2) O24→σ* (N8-H9)],和(σ(N38-H40)→σ* (O1群= C3)],以最大E(2)值为10.18,12.17,0.18,0.06,0.95,13.44,和0.36千卡摩尔−1分别揭示了强分子间氢键:(N15-H16···O1), (N38-H40···O14), (N15-H16···O1 = C3), (C22-H23···(N38-H39 O25 = C26)···O14), (N8-H9···O24)和(N38-H40···O1 = C3),因此形成的盐(我)(中描述补充表S9)。
4.7前沿分子轨道(fmo;HOMO-LUMO)分析
fmo被广泛用来解释分子系统的光学和电子性质Padmaja et al ., 2009)。最高占据分子轨道(人类)和最低未占据分子轨道(LUMO)表示捐赠和接受电子的能力,分别和他们的能量对应于电离势(IP)和电子亲和能(EA) (福井et al ., 1952)。这些轨道不仅解释电荷转移属性还分子化学反应和电子传输性能(Padmaja et al ., 2009)。分子有一个低HOMO-LUMO能源缺口是化学反应,更少的动力学稳定性,比那些柔软高能隙(Aihara 1999;Ruiz-Morales 2002;时et al ., 2021)。网卡的HOMO-LUMO情节,OXA,单体和二聚体盐所示的模型补充数据S11-S13;图8,分别。S11补充数据,向;图8显示的能隙值NIC (5.5356 eV) >单体盐(2.2044 eV) >二聚体(1.5896 eV);因此,盐(我)是化学反应(不稳定)比NIC (API)。
4.8化学反应描述符
化学反应描述符集中研究全球和当地的任何分子的反应,为下一部分要讨论的。
4.8.1全球反应描述符
提供洞察全球反应描述符是有用的化学反应活性和稳定性的化学实体。各种全球反应描述符,如电负性(χ),全球硬度(η),全球柔软(S),化学势(μ),与亲电性指数(ω),可以定义全球反应趋势,计算通过使用前沿分子轨道的能量(HOMO和LUMO)。χ的参数定义了电子吸引能力的分子向自己和是一个重要的参数提供卓有成效的亲电试剂的活跃的站点信息。根据皮尔森,分子的最低反应是由于更大的电阻值的要求努力改变它的电子分布和性质(皮尔森,2005)。对于任何分子,η和S值成反比。表3报道全球反应描述符NIC, OXA,单体和二聚体的模型NIC-OXA(盐)。
表3。计算E人类ELUMO,能量带隙(El- eH),化学势(μ)electro-negativity(χ),全球硬度(η),全球疲软(S)和全球亲电性指数(ω)在298.15 K为烟酰胺(NIC),草酸(OXA)和NIC-OXA(我)盐用B3LYP / 6 - 311 + + G (d, p)。
所示表3χ值表明,OXA最高倾向于吸引电子,所以它有利于亲核攻击。η值躺在的OXA >网卡>单体NIC-OXA(我)盐>二聚体NIC-OXA(我)盐。这意味着NIC-OXA的二聚体(我)盐比单体柔和,网卡,OXA,因此,OXA将与网卡容易反应。良好的亲电试剂的特点是μη低和高反之亦然一个好的亲核试剂。此外,稳定的分子一定是ω的最低价值。ω的价值在于NIC的范围< OXA <单体的NIC-OXA(盐I) <二聚体的NIC-OXA(盐)。所以,很明显,NIC OXA中最为稳定的单体、二聚体模型的盐(我)。
4.8.2当地反应描述符
福井函数(FF),当地的柔软,亲电性指数都是有用的工具来推断亲核/亲电反应的趋势(Chattaraj和罗伊,2007)。福井三种压缩功能(
从报告的值补充表S11为二聚体模型,观察盐(我),最亲核活性网站、亲电攻击O25和C41原子的网站(
4.9分子静电势(MESP, V (r))表面
MESP映射一直常用的预测分子亲核和亲电子网站,策划通过静电势(ESP)映射到电子密度(ED)表面,这是作为常数。蓝色(大多数阳性潜在表面)和红色(大多数的电负性的潜在地区)颜色代表一个分子的传导和捐赠能力(加戈et al ., 2021)。MESP拓扑分析决定了理解的化学键和分子识别的活跃团体分子系统(Gadre et al ., 2021)。NIC和OXA的MESP图绘制在一个先前的研究(时et al ., 2022所示)(补充表S14系列,S15),显示两个组分的活性结合位点。显示的网卡MESP负电位结合位点的N3原子吡啶环和积极的潜力在北半球2官能团。此外,MESP OXA显示积极的潜力是本地化的羟基(−哦)集团,其中一个参与一个质子的转移从OXA NIC和负电位两个羰基(C = O)团体的形成盐。MESP单体图(补充图S16)和二聚体模型(图9)NIC-OXA(我)盐来自最深最深的红色到蓝色区域。
在单体模型中,氧原子O6和O2 OXA参与形成氢键(N8-H9···O2和N8-H9···O6)中,其中显示一个中立的(绿色)地区(补充图S16)。另一个原子氧(O1和O5)被发现在红色区域,让自己的网站适合亲核攻击的地区。在北半球2组网卡出现在蓝色区域,导致可能的亲电攻击地区的网站。
在二聚体模型中,MESP表面表明(N31-H32)在蓝色区域(适合亲核攻击)和(C27 = O28)在红色区域(容易受到亲电攻击),所示图9。所以,这些反应网站进一步适合与邻近分子形成分子间氢键。
5的结论
分子结构、几何图形、氢键相互作用,反应活性NIC-OXA(形式I)研究了盐使用的实验技术(红外和拉曼)结合DFT计算。结构和振动光谱签名显示在二聚体模型中,吡啶环的(N8-H9)网卡参与分子间的氢键(C26 = O24)群OXA和拥有一个红色的转移在波数增加债券的长度(C26 = O24)与个人OXA二聚体。发现振动波数和盐的二聚体模型匹配的键长值与实验结果由于考虑最近的可能比单体氢键相互作用。获得NIC-OXA XRPD和DSC分析表明,(我)盐有单个水晶纯相。目的分析描述,氢键的强度(N8-H9···O24)是最大的相互作用能值−12.1千卡摩尔−1这是媒体和部分共价性质。此外,NBO结果表明,电荷转移相互作用从OXA→网卡;[LP (2) O24→σ* (N8-H9)]负责NIC-OXA的稳定(我)盐。小HOMO-LUMO能源缺口的二聚体(1.5896 eV) <单体(2.2044 eV) < NIC (5.5356 eV)确保网卡已得到增强后的化学反应形成NIC-OXA(我)盐。因此,盐比NIC柔软和更不稳定。当地的反应性参数定义的大多数活性亲核,亲电攻击网站O25和C41原子,分别在二聚体模型。单体,亲核,亲电反应主要是发现在北半球2NIC和(OXA O1和O5),分别从MESP地图。这些网站是中和在二聚体模型中,由于参与分子间氢键。因此,这项研究提供了有用的实验和理论结果定义structure-reactivity和氢键相互作用,这可能会进一步了解其他药物分子是有用的。
数据可用性声明
最初的贡献提出了研究中都包含在这篇文章/补充材料;进一步调查可以直接到相应的作者。
作者的贡献
PV概念化、执行计算和写了初稿;为提高工作质量;PT指导和监督工作;女士提供了实验的细节。所有作者的文章和批准提交的版本。
确认
金融支持从DST-SERB PT项目(格兰特没有:中国中铁/ 2019/006671)。PV和感谢中央设施计算研究(CFCR),勒克瑙大学印度勒克瑙。感谢克努特女士和爱丽丝•瓦伦堡基金会资助。2019.0124)对金融支持。
的利益冲突
作者声明,这项研究是在没有进行任何商业或财务关系可能被视为一个潜在的利益冲突。
出版商的注意
本文表达的所有索赔仅代表作者,不一定代表的附属组织,或那些出版商编辑和评论员。任何产品,可以评估在这篇文章中,或声称,可能是由其制造商,不保证或认可的出版商。
补充材料
本文的补充材料在网上可以找到:https://www.雷竞技rebatfrontiersin.org/articles/10.3389/fchem.2023.1203278/full补充材料
引用
Alkorta,我。Rozas,我。,和Elguero, J. (1998). Non-conventional hydrogen bonds.化学。Soc。牧师。27日,163 - 170。doi: 10.1039 / a827163z
Anugrah, d . s . B。、Darmalim l . V。,Polanen, M. R. I., Putro, P. A., Sasongko, N. A., Siahaan, P., et al. (2022). Quantum chemical calculation for intermolecular interactions of alginate dimer-water molecules.凝胶8,703。doi: 10.3390 / gels8110703
Athimoolam, S。,和Natarajan, S. (2007). Lamellar aggregation and hydrogen-bonding motifs in 3-(aminocarbonyl)pyridinium perchlorate and hydrogen oxalate.Acta Crystallogr。结晶的。结构体。Commun。63年,o263-o266。doi: 10.1107 / s0108270107010098
贝尼托·M。,Barceló-Oliver, M., Frontera, A., and Molins, E. (2022a). Oxalic acid, a versatile coformer for multicomponent forms with 9-ethyladenine.晶体12 (1),89。doi: 10.3390 / cryst12010089
贝尼托·M。,Barceló-Oliver, M., Frontera, A., and Molins, E. (2022b). Oxalic acid, a versatile coformer for multicomponent forms with 9-ethyladenine.晶体12日,89年。doi: 10.3390 / cryst12010089
由漫画家,v E。,Shephard, J. J., Ochiai, T., Price, A. N., Arnold, P. L., Parsons, S., et al. (2020). Quantum chemical topology and natural bond orbital analysis of M–O covalency in M(OC6H5)4(M =钛、锆、高频、Ce、Th, Pa, U, Np)。理论物理。化学。化学。理论物理。22日,16804 - 16812。doi: 10.1039 / d0cp02947e
Bhogala, b R。、Basavoju年代。,和Nangia, A. (2005). Tape and layer structures in cocrystals of some di- and tricarboxylic acids with 4,4′-bipyridines and isonicotinamide. From binary to ternary cocrystals.CrystEngComm7,551 - 562。doi: 10.1039 / b509162d
Buelna-Garcia, c, E。卡贝略港,j·L。,Quiroz-Castillo, J. M., Martinez-Guajardo, G., Castillo-Quevedo, C., de-Leon-Flores, A., et al. (2020). Exploration of free energy surface and thermal effects on relative population and infrared spectrum of the Be6B11− fluxional cluster.材料14 (1),112。doi: 10.3390 / ma14010112
Buelna-Garcia, c, E。,Robles-Chaparro, E., Parra-Arellano, T., Quiroz-Castillo, J. M., del-Castillo-Castro, T., Martínez-Guajardo, G., et al. (2021). Theoretical prediction of structures, vibrational circular dichroism, and infrared spectra of chiral Be4B8 cluster at different temperatures.分子26(13),3953年。doi: 10.3390 / molecules26133953
曹,H。,Qi, F., Liu, R., Wang, F., Zhang, C., Zhang, X., et al. (2017). The influence of hydrogen bonding on N-methyldiethanolamine-extended polyurethane solid–solid phase change materials for energy storage.RSC睡觉。7,11244 - 11252。doi: 10.1039 / c7ra00405b
Chattaraj p K。,和Roy, D. R. (2007). Update 1 of: Electrophilicity index.化学。牧师。107年,PR46-PR74。doi: 10.1021 / cr078014b
陈,Y。李,L。,Yao, J., Ma, Y. Y., Chen, J. M., and Lu, T. B. (2016). Improving the solubility and bioavailability of apixaban via apixaban–oxalic acid cocrystal.结晶的。增长Des。16,2923 - 2930。doi: 10.1021 / acs.cgd.6b00266
承德,L。Weitao, Y。,和Robert, G. P. (1988). Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density.理论物理。启B37岁,785 - 789。doi: 10.1103 / physrevb.37.785
孩子,s . L。,Stahly, G. P., and Park, A. (2007). The Salt−Cocrystal Continuum: the influence of crystal structure on ionization state.摩尔。制药。4,323 - 338。doi: 10.1021 / mp0601345
Chukwuebuk E。,和Chinenye, I. J. (2015). Antidiabetic potential of medicinal plants and their active components.j .制药。医学杂志。科学。10 (2),10 - 19。3008 - 10231019 . doi: 10.9790 /
Desiraju g r (1995)。超分子合成纤维在晶体工程新的有机合成。Angew。化学。Int,艾德。英格兰。34岁,2311 - 2327。doi: 10.1002 / anie.199523111
Duggirala: K。、木、g p F。费舍尔,。、WojtasŁ。,Perry, M. L., and Zaworotko, M. J. (2015). Hydrogen bond hierarchy: Persistent Phenol···Chloride hydrogen bonds in the presence of carboxylic acid moieties.结晶的。增长Des。15日,4341 - 4354。doi: 10.1021 / acs.cgd.5b00628
Fogarasi G。周,X。,Taylor, P. W., and Pulay, P. (1992). The calculation of从头开始分子的几何图形:高效优化自然内部坐标被抵消的力量和经验修正。j。化学。Soc。114年,8191 - 8201。doi: 10.1021 / ja00047a032
行为。,和Pilati, T. (2002). The 'halogen bond' I.O, I.N, Br.N as a temperature function.学报。A58 C321。doi: 10.1107 / s0108767302097829
弗里希,m . J。,Trucks, G. W., Schlegel, H. B., Scuseria, G. E., Cheeseman, J. R., Robb, M. A., et al. (2009).高斯09年修订。瓦林福德CT:高斯公司。
福井,K。,Yonezawa, T., and Shingu, H. (1952). A molecular orbital theory of reactivity in aromatic hydrocarbons.j .化学。理论物理。20岁,722 - 725。doi: 10.1063/1.1700523
Gadre, s R。,Suresh, C. H., and Mohan, N. (2021). Electrostatic potential topology for probing molecular structure, bonding and reactivity.分子26日,3289年。doi: 10.3390 / molecules26113289
加戈,U。,一个zim, Y., and Alam, M. (2021). In acid-aminopyrimidine continuum: Experimental and computational studies of furan tetracarboxylate-2-aminopyrimidinium salt.RSC睡觉。11日,21463 - 21474。doi: 10.1039 / d1ra01714d
Guirgis, g。Klaboe, P。沈,S。,Powell, D. L., Gruodis, A., Aleska, V., et al. (2003).j·拉曼Spectrosc。34岁,322 - 336。doi: 10.1002 / jrs.989
Hiendrawan, S。,Widjojokusumo, E., Veriansyah, B., and Tjandrawinata, R. R. (2017). Pharmaceutical salts of carvedilol: Polymorphism and physicochemical properties.aap。Pharmscitech。18日,1417 - 1425。doi: 10.1208 / s12249 - 016 - 0616 - x
老大J。,Inagaki, E., Fujita, M., Nakaya, H., Mitsuishi, M., Yamaguchi, S., et al. (2020). Effect of oral administration of nicotinamide mononucleotide on clinical parameters and nicotinamide metabolite levels in healthy Japanese men.Endocr。J。67年,153 - 160。doi: 10.1507 / endocrj.ej19 - 0313
Jabłoński, m (2019)。债券的无用路径连接遥远的原子和违反的特权交换通道的概念。ChemistryOpen8,497 - 507。doi: 10.1002 / open.201900109
Jabłoński, m (2015)。QTAIM-based比较agostic债券和分子内charge-inverted氢键。期刊。化学。一个119年,4993 - 5008。doi: 10.1021 / acs.jpca.5b02041
Jarzembska k . N。软管,A。Kamiński, R。,Madsen, A. Ø., Durka, K., and Wozniak, K. (2014). Combined experimental and computational studies of pyrazinamide and nicotinamide in the context of crystal engineering and thermodynamics.结晶的。增长Des。14日,3453 - 3465。doi: 10.1021 / cg500376z
Klitou, P。,Pask, C. M., Onoufriadi, L., Rosbottom, I., and Simone, E. (2020). Solid-state characterization and role of solvent molecules on the crystal structure, packing, and physiochemical properties of different quercetin solvates.结晶的。增长Des。20岁,6573 - 6584。doi: 10.1021 / acs.cgd.0c00751
科赫,U。,和Popelier, P. L. (1995). Characterization of C-H-O hydrogen bonds on the basis of the charge density.期刊。化学。99年,9747 - 9754。doi: 10.1021 / j100024a016
Koczon, P。,lewandowski, W., and Mazurek, A. (1999). Vibrational (FT-IR and FT-Raman) and NMR studies on selected metal (Ca, Mn, Zn) complexes with ortho-meta-and para-iodobenzoic acids.振动器Spectrosc。20岁,143 - 149。doi: 10.1016 / s0924 - 2031 (99) 00032 - 6
库马尔,。,Kumar, S., and Nanda, A. (2018). A review about regulatory status and recent patents of pharmaceutical Co-crystals.睡觉。制药。公牛。8,355 - 363。doi: 10.15171 / apb.2018.042
库马尔,美国(2018年)。制药cocrystals:概述。印第安纳州,j .制药。科学。79年,858 - 871。doi: 10.4172 / pharmaceutical-sciences.1000302
廖,B。,Zhao, Y., Wang, D., Zhang, X., Hao, X., Hu, M., et al. (2021). Nicotinamide mononucleotide supplementation enhances aerobic capacity in amateur runners: a randomized, double-blind study.j . Int, Soc。运动营养。18(1),54岁。doi: 10.1186 / s12970 - 021 - 00442 - 4
刘,L。,Zou, D., Zhang, Y., Zhang, Q., Feng, Y., Guo, Y., et al. (2020). Pharmaceutical salts/cocrystals of enoxacin with dicarboxylic acids: Enhancing在体外抗菌活性enoxacin通过改善溶解性和渗透性。欧元。j .制药。Biopharm。154年,62 - 73。doi: 10.1016 / j.ejpb.2020.06.018
陆,J。,Paci, I., and Leitch, D. (2022). A broadly applicable quantitative relative reactivity model for nucleophilic aromatic substitution (SN基于“增大化现实”技术)使用简单的描述符。化学。科学。43岁,12681 - 12695。doi: 10.1039 / d2sc04041g
马丁,f。、流行、M M。Borodi, G。菲利普,X。,和Kacso, I. (2013). Ketoconazole salt and Co-crystals with enhanced aqueous solubility.结晶的。增长Des。13日,4295 - 4304。doi: 10.1021 / cg400638g
米尔斯,k . F。,Yoshida, S., Stein, L. R., Grozio, A., Kubota, S., Sasaki, Y., et al. (2016). Long-term administration of nicotinamide mononucleotide mitigates age-associated physiological decline in mice.细胞金属底座。24岁,795 - 806。doi: 10.1016 / j.cmet.2016.09.013
默尔顿,B。,和Zaworotko, M. J. (2001). From molecules to crystal Engineering: supramolecular isomerism and polymorphism in network solids.化学。牧师。101 (6),1629 - 1658。doi: 10.1021 / cr9900432
Padmaja, L。,Ravikumar, C., Sajan, D., Hubert Joe, I., Jayakumar, V. S., Pettit, G. R., et al. (2009). Density functional study on the structural conformations and intramolecular charge transfer from the vibrational spectra of the anticancer drug combretastatin-A2.j·拉曼Spectrosc。40岁,419 - 428。doi: 10.1002 / jrs.2145
帕尔·r·G。,和Yang, W. (1995). Density-functional theory of the electronic structure of molecules.为基础。启。化学。46岁,701 - 728。doi: 10.1146 / annurev.pc.46.100195.003413
Ruiz-Morales, y (2002)。HOMO−LUMO差距作为一个索引的分子大小和结构多环芳烃(多环芳烃)和沥青质:理论研究。我。期刊。化学。一个106年,11283 - 11308。doi: 10.1021 / jp021152e
龙格,E。,和Gross, E. K. U. (1984). Density-functional theory for time-dependent systems.理论物理。启。52岁,997 - 1000。doi: 10.1103 / physrevlett.52.997
Sareen, S。,Mathew, G., and Joseph, L. (2012). Improvement in solubility of poor water-soluble drugs by solid dispersion.Int。j .制药。Investig。2 (1),12。2230 - 973 - x.96921 doi: 10.4103 /
山,N。,和Zaworotko, M. J. (2008). The role of cocrystals in pharmaceutical science.药物。今天13 (9),440 - 446。doi: 10.1016 / j.drudis.2008.03.004
Shimpi, m R。,Velaga, S. P., Shah, F. U., and Antzutkin, O. N. (2017). Pharmaceutical crystal engineering using ionic liquid anion–solute interactions.结晶的。增长Des。17日,1729 - 1734。doi: 10.1021 / acs.cgd.6b01698
施立夫、o D。,Heether, M. R., Knight, H. B., and Swern, D. (1951). Infrared absorption spectra of some epoxy compounds.肛交。化学。23日,277 - 282。doi: 10.1021 / ac60050a014
Soma, M。,和lalam, S. K. (2022). An international journal on molecular and cellular biology.摩尔。生物。代表。49岁的1 - 12。doi: 10.1007 / s11033 - 022 - 07459 - 1
Spiliopoulou, M。沃玛,。,Triandafillidis, D. P., Kosinas, C., Fitch, A., Karavassili, F., et al. (2020). Applications of X-ray powder diffraction in protein crystallography and drug screening.晶体10(2),54岁。doi: 10.3390 / cryst10020054
斯利瓦斯塔瓦,K。舒克拉,。,Karthick, T., Velaga, S. P., Tandon, P., Sinha, K., et al. (2019). Molecular structure, spectroscopic signature and reactivity analyses of paracetamol hydrochloride monohydrate salt using density functional theory calculations.CrystEngComm21日,857 - 865。doi: 10.1039 / C8CE01761A
Tognetti, V。,和Joubert, L. (2014). Density functional theory and bader's atoms-in-molecules theory: Towards a vivid dialogue.理论物理。化学。化学。理论物理。16,14539 - 14550。doi: 10.1039 / c3cp55526g
Verma P。,Srivastava, A., Shukla, A., Tandon, P., and Shimpi, M. R. (2019). Vibrational spectra, hydrogen bonding interactions and chemical reactivity analysis of nicotinamide–citric acid cocrystals by an experimental and theoretical approach.新的j .化学。43岁,15956 - 15967。doi: 10.1039 / c9nj03085a
Vioglio, p . C。,Chierotti, M. R., and Gobetto, R. (2017). Pharmaceutical aspects of salt and cocrystal forms of APIs and characterization challenges.药物Deliv放置。牧师。117年,86 - 110。doi: 10.1016 / j.addr.2017.07.001
王,j . R。Yu, X。,Zhou, C., Lin, Y., Chen, C., Pan, G., et al. (2015). Improving the dissolution and bioavailability of 6-mercaptopurine via co-crystallization with isonicotinamide.Bioorg。地中海,化学。列托人。25日,1036 - 1039。doi: 10.1016 / j.bmcl.2015.01.022
温格,M。,和Bernstein, J. (2008). An alternate crystal form of gabapentin: A cocrystal with oxalic acid.结晶的。增长Des。8,1595 - 1598。doi: 10.1021 / cg7008732
关键词:烟酰胺、草酸、NIC-OXA(我)盐、特征光谱,密度泛函理论,反应性
引用:时P,斯利瓦斯塔瓦,经脉P和Shimpi先生(2023)洞察结构、光谱和氢键相互作用模式nicotinamide-oxalic酸盐(我)用实验和理论方法。前面。化学。11:1203278。doi: 10.3389 / fchem.2023.1203278
收到:2023年4月10日;接受:2023年6月20日;
发表:2023年7月5日。
编辑:
尼诺Russo意大利卡拉布里亚大学版权©2023 Verma斯利瓦斯塔瓦,经脉和Shimpi。这是一个开放分布式根据文章知识共享归属许可(CC)。使用、分发或复制在其他论坛是允许的,提供了原始作者(年代)和著作权人(s)认为,最初发表在这个期刊引用,按照公认的学术实践。没有使用、分发或复制是不符合这些条件的允许。
*通信:Poonam经脉,poonam_tandon@yahoo.co.uk;Manishkumar r . Shimpimanishshimpi@gmail.com;Anubha斯利瓦斯塔瓦,anubhaphysics@gmail.com