评估抗菌肽抗菌力的系统及其使用方法与流程

本发明涉及一种评估抗菌肽之抗菌力之系统及其使用方法，其是通过规律移动抗菌肽上的氨基酸而产生的具有不同抗菌效果的多肽，并用新创的计算机平台来评估/预测抗菌肽的抗菌效果。

【先前技术】

抗生素可分为二种，一为是由天然微生物(包括细菌、真菌、放线菌属)所产生的一种具有抑制其它微生物生长、生存的次级代谢产物，另一类系通过化学方法合成或半合成的类似化合物。其可包括抗细菌抗生素、抗真菌抗生素以及抑制其他微小病原的抗生素；但临床实验中，抗生素通常是指抗细菌抗生素。自从抗生素被发现至今，已被大量使用乃至滥用，使得具有抗药性的病菌越来越多，人类用来对抗病菌的抗生素武器中持续有效的已愈来愈少，因此积极开发新的、更有效的“抗生素”药物为目前制药产业亟待发展的目标之一。

据美国数据库公司fiercebiotech发布的数据显示，全球抗生素市场，预计从2013年到2018年之间，以1.41％的年复合成长率扩大；虽然抗生素的抗药性问题引起了国际社会的高度关注，部分国家采取了一定的限抗措施，但是由于其特定的功效，全球抗生素市场规模仍将延续增长的态势，根据ims统计，抗感染药品2014年将全球市场规模403亿美元(约2.67千亿人民币)。

抗菌肽，普遍称为抗微生物多肽(antimicrobialpeptides,amps)，其长度大约是由10～80个不等的氨基酸所组成，是许多生物体中一种对抗病原的微小物质。它的功能主要就是帮助生物体对抗体内的有害物质，其作用大致有渗透细菌之细胞膜、穿孔、干扰细菌增殖与促进伤口愈合等功能，未来具有取代抗生素的潜力。

然而，现今大部分的研究皆是对单一抗菌肽的抗菌能力进行研究，并通过置换抗菌肽上的氨基酸提高抗菌能力，然而此种方法需耗费大量时间评估置换哪个氨基酸及位置可提高抗菌能力，且多肽抗菌能力是一个复杂的残基协同现象，而以往仅能通过疏水力矩、非键结作用位能以及接触表面积预测多肽的抗菌效果，且无法准确地预测多肽的抗菌效果，需进行额外实验才能证明多肽的抗菌效果。因此，若不通过分子动力学模拟是难以提供一种系统性并快速的评估设计抗菌肽的方法。

技术实现要素：

本发明提供一种通过“逐步循环重排”(stepwisecircularreshuffling)的方式依序移动已知有抗菌效果之抗菌肽上的氨基酸序列，来设计具有不同抗菌效果、溶血性及抗盐性之多肽，再通过计算机模拟及自由能计算的平台来评估这些抗菌肽的抗菌效果及溶血性。

本发明提供了一个新的分子动力学模拟系统来评估任意一种抗菌肽的抗菌活性及溶血性。此套分子动力学模拟系统，特色在于透过“下沉而上浮”的方式，将抗菌肽快速地拉入仿细菌的仿真生物膜内，并且相较于自然插入大幅缩短所需的平衡时间达10倍左右。根据其模拟轨迹可计算一个新设计出的分配自由能，δgp。此分配自由能δgp可作为一个更好的物理指标来预测某一种多肽的抗菌活性。

有鉴于此，本发明提供一种评估抗菌肽抗菌力的系统，其包括：

(a)第一输入单元，用于多肽的建构；

(b)第二输入单元，用于磷酸双脂层的建构；

(c)第一处理单元，其与所述第一输入单元以及所述第二输入单元连接，该第一处理单元用于将所述第一输入单元建构的多肽与第二输入单元建构的磷酸双脂层置于水溶液环境中进行分子动力学模拟；

(d)第二处理单元，其与该第一处理单元相连接，用于计算所述多肽在每格模拟快照(snapshot)中与膜内脂链尾有碰触的重原子数ni及与膜内脂链尾未碰触的重原子数no，对模拟平衡后4纳秒中所收集的快照做平均，得到<ni>与<no>，进而计算分配自由能δgp，其中所述重原子系指氢原子之外的所有原子；

(e)输出单元，其与该第二处理单元连接，用于将预测信息输出；

其中，所述第一输入单元、所述第二输入单元、所述第一处理单元、所述第二处理单元及所述输出单元通过计算机操作。

本发明又提供一种利用本发明系统，来评估/预测抗菌肽抗菌力的方法，其包含以下步骤：

i.利用第一输入单元建构多肽，将该多肽置于水溶液环境中持续第一时间，进行模拟平衡；

ii.利用第二输入单元建构磷酸双脂层，将该磷酸双脂层置于水溶液环境中持续第二时间，进行模拟平衡；

iii.利用第一处理单元，将所述第一输入单元所建构的多肽与第二输入单元建构的磷酸双脂层于水溶液环境中进行分子动力学模拟；

iv.利用第二处理单元，计算该多肽在每格模拟快照(snapshot)中与膜内脂链尾有碰触的重原子数ni及与膜内脂链尾未碰触的重原子数no，对模拟平衡后4纳秒中所收集的快照做平均，得到<ni>与<no>，进而计算分配自由能δgp，其中所述重原子是指氢原子之外的所有原子；

v.利用输出单元，该输出单元与该第二处理单元连接，将预测信息输出；

其中，所述第一输入单元、所述第二输入单元、所述第一处理单元、所述第二处理单元及该输出单元通过计算机操作。

较佳地，该多肽构型为α-螺旋。

较佳地，该第一输入单元可包含软件(如pymol或discoverystudiovisualizer)创造的模拟多肽，或由实验结果(如从x-raycrystallography或nmr)得知的多肽。模拟平衡的多肽包含利用模拟软件(如namd、vmd、amber、gromacs)以及常用力场(如charmm、amber、gromos、opls)来达到多肽在类生理环境下的平衡状态。

较佳地，所述第二输入单元系包含以“charmm-gui”服务器及相似软件来建膜(磷酸双脂层)。模拟平衡的磷酸双脂层包含利用模拟软件(如namd、vmd、amber、gromacs)以及常用力场(如charmm、amber、gromos、opls)来达到磷酸双脂层在类生理环境下的平衡状态。

较佳地，所述第一处理单元系包含利用模拟软件(如namd、vmd、amber、gromacs)以及常用力场(如charmm、amber、gromos、opls)来达到多肽稳定插入磷酸双脂层的平衡状态。

较佳地，所述第二处理单元的输入元素(多肽与膜内脂链尾有碰触及碰触之重原子数)是依据第一处理单元所输出的模拟轨迹得到的。

较佳地，所述磷酸双脂层是由十六酰油酰磷脂酰胆碱(palmitoyloleoylphosphatidylcholine,简称popc)及/或十六酰油酰磷脂酰甘油酰甘油(palmitoyloleoyl-phosphatidylglycerol,简称popg)构成。

较佳地，所述第一时间为数百皮秒(ps)。

其中，所述预测信息为所述多肽的抗菌效果。

其中，所述分子动力学模拟为，使所述多肽疏水面朝向预先平衡好的磷酸双脂层，而亲水面(即带正电端)朝向另一端，然后当两者的质心投影到与膜垂直的轴上相距为时，再展开“下沉而上浮表面”模拟。

其中，所述“下沉而上浮表面”模拟为所述螺旋构型的多肽被微弱的拉力拖至膜内磷酸双脂上层磷原子之平均坐标下方的位置并维持数纳秒，接着去除对螺旋构型多肽的限制力，再进行10～35纳秒的模拟，使其上浮至膜表面并且平衡。

较佳地，所述计算分配自由能δgp＝-kbtln(<ni>/<no>)或δgp＝-kbt(<lnni>/<lnno>)(当ni及no均不为零时)，其中kb为波兹曼常数、t为温度、ni与no分别为每格模拟快照(snapshot)中与膜内脂链尾有碰触(视为已插入)和未碰触(视为未插入)的多肽内重原子之数量，<>表示4纳秒的移动平均。

较佳地，在本发明抗菌肽中，经计算该δgp与其对应的抗菌效果(mic)彼此间的相关系数(r)大于0.84(“1”代表完美的相关)。在此系列中，mic＝13.49*δgp-5.79。

附图说明

附图1为分子模拟系统的简单示意图。

附图2为分子模拟系统的流程图。

附图3为分子模拟结果的示意图。

附图4为表示多肽溶血性的效果图

【实施方式】

利用本发明系统所设计出的具低溶血性的抗菌肽，该抗菌肽如下列式i所示：

xap1ybp2zcn(p3vd)；

其中x、y、z、v以及n均为正整数；

且x＝0～6、y＝0或4、z＝0～6、v＝0～4以及n＝0或1；

其中p1、p2以及p3系指-trp-leu-lys-；

其中a、b以及c系选自由trp、leu以及lys所组成之群组；

其中d系指lys。

于较佳实施例中，抗菌肽于分子动力学模拟系统的结果如表1所示：

表1本发明抗菌肽的模拟结果

于较佳实施例中，抗菌肽的抗菌效果如表2所示。

最低抑菌浓度测试

抗菌活性的最低抑菌浓度值被定义为能抑制90％的细菌生长的最低蛋白浓度，而大肠杆菌(atcc25922)被选作为此次wlk异构多肽系列(seqidno:1～13)的最低抑菌浓度测试的代表菌株。蛋白浓度由紫外光/可见光光谱仪(ultrospec1100profromamershambiosciences)在波长280nm的吸光值搭配消光系数计算(gills.c.&vonhippelp.h.,analyticalbiochemistry1989,182(2),319-326)得到，然后被配置为以下七个初始浓度：5、2.5、1.25、0.625、0.313、0.156、0.078mg/ml。将处于对数生长中期、浓度为每毫升5×10⁵菌落形成单位(cfu)的细菌培养液与上述浓度的蛋白溶液于96孔盘中均匀混和(每孔中含有1微升的蛋白溶液与99微升的菌液)。在37℃环境下震荡培养16小时后，由微量孔盘仪(thermomax,moleculardevices)在波长600nm环境下的吸光值来侦测该细菌生长抑制情形。最后，最低抑菌浓度值将通过三次上述重复实验结果来确定。

溶血活性测试

溶血活性是由受测多肽对于人类红血球的溶血性测定。新鲜的红血球和pbs缓冲溶液(ph7.4)均匀混和后在相对离心力(rcf)800g环境下采用离心10分钟进行，然后再由pbs缓冲溶液依体积百分比浓度(v/v)稀释至10％。各个受测多肽预先配置好不同浓度(800、400、200、100、50、25、12.5、6.25、3.13、1.56μg/ml)，分别与pbs缓冲溶液依1：1的体积比例混合(100μl：100μl)。而红血球溶液依体积比例1：1，分别与包含或不包含2％聚乙二醇辛基苯基醚x-100(tritonx-100)(v/v)的pbs缓冲液混合，用来作为正、负控制组样本。然后所有受测样品在37℃环境下静置1小时。接着将样品在rcf800g环境下再离心10分钟，取其上清液并测其波长450nm下的吸光值，溶血活性将依照以下公式而测定：

溶血性(％)＝[a样本-apbs]/[atritonx100-apbs]

表2本发明之抗菌肽的溶血性、抗盐性及抗菌效果

针对原先为螺旋构型的seqidno:1，经”序列逐步位移置换”过程中新设计的12段异构多肽(seqidno:2～13)，其抗菌效果、溶血性及其盐类耐受性如表2所示，进一步运用上述“下沉而上浮”模拟策略与δgp计算方式来衡量其抗菌活性，并分析该13段多肽(含原始序列)的分配自由能(δgp)与实验抗菌性(mic)之间的相关性，得到其相关系数大于0.8(mic＝13.50*δgp-5.79,r＝0.84)，表2的结果更说明该13段多肽在高盐浓度中仍然具有抑制效果，且在本发明提供的这12段异构多肽中有3段的抗菌活性比原先序列更佳，其中seqidno:12同时表现出对人类红血球更低的溶血性。本发明明显地提供了一个新的平台，运用已知的抗菌肽来合理地设计更有效的抗菌肽，并且利用计算机模拟来预测新抗菌肽的活性。

序列表

<110>杨立威

<120>评估抗菌肽抗菌力的系统及其使用方法

<160>13

<170>patentinversion3.5

<210>1

<211>13

<212>prt

<213>人工序列

<220>

<223>prt

<400>1

lyslystrpleulystrpleulystrpleulyslyslys

1510

<210>2

<211>13

<212>prt

<213>人工序列

<220>

<223>prt

<400>2

lystrpleulystrpleulystrpleulyslyslyslys

1510

<210>3

<211>13

<212>prt

<213>人工序列

<220>

<223>prt

<400>3

trpleulystrpleulystrpleulyslyslyslyslys

1510

<210>4

<211>13

<212>prt

<213>人工序列

<220>

<223>prt

<400>4

leulystrpleulystrpleulyslyslyslyslystrp

110

<210>5

<211>13

<212>prt

<213>人工序列

<220>

<223>prt

<400>5

lystrpleulystrpleulyslyslyslyslystrpleu

1510

<210>6

<211>13

<212>prt

<213>人工序列

<220>

<223>prt

<400>6

trpleulystrpleulyslyslyslyslystrpleulys

1510

<210>7

<211>13

<212>prt

<213>人工序列

<220>

<223>prt

<400>7

leulystrpleulyslyslyslyslystrpleulystrp

1510

<210>8

<211>13

<212>prt

<213>人工序列

<220>

<223>prt

<400>8

lystrpleulyslyslyslyslystrpleulystrpleu

1510

<210>9

<211>13

<212>prt

<213>人工序列

<220>

<223>prt

<400>9

trpleulyslyslyslyslystrpleulystrpleulys

1510

<210>10

<211>13

<212>prt

<213>人工序列

<220>

<223>prt

<400>10

leulyslyslyslyslystrpleulystrpleulystrp

1510

<210>11

<211>13

<212>prt

<213>人工序列

<220>

<223>prt

<400>11

lyslyslyslyslystrpleulystrpleulystrpleu

1510

<210>12

<211>13

<212>prt

<213>人工序列

<220>

<223>prt

<400>12

lyslyslyslystrpleulystrpleulystrpleulys

1510

<210>13

<211>13

<212>prt

<213>人工序列

<220>

<223>prt

<400>13

lyslyslystrpleulystrpleulystrpleulyslys

1510