【系统生物学】利用XPPAUT求解SOD系统-程序员宅基地

利用XPPAUT求解SOD系统

微分方程背景介绍

常微分方程是描述一个未知函数及其导数之间关系的方程，其中未知函数只有一个自变量。在生物系统建模中，常微分方程被用来描述生物系统中状态变量随时间变化的关系。

在生物系统建模中, 刻画某状态变量随时间变化的关系, 用常微分方程来描述。
对于确定性的系统而言, 系统的状态变量随时间的改变, 可以表达成如下的微分方程组:
$\frac{d x_i}{d t}=\dot{x}=f_i\left(x_1, \cdots, x_n, p_1, \cdots, p_l, t\right) \quad i=1, \cdots, n$
这里 $x_i$ 表示变量, 如浓度; $p_j$ 表示参数, 如酶浓度或动力学常数; $t$ 表示时间。
写成矢量形式: $\quad \frac{d \mathbf{x}}{d t}=\dot{\mathbf{x}}=\mathbf{f}(\mathbf{x}, \mathbf{p}, t)$
其中, $\mathbf{x}=\left(x_1, \cdots, x_n\right)^T, \mathbf{f}=\left(f_1, \cdots, f_n\right)^T, \mathbf{p}=\left(p_1, \cdots, p_l\right)^T$ 。
例如: 代谢建模中, 底物降解, 浓度变化正比于实际的浓度, 就可以表示成:
$\dot{c}(t)=k \cdot c(t)$

SOD系统

SOD简介

SOD是指超氧化物歧化酶（Superoxide Dismutase），是一种抗氧化酶，它可以将细胞内产生的超氧自由基转化为过氧化氢和氧分子。超氧自由基是一种高度反应性的分子，它们可以与细胞内的其他分子发生反应，导致细胞损伤和炎症等问题。因此，SOD在维持细胞健康方面起着重要作用。

SOD通常包括多个亚型，如Cu/Zn-SOD、Mn-SOD和Fe-SOD等。这些亚型在不同的组织和环境中发挥不同的作用。例如，Cu/Zn-SOD主要存在于细胞质中，而Mn-SOD则主要存在于线粒体中。

图1 CU/ZN superoxide dismutase from neisseria meningitidis E73A mutant

图2 人体中的SOD

反应体系和微分方程

反应体系: $\quad \stackrel{\mathrm{c}_1}{\longrightarrow} \mathrm{O}_2^{--} \stackrel{\mathrm{SOD}}{\stackrel{\mathrm{C}_2}{\longrightarrow}} \mathrm{H}_2 \mathrm{O}_2 \stackrel{\text { cat }}{\stackrel{\mathrm{C}_3}{\longrightarrow}} \mathrm{H}_2 \mathrm{O}$
系统动力学方程:
$\left\{\begin{array}{l} \frac{d O_2^{\bullet-}}{d t}=c_1-c_2 \cdot S O D \cdot O_2^{\bullet-} \\ \frac{d H_2 O_2}{d t}=c_2 \cdot S O D \cdot O_2^{\bullet-}-c_3 \cdot c a t \cdot \mathrm{H}_2 \mathrm{O}_2 \end{array}\right.$

j简写为

$\left\{\begin{array}{l} \frac{d x}{d t}=a-b x \\ \frac{d y}{d t}=b x-c y \end{array}\right.$
(其中， $\left.x=O_2^{\bullet-}, y=H_2 O_2, a=c_1, b=c_2 \cdot S O D, c=c_3 \cdot c a t\right)$
给定参数值:
$c_1=6.6 \times 10^{-7}, \quad c_2=1.6 \times 10^9, \quad c_3=3.4 \times 10^7, \quad S O D=10^{-5}, \quad c a t=10^{-5}$
和初始条件: $\quad O_2^{*-}(0)=0, \mathrm{H}_2 \mathrm{O}_2(0)=0$

求解常微分方程

XPP-AUT是一种用于求解微分方程的软件工具。它可以用于求解多种类型的微分方程，包括常微分方程、差分方程、延迟微分方程、边界值问题和随机微分方程等。

ode文件

x'=a-b*x
y'=b*x-c*y
par a=6.6e-7,b=1.6e4,c=3.4e2
init x=0,y=0
done

图3 导入ODE文件

求解并画图的步骤如下：

将ODE文件拖入xpp – Shortcut中
点击nUmerics将Dt改为0.0001并设置Method为Q
Esc返回并点击Initialconds中的(G)o
选择Viewwaxes中的2D X设置合适参数
点击Graphic stuff的exp(o)rt data导出数据
在matlab（或其他程序）中画图

结果


图3 02(XPPAUT)	图4 H202(XPPAUT)


图5 O2(XPPAUT)	图6 H202(XPPAUT)

MATLAB画图代码如下：

data = table2array(readtable("H2O2.dat"));
plot(data(:,1),data(:,2))
xlim([0, 0.01]);
ylim([0, 2e-9]);
xlabel('Time');
ylabel('Concentration');
title('H2O2');
data = table2array(readtable("O2.dat"));
plot(data(:,1),data(:,2))
xlim([0, 0.01]);
ylim([0, 4.3e-11]);
xlabel('Time');
ylabel('Concentration');
title('O2');

本文链接：https://blog.csdn.net/dream_of_grass/article/details/129979952

原作者删帖不实内容删帖广告或垃圾文章投诉

智能推荐

分布式光纤传感器的全球与中国市场2022-2028年：技术、参与者、趋势、市场规模及占有率研究报告_预计2026年中国分布式传感器市场规模有多大-程序员宅基地

文章浏览阅读3.2k次。本文研究全球与中国市场分布式光纤传感器的发展现状及未来发展趋势，分别从生产和消费的角度分析分布式光纤传感器的主要生产地区、主要消费地区以及主要的生产商。重点分析全球与中国市场的主要厂商产品特点、产品规格、不同规格产品的价格、产量、产值及全球和中国市场主要生产商的市场份额。主要生产商包括：FISO TechnologiesBrugg KabelSensor HighwayOmnisensAFL GlobalQinetiQ GroupLockheed MartinOSENSA Innovati_预计2026年中国分布式传感器市场规模有多大

07_08 常用组合逻辑电路结构——为IC设计的延时估计铺垫_基4布斯算法代码-程序员宅基地

文章浏览阅读1.1k次，点赞2次，收藏12次。常用组合逻辑电路结构——为IC设计的延时估计铺垫学习目的：估计模块间的delay，确保写的代码的timing 综合能给到多少HZ，以满足需求！_基4布斯算法代码

OpenAI Manager助手（基于SpringBoot和Vue）_chatgpt网页版-程序员宅基地

文章浏览阅读3.3k次，点赞3次，收藏5次。OpenAI Manager助手（基于SpringBoot和Vue）_chatgpt网页版

关于美国计算机奥赛USACO，你想知道的都在这_usaco可以多次提交吗-程序员宅基地

文章浏览阅读2.2k次。USACO自1992年举办，到目前为止已经举办了27届，目的是为了帮助美国信息学国家队选拔IOI的队员，目前逐渐发展为全球热门的线上赛事，成为美国大学申请条件下，含金量相当高的官方竞赛。USACO的比赛成绩可以助力计算机专业留学，越来越多的学生进入了康奈尔，麻省理工，普林斯顿，哈佛和耶鲁等大学，这些同学的共同点是他们都参加了美国计算机科学竞赛（USACO)，并且取得过非常好的成绩。适合参赛人群USACO适合国内在读学生有意向申请美国大学的或者想锻炼自己编程能力的同学，高三学生也可以参加12月的第_usaco可以多次提交吗

MySQL存储过程和自定义函数_mysql自定义函数和存储过程-程序员宅基地

文章浏览阅读394次。1.1 存储程序1.2 创建存储过程1.3 创建自定义函数1.3.1 示例1.4 自定义函数和存储过程的区别1.5 变量的使用1.6 定义条件和处理程序1.6.1 定义条件1.6.1.1 示例1.6.2 定义处理程序1.6.2.1 示例1.7 光标的使用1.7.1 声明光标1.7.2 打开光标1.7.3 使用光标1.7.4 关闭光标1.8 流程控制的使用1.8.1 IF语句1.8.2 CASE语句1.8.3 LOOP语句1.8.4 LEAVE语句1.8.5 ITERATE语句1.8.6 REPEAT语句。_mysql自定义函数和存储过程

半导体基础知识与PN结_本征半导体电流为0-程序员宅基地

文章浏览阅读188次。半导体二极管——集成电路最小组成单元。_本征半导体电流为0

随便推点

【Unity3d Shader】水面和岩浆效果_unity 岩浆shader-程序员宅基地

文章浏览阅读2.8k次，点赞3次，收藏18次。游戏水面特效实现方式太多。咱们这边介绍的是一最简单的UV动画（无顶点位移），整个mesh由4个顶点构成。实现了水面效果（左图），不动代码稍微修改下参数和贴图可以实现岩浆效果（右图）。有要思路是1，uv按时间去做正弦波移动2，在1的基础上加个凹凸图混合uv3，在1、2的基础上加个水流方向4，加上对雾效的支持，如没必要请自行删除雾效代码(把包含fog的几行代码删除)S..._unity 岩浆shader

广义线性模型——Logistic回归模型（1）_广义线性回归模型-程序员宅基地

文章浏览阅读5k次。广义线性模型是线性模型的扩展，它通过连接函数建立响应变量的数学期望值与线性组合的预测变量之间的关系。广义线性模型拟合的形式为：其中g(μY)是条件均值的函数（称为连接函数）。另外，你可放松Y为正态分布的假设，改为Y 服从指数分布族中的一种分布即可。设定好连接函数和概率分布后，便可以通过最大似然估计的多次迭代推导出各参数值。在大部分情况下，线性模型就可以通过一系列连续型或类别型预测变量来预测正态分布的响应变量的工作。但是，有时候我们要进行非正态因变量的分析，例如：（1）类别型.._广义线性回归模型

HTML+CSS大作业环境网页设计与实现(垃圾分类) web前端开发技术 web课程设计网页规划与设计_垃圾分类网页设计目标怎么写-程序员宅基地

文章浏览阅读69次。环境保护、保护地球、校园环保、垃圾分类、绿色家园、等网站的设计与制作。总结了一些学生网页制作的经验：一般的网页需要融入以下知识点：div+css布局、浮动、定位、高级css、表格、表单及验证、js轮播图、音频视频 Flash的应用、ul li、下拉导航栏、鼠标划过效果等知识点，网页的风格主题也很全面：如爱好、风景、校园、美食、动漫、游戏、咖啡、音乐、家乡、电影、名人、商城以及个人主页等主题，学生、新手可参考下方页面的布局和设计和HTML源码（有用点赞△）一套A+的网_垃圾分类网页设计目标怎么写

C# .Net 发布后,把dll全部放在一个文件夹中,让软件目录更整洁_.net dll 全局目录-程序员宅基地

文章浏览阅读614次，点赞7次，收藏11次。之前找到一个修改 exe 中 DLL地址的方法, 不太好使,虽然能正确启动, 但无法改变 exe 的工作目录,这就影响了.Net 中很多获取 exe 执行目录来拼接的地址 ( 相对路径 ),比如 wwwroot 和代码中相对目录还有一些复制到目录的普通文件等等,它们的地址都会指向原来 exe 的目录, 而不是自定义的 “lib” 目录,根本原因就是没有修改 exe 的工作目录这次来搞一个启动程序,把 .net 的所有东西都放在一个文件夹,在文件夹同级的目录制作一个 exe._.net dll 全局目录

BRIEF特征点描述算法_breif description calculation 特征点-程序员宅基地

文章浏览阅读1.5k次。本文为转载，原博客地址：http://blog.csdn.net/hujingshuang/article/details/46910259简介 BRIEF是2010年的一篇名为《BRIEF:Binary Robust Independent Elementary Features》的文章中提出，BRIEF是对已检测到的特征点进行描述，它是一种二进制编码的描述子，摈弃了利用区域灰度..._breif description calculation 特征点

房屋租赁管理系统的设计和实现，SpringBoot计算机毕业设计论文_基于spring boot的房屋租赁系统论文-程序员宅基地

文章浏览阅读4.1k次，点赞21次，收藏79次。本文是《基于SpringBoot的房屋租赁管理系统》的配套原创说明文档，可以给应届毕业生提供格式撰写参考，也可以给开发类似系统的朋友们提供功能业务设计思路。_基于spring boot的房屋租赁系统论文