基于伪苹果酸-天冬氨酸穿梭的计算建模揭示外源性酮体急性给药后人脑谷氨酸和GABA的动态变化

《Neurochemical Research》：Computational Modeling of Neurotransmitter Cycling Predicts Human Brain Glutamate and GABA Dynamics in Response to Acute Administration of Exogenous Ketones

　　

酮体的给药被用作多种疾病（包括癫痫、精神健康障碍和脑衰老）的治疗选择。其作用机制之一涉及调节谷氨酸（glutamate）和GABA——大脑主要的兴奋性和抑制性神经递质——二者共同调控兴奋-抑制平衡。然而，酮体影响这些神经递质的确切机制仍不清楚。在本研究中，研究人员假设酮体通过伪苹果酸-天冬氨酸穿梭（pseudo-malate-aspartate shuttle, PMAS）调节谷氨酸和GABA。为了验证这一点，研究人员开发了一个以PMAS为核心的神经递质循环计算模型，模拟了作为酮体代谢函数的谷氨酸和GABA的时间动态和稳态浓度。随后，研究人员将模型输出与来自酮体给药实验的磁共振波谱（MRS）数据进行比较，发现与模型预测一致，从而为该模型提供了定量支持。基于这一一致性，研究人员进行了代谢控制分析，确定了通过PMAS部分替代葡萄糖代谢是观察到的谷氨酸和GABA降低的主导机制，并揭示了选择性调节每种神经递质的关键酶。总体而言，该模型为研究人员和临床医生提供了一个用于假设检验和治疗优化的框架，同时也为未来的模型扩展奠定了基础。

研究背景、存在问题与研究目的

酮体是葡萄糖的替代能量底物，在多种脑部疾病中具有治疗潜力，已在难治性癫痫中确立临床地位，并在精神健康障碍、年龄相关的神经退行性疾病甚至健康老龄化中显示出潜在益处。正电子发射断层扫描（PET）和¹³C磁共振波谱（¹³C MRS）研究表明，无论是内源性产生还是外源性给予的酮体都能被人脑摄取和利用。除作为能量底物外，酮体还可能通过改变细胞内的谷氨酸（glutamate）和γ-氨基丁酸（GABA）水平发挥神经调节作用。最近一项磁共振波谱（MRS）研究发现，急性给予β-羟基丁酸（β-hydroxybutyrate, β-HB）可使谷氨酸和GABA分别降低11%和37%。尽管酮体具有治疗潜力，但其调节神经递质水平的机制缺乏共识，这对其治疗应用构成挑战。神经递质的维持与酮体氧化共享重叠的能量代谢通路，而整合不同实验观测结果（涵盖酮体药代动力学、神经递质代谢和既定能量约束）极具挑战性，因为存在三种高度特化的细胞类型：谷氨酸能神经元、GABA能神经元和胶质细胞。现有计算模型的参数多源于体外动物研究，可能无法准确反映活体人脑代谢。因此，本研究旨在开发一个结合酮体利用药代动力学与神经递质循环的计算代谢模型，以确定急性酮体暴露引起谷氨酸和GABA代谢变化的机制基础。

研究方法与技术路线

研究人员开发了一个以伪苹果酸-天冬氨酸穿梭（PMAS）为核心的计算代谢模型，整合了代数表达式和微分方程，参数主要来源于体内研究。模型通过实验发现的神经递质循环速率与氧化葡萄糖代谢之间接近化学计量学的1:1耦合关系进行约束。模型包含三个主要细胞区室：谷氨酸能神经元、GABA能神经元和星形胶质细胞。对于酮体代谢，研究人员采用可逆米氏动力学描述D-β-HB从血液到大脑的转运，并基于线性关系建立其氧化速率。代谢控制分析（MCA）用于评估不同酶和通路活性对谷氨酸、GABA和谷氨酰胺池的敏感性。

研究结果

神经递质循环与葡萄糖代谢的耦合

研究人员基于静息态人脑PET和MRSI数据的汇总结果，采用总灰质葡萄糖代谢率（CMR_glc）为0.41 mM/min。通过氧-葡萄糖指数（OGI）推导出氧化部分，再经神经元分数校正，得到神经元氧化葡萄糖代谢率（CMR_glc(ox)^N）为0.28 mM/min。根据先前工作，谷氨酸能和GABA能神经元的CMR_glc(ox)分别占83%和17%，由此计算出各自的PMAS速率（V_PMAS^e=0.19 mM/min, V_PMASⁱ=0.048 mM/min）。基于突触囊泡装载与胞质谷氨酸呈线性关系的原理，研究人员确定了谷氨酸释放速率常数。

GABA循环的动力学特性

研究人员区分了由两种谷氨酸脱羧酶（GAD65和GAD67）产生的两个不同的GABA池：GAD65供应突触囊泡GABA池（GABA_i^v），参与信号传导；GAD67产生胞质GABA池（GABA_i^c），通过GABA分流循环。在稳态下，两个GAD同工酶的流量均设为等于V_PMASⁱ。GAD65的低米氏常数（K_M^GAD65=0.01 mM）使其对谷氨酸浓度变化不敏感，而GAD67的K_M等于胞内谷氨酸浓度（0.13 mM）。GABA转氨酶（GABA-T）的K_M为0.16 mM，其对底物的低敏感性导致扰动后胞质GABA发生不成比例的大幅偏移。

星形胶质细胞反应的动力学特性

星形胶质细胞中的谷氨酸主要通过谷氨酰胺合成酶（GS）转化为谷氨酰胺，另有小部分经谷氨酸脱氢酶（GDH）氧化。GDH采用希尔动力学建模（希尔系数h_GDH=4.0），其流量设定为神经元谷氨酸和GABA释放总量的20%。回补反应（anaplerosis）被建模为恒定流量，以补充因氧化和谷氨酰胺外流而损失的谷氨酸。GS的动力学包含产物抑制项，其最大速率基于癫痫发作期间的最大氧化代谢率确定。

酮体代谢及其对神经递质循环的影响

研究人员采用可逆米氏动力学描述D-β-HB的血脑屏障转运，并从两项急性人体实验中推导出转运参数（T_max^ket=0.462 mM/min, K_T^ket=27.5 mM）。酮体氧化速率与脑组织酮体浓度呈线性关系。当使用实验中的血酮浓度驱动模型时，模型无需修改即可预测出与实验一致的谷氨酸和GABA降低幅度。

代谢控制分析与机制识别

通过代谢控制分析，研究人员确定β-羟基丁酸部分替代葡萄糖代谢是导致谷氨酸和GABA降低的主要机制。这种替代减少了通过PMAS的通量，从而减缓神经递质循环并降低谷氨酸和GABA的稳态浓度。GABA更大百分比的降低归因于PMAS抑制与GAD65、GAD67和GABA-T动力学特性的共同作用。

讨论与结论

该计算模型成功再现了急性酮体给药后实验观察到的谷氨酸和GABA降低，为假定的PMAS机制提供了定量支持。代谢控制分析表明，酮体部分取代葡萄糖代谢是主导机制，而GABA更大的降幅则源于其独特的酶动力学特征，特别是GAD67和GABA-T之间的放大效应。虽然本研究聚焦于急性酮体输注，但模型的代谢网络方程、化学计量关系和酶动力学参数均独立确定，因此有望用于解释和指导脑部药理和代谢干预。模型为研究人员和临床医生提供了一个用于假设检验和治疗优化的框架，也为未来模型扩展奠定了基础。