论文部分内容阅读
近十年来，系统生物学的概念和理论得到了迅速的发展，系统生物学的研究如火如荼地开展。系统生物学的核心是要整合各个层面，各个构成要素的知识。这些层面当中，除开耳熟能详的基因组，转录组，蛋白质组和代谢组外，通量组(Fluxome)成为越来越重要一个组成部分。
对这个体系的研究产生了新的组学，即代谢通量组学(Fluxomics)。其重点是对目标体系内的代谢通量的速度的方向和取值进行测定。所获得的信息对刻画生物系统的生理特性、评估遗传和环境对细胞的影响具有重要意义。
代谢通量组学的主要研究手段是代谢通量分析，它可以分为纯计量性代谢通量分析和13C代谢通量分析两类。其中13C代谢通量分析是近十多年来，国际上新发展并仍然处于发展当中的技术。它除开具有纯计量性代谢通量分析的能力外，还具有获取信息更多，适用范围更广的优点。它能获得代谢途径的交换通量，并且不依赖于除稳态假设外的其他生理假设，因而适用各种试验条件。但是对于13C代谢通量分析技术，国内还未见有博士论文专题研究。
本博士论文的内容主要分成二部分，一部分描述在中国科技大学建立13C代谢通量分析平台的工作，另一部分是对通量分析平台中的13C质量同位素异构体检测技术进行改进的工作。
第一章详细介绍了13C代谢通量分析的基本理论。
代谢通量分析的基础是代谢通量拟稳态假设。在此假设基础上，可以建立各个代谢物的物质平衡方程。13C代谢通量分析实验在培养底物里引入不完全的13C标记。13C会分布到所有代谢物中，而由于13C在代谢物中数量和位置的不同，可以得到同一种物质的不同的13C同位素异构体和质量同位素异构体(同位素异构体的一类特殊组合形式)。从同位素异构体概念出发，通过原子转移矩阵和同位素异构体转移矩阵，可以将物质平衡方程推广到同位素异构体守恒方程。通过对中间代谢物建立同位素异构体守恒方程组，可以得到代谢通量和代谢物的同位素异构体分布之间的确定但是隐式的关系。这个关系是一个较为复杂的多元二次方程组。当我们测量了同位素异构体分布的信息，依据此大规模的方程组和参数优化的数值方法，可以反推出代谢通量。同时，通过构建理论的测量值空间，利用Monte Carlo采样，仿照此求解过程，可以估计代谢通量的置信区间。
第二章的主要内容是13C通量分析平台的建立和实例研究。
第二章的第一部分详细介绍了13C代谢通量分析的检测技术。
首先，详细介绍了13C标记技术，这是进行13C代谢通量分析的基础。对于估计最后的代谢通量而言，我们还需要获取两方面的信息。一方面是关于同位素异构体比例的信息。通过核磁(Nuclear magnetic resonance,NMR)波谱，可以获得特定的同位素异构体组合的比例信息。这包括J-resovled HSQC脉冲序列，它主要用于测量单键相连接的13C异构体组合。通过质谱(Mass spectrometry,MS)，可以测量质量同位素异构体比例的信息。另一方面是关于代谢物流出代谢体系的量值信息。一部分中间代谢物的流出量值通过细胞的生物质的质量和成分构成；还有一部分中间代谢物的流出量值，通过测量这些中间代谢物在培养液中的残留浓度来获取。最后介绍了参数优化的目标函数和置信区间估计的方法。
第二章的第二部分，以连续培养的大肠杆菌JM101株的代谢通量分析为例子，描述了连续培养下进行13C代谢通量分析的完整流程，也就是13C代谢通量分析实验平台的建立和实现。
第三章阐述了对13C质量同位素异构体检测技术所作的改进。改进后的方法可以提高质量同位素异构体比例的测量确度，从而提高13C代谢通量分析的准确度。
第三章的第一部分首先粗略地介绍了质谱测量的原理和基本过程。然后按照测量发生的顺序，介绍了质谱仪的基本的结构和功能单元。这些单元包括进样器，离子源，质量分析仪和信号检测器。之后论文较为详细地描述了用于气质联用的衍生化方法，特别是氨基酸的特丁基二甲基硅烷衍生化方法。最后，仔细地讨论了用以解决重同位素引起的质量位移作用的数学框架修正举证。
第三章的第二部分首先详细介绍了用质谱进行13C质量同位素异构体检测的基本假设，并对该假设进行了数学化的描述。利用该数学描述，对之前的各种提高质谱测量确度的方法进行了分类和描述。之后则介绍了如何提高13C质量同位素异构体测量确度的方法。利用已知了同位素异构体比例的样品，作为标准样，制作任意质谱仪的校准曲线，是一种最具有广泛性的提高测量精度的方法。我们使用以甲醇作为碳源的嗜甲基菌Methylobacterium salsuginis sp.，将其在排除了天然CO2的环境下，以13C标记的甲醇作为唯一碳源来进行培养。产生的各种生成物的同位素异构体比例可以用二项式公式来计算。我们引进一个新的量，它是两个连续的质量同位素异构体(13C原子数目相差1)的比值，由此避开了对多维向量进行校准的麻烦。而且，我们后验地测定每个产物的整个分子的平均13C比例，从而使得生物同位素效应的对13C分布不均匀性的影响降至最低。
第三章的第三部分则是对第二部分理论的实验实现，包括实验程序，理论模拟和结果讨论。经过校准后的结果，其绝对误差都小于0.9mol％，比未校准的结果有明显进步。而且，若采用理想的实验步骤，则误差可预计下降到0.4mol％以下。
在13C MFA平台上，我们分析研究了大肠肝菌中心代谢通量在超氧化应激情况下的变化情况。这部分内容是由芮斌同学和我共同完成的，请参看芮斌同学的博士论文。
此外，我们还进行了“不同剂量的氧化剂引起代谢系统不同强度反应”的研究，正在成文中。
对这个体系的研究产生了新的组学，即代谢通量组学(Fluxomics)。其重点是对目标体系内的代谢通量的速度的方向和取值进行测定。所获得的信息对刻画生物系统的生理特性、评估遗传和环境对细胞的影响具有重要意义。
代谢通量组学的主要研究手段是代谢通量分析，它可以分为纯计量性代谢通量分析和13C代谢通量分析两类。其中13C代谢通量分析是近十多年来，国际上新发展并仍然处于发展当中的技术。它除开具有纯计量性代谢通量分析的能力外，还具有获取信息更多，适用范围更广的优点。它能获得代谢途径的交换通量，并且不依赖于除稳态假设外的其他生理假设，因而适用各种试验条件。但是对于13C代谢通量分析技术，国内还未见有博士论文专题研究。
本博士论文的内容主要分成二部分，一部分描述在中国科技大学建立13C代谢通量分析平台的工作，另一部分是对通量分析平台中的13C质量同位素异构体检测技术进行改进的工作。
第一章详细介绍了13C代谢通量分析的基本理论。
代谢通量分析的基础是代谢通量拟稳态假设。在此假设基础上，可以建立各个代谢物的物质平衡方程。13C代谢通量分析实验在培养底物里引入不完全的13C标记。13C会分布到所有代谢物中，而由于13C在代谢物中数量和位置的不同，可以得到同一种物质的不同的13C同位素异构体和质量同位素异构体(同位素异构体的一类特殊组合形式)。从同位素异构体概念出发，通过原子转移矩阵和同位素异构体转移矩阵，可以将物质平衡方程推广到同位素异构体守恒方程。通过对中间代谢物建立同位素异构体守恒方程组，可以得到代谢通量和代谢物的同位素异构体分布之间的确定但是隐式的关系。这个关系是一个较为复杂的多元二次方程组。当我们测量了同位素异构体分布的信息，依据此大规模的方程组和参数优化的数值方法，可以反推出代谢通量。同时，通过构建理论的测量值空间，利用Monte Carlo采样，仿照此求解过程，可以估计代谢通量的置信区间。
第二章的主要内容是13C通量分析平台的建立和实例研究。
第二章的第一部分详细介绍了13C代谢通量分析的检测技术。
首先，详细介绍了13C标记技术，这是进行13C代谢通量分析的基础。对于估计最后的代谢通量而言，我们还需要获取两方面的信息。一方面是关于同位素异构体比例的信息。通过核磁(Nuclear magnetic resonance,NMR)波谱，可以获得特定的同位素异构体组合的比例信息。这包括J-resovled HSQC脉冲序列，它主要用于测量单键相连接的13C异构体组合。通过质谱(Mass spectrometry,MS)，可以测量质量同位素异构体比例的信息。另一方面是关于代谢物流出代谢体系的量值信息。一部分中间代谢物的流出量值通过细胞的生物质的质量和成分构成；还有一部分中间代谢物的流出量值，通过测量这些中间代谢物在培养液中的残留浓度来获取。最后介绍了参数优化的目标函数和置信区间估计的方法。
第二章的第二部分，以连续培养的大肠杆菌JM101株的代谢通量分析为例子，描述了连续培养下进行13C代谢通量分析的完整流程，也就是13C代谢通量分析实验平台的建立和实现。
第三章阐述了对13C质量同位素异构体检测技术所作的改进。改进后的方法可以提高质量同位素异构体比例的测量确度，从而提高13C代谢通量分析的准确度。
第三章的第一部分首先粗略地介绍了质谱测量的原理和基本过程。然后按照测量发生的顺序，介绍了质谱仪的基本的结构和功能单元。这些单元包括进样器，离子源，质量分析仪和信号检测器。之后论文较为详细地描述了用于气质联用的衍生化方法，特别是氨基酸的特丁基二甲基硅烷衍生化方法。最后，仔细地讨论了用以解决重同位素引起的质量位移作用的数学框架修正举证。
第三章的第二部分首先详细介绍了用质谱进行13C质量同位素异构体检测的基本假设，并对该假设进行了数学化的描述。利用该数学描述，对之前的各种提高质谱测量确度的方法进行了分类和描述。之后则介绍了如何提高13C质量同位素异构体测量确度的方法。利用已知了同位素异构体比例的样品，作为标准样，制作任意质谱仪的校准曲线，是一种最具有广泛性的提高测量精度的方法。我们使用以甲醇作为碳源的嗜甲基菌Methylobacterium salsuginis sp.，将其在排除了天然CO2的环境下，以13C标记的甲醇作为唯一碳源来进行培养。产生的各种生成物的同位素异构体比例可以用二项式公式来计算。我们引进一个新的量，它是两个连续的质量同位素异构体(13C原子数目相差1)的比值，由此避开了对多维向量进行校准的麻烦。而且，我们后验地测定每个产物的整个分子的平均13C比例，从而使得生物同位素效应的对13C分布不均匀性的影响降至最低。
第三章的第三部分则是对第二部分理论的实验实现，包括实验程序，理论模拟和结果讨论。经过校准后的结果，其绝对误差都小于0.9mol％，比未校准的结果有明显进步。而且，若采用理想的实验步骤，则误差可预计下降到0.4mol％以下。
在13C MFA平台上，我们分析研究了大肠肝菌中心代谢通量在超氧化应激情况下的变化情况。这部分内容是由芮斌同学和我共同完成的，请参看芮斌同学的博士论文。
此外，我们还进行了“不同剂量的氧化剂引起代谢系统不同强度反应”的研究，正在成文中。