《Frontiers in Network Physiology》：The nephron arterial network and its interactions: nonlinear and information theoretic analyses
编辑推荐：
肾单位从树状动脉和小动脉网络接收血液，并通过影响其入球小动脉直径的动作来调节需求。网络结构是不对称的。血液流动需要静水压梯度，因此不同入球小动脉起点处的血压必然彼此不同。作用于入球小动脉的肾单位机制包括肌源性机制（平滑肌细胞全身普遍存在的特征）和管球反馈（tu
肾单位从树状动脉和小动脉网络接收血液，并通过影响其入球小动脉直径的动作来调节需求。网络结构是不对称的。血液流动需要静水压梯度，因此不同入球小动脉起点处的血压必然彼此不同。作用于入球小动脉的肾单位机制包括肌源性机制（平滑肌细胞全身普遍存在的特征）和管球反馈（tubuloglomerular feedback，TGF）（一种负反馈机制，当肾小管液从髓袢升支粗段进入远曲小管时，感知其中常见电解质的浓度）。肌源性机制在质膜的电位差中产生自主振荡。TGF信号也振荡，尽管频率较低，并调节肌源性机制电位差振荡的频率和幅度。TGF与肌源性机制的相互作用产生一种双峰电信号，该信号进入相邻的动脉段。网络的动脉段具有电传导性，为肾单位之间的相互作用提供了通道。此外，由肌源性-TGF相互作用引起的肾单位血流振荡导致局部动脉血压的变化，从而在肾单位之间提供血流动力学相互作用。研究人员开发了一个包含这些特征的数学模型。该模型展示了肾自身调节以及实验中发现的肾小管液流量和压力的振荡。结果说明了相互作用对网络的影响。本文中，研究人员估计了李雅普诺夫指数（Lyapunov exponents）以确定每个肾单位的动态状态，并利用互信息（mutual information）对各个信息流对肾单位动力学的重要性进行定量评估。结果表明，所有肾单位均处于弱混沌动力学状态，且通过动脉的电信号传导对肾单位动力学的影响大于血流动力学相互作用。混沌吸引子对初始条件敏感，因此研究人员认为肾单位簇对系统中的信息流敏感，从而使肾单位簇能够适应变化的外部条件。
本研究围绕肾单位动脉网络及其相互作用展开，旨在阐明网络结构如何通过非线性动力学和信息传递影响肾单位的整体行为。肾单位通过入球小动脉调节血流，其机制包括全身平滑肌共有的肌源性机制和肾脏特有的管球反馈（tubuloglomerular feedback，TGF）。肌源性机制产生质膜电位差振荡，TGF信号以较低频率调节该振荡的幅度和频率，二者共同形成双峰电信号。动脉网络具有电传导性，肾单位间通过电信号和血流动力学相互作用。然而，以往单肾单位模型未能包含完整动脉网络，难以解释实验观察到的肾单位间同步及血压波动的复杂动态。为此，Marsh和Holstein-Rathlou开发了包含10个肾单位及供应动脉的网络模型，并采用非线性时间序列分析工具定量研究网络中的信息流。该论文发表在《Frontiers in Network Physiology》。
研究人员基于大鼠肾血管铸型的显微计算机断层扫描数据构建了网络拓扑，将单肾单位模型与动脉模型耦合，模拟了85至115 mmHg七种动脉压水平下的急性高血压响应，并施加1/f噪声血压输入模拟清醒动物血压的生理波动。主要方法包括：单肾单位模型（肾小球滤过率由常微分方程计算，肾小管压力、流量和NaCl浓度由偏微分方程组描述）；入球小动脉模型采用Gonzalez-Fernandez-Ermentrout模型模拟肌源性机制；网络模型求解各血管节点的血流和电流守恒；数据采用最大李雅普诺夫指数（Lyapunov exponent）判定动力学状态，采用互信息（mutual information）量化信息流强度，并利用小扰动和雅可比矩阵预测检验混沌特征，使用小波变换分析时频特征。
结果部分：
1. 模型验证与时间序列特征：模型成功模拟了肾血流的自身调节以及实验中观察到的肾小管液流量和压力振荡。随着根动脉压升高，肾小管压力波动幅度增大，TGF振荡周期延长；血管内皮电位变化主要表现为肌源性振荡的相位移动，幅度变化较小。
2. 李雅普诺夫指数与混沌动力学：对10个肾单位在7个血压水平下计算最大李雅普诺夫指数，全部70个值均为正，表明所有肾单位均处于混沌状态。指数值较小，提示为弱混沌。通过施加10
-3 mmHg的微小肾小管压力扰动，系统轨迹随时间指数发散，到达预测水平的时间与李雅普诺夫指数大小成反比，进一步确认混沌的存在。
3. 拓扑变异性检验：改变网络拓扑中肾单位簇的起源位置和动脉长度，构建两种变体结构。三种模型中210个肾单位的李雅普诺夫指数均为正，平均值和标准差相似，表明拓扑变异不影响混沌动力学的结论，且自身调节功能保持不变。
4. 互信息分析：计算肾小管压力与动脉节点电位差、肾小管压力与动脉节点血压之间的互信息。结果表明，肾小管压力与动脉节点电位差的互信息值在70个数据点中有55个更高，说明电信号传导对肾小管压动力学的影响强于血流动力学压力相互作用。
5. 1/f血压输入下的网络行为：当根动脉压施加1/f噪声扰动，平均压为85和115 mmHg时，模型仍能维持自身调节（自身调节指数为-0.01），且肾小管压力在TGF频率处出现持续峰，提示网络起到类似时钟的作用，可协调肾单位活动。
讨论部分指出，单肾单位模型本可产生极限环振荡，但网络耦合引入了电信号和血流动力学相互作用的额外扰动，使系统从极限环分岔进入弱混沌状态。弱混沌既保留了一定规律性，又扩大了相空间访问范围，有利于肾单位簇对环境变化做出适应性响应。电信号传导占主导地位，可能源于入球小动脉动力学对膜电位变化的直接依赖。研究局限性包括：模型基于麻醉年轻雄性大鼠肾腹侧表面数据，肾单位类型为浅表肾单位，髓袢局限于外髓，且模拟时间约为1小时，因此结果应视为网络相互作用的示例。
研究结论是：肾单位动脉网络中的电信号和血流动力学相互作用使所有肾单位进入弱混沌动力学状态；这种状态对初始条件敏感，使肾单位簇能够响应系统中的信息流，从而实现功能的适应性整合，并维持稳定的肾自身调节。
