一种基于布尔网络模拟干预肿瘤细胞状态的方法技术

本发明专利技术实施例公开了一种基于布尔网络模拟干预肿瘤细胞状态的方法，该方法包括：利用基因芯片采集技术获取生物样本基因的连续表达谱数据，并根据所述获取到的连续表达谱数据，构建布尔网络模型的基因调控网络；确定基因调控网络干预前所有状态下的暂态数BOS及其包含的所有吸引子，并根据基因调控网络干预前所有状态下的暂态数BOS，在所有吸引子中，筛选出满足预定条件的干预位及其干预后基因调控网络的暂态数BOS；根据所筛选干预位干预后基因调控网络的暂态数BOS，调整基因调控网络结构，并根据调整后的基因调控网络，模拟肿瘤细胞状态。实施本发明专利技术，能够模拟干预肿瘤细胞状态，为肿瘤治疗研究提供一个有力的理论框架。

【技术实现步骤摘要】

本专利技术涉及系统生物学研究
，尤其涉及一种基于布尔网络模拟干预肿瘤细胞状态的方法。
技术介绍
在基因组时代，人们主要探讨的是静态碱基测序的问题，其具体任务是建立可以体现生物全基因组结构的图谱，到了后基因组时代，人们则开始重点探究非静态的功能注释领域的问题，它通过对各种高通量技术所获得的基因表达数据进行具体的分析、处理和建模，来帮助人们理解基因的功能，更好的认清基因与疾病的关系，进而设计适当的干预策略来影响、改变系统的动态行为。因此，人们发现除了直接分析表达数据，还能够更深入地理解基因在影响健康与疾病工程过程中的基本机制。随着计算机技术的发展，基因调控网络(Genetic regulatory networks,GRN)的研究成为21世纪生物学研究的一个很重要的领域，也成为系统生物学研究的一个热点问题，使得越来越多的科研人员开始关注基因调控网络。基因调控网络是由细胞中参与基因调控作用的DNA，RNA，蛋白质以及代谢中间物所形成的相互作用的网络，其节点会随着时间的变化而变化的，并且基因调控网络中细胞内的蛋白质浓度是基因调控网络动力学机制的最基本的驱动因素，它决定了细胞分化的时间和空间特性，还可以作为细胞的一种“记忆机制”。因此通过构造出的基因调控网络，人们能够更好的认识基因与基因之间的相互作用关系，从而进一步理解特定组织的遗传调控机制，开发合适的疾病(如癌症)治疗方法，这对于揭示生命现象的本质具有重要意义。因此，亟需一种以基因调控网络为模型来模拟干预肿瘤细胞状态的方法，能够模拟干预肿瘤细胞状态，为肿瘤治疗研究提供一个有力的理论框架。
技术实现思路
本专利技术实施例的目的在于提供一种基于布尔网络模拟干预肿瘤细胞状态的方法，能够模拟干预肿瘤细胞状态，为肿瘤治疗研究提供一个有力的理论框架。为了解决上述技术问题，本专利技术实施例提供了一种基于布尔网络模拟干预肿瘤细胞状态的方法，所述方法包括：a、利用基因芯片采集技术获取生物样本基因的连续表达谱数据，并根据所述获取到的连续表达谱数据，构建布尔网络模型的基因调控网络；b、确定所述基因调控网络干预前所有状态下的暂态数BOS及其包含的所有吸引子，并根据所述基因调控网络干预前所有状态下的暂态数BOS，在所述确定的所有吸引子中，筛选出满足预定条件的干预位，且进一步得到所筛选干预位干预后基因调控网络的暂态数BOS；c、根据所述得到的所筛选干预位干预后基因调控网络的暂态数BOS，调整所述基因调控网络结构，并根据所述调整后的基因调控网络，模拟肿瘤细胞状态。其中，所述步骤a具体包括：确定生物样本，并采用基因芯片采集技术对所述生物样本进行特定时间间隔取样，获取M个基因N个时间的连续表达谱数据；其中，M、N均为自然数；根据所述获取到的M个基因N个时间的连续表达谱数据，并通过预设的M个基因两两之间的配对关系赋值，得到M个基因两两之间的调控关系距离，且进一步得到M个基因两两之间的调控关系方向及其对应的调控关系相位；根据所述得到的M个基因两两之间的调控关系距离、调控关系方向和调控关系相位，构建布尔网络模型的基因调控网络。其中，所述步骤b具体包括：确定所述基因调控网络在干预前所有状态下的暂态数BOS，并进一步确定所述基因调控网络中所包含的所有吸引子；在所述基因调控网络在干预前所有状态下的暂态数BOS及所包含的所有吸引子中，查询出对所有吸引子产生影响的一个或多个干预位；确定每一个干预位干预后分别对应基因调控网络的暂态数BOS，并根据所述基因调控网络在干预前所有状态下的暂态数BOS、每一个干预位干预后分别对应基因调控网络的暂态数BOS以及预设的目标函数，得到每一个干预位对应的目标函数值；在所述得到的每一个干预位对应的目标函数值中，筛选出目标函数值最大时对应的干预位，并在所述每一个干预位干预后分别对应基因调控网络的暂态数BOS中，得到所筛选干预位干预后基因调控网络的暂态数BOS。实施本专利技术实施例，具有如下有益效果：本专利技术实施例基于布尔网络确定基因调控网络的最佳干预位，使基因调控网络能够通过干预位形成的干预策略逆向回到原来的状态或转移到另外一个期望的状态，最终改善基因调控网络使其往期望的方向发展，从而用于模拟干预肿瘤细胞状态，为肿瘤治疗研究提供一个有力的理论框架。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，根据这些附图获得其他的附图仍属于本专利技术的范畴。图1为本专利技术实施例提供的基于布尔网络模拟干预肿瘤细胞状态的方法的流程图；图2为本专利技术实施例提供的基于布尔网络模拟干预肿瘤细胞状态的方法应用场景中的布尔网络示意图；图3为本专利技术实施例提供的基于布尔网络模拟干预肿瘤细胞状态的方法应用场景中的布尔网络一位干预后状态变化示意图。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本专利技术，并不用于限定本专利技术。如图1所示，为本专利技术实施例中，提出的一种基于布尔网络模拟干预肿瘤细胞状态的方法，所述方法包括：步骤S101、利用基因芯片采集技术获取生物样本基因的连续表达谱数据，并根据所述获取到的连续表达谱数据，构建布尔网络模型的基因调控网络；具体过程为，步骤S11、确定生物样本，并采用基因芯片采集技术对生物样本进行特定时间间隔取样，获取M个基因N个时间的连续表达谱数据；其中，M、N均为自然数；步骤S12、根据获取到的M个基因N个时间的连续表达谱数据，并通过预设的M个基因两两之间的配对关系赋值，得到M个基因两两之间的调控关系距离，且进一步得到M个基因两两之间的调控关系方向及其对应的调控关系相位；其中，M个基因两两之间的调控关系方向由其对应的配对关系赋值的正负符号决定；M个基因两两之间的调控关系相位由其对应的调控关系距离与其对应的配对关系赋值的绝对值决定；步骤S13、根据得到的M个基因两两之间的调控关系距离、调控关系方向和调控关系相位，且按照M个基因两两之间的调控关系距离的绝对值大小顺序进行循环选择，构建布尔网络模型的基因调控网络.应当说明的是，在步骤S12中，当M个基因两两之间的调控关系距离采用N个时间点对齐计算时，则预设的M个基因两两之间的配对关系赋值为0；当M个基因两两之间的调控关系距离采用小于或等于k个错位计算时，则预设的M个基因两两之间的配对关系赋值可从-k至k，并将计算得到的2k+1个调控关系距离中的最大值作为M个基因两两之间的最终调控关系距离；其中，k为自然数。应当说明的是，在步骤S13中，当检测到当前循环增加的两个基因与初始基因调控网络中已有的基因具有直接或间接链接关系时，则在后续的选择过程中忽略当前循环增加的两个基因之间的调控关系。作为一个例子，以M个基因N个时间点的连续表达谱数据为例：取M个基因中任一个和剩余的M-1个基因两两配对。对于调控关系距离Dp的计算，可以使两个基因的相应的N个时间点对齐，也可以前后进行小于或等于k个错位，得到2k+1个相位调控关系距离：D-k、D-k+1、D-k+2、...、D-1、D0、D1本文档来自技高网...

【技术保护点】
一种基于布尔网络模拟干预肿瘤细胞状态的方法，其特征在于，所述方法包括：a、利用基因芯片采集技术获取生物样本基因的连续表达谱数据，并根据所述获取到的连续表达谱数据，构建布尔网络模型的基因调控网络；b、确定所述基因调控网络干预前所有状态下的暂态数BOS及其包含的所有吸引子，并根据所述基因调控网络干预前所有状态下的暂态数BOS，在所述确定的所有吸引子中，筛选出满足预定条件的干预位，且进一步得到所筛选干预位干预后基因调控网络的暂态数BOS；c、根据所述得到的所筛选干预位干预后基因调控网络的暂态数BOS，调整所述基因调控网络结构，并根据所述调整后的基因调控网络，模拟肿瘤细胞状态。

【技术特征摘要】
1.一种基于布尔网络模拟干预肿瘤细胞状态的方法，其特征在于，所述方法包括：a、利用基因芯片采集技术获取生物样本基因的连续表达谱数据，并根据所述获取到的连续表达谱数据，构建布尔网络模型的基因调控网络；b、确定所述基因调控网络干预前所有状态下的暂态数BOS及其包含的所有吸引子，并根据所述基因调控网络干预前所有状态下的暂态数BOS，在所述确定的所有吸引子中，筛选出满足预定条件的干预位，且进一步得到所筛选干预位干预后基因调控网络的暂态数BOS；c、根据所述得到的所筛选干预位干预后基因调控网络的暂态数BOS，调整所述基因调控网络结构，并根据所述调整后的基因调控网络，模拟肿瘤细胞状态。2.如权利要求1所述的方法，其特征在于，所述步骤a具体包括：确定生物样本，并采用基因芯片采集技术对所述生物样本进行特定时间间隔取样，获取M个基因N个时间的连续表达谱数据；其中，M、N均为自然数；根据所述获取到的M个基因N个时间的连续表达谱数据，并通过预设的M个基因两两之间的配对关系赋值，得到M个基因两两之...

【专利技术属性】
技术研发人员：沈良忠，
申请(专利权)人：温州大学，
类型：发明
国别省市：浙江;33