L10 · 30 دقيقة
FEP وTI
تحولات كيميائية في خريطة ترموديناميكية.
أهداف التعلّم
- صياغة FEP (Zwanzig) وTI كمسارين كيميائيين تخيليين لنفس ΔG.
- تمييز ΔG المطلق عن ΔΔG النسبي ومتى يستحق كل منهما العنقود.
كيمياء تخيلية على خريطة ترموديناميكية
التحول الكيميائي (alchemical) يغيّر جزيئاً إلى آخر — أو إلى عدم — عبر معامل λ بين هاميلتونيين. لأن غيبس دالة حالة، يمكنك أن تذهب من A-bound إلى B-bound عبر مسار لا يحدث في المختبر، ما دامت النهايات فيزيائية والمعاينة على المسار كافية. FEP يقدّر الفرق كمتوسط أسي لفروق الكمون. TI يُكامل ⟨∂U/∂λ⟩ على λ. الاثنان يتفقان في حد المعاينة الكاملة؛ يختلفان في الاستقرار العددي حين يتداخل الهاميلتونيان قليلاً. النسبية (A→B في الجيب وفي الماء، دورة مغلقة) تلغي كثيراً من خطأ المجال إن كان A وB قريبين. المطلقة (ليجند → عدم) أصعب: نهاية الاختفاء تتطلب قيوداً وهيئات decoupling بعناية.
Zwanzig FEP
فرق الطاقة الحرة من متوسط أسي لفروق الكمون عيّنت في الحالة المرجعية 0.
- \Delta UU₁ − U₀ عند نفس الهيئة · kcal mol⁻¹
- \beta1/k_B T · kcal⁻¹ mol
- \langle \cdot \rangle_{0}متوسط على ensemble الحالة 0
المتوسط الأسي تهيمن عليه ذيول نادرة. إن لم تتداخل هيئات 0 مع 1، فالرقم متحيّز ولو طال المسار. لذلك نُقسّم λ إلى نوافذ، أو نستخدم BAR/MBAR لا اتجاهاً واحداً ساذجاً.
تطبيق في التصميم الدوائي: نافذة λ حيث الطاقة تتفجر تعني أن الخطوة في الفضاء الكيميائي أكبر من المعاينة. زد النوافذ أو ليّن soft-core لـ LJ.
التكامل الثرموديناميكي TI
ΔG تكامل متوسط مشتق الهاميلتوني على مسار λ.
- \lambdaمعامل التحول بين الهاميلتونيين
- U(\lambda)الكمون الممزوج · kcal mol⁻¹
المنحنى ⟨∂U/∂λ⟩ يجب أن يكون سلساً بما يكفي للتكامل. قمة حادة عند اختفاء vdW شائعة دون soft-core. TI لا يُلغي حاجة المعاينة؛ يُغيّر شكل المقدِّر.
تطبيق في التصميم الدوائي: ارسم ⟨∂U/∂λ⟩ مقابل λ قبل أن تثق بالتكامل. الفجوات بين النقاط ليست ديكوراً.
| FEP / BAR / MBAR | TI | |
|---|---|---|
| المقدر | متوسطات أسية أو نسب احتمال بين نوافذ. | تكامل مشتق. |
| حين يزعجك | تداخل ضعيف بين نوافذ λ. | منحنى ∂U/∂λ غير أملس أو نقاط قليلة. |
| في الممارسة | MBAR شائع كتحليل لاحق لنفس المعاينة. | ما زال قياساً واضحاً إن كانت النوافذ كثيفة. |
لماذا؟ — لماذا يستحق FEP النسبي أسبوع عنقود لسلسة متقاربة بينما لا يستحقه docking؟
لأن السؤال أصبح ΔΔG بحجم 1 kcal/mol داخل هيكل مشترك وهيئة مشتركة — وهو نطاق يمكن أن تصل إليه دورة كيميائية مع مجال حديث ومعاينة نوافذ. docking لا يملك حالة unbound ولا متوسطاً حرارياً. التكلفة تُشترى فقط إن كانت الكيمياء تسمح بالإلغاء (تشابه عالٍ، لا تغيّر شحنة إن أمكن، لا قلب هيئة).
ماذا لو؟ — ماذا لو حوّلت حلقة فينيل إلى بيريدين في الجيب بينما الهيئة تدور 180° في منتصف λ؟
المسار لم يعد تحويلاً بين حالتين متماثلتين البنى؛ هو مسار يخلط كيمياء بانقلاب هيئة. ΔΔG يخلط السؤالين. عالج الانقلاب أولاً (قيود، نوافذ إضافية، أو هيئتان نهائيتان منفصلتان) أو اعترف أن التحويل غير صالح.
تمرين
تريد معرفة أثر إضافة ميثيل في موقع ميتا. صف دورة نسبية واذكر علامتين أن المعاينة فشلت.
مواضيع مرتبطة