计算机如何辅助抗病毒药物研发？为特效药物研发赢得宝贵时间

很多人以为，药物的研发只有在穿着白大褂的实验室里可以做。近年来，随着计算机硬件、专业软件的发展，计算机辅助药物设计的方法已日趋成熟，其应用极大地加快了新药研发的速度与效率，已成为现代药物研发的常规方法之一。依赖于这样的技术，研发人员即便足不出户，也能为寻找治疗此次病毒感染的潜在药物上贡献自己的力量。

寻找关键蛋白质

蛋白质是生命的物质基础，是组成细菌、病毒及动植物的一切细胞、组织的重要成分。各类功能性蛋白质在机体中各司其职，维系着整个机体的正常运转。

以病毒为例，病毒是由核酸（DNA或RNA）与蛋白质构成的非细胞形态，介于生命体与非生命体之间，无法自我复制与繁殖，需要寄生在活的宿主细胞内，依赖于宿主细胞的原料、能量供给与场所，完成自我的复制与释放。

病毒的生命周期需要经历吸附、侵入、脱壳、生物合成、组装和释放等六大步骤，病毒的各类功能性蛋白在这些步骤中分工明确、高度协作，才能完成从感染宿主细胞到复制病毒的整个周期。目前，研究人员已从新型冠状病毒（2019-nCoV）中分离出orf1ab，S，E，M，N等10条基因组序列，各自编码相应的病毒蛋白。

通过与同为冠状病毒的SARS病毒类比，我们可以合理地推测出新型冠状病毒基因编码的各类蛋白质的功能。举例来讲，如，orf1ab基因编码orf1ab多聚蛋白，参与病毒RNA的转录与复制，并具有蛋白酶、甲基转移酶等多个功能；S基因编码冠状病毒的表面糖蛋白，也称为棘突蛋白，通过与人体内的ACE2蛋白质结合，直接介导病毒对宿主细胞的感染及融合，这类蛋白如同日冕般分布在病毒的包膜上，“冠状病毒”由此得名。这些功能性蛋白质对病毒的感染与复制发挥着重要作用，单独或同时干扰其中一个或多个蛋白质的功能，抑制其活性，便能阻断病毒感染宿主细胞或在宿主细胞内自我复制的进程，从而起到治疗的效果。

计算机辅助药物设计的技术

因此，一旦获取了病毒蛋白质的晶体结构，基于前期对其功能的了解，我们便能利用计算机辅助药物设计的技术，针对性地寻找潜在的有效药物。那么，这一过程是如何进行的呢？

若是“以貌取物”的话，在以“光滑弹嫩”为美的今天，蛋白质表面坑坑洼洼，“长得”可实在不算漂亮。可恰是这些坑坑洼洼的空腔，才是真正暗藏玄机的地方。

以SARS病毒orf1ab多聚蛋白中的一段，3C-like蛋白酶，也称Mpro蛋白为例。这一蛋白主要负责将多聚蛋白水解为功能性的多肽，以发挥其各自的功能。就好比做鱼的时候，用一把锋利的刀把整条鱼除鳞去脏，并分为鱼头、鱼身、鱼尾，分别预备炖汤、红烧、清蒸一样。这一蛋白酶的活性位点，就好比是这把菜刀的刀刃，藏身于蛋白质表面的空腔（口袋）中。

筛选或合成得到的具有生物活性的小分子可以很好地结合在这一活性位点的口袋中，从而抑制蛋白酶的活性，阻止其将多聚蛋白切割为功能性多肽，从而阻止这些多肽在后续病毒的复制和感染中发挥功能。就好比没收了这把切鱼的刀，或者在刀刃上套了一层使其钝化的保护套，使其无法处理这条鱼，也就无法进行后续的烹饪。

对每一个蛋白质，寻找具有令人满意的生物活性的分子的过程，就好比是面对一把精美绝伦的锁，需要找到一把同样精美绝伦的钥匙，来与之紧密契合。而计算机辅助药物设计所需要做的，就是以高效和低成本的方式，去找到这样一把合理的钥匙，即，基于对蛋白质功能的前期了解，及对蛋白质结构的系统分析，通过计算的方式，评估各类分子在蛋白质口袋中的结合强度和作用模式，从而筛选或设计出最有可能成为蛋白功能抑制剂（或激动剂）的分子。

计算机辅助药物发现

不幸的是，尽管近年来相关技术飞速发展，但找到这样一把精美绝伦的钥匙绝非易事。一个新药从研发到上市往往需要消耗数十年的时间和数以十亿计的美金。但好消息是，总有这样几把锁长得有点像，其钥匙可以通用，就像艾滋病蛋白酶的抑制剂可能可以作用于冠状病毒的蛋白酶，负责埃博拉或流感病毒RNA复制的RNA聚合酶抑制剂可能可以作用于冠状病毒的RdRp一样。

面对急性暴发的疫情，从已上市或已在临床上的“老药物”中寻找合适的分子显然比从头研发新的分子更具时间优势。运用计算机辅助药物设计，基于一种叫作分子对接的技术进行虚拟筛选，我们可以模拟出每个“老药”分子在病毒蛋白质口袋中的结合构象，通过打分函数、自由能计算等方式评估其理论上的结合强度，从而分析该分子成为潜在抑制剂的可能性。就好比当我们已知“锁”的构造时，运用这样的技术，可以不必逐自把每把钥匙都在锁孔里插一遍来寻找能开锁的那把。通过计算机模拟分析的方式，筛选出最可能打开锁的几把钥匙，而后只对这几把钥匙进行测试就可以了。

利用同源模建的技术，我们甚至可以不用知道当前病毒蛋白这把“锁”的构造，在仅有蛋白质氨基酸序列的情况下，构建出病毒蛋白质结构可能的模型，从而提交虚拟筛选作业。尽管据早期的新闻报道，2019-nCoV病毒和SARS病毒在基因组水平的相似度只有70%，但事实上，通过对病毒蛋白的序列比对可以发现，此次新型冠状病毒的某些关键蛋白和SRAS病毒的氨基酸同源性能达到95%以上。因此，依赖于蛋白质晶体结构数据库中，早年研究SRAS病毒时获得的SARS病毒蛋白的晶体结构，我们便可以构建出合理的2019-nCoV的相应蛋白质的结构。