⌇

Hill Equation

Math

HillEquation 模块实现了药效动力学（Pharmacodynamics, PD）中经典的 S 型 $E_{\text{max}}$ 曲线模型：用于将效应部位的药物浓度非线性映射为生物体产生的实际生理药效。在现代临床精准麻醉控制中，该模型常用于将大脑效应室丙泊酚（异丙酚）浓度转换为双频指数（Bispectral Index, BIS）——清醒状态下 BIS 约为 93，外科手术维持期则处于 40 至 60 之间。

响应特性在极低浓度下十分平坦，在半数有效浓度 $EC_{50}$ 附近呈现极其陡峭的响应斜率，而在高浓度区域因受体全部饱和而趋于定值。希尔系数 $\gamma$ 决定了患者从清醒跌入深度麻醉的阶跃突变敏感程度。该模块为无状态静态映射；浓度随时间的动态演化过程（药代动力学 PK）应在输入端通过传递函数建模。

数学模型

\text{occ} = \frac{C_e^{\gamma}}{EC_{50}^{\gamma} + C_e^{\gamma}}, \qquad E = E_0 - E_{max} \cdot \text{occ}

其中：

Ce（输入端口）为效应部位药物浓度（负值自动截断为 0）。
occ（输出端口 Occ.）为生物受体占有率，范围在 0 至 1 之间。
E（输出端口 Effect）为产生的生理效应指标（如 BIS 脑电麻醉深度评分）。

受体占有率在对数空间中计算为 1 / (1 + (EC50/Ce)^γ)，确保在超高陡度或极大浓度下数值平滑饱和而不产生浮点溢出。

输入与输出

Direction	ID	Label	类型	Status
→ In	`in1`	Ce	number	Required
← Out	`effect`	Effect	number	Output
← Out	`occupancy`	Occ.	number	Output

参数

参数	Label	类型	默认值	描述
`e0`	Baseline (E0)	number	`93`	Effect at zero concentration, for example an awake BIS of 93
`emax`	Maximum Effect (Emax)	number	`93`	Largest possible drop from the baseline
`ec50`	EC50	number	`3.4`	Concentration giving half the maximum effect, in the unit of Ce
`gamma`	Hill Slope (γ)	number	`3`	Steepness of the curve: 1 is hyperbolic (Michaelis-Menten), above 1 sigmoidal

使用示例

丙泊酚浓度至 BIS 脑电麻醉深度映射

3.4 µg/ml 的血浆浓度阶跃通过一阶惯性滞后环节 1/(10s+1)（表征血脑转运平衡速率 $k_{e0}$ ）到达脑组织效应室。BIS 自 93 出发，短暂停留后陡峭下沉并稳定在 46.5。

受体配体结合与剂量响应关系

定量分析药物分子在不同浓度梯度下与受体靶点的饱和结合动力学。

说明与最佳实践

模型个体差异：默认参数（ $E_0 = 93, E_{\text{max}} = 93, EC_{50} = 3.4, \gamma = 3$ ）用于教学展示；临床实践中必须依据群体药代药效模型按患者进行参数辨识。
闭环增益剧烈变化：由于 $EC_{50}$ 邻域的一阶导数极高，在此工作点附近的反馈控制回路等效开环增益会显著升高。