续随子醇Ｌ７１３，２８３二萜类化合物溶液构象的分子模拟

胡建平 刘嵬 唐典勇 成英 陈封政 常珊

摘要：将续随子醇Ｌ７１３，２８３二萜类化合物放入甲醇和氯仿两种不同极性的溶液中，分别进行长时间的分子动力学模拟。结果表明，续随子醇Ｌ７１３，２８３二萜类化合物在氯仿溶液中更稳定，构象波动幅度更小；续随子醇Ｌ７１３，２８３二萜类化合物五元环和十一元环侧链上的Ｏ７１及Ｏ１５与甲醇羟基所形成的较稳定氢键是其易溶于含氧有机溶剂的重要驱动力。基于续随子醇Ｌ７１３，２８３二萜类化合物在氯仿溶液中的３个代表性构象，用密度泛函理论优化获得它们的稳定构象、静电势电荷及键参数。模拟结果与试验吻合较好，为基于配体结构的药物分子设计提供了理论基础。

关键词：续随子醇；构象；药物设计；分子动力学模拟；密度泛函理论

中图分类号：Ｑ９４９．７５３．５；Ｏ６２９．６＋１；Ｏ６２１．１４＋２文献标识码：Ａ文章编号：0439－８114（２０12）16－3545-05

Molecular Modeling of the Solution Conformation of Euphorbia lathyris Alcohol

L713， 283 Diterpenoid Compound

ＨＵ Ｊｉａｎ－ｐｉｎｇ１，ＬＩＵ Ｗｅｉ２，ＴＡＮＧ Ｄｉａｎ－ｙｏｎｇ１，ＣＨＥＮG Ｙｉｎｇ１，ＣＨＥＮ Ｆｅｎｇ－ｚｈｅｎｇ１，ＣＨＡＮＧ Ｓｈａｎ３

（１． Ｍｏｌｅｃｕｌａｒ Ｄｅｓｉｇｎ Ｃｅｎｔｅｒ，Ｃｏｌｌｅｇｅ of Ｌｉｆｅ ａｎｄ Ｃｈｅｍｉｓｔｒｙ， Ｌｅｓｈａｎ Ｎｏｒｍａｌ Ｕｎｉｖｅｒｓｉｔｙ， Ｌｅｓｈａｎ ６１４００４， Ｓｉｃｈｕaｎ， Ｃｈｉｎａ；

２． Ｋｅｙ Ｌａｂｏｒａｔｏｒｙ ｏｆ Ｍｅｄｉｃｉｎａｌ ａｎｄ Ｅｄｉｂｌｅ Ｐｌａｎｔｓ Ｒｅｓｏｕｒｃｅｓ Ｄｅｖｅｌｏｐｍｅｎｔ， Faculty of Biotechnology Industry， Ｃｈｅｎｇｄｕ Ｕｎｉｖｅｒｓｉｔｙ， Ｃｈｅｎｇｄｕ ６１０１０６， Ｃｈｉｎａ； ３． Ｃｏｌｌｅｇｅ ｏｆ Ｉｎｆｏｒｍａｔｉｃｓ， Ｓｏｕｔｈ Ｃｈｉｎａ Ａｇｒｉｃｕｌｔｕｒａｌ Ｕｎｉｖｅｒｓｉｔｙ， Ｇｕａｎｇｚｈｏｕ ５１０６４２， Ｃｈｉｎａ）

Ａｂｓｔｒａｃｔ： Ｔｗｏ ｌｏｎｇ－ｔｉｍｅ ｍｏｌｅｃｕｌａｒ ｄｙｎａｍｉｃｓ ｓｉｍｕｌａｔｉｏｎｓ ｗｅｒｅ ｃａｒｒｉｅｄ ｏｕｔ ｆｏｒ Ｅｕｐｈｏｒｂｉａ （Ｅｕｐｈｏｒｂｉａ ｌａｔｈｙｒｉｓ Ｌ．） ａｌｃｏｈｏｌ Ｌ７１３， ２８３ ｄｉｔｅｒｐｅｎｏｉｄ ｃｏｍｐｏｕｎｄ ｉｎ ｍｅｔｈａｎｏｌ ａｎｄ ｃｈｌｏｒｏｆｏｒｍ ｓｏｌｖｅｎｔｓ， ｒｅｓｐｅｃｔｉｖｅｌｙ． Ｔｈｅ ｒｅｓｕｌｔｓ ｓｈｏｗｅｄ ｔｈａｔ Ｅｕｐｈｏｒｂｉａ ａｌｃｏｈｏｌ Ｌ７１３， ２８３ ｄｉｔｅｒｐｅｎｏｉｄ ｃｏｍｐｏｕｎｄ ｗａｓ ｍｏｒｅ ｓｔａｂｌｅ ｉｎ ｃｈｌｏｒｏｆｏｒｍ ｓｏｌｖｅｎｔ； ａｎｄ ｉｔｓ ｃｏｎｆｏｒｍａｔｉｏｎ ｆｌｕｃｔｕａｔｉｏｎ ａｍｐｌｉｔｕｄｅ ｗａｓ ｓｍａｌｌｅｒ． Ｔｈｅ ｍａｉｎ ｄｒｉｖｉｎｇ ｆｏｒｃｅ ｆｏｒ Ｅｕｐｈｏｒｂｉａ ａｌｃｏｈｏｌ Ｌ７１３， ２８３ ｄｉｔｅｒｐｅｎｏｉｄ ｃｏｍｐｏｕｎｄ ｂｅｉｎｇ ｅａｓｉｌｙ－ｓｏｌｕｂｌｅ ｉｎ ｏｘｙｇｅｎ－ｃｏｎｔａｉｎｉｎｇ ｏｒｇａｎｉｃ ｓｏｌｖｅｎｔ ｗａｓ ｔｈａｔ ｔｈｅ Ｏ７１ ａｎｄ Ｏ１５ ｏｎ ｔｈｅ ｓｉｄｅ－ｃｈａｉｎｓ ｏｆ ｔｈｅ ｆｉｖｅ－ｃｙｃｌｅ ａｎｄ ｅｌｅｖｅｎ－ｃｙｃｌｅ ｆｏｒｍｅｄ ｓｔａｂｌｅ ｈｙｄｒｏｇｅｎ－ｂｏｎｄｓ ｗｉｔｈ ｔｈｅ ｈｙｄｒｏｘｙｌ ｏｆ ｍｅｔｈａｎｏｌ． Ｂａｓｅｄ ｏｎ ｔｈｅ ｔｈｒｅｅ ｒｅｓｐｅｃｔｉｖｅ ｃｏｎｆｏｒｍａｔｉｏｎｓ ｏｆ Ｅｕｐｈｏｒｂｉａ ａｌｃｏｈｏｌ Ｌ７１３， ２８３ ｄｉｔｅｒｐｅｎｏｉｄ ｃｏｍｐｏｕｎｄ ｉｎ ｃｈｌｏｒｏｆｏｒｍ ｓｏｌｖｅｎｔ， ｔｈｅ ｓｔａｂｌｅ ｃｏｎｆｏｒｍａｔｉｏｎｓ， ｅｌｅｃｔｒｏｓｔａｔｉｃ ｐｏｔｅｎｔｉａｌ （ＥＳＰ） ｃｈａｒｇｅｓ ａｎｄ ｂｏｎｄ ｐａｒａｍｅｔｅｒｓ ｗｅｒｅ ｏｂｔａｉｎｅｄ ｔｈｒｏｕｇｈ ｄｅｎｓｉｔｙ ｆｕｎｃｔｉｏｎａｌ ｔｈｅｏｒｙ． Ｔｈｅ ｓｉｍｕｌａｔｉｏｎ ｒｅｓｕｌｔ ｗａｓ ａｃｃｏｒｄａｎｔ ｗｉｔｈ ｔｈｅ ｅｘｐｅｒｉｍｅｎｔａｌ ｄａｔａ， ｗｈｉｃｈ ｐｒｏｖｉｄｅｄ ｔｈｅ ｔｈｅｏｒｅｔｉｃａｌ ｂａｓｉｓ ｆｏｒ ｔｈｅ ｄｒｕｇ ｄｅｓｉｇｎ ｂａｓｅｄ ｏｎ ｔｈｅ ｓｔｒｕｃｔｕｒｅ ｏｆ ｌｉｇａｎｄ．

Ｋｅｙ ｗｏｒｄｓ： Ｅｕｐｈｏｒｂｉａ ａｌｃｏｈｏｌ； ｃｏｎｆｏｒｍａｔｉｏｎ； ｄｒｕｇ ｄｅｓｉｇｎ； ｍｏｌｅｃｕｌａｒ ｄｙｎａｍｉｃｓ ｓｉｍｕｌａｔｉｏｎ； ｄｅｎｓｉｔｙ ｆｕｎｃｔｉｏｎａｌ ｔｈｅｏｒｙ

分子模拟作为理论化学的一个重要领域，已经成为继试验、理论之后的第三种研究手段。常用的计算方法有分子动力学（Ｍｏｌｅｃｕｌａｒ ｄｙｎａｍｉｃｓ，ＭＤ）模拟以及密度泛函理论（Ｄｅｎｓｉｔｙ ｆｕｎｃｔｉｏｎａｌ ｔｈｅｏｒｙ， ＤＦＴ）等。前者不但能给出体系在原子水平上的运动细节，而且还可提供位置涨落和构象变化等详细信息；ＤＦＴ则可以解决原子、分子中的许多问题，如单点能的计算、振动光谱的研究、催化活性部位的预测以及分子结构优化等。综合使用ＭＤ模拟和ＤＦＴ来研究药物分子的优势构象及电荷分布等技术一直是国际上药物化学理论研究的主流［１，２］。续随子（Ｅｕｐｈｏｒｂｉａ ｌａｔｈｙｒｉｓ Ｌ．）为大戟科（Ｅｕｐｈｃｒｂｉａｃｅａｅ）大戟属（Ｅｕｐｈｏｒｂｉａ Ｌ．）二年生草本植物，广泛分布于中国各地。续随子的种子又名千金子，形状椭圆，含有二萜类、三萜类、香豆素和挥发油等成分，在中医临床上主要用于治疗水肿、痰饮、积聚胀满、血瘀经闭等症，还可用于治疗血吸虫、肝硬变、肾性水肿等疾病［３，４］。细胞学研究表明，续随子二萜类化合物有刺激皮肤和抑制肿瘤生长的作用，在体外表现出明显的抗肿瘤和抗ＨＩＶ的生物活性［５］。李忌等［６］以长春碱（Ｖｉｎｃａｌｅｕｋｏｂｌａｓｔｉｎｅ ｓｕｌｆａｔｅ）为阳性对照，从大戟属植物中分离得到了５个对肝癌ＳＭＭＣ－７２２１、肺腺癌以及胃腺癌ＭＣ８０－３细胞的生长均具有较强抑制活性的天然二萜类化合物。Ｄｕａｒｔｅ等［７］和Ｃｏｒｅａ等［８］的体外化疗试验表明，续随子醇二萜类化合物的抗癌活性强度与其化学结构密切相关，在分子五元环和十一元环均有影响抗癌活性的关键碳原子及官能团存在。细胞膜Ｐ－糖蛋白（Ｐ－ｇｌｙｃｏｐｒｏｔｅｉｎ，Ｐ－ｇｐ）表达增加是癌症治疗产生多药抗药性（Ｍｕｌｔｉｄｒｕｇ ｒｅｓｉｓｔａｎｃｅ，ＭＤＲ）的主要原因，而续随子二萜类化合物可有效调节Ｐ－ｇｐ的表达，进而可能避免ＭＤＲ的发生［９－１１］。另外，续随子拥有独特的千金二萜烷型结构，进而可通过环合作用而产生复杂的多环结构，这是生物体中合成多种二萜类化合物的前体物［１２］。最近，Ｊｉａｏ等［１３］用包括二维ＮＭＲ在内的系列光谱方法确定了２个续随子二萜类化合物的结构，并用Ｘ－ｒａｙ方法给出了它们的结构及立体构象。

尽管之前有较多对续随子二萜类化合物结构以及药理学方面的研究，但是结构－功能关系、活性官能团、构象变化和极性条件对分子构象的影响等一系列重要问题尚未见报道。试验以目前结构研究最为彻底的续随子二萜类化合物续随子醇 Ｌ７１３，２８３（剑桥晶体数据库Ｎｏｓ． ＣＣＤＣ７１３２８３）［１３］为例，首先将其置于甲醇和氯仿溶液中进行ＭＤ模拟，研究其构象变化；接着对代表性构象进行ＤＦＴ优化，获得了续随子醇Ｌ７１３，２８３二萜类化合物的收敛精细构象及其结构参数，并对体系稳定性、极性对构象的影响以及分子结构参数进行了较为详细的研究。图１给出了续随子醇Ｌ７１３，２８３二萜类化合物的分子结构，为了后续对其电荷及氢键进行分析，所有原子序号均给出。从图１可以看出，续随子醇Ｌ７１３，２８３二萜类化合物的骨架是由五元环与十一元环骈合形成的，同时在其Ｃ２６和Ｃ３７位顺式骈合了一个环丙烷。

１计算方法

综合使用ＭＤ模拟和ＤＦＴ对续随子醇Ｌ７１３，２８３二萜类化合物的溶液构象进行较为详细的研究。将续随子醇 Ｌ７１３，２８３二萜类化合物的剑桥晶体数据库结构［１３］作为ＭＤ模拟的初始构象。ＭＤ模拟采用ＡＭＢＥＲ ９程序和ＧＡＦＦ力场［１４，１５］，具体操作是首先在溶质外围分别加上甲醇和氯仿１．０ ｎｍ的去头八面体溶剂盒，总共加了４０７个甲醇和１８６个氯仿分子，体系总原子分别为２ ５１８和１ ０１２个；在ＭＤ模拟之前，对体系进行２次能量优化，先约束溶质（约束力常数为 ２．０９×１０５ ｋＪ／（ｍｏｌ·ｎｍ２），用最陡下降法优化５ ０００步，再用共轭梯度法优化５ ０００步，然后去约束后再进行５ ０００步最陡下降法优化和５ ０００步共轭梯度法优化，收敛条件为能量梯度小于４．１８×１０－４ ｋＪ／（ｍｏｌ·ｎｍ２）。ＭＤ模拟分为两步，第一步进行了２ nｓ约束溶质的ＭＤ模拟，约束力常数为４．１８×１０３ ｋＪ／（ｍｏｌ·ｎｍ２），期间温度从０ Ｋ逐步升高到３００ Ｋ，第二步进行１０ ｎｓ的无约束恒温ＭＤ模拟。在模拟中，用ＶＭＤ图像显示软件实时跟踪体系构象。采用ＳＨＡＫＥ算法［１６］约束键长，ＭＤ模拟的积分步长设为２ ｆｓ，非键相互作用截断半径设为１．２ ｎｍ，每隔１ ｐｓ采集１次构象，共采集了１０ ２００个构象。基于ＭＤ模拟方法获得了续随子醇Ｌ７１３，２８３二萜类化合物的３个代表性构象（ＣＬ３＿Ｓｎ，ｎ＝１～３），再使用Ｇａｕｓｓｉａｎ ０３程序包中的Ｂ３ＬＹＰ／６－３１１Ｇ（ｄ，ｐ）计算方法［２，１７］对这些构象进行几何构型优化和振动频率计算，得到它们的能量和精细构型，并计算了静电势（Ｅｌｅｃｔｒｏｓｔａｔｉｃ ｐｏｔｅｎｔｉａｌ，ＥＳＰ）电荷信息。频率计算表明，ＣＬ３＿Ｓｎ（ｎ＝１～３）的所有优化构型均为稳定结构（无虚频）。

２结果与分析

２．１ＭＤ模拟

续随子醇Ｌ７１３，２８３二萜类化合物在甲醇及氯仿溶液中势能和方均根偏差（Ｒｏｏｔ ｍｅａｎ ｓｑｕａｒｅ ｄｅｖｉａｔｉｏｎ，ＲＭＳＤ，反映了模拟体系中原子的空间位置与其初始位置的偏移程度）随时间的动态变化情况见图２。从图２ａ可知，在２ nｓ之后，体系能量基本稳定，甲醇溶液中体系势能是（－２６０±３９．９７） ｋｃａｌ／ｍｏｌ，氯仿溶液中则为（－９６７．９４±２７．７２） ｋｃａｌ／ｍｏｌ，波动率分别为１５．４０％和２．８６％。从图２ｂ可以看出，体系在６ ｎｓ前基本稳定，之后有一定的波动。这可能是因为随着ＭＤ模拟时间的延长，模拟采样更加充分，会在不同的采样空间收集到构象。续随子醇Ｌ７１３，２８３二萜类化合物在甲醇和氯仿溶液中的整体ＲＭＳＤ分别为（０．１０７±０．０２１） ｎｍ和（０．０９３±０．０１９） ｎｍ，体系在非极性氯仿溶液中的构象与初始晶体结构的差别以及构象涨落幅度均小于在极性甲醇溶液中的。势能和ＲＭＳＤ分析结果均表明，非极性环境有利于续随子醇Ｌ７１３，２８３二萜类化合物的稳定，这可能与续随子醇Ｌ７１３，２８３二萜类化合物是非极性的化合物相关。

随后分析了续随子醇Ｌ７１３，２８３二萜类化合物与甲醇溶液所形成的氢键，氢键采用几何方法判定［１８］，即氢供体－氢原子－氢受体角度大于１３５°，氢供体－氢受体间距离小于０．３５ ｎｍ，此处氢供体是甲醇羟基上的氧原子，氢受体为续随子醇Ｌ７１３，２８３二萜类化合物上的９个氧原子（图１）。计算结果表明，除了Ｏ１３、Ｏ４９和Ｏ６９外，其他氧原子（包括Ｏ７）都能与甲醇形成中等强度的氢键，不过参与形成氢键的甲醇分子呈动态变化，可见续随子醇Ｌ７１３，２８３二萜类化合物分子在甲醇溶液中整体运动较为剧烈。通过计算氢键占有率（即氢键构象在整个模拟轨迹中所出现的几率）发现，续随子醇五元环和十一元环上的Ｏ７１和Ｏ１５与甲醇羟基所形成的氢键占有率最高，均超过了２５％，推测这2个氧原子所形成的氢键是续随子二萜类化合物易溶于甲醇等含氧有机溶液［７，１９－２１］的驱动力之一。鉴于续随子醇Ｌ７１３，２８３二萜类化合物与甲醇无法形成很稳定的氢键，不难理解续随子醇Ｌ７１３，２８３二萜类化合物在非极性溶液氯仿中依赖很强的非极性相互作用，其溶解性和稳定性均比在极性溶液甲醇中的大。

２．２结构叠落

从图２ｂ可知，甲醇溶液中的续随子醇Ｌ７１３，２８３二萜类化合物在６．０～６．３ ｎｓ、７．２～７．５ ｎｓ以及９．４～１０．０ ｎｓ段的ＲＭＳＤ分别处于较高的平台，将占大部分时间以及处于这３个时间段的平均结构分别定义为ＭＯＨ＿Ｓ１、ＭＯＨ＿Ｓ２、ＭＯＨ＿Ｓ３和ＭＯＨ＿Ｓ４；而将氯仿溶液中续随子醇Ｌ７１３，２８３二萜类化合物的３个代表构象（即大部区域，７．２～７．５、８．６～８．９ ｎｓ段）分别命名为ＣＬ３＿Ｓ１，ＣＬ３＿Ｓ２和ＣＬ３＿Ｓ３。为了分析它们的构象差别，图３给出了续随子醇Ｌ７１３，２８３二萜类化合物在２种溶液中的代表性构象的叠落结果，图中用灰色Ｓｔｉｃｋ表示分子骨架，而体系中的氧原子用ＣＰＫ模型表示。从图３可以看出，体系在２种溶液中的代表性构象骨架位置和３个环方向变化较小，而涨落较大区域主要是较长的Ｃ６７和Ｃ１１侧链，尤以Ｃ１１侧链明显，侧链围绕Ｃ８－Ｃ１１单键发生了一定的转动。考虑到单键旋转导致的能差不大以及分子在２种溶液中的分子骨架位置接近，所有后面分析仅采用了氯仿溶液中的３个代表性构象作为初始结构，用ＤＦＴ优化获得它们的收敛构象，进而获得体系的力场参数。

２．３ＤＦＴ计算

表１给出了续随子醇Ｌ７１３，２８３二萜类化合物在氯仿溶液中３个代表性构象优化后的相关参数，从表１可以看出，体系电子总能量、焓及吉布斯自由能排序为ＣＬ３＿Ｓ２＞ＣＬ３＿Ｓ３＞ＣＬ３＿Ｓ１，能差都小于３ ｋｃａｌ／ｍｏｌ，这属于单键旋转所导致的能差数量级范围以内。ＨＯＭＯ轨道表示已占有电子能级最高的轨道，是主要的电子给予体，而ＬＵＭＯ轨道表示未占有电子能级最低的轨道，是主要的电子接受体。这２个前线轨道的能差（ＥＬＵＭＯ－ＨＯＭＯ）反映了体系的稳定性，能差越大，电子转移反应概率越低。前线轨道能差分析显示，３个最优体系的反应活性排序为ＣＬ３＿Ｓ３＞ＣＬ３＿Ｓ１＞ＣＬ３＿Ｓ２。从表１还可以看出，分子体积排序为ＣＬ３＿Ｓ２＞ＣＬ３＿Ｓ３＞ＣＬ３＿Ｓ１，大部分时间段所出现的ＣＬ３＿Ｓ１优化构象发生了一定收缩，代表着有利于进入靶点蛋白结合位点的构象。偶极矩计算结果大小顺序为ＣＬ３＿Ｓ２＞ＣＬ３＿Ｓ３＞ＣＬ３＿Ｓ１，绝对值差别在１０％左右，方向均为Ｃ６５指向Ｏ１５和Ｏ４９连线的中点。另外，３个ＤＦＴ最终优化构象叠落结果与ＭＤ模拟后的代表性构象叠落结果（图３）类似，只是骨架叠落（包括Ｃ６７上的苯甲酰氧基和Ｃ８－Ｃ１１键）更好，惟一例外的是Ｃ８－Ｏ７单键所连的乙酰氧基（含Ｏ１和Ｏ７）发生了明显的旋转，这可能也是导致构象间有低于３ ｋｃａｌ／ｍｏｌ能差的主要原因。综上所述，除了偶极矩稍有变化外，其他参数变化较小。试验将ＣＬ３＿Ｓ１的收敛构象确定为续随子醇Ｌ７１３，２８３二萜类化合物的稳定结构，后面将分析最终稳定结构的ＥＳＰ电荷及分子构象模型，这将为后续精细ＭＤ模拟以及抗癌靶点确定的分子对接研究提供有价值的参考。

表２给出了ＤＦＴ优化后３个构象的ＥＳＰ电荷分布。从表２可知，除Ｃ１１、Ｃ１２和Ｃ６５等几个原子的静电势电荷差别较大外，其他原子的类似。观察之前的ＭＤ模拟轨迹发现，Ｃ１１、Ｃ１２和Ｃ６５原子的构象波动较明显。柔性和电荷分布分析均表明这３个碳原子有可能参与了分子的化学反应，该结果与试验预测［１３，２２］一致。基于表２构象１的ＥＳＰ电荷数据，正电性由强到弱的原子排序是Ｃ２、Ｃ１４、Ｃ４６、Ｃ５０、Ｃ７０、Ｃ１２、Ｃ５７、Ｃ２７，其中Ｃ２、Ｃ１４、Ｃ４６、Ｃ５０、Ｃ７０都是羰基碳原子，Ｃ４６还是续随子醇Ｌ７１３，２８３二萜类化合物主体结构（非侧链）中所固有的羰基，不容易发生酯化水解，可能是一个较为重要的活性位点。Ｃ１２，Ｃ５７和Ｃ２７也是易受到亲核进攻的位点，尤以Ｃ１２明显。在构象２中，Ｃ１２的正电性最高；前面氢键分析曾指出，Ｃ１２所连酯键上的Ｏ１５是体系中最容易与周围环境形成氢键的氢受体。因此在续随子醇Ｌ７１３，２８３二萜类化合物活性改造时，可考虑加强Ｃ１２的正电性，比如用羧酸根替代Ｃ１２上的Ｈ２０。

图４给出了构象１分别与构象２及构象３的所有非氢原子的电荷相关性。从图４可以看出，构象１与构象２及构象３的ＥＳＰ电荷相关性分别达到了０．９４４和０．９６７，可见续随子醇Ｌ７１３，２８３二萜类化合物的电荷计算结果可信。基于ＣＬ３＿Ｓ１的最终稳定构象的原子坐标，可获得体系键长、键角、二面角以及非正常二面角的信息，这些均为续随子醇Ｌ７１３，２８３二萜类化合物的精细ＭＤ模拟及改造提供了力场支持。

３小结

用ＭＤ模拟和ＤＦＴ相结合的方法研究了续随子醇Ｌ７１３，２８３二萜类化合物的溶液构象。将续随子醇Ｌ７１３，２８３二萜类化合物放在甲醇和氯仿２个不同极性的溶液中进行了１０ ｎｓ的ＭＤ模拟，结果显示，非极性环境有利于体系稳定。随后用ＤＦＴ优化了ＭＤ模拟所获得的氯仿溶液中的３个代表性构象，进而给出续随子醇Ｌ７１３，２８３二萜类化合物最终稳定构象的电荷参数及键参数，这个试验结果对后续的续随子醇Ｌ７１３，２８３二萜类化合物的精细ＭＤ模拟、分子对接及基于配体结构的抗癌药物分子设计具有一定的指导意义。

参考文献：

［１］ ＳＭＩＴＨ Ｌ Ｊ， ＤＯＢＳＯＮ Ｃ Ｍ， ＶＡＮ ＧＵＮＳＴＥＲＥＮ Ｗ Ｆ． Ｓｉｄｅ ｃｈａｉｎ ｃｏｎｆｏｒｍａｔｉｏｎａｌ ｄｉｓｏｒｄｅｒ ｉｎ ａ ｍｏｌｔｅｎ ｇｌｏｂｕｌｅ： Ｍｏｌｅｃｕｌａｒ ｄｙｎａｍｉｃｓ ｓｉｍｕｌａｔｉｏｎｓ ｏｆ ｈｕｍａｎ ａｌｐｈａ－ｌａｃｔａｌｂｕｍｉｎ［Ｊ］． Ｊ Ｍｏｌ Ｂｉｏｌ，１９９９，２８６（５）：１５６７－１５８０．

［２］ ＨＥＨＲＥ Ｗ Ｊ， ＲＡＤＯＭ Ｌ， ＳＣＨＬＥＹＥＲ Ｐ Ｖ Ｒ， ｅｔ ａｌ． Ａｂ Ｉｎｉｔｉｏ Ｍｏｌｅｃｕｌａｒ Ｏｒｂｉｔａｌ Ｔｈｅｏｒｙ［Ｍ］． Ｊｏｈｎ Ｎｅｗ Ｙｏｒｋ： Ｗｉｌｅｙ ＆ Ｓｏｎｓ，１９８６．

［３］ ＪＡＳＳＢＩ Ａ Ｒ． Ｃｈｅｍｉｓｔｒｙ ａｎｄ ｂｉｏｌｏｇｉｃａｌ ａｃｔｉｖｉｔｙ ｏｆ ｓｅｃｏｎｄａｒｙ ｍｅｔａｂｏｌｉｔｅｓ ｉｎ Ｅｕｐｈｏｒｂｉａ ｆｒｏｍ Ｉｒａｎ［Ｊ］． Ｐｈｙｔｏｃｈｅｍｉｓｔｒｙ，２００６，

６７（１８）：１９７７－１９８４．

［４］ ＳＩＮＧＬＡ Ａ Ｋ， ＫＡＭＬＡ Ｐ． Ｐｈｙｔｏｃｏｎｓｔｉｔｕｅｎｔｓ ｏｆ Ｅｕｐｈｏｒｂｉａ ｓｐｅｃｉｅｓ［Ｊ］． Ｆｉｔｏｔｅｒａｐｉａ，１９９０，４１（６）：４８３－５１６．

［５］ ＨＡＳＬＥＲ Ｃ Ｍ， ＡＣＳ Ｇ， ＢＬＵＭＢＥＲＧ Ｐ Ｍ． Ｓｐｅｃｉｆｉｃ ｂｉｎｄｉｎｇ ｔｏ ｐｒｏｔｅｉｎ－ｋｉｎａｓｅ－Ｃ ｂｙ ｉｎｇｅｎｏｌ ａｎｄ ｉｔｓ ｉｎｄｕｃｔｉｏｎ ｏｆ ｂｉｏｌｏｇｉｃａｌ ｒｅｓｐｏｎｓｅｓ［Ｊ］． Ｃａｎｃｅｒ Ｒｅｓ，１９９２，５２（１）：２０２－２０８．

［６］ 李忌，郑耘，郑荣梁，等．天然二萜类化合物的抗肿瘤活性［Ｊ］． 肿瘤防治研究，１９９５，２２（５）：２７1－２７２．

［７］ ＤＵＡＲＴＥ Ｎ，ＶＡＲＧＡ Ａ，ＣＨＥＲＥＰＮＥＶ Ｇ，ｅｔ ａｌ． Ａｐｏｐｔｏｓｉｓ ｉｎｄｕｃｔｉｏｎ ａｎｄ ｍｏｄｕｌａｔｉｏｎ ｏｆ Ｐ－ｇｌｙｃｏｐｒｏｔｅｉｎ ｍｅｄｉａｔｅｄ ｍｕｌｔｉｄｒｕｇ ｒｅｓｉｓｔａｎｃｅ ｂｙ ｎｅｗ ｍａｃｒｏｃｙｃｌｉｃ ｌａｔｈｙｒａｎｅ－ｔｙｐｅ ｄｉｔｅｒｐｅｎｏｉｄｓ ［Ｊ］． Ｂｉｏｏｒｇ Ｍｅｄ Ｃｈｅｍ，２００７，１５（１）：５４６－５５４．

［８］ ＣＯＲＥＡ Ｇ， ＦＡＴＴＯＲＵＳＳＯ Ｅ， ＬＡＮＺＯＴＴＩ Ｖ， ｅｔ ａｌ． Ｊａｔｒｏｐｈａｎｅ ｄｉｔｅｒｐｅｎｅｓ ａｓ ｍｏｄｕｌａｔｏｒｓ ｏｆ ｍｕｌｔｉｄｒｕｇ ｒｅｓｉｓｔａｎｃｅ． Ａｄｖａｎｃｅｓ ｏｆ ｓｔｒｕｃｔｕｒｅ－ａｃｔｉｖｉｔｙ ｒｅｌａｔｉｏｎｓｈｉｐｓ ａｎｄ ｄｉｓｃｏｖｅｒｙ ｏｆ ｔｈｅ ｐｏｔｅｎｔ ｌｅａｄ ｐｅｐｌｕａｎｉｎ Ａ［Ｊ］． Ｊ Ｍｅｄ Ｃｈｅｍ，２００４，４７（４）：９８８－９９２．

［９］ ＨＯＨＭＡＮＮ Ｊ， ＭＯＬＮＡＲ Ｊ， ＲＥＤＥＩ Ｄ， ｅｔ ａｌ． Ｄｉｓｃｏｖｅｒｙ ａｎｄ ｂｉｏｌｏｇｉｃａｌ ｅｖａｌｕａｔｉｏｎ ｏｆ ａ ｎｅｗ ｆａｍｉｌｙ ｏｆ ｐｏｔｅｎｔ ｍｏｄｕｌａｔｏｒｓ ｏｆ ｍｕｌｔｉｄｒｕｇ ｒｅｓｉｓｔａｎｃｅ： Ｒｅｖｅｒｓａｌ ｏｆ ｍｕｌｔｉｄｒｕｇ ｒｅｓｉｓｔａｎｃｅ ｏｆ ｍｏｕｓｅ ｌｙｍｐｈｏｍａ ｃｅｌｌｓ ｂｙ ｎｅｗ ｎａｔｕｒａｌ ｊａｔｒｏｐｈａｎｅ ｄｉｔｅｒｐｅｎｏｉｄｓ ｉｓｏｌａｔｅｄ ｆｒｏｍ Ｅｕｐｈｏｒｂｉａ ｓｐｅｃｉｅｓ［Ｊ］． Ｊ Ｍｅｄ Ｃｈｅｍ，２００２，４５（１２）：２４２５－２４３１．

［１０］ ＢＡＲＩＬＥ Ｅ， ＣＯＲＥＡ Ｇ， ＬＡＮＺＯＴＴＩ Ｖ． Ｄｉｔｅｒｐｅｎｅｓ ｆｒｏｍ Ｅｕｐｈｏｒｂｉａ ａｓ ｐｏｔｅｎｔｉａｌ ｌｅａｄｓ ｆｏｒ ｄｒｕｇ ｄｅｓｉｇｎ［Ｊ］． Ｎａｔ Ｐｒｏｄ Ｃｏｍｍｕｎ，２００８，３（６）：１００３－１０２０．

［１１］ ＰＵＳＺＴＡＩ Ｒ，ＦＥＲＲＥＩＲＡ Ｍ Ｊ Ｕ，ＤＵＡＲＴＥ Ｎ，ｅｔ ａｌ． Ｍａｃｒｏｃｙｃｌｉｃ ｌａｔｈｙｒａｎｅ ｄｉｔｅｒｐｅｎｅｓ ａｓ ａｎｔｉｔｕｍｏｒ ｐｒｏｍｏｔｅｒｓ［Ｊ］． Ａｎｔｉｃａｎｃｅｒ Ｒｅｓ，２００７，２７（１Ａ）：２０１－２０５．

［１２］ ＪＵＲＹ Ｓ Ｌ，ＲＥＹＮＯＬＤＳ Ｔ，ＣＵＴＬＥＲ Ｄ Ｆ，ｅｔ ａｌ． Ｔｈｅ Ｅｕｐｈｏｒｂｉａｌｅｓ：Ｃｈｅｍｉｓｔｒｙ，Ｔａｘｏｎｏｍｙ ａｎｄ Ｅｃｏｎｏｍｉｃ Ｂｏｔａｎｙ［Ｍ］． Ｌｏｎｄｏｎ：Ａｃａｄｅｍｉｃ Ｐｒｅｓｓ，１９８７．

［１３］ ＪＩＡＯ Ｗ， ＤＯＮＧ Ｗ Ｗ， ＬＩ Ｚ Ｆ， ｅｔ ａｌ． Ｌａｔｈｙｒａｎｅ ｄｉｔｅｒｐｅｎｅｓ ｆｒｏｍ Ｅｕｐｈｏｒｂｉａ ｌａｔｈｙｒｉｓ ａｓ ｍｏｄｕｌａｔｏｒｓ ｏｆ ｍｕｌｔｉｄｒｕｇ ｒｅｓｉｓｔａｎｃｅ ａｎｄ ｔｈｅｉｒ ｃｒｙｓｔａｌ ｓｔｒｕｃｔｕｒｅｓ［Ｊ］． Ｂｉｏｏｒｇ Ｍｅｄ Ｃｈｅｍ，２００９，１７：４７８６－４７９２．

［１４］ ＷＡＮＧ Ｊ Ｍ， ＣＩＥＰＬＡＫ Ｐ， ＫＯＬＬＭＡＮ Ｐ Ａ． Ｈｏｗ ｗｅｌｌ ｄｏｅｓ ａ ｒｅｓｔｒａｉｎｅｄ ｅｌｅｃｔｒｏｓｔａｔｉｃ ｐｏｔｅｎｔｉａｌ （ＲＥＳＰ） ｍｏｄｅｌ ｐｅｒｆｏｒｍ ｉｎ ｃａｌｃｕｌａｔｉｎｇ ｃｏｎｆｏｒｍａｔｉｏｎａｌ ｅｎｅｒｇｉｅｓ ｏｆ ｏｒｇａｎｉｃ ａｎｄ ｂｉｏｌｏｇｉｃａｌ ｍｏｌｅｃｕｌｅｓ［Ｊ］． Ｊ Ｃｏｍｐｕｔ Ｃｈｅｍ，２０００，２１（１２）：１０４９－１０７４．

［１５］ ＷＡＮＧ Ｊ Ｍ，ＷＯＬＦ Ｒ Ｍ，ＣＡＬＤＷＥＬＬ Ｊ Ｗ，ｅｔ ａｌ． Ｄｅｖｅｌｏｐｍｅｎｔ ａｎｄ ｔｅｓｔｉｎｇ ｏｆ ａ ｇｅｎｅｒａｌ ａｍｂｅｒ ｆｏｒｃｅ ｆｉｅｌｄ［Ｊ］． Ｊ Ｃｏｍｐｕｔ Ｃｈｅｍ，２００４，２５（９）：１１５７－１１７４．

［１６］ ＲＹＣＫＡＥＲＴ Ｊ Ｐ， ＣＩＣＣＯＴＴＩ Ｇ， ＢＥＲＥＮＤＳＥＮ Ｈ Ｊ Ｃ． Ｎｕｍｅｒｉｃａｌ ｉｎｔｅｇｒａｔｉｏｎ ｏｆ ｔｈｅ ａｒｔｅｓｉａｎ ｅｑｕａｔｉｏｎｓ ｏｆ ｄｙｎａｍｉｃｓ ｏｆ ｎａｌｋａｎｅｓ［Ｊ］． Ｊ Ｃｏｍｐｕｔ Ｐｈｙｓ，１９７７，２３（３）：３２７－３４１．

［１７］ ＢＥＣＫＥ Ａ Ｄ． Ｄｅｎｓｉｔｙ－ｆｕｎｃｔｉｏｎａｌ ｅｘｃｈａｎｇｅ－ｅｎｅｒｇｙ ａｐｐｒｏｘｉｍａｔｉｏｎ ｗｉｔｈ ｃｏｒｒｅｃｔ ａｓｙｍｐｔｏｔｉｃ－ｂｅｈａｖｉｏｒ［Ｊ］． Ｐｈｙｓ Ｒｅｖ Ａ，１９８８， ３８（６）：３０９８－３１００．

［１８］ ＨＵ Ｊ Ｐ， ＧＯＮＧ Ｘ Ｑ， ＳＵ Ｊ Ｇ， ｅｔ ａｌ． Ｓｔｕｄｙ ｏｎ ｔｈｅ ｍｏｌｅｃｕｌａｒ ｍｅｃｈａｎｉｓｍ ｏｆ ｉｎｈｉｂｉｔｉｎｇ ＨＩＶ－１ ｉｎｔｅｇｒａｓｅ ｂｙ ＥＢＲ２８ ｐｅｐｔｉｄｅ ｖｉａ ｍｏｌｅｃｕｌａｒ ｍｏｄｅｌｉｎｇ ａｐｐｒｏａｃｈ［Ｊ］． Ｂｉｏｐｈｙｓ Ｃｈｅｍ，２００８，１３２（２／３）：６９－８０．

［１９］ ＬＵ Ｚ Ｑ， ＧＵＡＮ Ｓ Ｈ， ＬＩ Ｘ Ｎ， ｅｔ ａｌ． Ｃｙｔｏｔｏｘｉｃ ｄｉｔｅｒｐｅｎｏｉｄｓ ｆｒｏｍ Ｅｕｐｈｏｒｂｉａ ｈｅｌｉｏｓｃｏｐｉａ［Ｊ］． Ｊ Ｎａｔ Ｐｒｏｄ，２００８，７１（５）：８７３－８７６．

［２０］ ＤＵＡＴＲＥ Ｎ， ＭＡＲＩＡ Ｕ Ｆ． Ｌａｇａｓｐｈｏｌｏｎｅｓ Ａ ａｎｄ Ｂ： Ｔｗｏ ｎｅｗ ｊａｔｒｏｐｈｏｌａｎｅ－ｔｙｐｅ ｄｉｔｅｒｐｅｎｅｓ ｆｒｏｍ Ｅｕｐｈｏｒｂｉａ ｌａｇａｓｃａｅ［Ｊ］． Ｏｒｇ Ｌｅｔｔ，２００７，９（３）：４８９－４９２．

［２１］ 焦威，鲁璐，邓美彩，等．千金子化学成分的研究［Ｊ］．中草药，２０１０，４１（２）：１８１－１８７．

［２２］ ＡＬＹＥＬＡＡＧＢＥ Ｏ Ｏ， ＡＤＥＳＯＧＡＮ Ｋ， ＥＫＵＮＤＡＹＯ Ｏ， ｅｔ ａｌ． Ａｎｔｉｂａｃｔｅｒｉａｌ ｄｉｔｅｒｐｅｎｏｉｄｓ ｆｒｏｍ Ｊａｔｒｏｐｈａ ｐｏｄａｇｒｉｃａ Ｈｏｏｋ ［Ｊ］． Ｐｈｙｔｏｃｈｅｍｉｓｔｒｙ，２００７，６８（９）：２４２０－２４２５．