DOI: https://doi.org/10.1098/rsos.250910
تاريخ النشر: 2026-03-25
المؤلف: Xabier Rey Barreiro وآخرون
الموضوع الرئيسي: علم المناعة للخلايا التائية والخلايا البائية
نظرة عامة
تناقش هذه الفقرة تطوير نماذج رياضية لتنشيط خلايا T، مع التركيز على التفاعل بين المستضدات ومستقبلات خلايا T (TCRs) التي تبدأ الاستجابات المناعية. يبرز المؤلفون أهمية كميتين رئيسيتين قابلتين للقياس: $E_{\text{max}}$، المستوى الأقصى الممكن لتنشيط خلايا T، و$EC_{50}$، تركيز المستضد الذي يحفز 50% من $E_{\text{max}}$. بينما تعتبر هذه الكميات حاسمة لتحليل خصائص النموذج مثل القابلية للتحديد، والملاحظة، والحساسية، فإن تعبيراتها تم اشتقاقها فقط لعدد محدود من النماذج الموجودة. يقترح المؤلفون طريقة عامة لاشتقاق هذه التعبيرات لنماذج إضافية، مما يعزز القدرات التحليلية المتعلقة بتنشيط خلايا T وآثارها على التحكم في الاستجابة المناعية.
في استنتاجاتهم، يؤكد المؤلفون مساهمتهم في الأدبيات من خلال تقديم التعبيرات غير المعروفة سابقًا لـ $E_{\text{max}}$ و$EC_{50}$ عبر نماذج مختلفة لتنشيط خلايا T. يلاحظون أنه بينما تنتج النماذج الأبسط تعبيرات مباشرة، فإن النماذج الأكثر تعقيدًا تؤدي إلى صيغ معقدة، مما يعقد العلاقة بين معلمات النموذج وهذه الكميات. قد تؤدي هذه التعقيدات إلى تحديات في استنتاج قيم المعلمات بدقة من قياسات $E_{\text{max}}$ و$EC_{50}$، مما يشير إلى مشاكل محتملة في القابلية للتحديد داخل هذه النماذج. يشير المؤلفون إلى التحقيقات الجارية في هذه المسألة، مع نتائج سيتم الإبلاغ عنها في الأعمال المستقبلية.
مقدمة
في الفقاريات، تتميز الاستجابة المناعية باستراتيجيتين رئيسيتين: المناعة الفطرية والمناعة التكيفية. يعمل النظام المناعي الفطري كآلية دفاع سريعة وغير محددة، بينما ينتج النظام المناعي التكيفي استجابات محددة للغاية لمستضدات معينة (Ags). في صميم الاستجابة التكيفية توجد الخلايا اللمفاوية T، التي تتفاعل مع خلايا تقديم المستضد (APCs) من خلال مستقبلات خلايا T (TCRs). يتم تنشيط خلايا T من خلال ارتباط TCRs بمجمعات معقد التوافق النسيجي الرئيسي (pMHC) على APCs، مما يؤدي إلى مسارات إشارات داخلية معقدة تؤثر على قرارات مصير الخلايا. على الرغم من أهمية هذه العمليات، إلا أن الآليات الأساسية لا تزال مفهومة جزئيًا فقط، مما يحفز تطوير نماذج رياضية لوصف تنشيط خلايا T بشكل كمي.
تتطلب هذه النماذج، التي تمثل غالبًا كمجموعات من المعادلات التفاضلية العادية (ODEs)، المعايرة مقابل البيانات التجريبية للحصول على تقديرات دقيقة للمعلمات. الخصائص الهيكلية الرئيسية لهذه النماذج، مثل القابلية للتحديد والملاحظة، حاسمة لفهم العلاقة بين معلمات النموذج والمخرجات القابلة للقياس، وبشكل خاص المستوى الأقصى لاستجابة خلايا T ($E_{max}$) وتركيز الفعالية نصف القصوى ($EC_{50}$). ومع ذلك، فإن التعبيرات الرياضية التي تربط هذه المخرجات بمعلمات النموذج محدودة في الأدبيات. تهدف هذه الورقة إلى معالجة هذه الفجوة من خلال تقديم منهجية عامة لاشتقاق التعبيرات لـ $E_{max}$ و$EC_{50}$ القابلة للتطبيق على نماذج مختلفة لتنشيط خلايا T. ستوفر الأقسام اللاحقة تعبيرات محددة لمجموعة من النماذج، بما في ذلك اشتقاقات مفصلة لتلك التي لم يتم الإبلاغ عنها سابقًا، وتقترح تعديلات على النماذج الموجودة حيث تم تحديد التباينات.
طرق
في هذا القسم، يحدد المؤلفون المنهجية المستخدمة لتحليل نموذج تصحيح كينتيك مكيثان (KPR)، كما تم الإشارة إليه في الدراسات السابقة [15، 3]. يتميز النموذج بنظام من المعادلات التفاضلية العادية (ODEs)، موصوفة بالتفصيل في القسم الفرعي 3.2، مع تقديم المعادلات ذات الصلة في المعادلة (19). يتم تمثيل مستوى تنشيط خلايا T رياضيًا كـ \( R(t) = C N(t) \)، حيث \( C \) هو ثابت و \( N(t) \) يدل على المتغير المعتمد على الزمن المتعلق بديناميات خلايا T. تشكل هذه الإطار الأساس لاستكشاف المزيد من آثار نموذج KPR في عمليات تنشيط خلايا T.
مناقشة
في هذا القسم، يناقش المؤلفون نماذج ظاهرة مختلفة لتنشيط خلايا T، ملخصين معلماتها الرئيسية وصيغها الرياضية. تشمل النماذج، الموضحة في الجدول 1 والمصورة في الشكل 1، ثوابت معدل حركية شائعة ومعلمات محددة فريدة لكل نموذج. يستخدم المؤلفون نهجًا منهجيًا باستخدام المعادلات التفاضلية العادية (ODEs) لوصف ديناميات تنشيط خلايا T، مع التركيز على المخرجات $E_{\text{max}}$ و$EC_{50}$ كقياسات لتنشيط خلايا T استجابةً لإجمالي الربط لكل خلية تقديم مستضد (APC).
يتضمن اشتقاق $E_{\text{max}}$ و$EC_{50}$ عدة خطوات، بدءًا من افتراض ظروف الحالة المستقرة، حيث يتم تعيين المعادلات التفاضلية إلى الصفر. يعرف المؤلفون مستوى استجابة خلايا T $R$ من حيث مستقبلات خلايا T المحفزة ويشتقون تعبيرات لـ $E_{\text{max}}$ كأقصى استجابة عند تشبع الربط و$EC_{50}$ كتركيز الربط الذي ينتج استجابة نصف قصوى. يوضح القسم الإجراءات الرياضية لاشتقاق هذه التعبيرات، مع التأكيد على الاختلافات في السلوك بين النماذج المختلفة، مثل نموذج الإشغال، ونماذج تصحيح كينتيك، وتلك التي تتضمن آليات التغذية الراجعة. يتم تقديم النتائج بالتفصيل، مع توفير تعبيرات مغلقة للعديد من النماذج، مما يساهم في فهم شامل لديناميات تنشيط خلايا T.
DOI: https://doi.org/10.1098/rsos.250910
Publication Date: 2026-03-25
Author(s): Xabier Rey Barreiro et al.
Primary Topic: T-cell and B-cell Immunology
Overview
This section discusses the development of mathematical models for T-cell activation, focusing on the interaction between antigens and T-cell receptors (TCRs) that initiates immune responses. The authors highlight the significance of two key measurable quantities: $E_{\text{max}}$, the maximum achievable level of T-cell activation, and $EC_{50}$, the antigen concentration that induces 50% of $E_{\text{max}}$. While these quantities are crucial for analyzing model properties such as identifiability, observability, and sensitivity, their expressions have only been derived for a limited number of existing models. The authors propose a general method to derive these expressions for additional models, thereby enhancing the analytical capabilities related to T-cell activation and its implications for immune response control.
In their conclusions, the authors emphasize their contribution to the literature by providing the previously unknown expressions for $E_{\text{max}}$ and $EC_{50}$ across various T-cell activation models. They note that while simpler models yield straightforward expressions, more complex models result in intricate formulas, complicating the relationship between model parameters and these quantities. This complexity may lead to challenges in accurately inferring parameter values from measurements of $E_{\text{max}}$ and $EC_{50}$, suggesting potential identifiability issues within these models. The authors indicate ongoing investigations into this matter, with results to be reported in future work.
Introduction
In vertebrates, the immune response is characterized by two main strategies: innate and adaptive immunity. The innate immune system serves as a rapid, non-specific defense mechanism, while the adaptive immune system generates highly specific responses to particular antigens (Ags). Central to the adaptive response are T lymphocytes, which interact with antigen-presenting cells (APCs) through T-cell receptors (TCRs). The activation of T cells is initiated by the binding of TCRs to peptide-major histocompatibility complex (pMHC) complexes on APCs, leading to complex intracellular signaling pathways that influence cell fate decisions. Despite the significance of these processes, the underlying mechanisms remain only partially understood, prompting the development of mathematical models to quantitatively describe T-cell activation.
These models, often represented as sets of ordinary differential equations (ODEs), require calibration against experimental data to yield accurate parameter estimates. Key structural properties of these models, such as identifiability and observability, are crucial for understanding the relationship between model parameters and measurable outputs, specifically the maximum T-cell response level ($E_{max}$) and the half-maximal effective concentration ($EC_{50}$). However, mathematical expressions linking these outputs to model parameters are limited in the literature. This paper aims to address this gap by presenting a general methodology for deriving expressions for $E_{max}$ and $EC_{50}$ applicable to various T-cell activation models. Subsequent sections will provide specific expressions for a range of models, including detailed derivations for those previously unreported, and propose modifications to existing models where discrepancies have been identified.
Methods
In this section, the authors outline the methodology employed to analyze McKeithan’s kinetic proofreading (KPR) model, as referenced in prior studies [15, 3]. The model is characterized by a system of ordinary differential equations (ODEs), specifically detailed in subsection 3.2, with the relevant equations presented in Equation (19). The T-cell activation level is mathematically represented as \( R(t) = C N(t) \), where \( C \) is a constant and \( N(t) \) denotes the time-dependent variable related to T-cell dynamics. This framework serves as the foundation for further exploration of the KPR model’s implications in T-cell activation processes.
Discussion
In this section, the authors discuss various phenotypic models of T-cell activation, summarizing their key parameters and mathematical formulations. The models, detailed in Table 1 and illustrated in Fig. 1, include common kinetic rate constants and specific parameters unique to each model. The authors employ a systematic approach using ordinary differential equations (ODEs) to describe T-cell activation dynamics, focusing on the outputs $E_{\text{max}}$ and $EC_{50}$ as measures of T-cell activation in response to total ligand per antigen-presenting cell (APC).
The derivation of $E_{\text{max}}$ and $EC_{50}$ involves several steps, starting with the assumption of steady-state conditions, where the differential equations are set to zero. The authors define the T-cell response level $R$ in terms of triggered T-cell receptors and derive expressions for $E_{\text{max}}$ as the maximal response at ligand saturation and for $EC_{50}$ as the ligand concentration that produces half-maximal response. The section outlines the mathematical procedures for deriving these expressions, emphasizing the differences in behavior among the various models, such as the occupancy model, kinetic proofreading models, and those incorporating feedback mechanisms. The results are presented in detail, with closed-form expressions provided for many models, contributing to a comprehensive understanding of T-cell activation dynamics.