系统生物学方法第三篇.pptx

1细胞网络细胞网络蛋白质相互作用蛋白质相互作用基因网络基因网络基因网络基因网络二维网络二维网络表达与调控表达与调控代谢代谢功能与行为功能与行为单个功能单个功能三三维维系系统统生生物物学学代谢、代谢组、代谢网络，代谢、代谢组、代谢网络，处于生命活动的末端，直处于生命活动的末端，直接导致生物功能与行为。接导致生物功能与行为。2 代谢处于生命活动的末端，是驱动生命活动的代谢处于生命活动的末端，是驱动生命活动的化学引擎。活细胞需要能量和物质来构建膜、贮存化学引擎。活细胞需要能量和物质来构建膜、贮存分子、补充酶、复制和修复分子、补充酶、复制和修复DNA、运动以及完成许、运动以及完成许多其他生理过程。细胞通过代谢获得这些能力并将多其他生理过程。细胞通过代谢获得这些能力并将其用于构建新的细胞。代谢是细胞生存和繁衍的手其用于构建新的细胞。代谢是细胞生存和繁衍的手段，它大体可以分为两大类：段，它大体可以分为两大类：分解代谢分解代谢（分解复杂化合物以获取（分解复杂化合物以获取能量能量和构建细和构建细胞所需胞所需模块模块）。）。合成代谢合成代谢（构建细胞功能所需的复杂化合物）。（构建细胞功能所需的复杂化合物）。代谢是一种高度有组织的过程，它涉及了由酶催化代谢是一种高度有组织的过程，它涉及了由酶催化的数以千计的反应。的数以千计的反应。第一章代谢系统建模第一章代谢系统建模3 代谢网络把细胞内所有生化反应表示为一个网代谢网络把细胞内所有生化反应表示为一个网络，反映了所有参与代谢过程的化合物之间以及所络，反映了所有参与代谢过程的化合物之间以及所有催化酶之间的相互作用，是细胞代谢的数学模型。有催化酶之间的相互作用，是细胞代谢的数学模型。代谢网络是由将分子从一种形态转换为另一种代谢网络是由将分子从一种形态转换为另一种形态的反应所组成的。在建模关系中，分子的形态的反应所组成的。在建模关系中，分子的浓度浓度及其及其转化速率转化速率尤为重要。尤为重要。代谢反应网络的基本概念和方法同样可以应用代谢反应网络的基本概念和方法同样可以应用于其他类型的细胞反应网络，如信号转导通路和转于其他类型的细胞反应网络，如信号转导通路和转录因子对录因子对DNA的结合的结合等。4本章将在本章将在3个层次上对代谢进行介绍：个层次上对代谢进行介绍：（1）酶动力学研究）酶动力学研究孤立体系孤立体系中单个反应的动力学中单个反应的动力学性质。性质。（2）代谢的网络特性则需要用到对化合物进行生）代谢的网络特性则需要用到对化合物进行生成和降解平衡的成和降解平衡的化学计量分析化学计量分析。（3）代谢控制分析则通过分析个体浓度变化动力）代谢控制分析则通过分析个体浓度变化动力学并将其整合到学并将其整合到网络网络中，定量描述扰动对网络的影中，定量描述扰动对网络的影响。响。5一个例子一个例子 糖酵解代谢糖酵解代谢糖酵解中的最初糖酵解中的最初4步步，以及，以及平衡能量流平衡能量流ATP和和ADP的反应的反应：表示糖酵解的几步上游反应，糖酵解即葡萄糖的降解用以产生能量和细胞表示糖酵解的几步上游反应，糖酵解即葡萄糖的降解用以产生能量和细胞进程所需的构建模块。缩写：进程所需的构建模块。缩写：Gluc6P(葡萄糖葡萄糖-6-磷酸磷酸)，Fruc6P（果糖（果糖-6-磷酸），磷酸），Fruc1，6P2（果糖（果糖-1，6-二磷酸），二磷酸），ATP（三磷酸腺苷），（三磷酸腺苷），ADP（二磷酸腺苷），（二磷酸腺苷），AMP(单磷酸腺苷单磷酸腺苷)。反应。反应v1(己糖激酶己糖激酶)，v2(葡萄糖葡萄糖-6-磷酸在其他通路中的消耗磷酸在其他通路中的消耗)，v3(磷酸葡萄糖异构酶磷酸葡萄糖异构酶)，v4(磷酸果糖激酶磷酸果糖激酶)，v5(醛缩酶醛缩酶)，v6(后继糖酵解中的后继糖酵解中的ATP生成生成)，v7(ATP在其他通路中的生在其他通路中的生成成)，v8(腺苷酸激酶腺苷酸激酶)6这个反应体系的常微分方程（这个反应体系的常微分方程（ODE）组如下：）组如下：7速率表达式如下：速率表达式如下：8如何如何建模建模，速率常数速率常数如何确定是我们要研究的。如何确定是我们要研究的。动力系统会讨论微分方程组的时间发展与参数动力系统会讨论微分方程组的时间发展与参数的依赖性，或敏感性的依赖性，或敏感性结构稳定性。包括稳态吸结构稳定性。包括稳态吸引子、周期吸引子和混沌吸引子，以及对应的吸引引子、周期吸引子和混沌吸引子，以及对应的吸引域。吸引域的边界往往是分形几何。域。吸引域的边界往往是分形几何。9 第一节第一节 酶动力学和热力学酶动力学和热力学 本章将讲解象糖酵解本章将讲解象糖酵解上游模型上游模型这样的代这样的代谢模型如何被描述和分析的。将介绍对速率谢模型如何被描述和分析的。将介绍对速率方程各种可能的假设，并给出对这些方程各种可能的假设，并给出对这些假设到假设到速率的推理速率的推理。特别值得注意的是：。特别值得注意的是：信息包含于网络中是否有几条主线？哪些流信息包含于网络中是否有几条主线？哪些流可能出于稳态？底物是否被产生或被消耗，可能出于稳态？底物是否被产生或被消耗，或是否有保持浓度守恒关系的代谢物？代谢或是否有保持浓度守恒关系的代谢物？代谢控制分析使我们能够估算网络中一部分的参控制分析使我们能够估算网络中一部分的参量改变将对网络任何部分的变量（这里指稳量改变将对网络任何部分的变量（这里指稳态流或浓度）变化的影响？等等问题。态流或浓度）变化的影响？等等问题。10 本节讨论单个生化反应的确定性动力学建模。本节讨论单个生化反应的确定性动力学建模。最基本的量是底物最基本的量是底物S的浓度的浓度S（即单位体积（即单位体积V中分子中分子的个数的个数n）和反应速率）和反应速率v（即单位时间（即单位时间t中浓度中浓度S的改的改变）。相对于考察单分子及其相互作用的微观方变）。相对于考察单分子及其相互作用的微观方法而言，这种建模方法是宏观和表象的。法而言，这种建模方法是宏观和表象的。化学和生化动力学依赖于如下假设：化学和生化动力学依赖于如下假设：在特定时在特定时空点的反应速率空点的反应速率v能够被该时空点中所有底物浓度的能够被该时空点中所有底物浓度的唯一函数所表达唯一函数所表达。经典酶动力学为了简化处理，假。经典酶动力学为了简化处理，假设反应速率的空间分布式均一的（设反应速率的空间分布式均一的（“一室一室”试管），试管），而且不直接依赖于时间：而且不直接依赖于时间：11 催化生化反应的催化生化反应的酶是蛋白质酶是蛋白质，通常会和，通常会和辅因子辅因子复合。酶具有一个通常是高度特异的催化中心，并复合。酶具有一个通常是高度特异的催化中心，并在反应中保持不变。一个酶分子每秒大约催化一千在反应中保持不变。一个酶分子每秒大约催化一千个反应（所谓的转换数范围从）。这导致反应速率个反应（所谓的转换数范围从）。这导致反应速率比在非催化的自发反应下提高了比在非催化的自发反应下提高了1061012倍。倍。在朝全细胞建模方向发展的更高级建模方法中，在朝全细胞建模方向发展的更高级建模方法中，空间的不均一性会被考虑进来，如事实上存在这样空间的不均一性会被考虑进来，如事实上存在这样的事实：许多组分是与膜结合的，或细胞结构阻止的事实：许多组分是与膜结合的，或细胞结构阻止了分子的自由移动。然而，在绝大多数情况下可以了分子的自由移动。然而，在绝大多数情况下可以假设扩散足够迅速，使得所有底物能够在空间中均假设扩散足够迅速，使得所有底物能够在空间中均匀分布。匀分布。121、质量作用定律、质量作用定律 生化动力学是基于质量作用定律的，该定律由生化动力学是基于质量作用定律的，该定律由Guldberg和和Waage于于19世纪引入。该定律认为世纪引入。该定律认为反应反应速率与反应物的碰撞概率成比例速率与反应物的碰撞概率成比例，而碰撞概率又与，而碰撞概率又与以相应反应分子数为幂的反应浓度成比例。反应分以相应反应分子数为幂的反应浓度成比例。反应分子数即该反应物参加特定反应的（最简约化）数目。子数即该反应物参加特定反应的（最简约化）数目。例如，对于一个简单的反应如：例如，对于一个简单的反应如：反应速率为：反应速率为：V是净速率，是净速率，是正反应速率，是正反应速率，是逆反应速率，是逆反应速率，和和 分别是相应的比例因子，即所谓的分别是相应的比例因子，即所谓的动力学或动力学或速率常数速率常数。13底物的反应分子数对于正反应而言均为底物的反应分子数对于正反应而言均为1，对于逆反，对于逆反应则为应则为2.如果我们测量每升的摩尔数（如果我们测量每升的摩尔数（或或M）作为浓度，时间以秒（）作为浓度，时间以秒（s）计，这样速率单位为）计，这样速率单位为 。因此，双分子反应的速率常数单位应该是。因此，双分子反应的速率常数单位应该是 ，而单分子反应的速率常数量纲则为，而单分子反应的速率常数量纲则为 。对于底。对于底物浓度为物浓度为 Si，产物浓度为，产物浓度为 的反应，基于质量作的反应，基于质量作用定律的一般速率表达式为：用定律的一般速率表达式为：这里这里 和和 分别代表该反应中分别代表该反应中 和和 的反应的反应分子数（分子数（Heinrich和和Schuster 1966）。）。14平衡常数平衡常数 （我们还将使用更简单的符号（我们还将使用更简单的符号q）表）表征平衡态下产物和底物的浓度比（征平衡态下产物和底物的浓度比（和和 ），），平衡态即反应与逆反应速率相等时的状态。速率常平衡态即反应与逆反应速率相等时的状态。速率常数以下述方式与数以下述方式与 的关联：的关联：上述反应式的浓度动力学可由常微分方程组描述：上述反应式的浓度动力学可由常微分方程组描述：、和和P的时间进程可通过对上述常微分方程组进行的时间进程可通过对上述常微分方程组进行积分获得。积分获得。底物减少，产物底物减少，产物增加！增加！15例、简单分解的动力学如：例、简单分解的动力学如：可以用：可以用：描述。描述。对此常微分方程从时间对此常微分方程从时间 开始以初始浓度开始以初始浓度 积分到任何时刻积分到任何时刻 t 的浓度为的浓度为：得到随时间的表达式得到随时间的表达式：162、反应的动力学和热力学反应的动力学和热力学 代谢的一个重要目的就是从营养物中取得用于代谢的一个重要目的就是从营养物中取得用于合成生长和增殖所需分子的能量。我们将反应区分合成生长和增殖所需分子的能量。我们将反应区分为为供能反应供能反应、需能反应需能反应和和近零反应近零反应3种。可逆热力种。可逆热力学原理及其在化学反应上的应用使我们能够理解细学原理及其在化学反应上的应用使我们能够理解细胞中的能量循环。生物反应通常在胞中的能量循环。生物反应通常在恒温恒压和近恒恒温恒压和近恒容容水溶液环境中进行这一假设将使问题简化。水溶液环境中进行这一假设将使问题简化。无论反应是否自发的，反应进行的方向和平衡的无论反应是否自发的，反应进行的方向和平衡的位置对于生命过程而言都是重要的特性。热力学第一位置对于生命过程而言都是重要的特性。热力学第一定律即能量守恒定律告诉我们在任何孤立系统中的总定律即能量守恒定律告诉我们在任何孤立系统中的总能量保持为常数。热力学第二定律则宣称只有当系统能量保持为常数。热力学第二定律则宣称只有当系统的总熵增加时，过程才是自发的。然而熵通常是无法的总熵增加时，过程才是自发的。然而熵通常是无法直接测量到的，一个更适合的度量是吉布斯自由能直接测量到的，一个更适合的度量是吉布斯自由能G，它是指体，它是指体系统在恒温恒压条件下对外做功的能力系统在恒温恒压条件下对外做功的能力。17这个条件下，吉布斯自由能变为：这个条件下，吉布斯自由能变为：这里这里 为焓变，为焓变，为熵变，为熵变，T为开氏绝对温度。为开氏绝对温度。是对化学反应自发性和驱动力的度量是对化学反应自发性和驱动力的度量。如果。如果 则反应自发进行并伴随释放能量（放能过程）则反应自发进行并伴随释放能量（放能过程）；如果；如果 则反应在能量上是不利的因而不会则反应在能量上是不利的因而不会自发进行（吸能过程）；自发进行（吸能过程）；意味着系统达到了意味着系统达到了平衡。如果可以通过能量耦合从严格的放能反应中平衡。如果可以通过能量耦合从严格的放能反应中获取到能量，吸能反应也是可以发生的。获取到能量，吸能反应也是可以发生的。反应的自由能变可以通过产物自由能和与底物自由反应的自由能变可以通过产物自由能和与底物自由能之间的差值计算出来：能之间的差值计算出来：18酶不能改变反应中产物和底物的自由能，同样不能酶不能改变反应中产物和底物的自由能，同样不能改变其差，但是酶可以改变反应路径，从而降低反改变其差，但是酶可以改变反应路径，从而降低反应需要的活化能。过度态理论（应需要的活化能。过度态理论（2001，Haynie）阐）阐明：由于明：由于G0，许多物质或混合物是热力学不，许多物质或混合物是热力学不稳定的。然而它们在正常状态下可以储存很长时间，稳定的。然而它们在正常状态下可以储存很长时间，原因是在反应过程中，代谢物必须经过一个或更多原因是在反应过程中，代谢物必须经过一个或更多的具有极大自由能的过度态，在这种状态下，化学的具有极大自由能的过度态，在这种状态下，化学键会被解开或重新生成。过度态是不稳定的，其对键会被解开或重新生成。过度态是不稳定的，其对应的分子构型被称为活化络合物。反应物与活化络应的分子构型被称为活化络合物。反应物与活化络合物之间的自由能差合物之间的自由能差G决定反应的动力学：这个决定反应的动力学：这个能差越高，分子越过该能垒概率越低，反应速率也能差越高，分子越过该能垒概率越低，反应速率也就越小。就越小。（表给出一些反应的自由能差）（表给出一些反应的自由能差）19图展示了反应过程的简化图景。底物和产物均处于的局图展示了反应过程的简化图景。底物和产物均处于的局部极小自由能状态，活化络合物则被赋予了局部最大自部极小自由能状态，活化络合物则被赋予了局部最大自由能。自由能变化与各反应中平衡常数的对数成比例：由能。自由能变化与各反应中平衡常数的对数成比例：20G值对应于正反应动力学常数值对应于正反应动力学常数 k+，G=-RT*ln k+。而。而G+G则与逆反则与逆反应速率应速率k_相关。相关。反应物与酶的相互作用可能会改变反应路径，反应物与酶的相互作用可能会改变反应路径，从而导致更低的从而导致更低的G值。更进一步的，在值。更进一步的，在反应路径中可能呈现出更多的局部极小和极反应路径中可能呈现出更多的局部极小和极大自由能状态，分别对应不稳定的中间体复大自由能状态，分别对应不稳定的中间体复合物。一些生物学上重要反应的自由能变的合物。一些生物学上重要反应的自由能变的值在书上表值在书上表5.1中给出。中给出。213、米氏动力学、米氏动力学 Brown（1920）为转移酶反应提出了首个酶促）为转移酶反应提出了首个酶促反应机理，该反应全部是单底物的，没有逆反应和反应机理，该反应全部是单底物的，没有逆反应和效应物：效应物：上式包含从游离酶上式包含从游离酶E与底物与底物S结合形成酶底复合物结合形成酶底复合物ES的可逆过程，以及产物的可逆过程，以及产物P从酶从酶E上不可逆释放过程。上不可逆释放过程。这一反应的动力学的相应常微分方程组这一反应的动力学的相应常微分方程组如下：如下：22总反应速率与底物的消耗速率的负值以及产物总反应速率与底物的消耗速率的负值以及产物的生成速率相同：的生成速率相同：这个常微分方程组无法得到解析解。通过假设达到这个常微分方程组无法得到解析解。通过假设达到令人满意的简化。令人满意的简化。Michaelis、Menten(1913)假定假定E和和S转换为转换为ES以及其逆过程远远快于以及其逆过程远远快于ES分解为分解为E和和P的过程（所谓的游离酶与酶的过程（所谓的游离酶与酶-底复合物间的准平衡），底复合物间的准平衡），或者这或者这3个常数间有关系：个常数间有关系：23 Briggs和和Haldane（1925）则假定在反应过程）则假定在反应过程中会达到一个中会达到一个ES复合物浓度保持不变的状态。这复合物浓度保持不变的状态。这一假设在底物初始浓度远高于酶浓度时是合理的，一假设在底物初始浓度远高于酶浓度时是合理的，即：如果不满足这个条件，这种假设中的稳态是即：如果不满足这个条件，这种假设中的稳态是永远达不到的。他们提出的永远达不到的。他们提出的ES复合物准稳态一般复合物准稳态一般性假设可表示为：性假设可表示为：24 反应速率的表达式反应速率的表达式可以从方程可以从方程的处理的处理以及以及ES的准稳态假设中导出：的准稳态假设中导出：在该反应中，酶即不会被产生也不会被消耗，它们可在该反应中，酶即不会被产生也不会被消耗，它们可能以游离态或复合物形式存在，但总浓度保持不变。能以游离态或复合物形式存在，但总浓度保持不变。将上面稳态假设方程代入方程组将上面稳态假设方程代入方程组得：得：25求出反应速率，得：求出反应速率，得：在酶动力学中通常是将在酶动力学中通常是将方程以更简化的形式表方程以更简化的形式表示，该形式在理论和实示，该形式在理论和实际上都是很重要的：际上都是很重要的：此方程即为此方程即为Michaelis-Menten动力学表达式。各参量动力学表达式。各参量解释如下。最大速率是酶在被底物饱和时可以达到的解释如下。最大速率是酶在被底物饱和时可以达到的最大酶促反应速率：最大酶促反应速率：米氏常数等于酶促反应速率达到最大速率一半时所米氏常数等于酶促反应速率达到最大速率一半时所需的底物浓度：需的底物浓度：26对于准平衡假设，有对于准平衡假设，有 。各参量的意义。各参量的意义见速率对底物浓度图，该图像具有双曲线形状。见速率对底物浓度图，该图像具有双曲线形状。M-M动力学中反动力学中反应速率应速率v对底物浓对底物浓度度S的关系图。的关系图。表示高底物浓表示高底物浓度下能够达到的最度下能够达到的最大反应速率。大反应速率。当底物浓度很低时，当底物浓度很低时，v几乎随着几乎随着S线性地增加；而当线性地增加；而当底物浓度很高时，底物浓度很高时，v几乎与几乎与S无关无关。27推导速率方程的一般方法小结：推导速率方程的一般方法小结：1）将所有需要考虑的步骤都绘制出来。包括所有的）将所有需要考虑的步骤都绘制出来。包括所有的底物和产物（底物和产物（S和和P）以及）以及n种游离的和结合态的种游离的和结合态的酶（酶（E和和ES）。）。2）将所有生成和消耗该物质步骤的速率都加合起来，）将所有生成和消耗该物质步骤的速率都加合起来，置于浓度变化常微分方程的右边。反应速率遵循置于浓度变化常微分方程的右边。反应速率遵循质量作用定律。质量作用定律。3）所有含有酶的物质的量加起来应该等于总酶浓度）所有含有酶的物质的量加起来应该等于总酶浓度Etotal（含有酶的物质的所有微分方程右边加合起（含有酶的物质的所有微分方程右边加合起来应该等于来应该等于0）。这个关系会构成一个方程。）。这个关系会构成一个方程。4）对）对n-1种含酶物质准稳态假设（即各自常微分方种含酶物质准稳态假设（即各自常微分方程右边都为程右边都为0）与第）与第3步的结果合起来可以构成关步的结果合起来可以构成关于这于这n种含酶物质浓度的种含酶物质浓度的n个代数方程。个代数方程。5）反应速率等于产物生成的速率，从第）反应速率等于产物生成的速率，从第4步的结果步的结果中将各种含酶物质浓度代入即可。中将各种含酶物质浓度代入即可。28米氏方程的参数估计和线性化米氏方程的参数估计和线性化：为了估计一种单独的酶的参量为了估计一种单独的酶的参量 和和 值，值，需要测定不同初始底物浓度下的初始反应速率。由于需要测定不同初始底物浓度下的初始反应速率。由于速率是底物浓度的非线性函数，一种方法是通过非线速率是底物浓度的非线性函数，一种方法是通过非线性回归来求解参量。另一种方法则是将速率方程变换性回归来求解参量。另一种方法则是将速率方程变换为线性变量关系，然后进行线性回归。为线性变量关系，然后进行线性回归。变换后的方程的一个优点是或多或少能够从测量数变换后的方程的一个优点是或多或少能够从测量数据线性回归得到的图中直接读出参量值，在据线性回归得到的图中直接读出参量值，在Lineweaver-Burk图中，图中，和和 的值能够分别从的值能够分别从图像与纵、横坐标的交点获得。图像与纵、横坐标的交点获得。Linewaver-Burk图图同样有助于区分不同类型的抑制同样有助于区分不同类型的抑制。变换后的方程的。变换后的方程的不足在于当底物浓度或反应速率较低或较高时对误不足在于当底物浓度或反应速率较低或较高时对误差十分敏感。差十分敏感。Eadie和和Hofstee 以及以及Hanes 和和 Woolf 介绍了其他种类的线性化方法以克服这一局限。介绍了其他种类的线性化方法以克服这一局限。29Lineweaver-BurkEadie-HofsteeHanes-Woolf变换后的方程新变量图解展示30 可逆反应的米氏方程可逆反应的米氏方程 在现实中许多反应式可逆的，酶可以在双向在现实中许多反应式可逆的，酶可以在双向上催化反应。考察下列机理：上催化反应。考察下列机理：产物生成速率：产物生成速率：相应的速率方程相应的速率方程：31参量参量 和和 是各步反应的动力学常数；表观参是各步反应的动力学常数；表观参量量 和和 分别表示零产物浓度下正方向和零底分别表示零产物浓度下正方向和零底物下逆方向的最大反应速率；表观参量物下逆方向的最大反应速率；表观参量 和和 分别表示使得正反应或逆反应为各自最大速率一分别表示使得正反应或逆反应为各自最大速率一半时，底物或产物的浓度。它们之间的关系有半时，底物或产物的浓度。它们之间的关系有：例一中例一中 所描述的反应即是一个所描述的反应即是一个Michaelis-Menten型的可逆反应。（试着推导一下！）型的可逆反应。（试着推导一下！）324、基于蛋白质相互作用的酶活性调节、基于蛋白质相互作用的酶活性调节 酶可以极大地提高反应速率，但这并非它们的酶可以极大地提高反应速率，但这并非它们的全部功能。酶通过多种方式参与代谢调节，它们的全部功能。酶通过多种方式参与代谢调节，它们的生成与降解通常是与当前的细胞需求相适应的，更生成与降解通常是与当前的细胞需求相适应的，更进一步的，它们可能成为一些效应物的靶标，这些进一步的，它们可能成为一些效应物的靶标，这些效应物既可能是抑制剂也可能是活化剂。效应物既可能是抑制剂也可能是活化剂。这些效应物是能够影响酶促反应性能的蛋白质这些效应物是能够影响酶促反应性能的蛋白质或其他分子。效应物与酶之间的相互作用会改变反或其他分子。效应物与酶之间的相互作用会改变反应速率。我们将在这里考察这类对于代谢的精细调应速率。我们将在这里考察这类对于代谢的精细调控至关重要的调节性的相互作用。控至关重要的调节性的相互作用。33 依据效应物与之结合的酶的状态（即与游离依据效应物与之结合的酶的状态（即与游离酶酶E或酶底复合物或酶底复合物ES结合，还是与两者都结合），结合，还是与两者都结合），以及不同复合物释放产物的能力，可以将抑制作以及不同复合物释放产物的能力，可以将抑制作用分为几个基本类型。通常的抑制作用模式图解用分为几个基本类型。通常的抑制作用模式图解展示在下图中。如果其中一些反应不发生的话，展示在下图中。如果其中一些反应不发生的话，就会导致不同的抑制类型。就会导致不同的抑制类型。E+Sk1k-1ESk2E+Pk3k-3+IEI+Sk5k-5k4k-4+IESIk6E+P+I34速率方程的推导依照下列步骤进行。速率方程的推导依照下列步骤进行。（1）考虑化合物及其复合物间的结合平衡）考虑化合物及其复合物间的结合平衡:注意如果所有反应均能发生的话，则存在如下形式注意如果所有反应均能发生的话，则存在如下形式的的Wegscheider规则规则:这意味着这意味着E和和ESI之间的自由能差取决于反应选择之间的自由能差取决于反应选择的路径（通过的路径（通过ES或或EI）。）。两条路径的产物两条路径的产物ESI相同！相同！35（2）考虑到总酶量的部分守恒（只包括在反应中）考虑到总酶量的部分守恒（只包括在反应中出现的那些复合物）：出现的那些复合物）：(3)反应速率等于产物生成速率：反应速率等于产物生成速率：上述方程包含了关于上述方程包含了关于E、ES、EI和和ESI这这4种未知浓种未知浓度的度的4个独立的方程。将它们的解嵌入到速率方程个独立的方程。将它们的解嵌入到速率方程中。抑制剂的作用依赖于底物和抑制剂的浓度，以中。抑制剂的作用依赖于底物和抑制剂的浓度，以及各自与酶的亲合力。及各自与酶的亲合力。下表列出了下表列出了Michaelis-Menten动力学中各种动力学中各种可逆和不可逆抑制类型及其相应的速率方程。可逆和不可逆抑制类型及其相应的速率方程。36实现实现不可逆速率不可逆速率可逆速率方程可逆速率方程特性特性竞竞争争性性抑抑制制I只与游离的E结合；P只从ES复合物中释放；改变，保持；S与I竞争结合位点；高S可竞争过I 反反竞竞争争性性抑抑制制 I只与只与ES复合物结合；复合物结合；P只从只从ES复合物中释复合物中释放；放；和 均改变，但二者比率保持；S不会竞争过I 非竞争性抑制I与E和ES均结合；P只从ES中释放；保持，改变；S不会竞争过I 37实现实现不可逆速率不可逆速率可逆速率方程可逆速率方程特性特性混混合合抑抑制制 I与E和ES均结合；P只从ES中释放；和均改变；：竞争非竞争性抑制；：反竞争非竞争性抑制。部部分分抑抑制制 I可能与E和ES均结合；P从ES和ESI中释放；和均改变；则表现活化作用而非抑制作用。缩写：缩写：385、抑制剂对酶的不可逆结合抑制、抑制剂对酶的不可逆结合抑制不可逆抑制剂不可逆地与酶活性位点结合：不可逆抑制剂不可逆地与酶活性位点结合：这样就阻止了正常底物的结合并破坏了催化中心，这样就阻止了正常底物的结合并破坏了催化中心，还会导致酶的变性。在任何情况下，初始浓度为还会导致酶的变性。在任何情况下，初始浓度为I0 的的抑制剂会将酶的有效浓度从初始的浓度抑制剂会将酶的有效浓度从初始的浓度E0降低到降低到E0-I0。抑制剂的摩尔数过量会导致催化活性的完全丧。抑制剂的摩尔数过量会导致催化活性的完全丧失；而酶的摩尔数过量时，酶促反应最大速率减为失；而酶的摩尔数过量时，酶促反应最大速率减为:396、底物抑制、底物抑制 酶促反应的一个普遍特性是反映速率随着底物酶促反应的一个普遍特性是反映速率随着底物浓度浓度S增加，直至逼近最大反应速率。但在一些情增加，直至逼近最大反应速率。但在一些情况下，在超过一定况下，在超过一定S值后反应速率反而会降低。这值后反应速率反而会降低。这可能的一个原因是：酶可能的一个原因是：酶-底物复合物会进一步与底底物复合物会进一步与底物分子结合形成物分子结合形成ESS复合物，而复合物不能生成产复合物，而复合物不能生成产物。如果第二个底物能够释放，则这种抑制作用是物。如果第二个底物能够释放，则这种抑制作用是可逆的。将抑制剂用第二个结合的底物替换，速率可逆的。将抑制剂用第二个结合的底物替换，速率方程则可以从反竞争性抑制模式中推导出来，得：方程则可以从反竞争性抑制模式中推导出来，得：40这个表达式有最大值，当：这个表达式有最大值，当：V对对S的依赖现实在图中。底物抑制一个典型的例子的依赖现实在图中。底物抑制一个典型的例子就是两个琥珀酸分子结合于丙二酸脱氢酶上，该酶具就是两个琥珀酸分子结合于丙二酸脱氢酶上，该酶具有两个对羟基的结合袋。这个也图解展示于图中。有两个对羟基的结合袋。这个也图解展示于图中。41在底物抑制情况下反应速率在底物抑制情况下反应速率v对底物浓度对底物浓度S的图像。上面的曲的图像。上面的曲线显示了没有抑制时的线显示了没有抑制时的Michaelis-Menten动力学，下面的两动力学，下面的两条曲线显示了在结合常数下的动力学特性。参量值：。右边条曲线显示了在结合常数下的动力学特性。参量值：。右边的图显示了底物抑制的可能机理。酶（灰色物）具有对底物的图显示了底物抑制的可能机理。酶（灰色物）具有对底物分子不同部分结合的两个结合口袋（上图）。在高底物浓度分子不同部分结合的两个结合口袋（上图）。在高底物浓度时，可能会有两个不同的分子进入结合袋，因而阻止了特定时，可能会有两个不同的分子进入结合袋，因而阻止了特定的反应（下图）。的反应（下图）。CO2CH2CO2-CH2CO2CH2CO2-CH2CO2CH2=CH2CO2CO2CH2=CH2CO2427、抑制剂与底物结合导致的抑制、抑制剂与底物结合导致的抑制 如果抑制剂与底物形成紧密复合物，反应速如果抑制剂与底物形成紧密复合物，反应速率同样会降低：率同样会降低：有效的底物浓度会下降至：有效的底物浓度会下降至：依据质量作用定律，实际平衡浓度有关系：依据质量作用定律，实际平衡浓度有关系：有效酶浓度的计算：有效酶浓度的计算：在高底物浓度下。存在抑制剂时反应的最大速率可在高底物浓度下。存在抑制剂时反应的最大速率可以达到与无抑制剂时一样。该情况下的以达到与无抑制剂时一样。该情况下的Lineweaver-Burk图是非线性的。图是非线性的。438、配基于蛋白质的结合、配基于蛋白质的结合 不考虑是否是反应中的主体，凡是与蛋白质结合的不考虑是否是反应中的主体，凡是与蛋白质结合的分子都被称为分子都被称为配基配基。下面我们考虑配基与单体蛋白。下面我们考虑配基与单体蛋白质和寡聚蛋白质相结合的情形。在寡聚蛋白质中，质和寡聚蛋白质相结合的情形。在寡聚蛋白质中，可能存在着不同亚基中结合位点间的相互作用。可能存在着不同亚基中结合位点间的相互作用。考虑到一个配基于蛋白质只会有一个结合位点：考虑到一个配基于蛋白质只会有一个结合位点：结合常数结合常数由下式给出：由下式给出：的倒数的倒数 为离解常数为离解常数，蛋白质的饱和度分数，蛋白质的饱和度分数Y为结合了配基的亚基数除以总的亚基数。单亚基蛋为结合了配基的亚基数除以总的亚基数。单亚基蛋白质的饱和分数为：白质的饱和分数为：44类似于类似于Michaelis-Menten动力学中动力学中v对对S的图像，在的图像，在总酶浓度保持不变的情况下总酶浓度保持不变的情况下Y对对S的图像为一双曲线。的图像为一双曲线。在整个过程中，在整个过程中，S结合到结合到E上是第一步，接下来产物上是第一步，接下来产物释放。因此当释放。因此当S初始浓度远高于初始浓度远高于E初始浓度时，反应初始浓度时，反应速率与速率与ES浓度成正比，有：浓度成正比，有：如果蛋白质有几个结合位点，这些位点间可能发生如果蛋白质有几个结合位点，这些位点间可能发生相互作用，即在结合了一个或多个配基之后，随后相互作用，即在结合了一个或多个配基之后，随后上去的配基的亲合力的可能会有变化。这个现象被上去的配基的亲合力的可能会有变化。这个现象被称为协同。正、负协同分别表示蛋白质对后结合上称为协同。正、负协同分别表示蛋白质对后结合上的配基亲合力的上升和下降。同位协同和异位协同的配基亲合力的上升和下降。同位协同和异位协同分别表示施加影响的先结合配基与被影响的后结合分别表示施加影响的先结合配基与被影响的后结合配基是相同和不同的。配基是相同和不同的。459、同位正协同与希尔方程、同位正协同与希尔方程 考察一个具有两个相同结合位点的二聚蛋白质，考察一个具有两个相同结合位点的二聚蛋白质，第一个配基的结合促进第二个配基的结合：第一个配基的结合促进第二个配基的结合：这里这里E是单聚体是单聚体而E2是二聚体。饱和度分数为：二聚体。饱和度分数为：如果第二个配基与蛋白质的亲合力随第一个配基的如果第二个配基与蛋白质的亲合力随第一个配基的结合强烈增加，这时结合强烈增加，这时 一旦生成将与一旦生成将与S马上结合，马上结合，因而因而 的浓度可以被忽略的浓度可以被忽略.46在完全协同中，所有的蛋白质不是出于无配基结合在完全协同中，所有的蛋白质不是出于无配基结合就是处于与配基完全结合状态，反应式化简为：就是处于与配基完全结合状态，反应式化简为：结合常数为：结合常数为：因此饱和度分数为：因此饱和度分数为：通常对于具有通常对于具有n个亚基的蛋白质有：个亚基的蛋白质有：此为此为Hill方程的一般形式，该式意味着完全的同位协方程的一般形式，该式意味着完全的同位协同。饱和度分数同。饱和度分数Y对底物浓度对底物浓度S作图可获得一作图可获得一S形曲形曲线，该曲线的拐点在线，该曲线的拐点在 处。量处。量n（通常用（通常用“h”代替）被称为代替）被称为Hill系数系数.47该表达式的推导是基于氧与血红蛋白（该表达式的推导是基于氧与血红蛋白（Hb）结合）结合实验的结果（实验的结果（Hill1910、1913）。）。Bohr及其合作者及其合作者在在1904年发现携氧年发现携氧Hb的饱和度分数对氧分压作图的饱和度分数对氧分压作图呈现出呈现出S形。形。Hill（1909）用位于不同亚基上的结）用位于不同亚基上的结合位点间的相互作用解释了这一现象。这时已经合位点间的相互作用解释了这一现象。这时已经知道每个亚基只结合一个氧分子。知道每个亚基只结合一个氧分子。Hill假定过程为假定过程为完全协同的，因此预测出实验完全协同的，因此预测出实验Hill系数为系数为2.8。今天。今天我们知道血红蛋白有我们知道血红蛋白有4个结合位点，但它们间的协个结合位点，但它们间的协同是非完全的。同是非完全的。S形结合特征具有优势：在很高的形结合特征具有优势：在很高的养分压下，血红蛋白与肺中的氧结合；当养分压养分压下，血红蛋白与肺中的氧结合；当养分压很低时，则很容易将释放到体内。很低时，则很容易将释放到体内。4810、具有、具有S形动力学的酶的形动力学的酶的Monod-Wynma-Changeux速率表达式速率表达式 1965年年Monod及其同事提出了一个解释及其同事提出了一个解释S形酶动力形酶动力学的模型，该模型考虑到了酶亚基间的相互作用学的模型，该模型考虑到了酶亚基间的相互作用（Monld等等1965）。更全面的模型由）。更全面的模型由Koshland等等（1966）提出。）提出。Monod等的模型采用了下列的假设：等的模型采用了下列的假设：酶包含酶包含n个个相同的亚基；相同的亚基；每个亚基可以假定具有活性（每个亚基可以假定具有活性（R）和非活性（和非活性（T）两种构象；）两种构象；所有的亚基在同一时所有的亚基在同一时刻变换它们的构象（一致变换）；刻变换它们的构象（一致变换）；R和和T构象间的构象间的平衡由变构常数给出：平衡由变构常数给出：49Ti和和Ri的指数的指数i表示底物分子的范围数。活性和非活性表示底物分子的范围数。活性和非活性构象下的结合常数分别由构象下的结合常数分别由KR和和KT给出。如果底物分给出。如果底物分子只能与活性形式的亚基结合，即子只能与活性形式的亚基结合，即 ，那么反应，那么反应速率可表示为：速率可表示为：这里因子这里因子对应于对应于M-M速率表达式，而速率表达式，而 则为调节因子。则为调节因子。50对于对于L=0，v对对S的图形是一条的图形是一条Michaelis-Menten动动力学的双曲线。对于力学的双曲线。对于L0，曲线向右漂移从而得，曲线向右漂移从而得到到S形曲线。典型的边构常数值形曲线。典型的边构常数值:当底物也能与非活性状态的亚基结合时（当底物也能与非活性状态的亚基结合时（），反应速率为：），反应速率为：51已经考虑了已经考虑了Monod-Wyman-Changeux模型中的同模型中的同位正效应。但这一模型同样很适于解释反应速率与位正效应。但这一模型同样很适于解释反应速率与活化剂和抑制剂间的相关性。活化剂活化剂和抑制剂间的相关性。活化剂A仅仅与活性仅仅与活性构象结合，而抑制剂构象结合，而抑制剂I仅仅与肺活性构象结合。这仅仅与肺活性构象结合。这使得平衡向相应的构象移动。实际上，效应物的结使得平衡向相应的构象移动。实际上，效应物的结合会改变合会改变L：KI和和KA表示结合常数。与效应物的相互作用是同位表示结合常数。与效应物的相互作用是同位效应。如图效应。如图5.7所示，活化剂会减弱曲线的所示，