结合自由能的计算

格式：pdf
大小：176.63 KB
文档页数：4

结合自由能的计算

一:自由能计算简介

基于结构的药物设计方法有多种多样，包括分子对接、数据库搜索以及从头设计等等，但是它们最终往往都面临一个同样的问题就是药物分子活性的评

价。许多药物和其它生物分子的活性都是通过与受体大分子之间的相互作用表

现出来的，所以受体和配体之间结合自由能(bindingafinity)评价是基于结构的计算机辅助药物分子设计的核心问题。精确的自由能预测方法能够大大提高

药物设计的效率。在过去的20年中，随着受体一配体相互作用的理论研究以及

计算机辅助药物分子设计方法的快速发展，自由能预测方法的研究受到了越来

越多的关注。

在配体(包括药物)和受体的相互作用过程中，牵涉到两类相互作用，即非键相互作用和共价相互作用。非键相互作用包括静电相互作用、范德华相互作

用以及氢键相互作用等等，这些相互作用都可以通过力场计算进行比较简单的

表达;而共价相互作用则牵涉到化学键的断裂和生成，需要用量子力学的方法进

行考察。在药物分子和受体的相互作用时，一般只牵涉到非键相互作用.而共价相互作用只存在少数体系中。根据热力学原理，结合自由能和结合常数K,存在

定量的关系:ΔG=-RTlnK

自由能的计算方法主要可以分为三类:第一类方法包括自由能微扰(FreeEnergyPerturbationFEN和热力学积分(ThermodynamicIntegration.TI)

方法。这类方法是最为经典的方法，为大家所广泛接受。这类方法从原理上比

较严格，计算结果也较为精确.但需要长时间的数据采集，对计算体系有严格的

限制，只能适合较为简单的情况，因此在药物设计中还没有达到实用化的阶段。

第二类方法包括一系列基于经验方程的计算方法，这类方法把结合自由能分解为不同的相互作用能量项，通过一组训练集并利用统计方法来得到自由能计算

的经验公式。这类方法取样简单，计算量小，在药物设计中应用较多。但这类

方法得到的经验公式很依赖于训练集的选择，对不同的体系具有不同的预测能

力。而且这类方法不能很好的考虑体系的柔性以及溶剂效应，因此它一般只能

对一些体系具有有限的预测能力，也只能作为初筛的手段。第三类方法就是近

几年发展起来的基于分子动力学采样的自由能预测方法，主要包括LIE方法和

MM/PBSA方法。

二:常用的自由能计算方法

2.1自由能微扰(FEP)方法

FEP的基本思想是从一个己知体系出发，通过一系列微小的变化变到另一

个体系，在每一个变化步骤做分子动力学模拟，把每一步的体系势能代入相应

的公式中，就可以得到两步之间的自由能变化，把所有的自由能变化加起来，就得到两个体系之间的自由能变化。如果有两个配体分子A和B，它们和受体S

形成了复合物AS和BS。为了求算这两个配体分子和受体之间结合自由能的差

值，则需要通过热力学循环。先求算配体分子A和B之间的自由能变AG,，然

后再求算出复合物AS和BS之间的自由能变△G,二者之差就是我们需要得到的

自由能变化△AG.

在一般的自由能微扰的计算中，求算的都是两个配体分子和受体之间结合

能的差值。在有些应用中，也可以采用这种方法来计算单个配体和受体之间的

结合自由能。在这种情况下，突变的过程不是一个配体到另外一个配体，而是一个配体到相应的溶剂分子(一般为水分子).

FEP方法有很多变种，比如定长窗口增长(FixedWideWindowGrowth)、慢增

长(slowgrowth)以及动态窗口增长(dynamicmodifiedwindowgrowth)等，这些方

法的原理都是一样的，不同在于计算过程中的步长的选取，所谓定长窗口增长，它的步长是一个固定的数值，由变化步数i决定，若步数为20，则步长为0.05;

慢增长方法中，步长是趋于零的一个数，变化步数则趋于无穷，因为步数多，

每步之间的变化小，所以在每一步只做一步或几步动力学模拟;动态窗日增长方

法中，步长是根据上一步计算自由能的结果而变化的，若自由能变化太大则步

长小一点.若变化太小则步长大一点，在保证了计算精度的同时又大大减少了计算的步骤，可谓集前二者所长

2.2热力学积分(TI)方法

计算的步骤和FEP方法相似，也是通过两次计算先算出两底物之间的自由

能变，再算出复合物之间的自由能变，二者之差就是我们关心的相对自由能。

只是自由能差的定义上与FEP方法有所不同。TI也分普通的热力学积分和动态热力学积分(dynamicmodifiedthermodynamicintegration)两种。

FEP和TI是比较经典的自由能计算的方法，它们的优点在于理论严格，逻

辑清晰，具有普适性。但它的缺陷也非常明显:首先这两种方法计算量大，耗时;其次，它们只能计算差别较小的两态之间的相对结合自由能，当两态差别较大

时，我们很难指定变化的路径。由于这些缺陷，这两种方法在药物设计中的应

用受到了很大的局限。

2.3线性相互作用能(LIE)方法

线性相互作用能(linearinteractionenergy,LIE)的自由能计算方法,是近

几年发展起来的一种基于分子动力学采样的自由能计算方法。这种方法理论基

础来源于非平衡态统计物理学中的线性响应近似(linearresponseapproximation)理论，把自由能分解为极性和非极性的贡献。

它应该是介于前两类方法之间的一类方法。它也需要进行分子动力学的采样，

但和FEP比较，这个方法仅仅需要对体系的始态和终态做动力学采样，因此计

算量比FEP大大减少。和半经验方程方法比较，它的优势也非常明显，首先它

用分子动力学和蒙特卡罗模拟考虑了体系的柔性，在正确预测复合物结构的基础上对体系的构象空间进行采样;其次它通过采用在溶剂环境以及蛋白环境两次

分子动力学采样，考察了去溶剂化效应对结合自由能的影响。

LIE方法的理论基础来源于非平衡态统计物理学中的线性响应理论，这是一个在各个学科领域中都有广泛应用的重要理论。

线性响应理论用于简单离子的水合作用取得了一些成功，例如1990年

Aqvist用自由能微扰的方法模拟了Na+和Ca2+的水溶液体系，计算结果比较好。

1996年，qvist又对多个体系进行了进一步的系统研究，研究表明大部分体系都能较好的符合所给出的公式，但是对一些体系尤其是包含较多氢键给体或受

体的体系产生了一些偏离。

对药物分子和受体之间的结合自由能中的极性部分，可以通过线性响应近

似来处理。但对于非极性相互作用和疏水效应对自由能的影响，处理起来就困

难一些。对于一些长链烷烃实验已经证明其溶剂化自由能和烷烃的链长呈线性的关系。通过计算也已经证明对长链烷烃，其在水中和溶剂之间的平均非键相

互作用能和烷烃的链长有很好的线性关系。因此，我们不难想到自由能中的非

极性贡献和溶质与溶剂之间平均范得华相互作用能存在一定的线性关系

本文引用地址:空间完美搬家到新浪博客!

薛定谔计算结合自由能方法学

薛定谔计算结合自由能的方法学是一种基于量子力学原理的分子模拟方法，用于描述分子与固体物质中的原子和分子的运动。该方法首先对物质进行能级结构的计算，然后通过计算自由能的变化来描述物质中原子和分子的运动。这种方法能够给出非常精确的预测结果，因此在药物设计、材料科学、环境科学等领域有广泛的应用。

薛定谔计算结合自由能方法学的计算过程包括以下几个步骤：

1. 分子建模：首先需要建立分子模型，包括分子的几何结构和电子结构。这可以通过量子力学软件来完成，如Gaussian、Schrödinger等。

2. 薛定谔方程求解：根据建立的分子模型，使用薛定谔方程进行求解。薛定谔方程是描述微观粒子运动规律的偏微分方程，通过数值方法求解该方程可以得到粒子的波函数和能量。

3. 自由能计算：根据求解的薛定谔方程，可以计算分子的自由能。自由能是描述分子能量和熵等性质的物理量，对于化学反应和相变等过程具有重要的意义。

4. 结果分析：通过对计算结果进行分析，可以得到分子的性质和行为，如反应活性、扩散系数、热力学性质等。这些结果可以为实际应用提供重要的理论支持。

薛定谔计算结合自由能方法学的主要优点在于其精度高、适用范围广，可以应用于不同尺度、不同种类的分子模拟。此外，该方法还可以与其他计算方法相结合，如分子动力学模拟、量子化学计算等，进一步提高模拟的精度和可靠性。然而，该方法也存在一些局限性，如计算量大、计算时间长等，需要结合具体应用进行选择和使用。

算结合自由能能量

自由能是热力学中一个重要的概念，它描述了一个系统在特定条件下可以进行的最大有用功。自由能可以分为两种形式，分别是Helmholtz自由能和Gibbs自由能。本文将从这两个方面来探讨自由能的含义和应用。

一、Helmholtz自由能

Helmholtz自由能是描述一个恒温恒容系统的自由能，通常用符号F表示。它由内能U和熵S两个热力学函数组成，可以用以下公式表示：

F = U - TS

其中T是系统的温度，S是系统的熵。Helmholtz自由能描述了系统在恒温下可用的能量，也可以理解为系统可以做的最大有用功。

Helmholtz自由能在实际应用中有着广泛的用途。首先，它可以用来计算系统的平衡条件。在恒温恒容条件下，系统的平衡态对应着Helmholtz自由能取最小值的状态。通过最小化Helmholtz自由能，可以确定系统的平衡态以及各个平衡态之间的相变条件。

Helmholtz自由能还可以用来预测化学反应的进行方向。在化学反应中，如果反应物的Helmholtz自由能大于生成物的Helmholtz自由能，那么反应是自发进行的。反之，如果反应物的Helmholtz自由能小于生成物的Helmholtz自由能，则反应不会自发进行。通过计算反应物和生成物的Helmholtz自由能，我们可以预测反应的进行方向和是否需要外界能量的输入。

二、Gibbs自由能

Gibbs自由能是描述一个恒温恒压系统的自由能，通常用符号G表示。它由焓H、熵S和温度T三个热力学函数组成，可以用以下公式表示：

G = H - TS

其中H是系统的焓，S是系统的熵。Gibbs自由能描述了系统在恒温恒压条件下可用的能量，也可以理解为系统可以做的最大有用功。

Gibbs自由能在化学和生物学等领域有着广泛的应用。首先，它可以用来预测化学反应的进行性质。根据Gibbs自由能的定义，如果反应物的Gibbs自由能大于生成物的Gibbs自由能，那么反应是自发进行的。反之，如果反应物的Gibbs自由能小于生成物的Gibbs自由能，则反应不会自发进行。通过计算反应物和生成物的Gibbs自由能，我们可以预测反应的进行性质和是否需要外界能量的输入。

mmpbsa结合自由能数值

MMPBSA（分子力学泛函面积）是一种计算自由能的方法，它通常用于研究生物分子的结合自由能。MMPBSA方法结合了分子力学模拟和溶剂化自由能计算，通过对分子间相互作用、溶剂化和构象变化等因素进行全面考虑，可以较为准确地预测分子结合自由能。

MMPBSA方法的数值计算通常包括以下步骤，首先，通过分子动力学模拟或分子对接等方法得到一系列构象；然后，对每个构象进行能量计算，包括分子内能、分子间相互作用能以及溶剂化能；最后，根据这些能量值计算结合自由能及其各项贡献。

在MMPBSA方法中，结合自由能的数值通常包括总的结合自由能以及其各项贡献，如范德华相互作用能、电荷-静电相互作用能、溶剂化自由能等。这些数值可以通过数值计算得到，通常以千焦耳/摩尔（kJ/mol）为单位。

需要注意的是，MMPBSA方法的计算结果受到多个因素的影响，包括模拟系统的选择、力场参数的准确性、计算方法的精度等。因此，在使用MMPBSA方法计算结合自由能时，需要对计算结果进行适当的验证和分析，以确保结果的可靠性和准确性。

总之，MMPBSA方法是一种用于计算生物分子结合自由能的有效方法，通过对分子间相互作用、溶剂化等因素进行全面考虑，可以得到较为准确的结合自由能数值。在实际应用中，需要综合考虑多个因素，并对计算结果进行适当的验证和分析。