PhysioEng 3D: First-Principles Physiology for Engineers | Chapter 1: Biochemical Kinetics & Mass Action

1.1 Elementary Chemical Reactions & Activation Energy Barrier

First-Principles Reaction Kinetics

रासायनिक अभिक्रियाओं (Chemical Reactions) को मैकेनिकल इंजीनियरिंग के Potential Energy Landscape और Thermal Collisions के रूप में समझें। किसी भी बाइंडिंग या बॉन्ड टूटने के लिए अणुओं (Molecules) को एक थ्रेशोल्ड Activation Energy ($\Delta H^\ddagger$) बैरियर पार करना होता है। एंजाइम (Enzyme) इस बैरियर की ऊँचाई को घटाकर एक वैकल्पिक मार्ग प्रदान करता है, जिससे कम ऊर्जा पर भी टकराव सफल हो जाते हैं।

3D Molecular Collision Chamber

Substrate [A]

Substrate [B]

Product [C]

Enzyme [E] (Catalytic Site)

मैकेनिकल एनालॉजी: चैंबर में गति कर रहे अणु पिस्टन-सिलेंडर में गैस के कणों जैसे हैं। जब पर्याप्त गतिज ऊर्जा (Kinetic Energy) के साथ $[A]$ और $[B]$ टकराते हैं, तो वे ऊर्जा बैरियर को पार कर उत्पाद $[C]$ में परिवर्तित हो जाते हैं। एंजाइम की उपस्थिति में यह प्रभावकारी टकराव कम गति पर भी घटित होता है।

Potential Energy Barrier ($\Delta H^\ddagger$ Hill)

Thermodynamic State

Activation Energy ($\Delta H^\ddagger$)

68.5 kJ/mol

Catalytic Barrier

Free Energy Change ($\Delta G^\circ$)

-14.2 kJ/mol

Thermodynamic Affinity

Collision Efficiency

8.4%

Fraction Over Barrier

Active Product Count [C]

Synthesized Molecules

Arrhenius Rate Law (First-Principles Equation)

k = A \exp\left(-\frac{\Delta H^\ddagger}{R T}\right)

एंजाइम केवल एक्टिवेशन बैरियर $\Delta H^\ddagger$ को कम करता है, किंतु शुरुआत और अंतिम उत्पादों के बीच कुल फ्री एनर्जी $\Delta G$ अपरिवर्तित रहती है।

Law of Mass Action Dynamic Time-Series Simulator

$A + B \underset{k_{-}}{\stackrel{k_{+}}{\rightleftharpoons}} C$

Kinetics Control Inputs

Forward Rate Constant ($k_+$) 0.40 mM⁻¹s⁻¹

Reverse Rate Constant ($k_-$) 0.10 s⁻¹

Initial Substrate $[A]_0$ 2.00 mM

Initial Substrate $[B]_0$ 1.50 mM

Equilibrium Constant ($K = \frac{k_-}{k_+}$): 0.250 mM

Equilibrium Product $[\bar{c}]$: 1.186 mM

Mass Conservation ($[A] + [C] = [A]_0$): Verified

क्षणिक प्रतिक्रिया (Transient Response): प्रारंभिक अवस्था से सिस्टम समय के साथ गतिशील संतुलन (Dynamic Equilibrium) की ओर अग्रसर होता है जहाँ $\frac{d[A]}{dt} = 0$, और अनुपात $\frac{\bar{a}\bar{b}}{\bar{c}} = \frac{k_-}{k_+}$ स्थापित हो जाता है।

1.2 Michaelis-Menten Kinetics: Equilibrium vs. Briggs-Haldane QSSA

Enzymatic Rate Formulations

एंजाइम-सब्सट्रेट बाइंडिंग $E + S \underset{k_{-1}}{\stackrel{k_{+1}}{\rightleftharpoons}} ES \xrightarrow{k_{+2}} E + P$ के दो प्रमुख शास्त्रीय समाधान हैं: Equilibrium Approximation (यह मानकर कि $ES$ कॉम्प्लेक्स बहुत तेज़ी से संतुलन में आ जाता है और $k_{+2} \ll k_{-1}$) तथा Briggs-Haldane Quasi-Steady-State Approximation (QSSA) (जहाँ $\frac{d[ES]}{dt} \approx 0$).

1. Equilibrium Approximation ($K_s$)

Henri-Michaelis-Menten

V = \frac{V_{\max} s}{K_s + s}, \quad \text{where } K_s = \frac{k_{-1}}{k_{+1}}

यहाँ माना जाता है कि सब्सट्रेट का जुड़ना और अलग होना ($E + S \rightleftharpoons ES$) उत्पाद बनने ($ES \to E + P$) की तुलना में अत्यधिक तीव्र गति से होता है, अतः $K_s$ शुद्ध पृथक्करण स्थिरांक (Dissociation Constant) है।

2. Briggs-Haldane QSSA ($K_m$)

Rigorous Engineering Formulation

V = \frac{V_{\max} s}{K_m + s}, \quad \text{where } K_m = \frac{k_{-1} + k_{+2}}{k_{+1}} = K_s + \frac{k_{+2}}{k_{+1}}

इंजीनियरिंग में यह अधिक सटीक है क्योंकि यह कॉम्प्लेक्स के टूटने के दोनों रास्तों ($k_{-1}$ वापस और $k_{+2}$ उत्पाद की ओर) का हिसाब रखता है। चूँकि $k_{+2} > 0$, इसलिए हमेशा $K_m > K_s$ होता है।

Dual Velocity vs. Substrate Concentration Curve

Live Comparison

ग्राफ अवलोकन: नीली वक्र रेखा (Equilibrium) हमेशा एम्बर वक्र रेखा (QSSA) से ऊपर रहती है, क्योंकि $K_m > K_s$ होने से हर परिमित $s$ पर QSSA का हर (denominator) बड़ा होता है।

Kinetic Parameter Controls

Live Parameter Identification

Maximum Velocity ($V_{\max}$) 1.00 mM/s

Equilibrium Constant ($K_s$) 1.20 mM

Turnover Extra Ratio ($\frac{k_{+2}}{k_{+1}}$) 0.80 mM

गणितीय प्रमाण: $V_{\text{QSSA}} < V_{\text{Equil}}$

K_m = K_s + \frac{k_{+2}}{k_{+1}} > K_s \implies K_m + s > K_s + s$$ $$\implies \frac{V_{\max}s}{K_m + s} < \frac{V_{\max}s}{K_s + s} \quad (\forall s > 0)

अतः उत्पाद बनने की दर ($k_{+2}$) जितना अधिक होगी, कॉम्प्लेक्स उतनी ही तेज़ी से खपेगा, और वास्तविक प्रतिक्रिया वेग (Reaction Velocity) संतुलन अनुमान से उतना ही कम होगा।

Lineweaver-Burk Double-Reciprocal Transformation

Linear Parameter Extraction

Governing Linearized Equation

\frac{1}{V} = \frac{K_m}{V_{\max}}\left(\frac{1}{s}\right) + \frac{1}{V_{\max}}

y-intercept ($1/s = 0$)

1.000 (mM/s)⁻¹

$\frac{1}{V_{\max}}$ (Yields $V_{\max}$)

x-intercept ($1/V = 0$)

-0.500 mM⁻¹

$-\frac{1}{K_m}$ (Yields $K_m$)

Line Slope

2.000 s

$\frac{K_m}{V_{\max}}$

Identified Parameters

$V_{\max}=1.0, K_m=2.0$

Direct Extraction

इंजीनियरिंग उपयोग: प्रयोगशाला में गैर-रैखिक (Non-linear) अतिपरवलयिक वक्र से सीधे $V_{\max}$ और $K_m$ निकालना कठिन होता है। Lineweaver-Burk प्लॉट डेटा को एक सीधी रेखा में बदलकर एक्सिस इंटरसेप्ट्स से सीधे मापदंडों का निर्धारण संभव बनाता है।

1.3 Hill Cooperativity & Oxygen-Hemoglobin Binding Switch

Allosteric Transition & Sigmoidal Dynamics

जब एक एंजाइम या ट्रांसपोर्ट प्रोटीन में एक से अधिक बाइंडिंग साइट्स होती हैं और एक अणु के जुड़ने से बाकी साइट्स की आकर्षण क्षमता (Affinity) बढ़ जाती है, तो इसे Positive Cooperativity कहा जाता है। हीमोग्लोबिन (Hemoglobin) में चार सबयूनिट्स ($n \approx 2.5$ से $2.8$) फेफड़ों ($pO_2 \approx 100\text{ mmHg}$) में तेजी से ऑक्सीजन लोड करने और मांसपेशियों के ऊतकों ($pO_2 \approx 20\text{ mmHg}$) में तुरंत अनलोड करने का एक सटीक ऑल-ऑर-नन स्विच बनाते हैं।

The Hill Cooperativity Simulator

Physiological Switch

Hill Coefficient ($n$) 2.50

Half-Saturation Constant ($K$) 1.50 mM

Physiological Transport Dynamics

Pulmonary vs Peripheral Delivery

V = \frac{V_{\max} s^n}{K^n + s^n} \quad \implies \quad \theta = \frac{s^n}{K^n + s^n}

Alveolar Loading ($s = 3.5\text{ mM}$)

91.8%

Oxygen Saturation ($\theta_{\text{lungs}}$)

Tissue Capillary Unloading ($s = 0.8\text{ mM}$)

20.4%

Oxygen Saturation ($\theta_{\text{tissue}}$)

Delivered Oxygen Differential ($\Delta \theta$)

71.4% Delivered

Net Physiological Oxygen Release to Tissues

यदि हीमोग्लोबिन में सहकारिता (Cooperativity) न होती ($n=1.0$), तो ऊतकों में पर्याप्त ऑक्सीजन रिलीज करने के लिए शरीर को अत्यधिक निम्न पार्शियल प्रेशर की आवश्यकता होती, जिससे ऊतकों में हाइपोक्सिया (Hypoxia) हो जाता। $n \approx 2.5$ होने के कारण थोड़ा सा प्रेशर ड्रॉप होते ही हीमोग्लोबिन भारी मात्रा में ऑक्सीजन छोड़ देता है।

Exercise E: Hill Linearization & Experimental Dataset Regression

$\ln\left(\frac{V}{V_{\max}-V}\right) = n \ln s - n \ln K$

Experimental Data Matrix (Chapter 1, Exercise E)

Point	Substrate $s$ (mM)	Velocity $V$ (mM/s)	$\ln(s)$	$\ln[V / (V_{\max}-V)]$

Linear Slope ($n$): 2.394 (Textbook Fit: $n = 2.394$)

Intercept ($-n \ln K$): -0.971 ($\ln K = 0.405$)

Computed $K$: 1.500 mM (Textbook Fit: $K = 1.50\text{ mM}$)

Regression Fit Quality ($R^2$): 0.9998 High-precision fit

1.4 Enzyme Inhibition Studio: Competitive vs. Allosteric/Uncompetitive

Pharmacodynamics & Drug Reversal Mechanics

दवाओं (Pharmaceutical Inhibitors) का एंजाइम्स पर प्रभाव दो मौलिक तरीकों से होता है: Competitive Inhibition: अवरोधक (Inhibitor) सीधे एक्टिव साइट पर सब्सट्रेट से मुकाबला करता है। अत्यधिक सब्सट्रेट डालकर इसे रिवर्स किया जा सकता है (Apparent $K_m$ बढ़ता है, $V_{\max}$ स्थिर रहता है)। Allosteric / Uncompetitive Inhibition: अवरोधक एंजाइम-सब्सट्रेट कॉम्प्लेक्स ($ES$) से जुड़ता है, जिससे उत्प्रेरक क्षमता घट जाती है (Apparent $V_{\max}$ घट जाता है)।

Competitive Inhibition Pathway

Surmountable with Substrate

$E + S \rightleftharpoons ES \to E + P$

parallel with

$E + I \rightleftharpoons EI \quad (\text{Dead-End})$

V = \frac{V_{\max} s}{K_m\left(1 + \frac{[I]}{K_i}\right) + s} = \frac{V_{\max} s}{K_m^{\text{app}} + s}

प्रभाव: $K_m^{\text{app}} = K_m(1 + [I]/K_i)$ बढ़ जाता है। सब्सट्रेट सांद्रता $s \to \infty$ बढ़ाने पर अधिकतम गति $V_{\max}$ प्राप्त की जा सकती है।

Allosteric / Uncompetitive Pathway

Insurmountable by Substrate

$E + S \rightleftharpoons ES \to E + P$

coupled with

$ES + I \rightleftharpoons ESI \quad (\text{Inactivated})$

V = \frac{\left(\frac{V_{\max}}{1 + [I]/K_i}\right) s}{\left(\frac{K_m}{1 + [I]/K_i}\right) + s} = \frac{V_{\max}^{\text{app}} s}{K_m^{\text{app}} + s}

प्रभाव: $V_{\max}^{\text{app}}$ घट जाता है। आप चाहे कितना भी अधिक सब्सट्रेट डालें, अधिकतम गति पूर्व स्तर पर नहीं पहुँच सकती।

Interactive Drug Inhibition Sandbox

Live Parameter Dynamics

Select Inhibition Mechanism:

Competitive Allosteric / Uncompetitive

Inhibitor Concentration ($[I]$) 1.50 mM

Inhibition Constant ($K_i$) 1.00 mM

Inhibition Factor $(1 + [I]/K_i)$: 2.50

Apparent $K_m^{\text{app}}$: 3.75 mM

Apparent $V_{\max}^{\text{app}}$: 1.00 mM/s

Michaelis-Menten Velocity Profile

Lineweaver-Burk Double-Reciprocal Transformation

मैकेनिकल डायग्नोस्टिक संकेत: Competitive Inhibition में Lineweaver-Burk रेखा y-अक्ष ($\frac{1}{V_{\max}}$) के चारों ओर घूमती है (Slope बढ़ता है किंतु y-intercept स्थिर रहता है)। इसके विपरीत, Uncompetitive Inhibition में रेखा समानांतर (Parallel shift) खिसकती है क्योंकि Slope ($\frac{K_m}{V_{\max}}$) अपरिवर्तित रहता है जबकि दोनों अक्षीय इंटरसेप्ट्स बदल जाते हैं।

1.5 Complete Engineering Lecture Notes & Rigorous Exercise Derivations

Textbook-Grade First-Principles Reference

यह खंड Michael Chappell & Stephen Payne की टेक्स्टबुक के Chapter 1 के मुख्य वैचारिक नोट्स और सभी छह अभ्यासों (Exercises A through F) के संपूर्ण, चरण-दर-चरण गणितीय समाधान प्रदान करता है।

1. The Cell as a Thermodynamic Engineering System

Biophysical Foundations

एक बायोलॉजिकल सेल को इलेक्ट्रिकल और केमिकल इंजीनियरिंग के बुनियादी सिद्धांतों के संयोजन के रूप में देखा जा सकता है:

1. Lipid Dielectric Membrane

फॉस्फोलिपिड बाइलेयर (Phospholipid Bilayer) लगभग 3 से 4 नैनोमीटर मोटाई की एक गैर-चालक (Insulating) परत है। यह एक उच्च-गुणवत्ता वाले Dielectric Capacitor की तरह कार्य करती है ($C_m \approx 1\,\mu\text{F/cm}^2$), जो बाह्य और आंतरिक आयनिक घोल को अलग रखती है।

2. Cytoplasmic Electrolyte Solution

कोशिका का आंतरिक द्रव (Cytoplasm) उच्च $[K^+]$ और कम $[Na^+]$, $[Ca^{2+}]$ वाला एक कंडक्टिव आयनिक माध्यम है। यह सर्किट में रेसिस्टिव पाथ (Resistive Fluid Conduit) का निर्माण करता है।

3. Transmembrane Channels & Pumps

आयन चैनल्स वोल्टेज या लिगैंड-नियंत्रित Variable Conductances (Resistors) हैं। वहीं $Na^+/K^+$ ATP-ase पंप एक निरंतर चालू Active Current Source है जो केमिकल ऊर्जा को खर्च करके विद्युत-रासायनिक प्रवणता (Electrochemical Gradient) को बनाए रखता है।

4. ATP Energy Storage & Reaction Kinetics

सेलुलर ईंधन का दहन एरोबिक मेटाबॉलिज्म द्वारा होता है:

C_6H_{12}O_6 + 6O_2 \longrightarrow 6CO_2 + 6H_2O + 36\,\text{ATP}

ATP का हाइड्रोलिसिस ($\text{ATP} \to \text{ADP} + P_i$) लगभग $-30.5\text{ kJ/mol}$ की फ्री ऊर्जा छोड़ता है, जो जैविक मोटरों और पंपों को ऊर्जा प्रदान करता है।

Exercise A

Stoichiometric Dynamic Mass Action: $A + 2B \underset{k_-}{\stackrel{k_+}{\rightleftharpoons}} C + D$

Problem Statement: Write the four Ordinary Differential Equations (ODEs) describing the rates of change of species concentrations for the elementary reaction $A + 2B \underset{k_-}{\stackrel{k_+}{\rightleftharpoons}} C + D$. At dynamic equilibrium, show that the thermodynamic equilibrium constant is given by $K = \frac{\bar{a}\bar{b}^2}{\bar{c}\bar{d}} = \frac{k_-}{k_+}$.

Step 1: Stoichiometric Reaction Rates

द्रव्य-क्रिया के नियम (Law of Mass Action) के अनुसार, अग्रगामी दर (Forward Rate) अभिकारकों की सांद्रता के स्टॉइकियोमेट्रिक घात (Power) के समानुपाती होती है:

r_+ = k_+ a b^2, \qquad r_- = k_- c d

Step 2: 4 Coupled Governing ODEs

\frac{da}{dt} = -k_+ a b^2 + k_- c d$$ $$\frac{db}{dt} = -2\left(k_+ a b^2 - k_- c d\right) = -2 k_+ a b^2 + 2 k_- c d$$ $$\frac{dc}{dt} = k_+ a b^2 - k_- c d$$ $$\frac{dd}{dt} = k_+ a b^2 - k_- c d

Step 3: Dynamic Equilibrium Condition

संतुलन की स्थिति में सांद्रता में परिवर्तन की दर शून्य हो जाती है ($\frac{da}{dt} = \frac{db}{dt} = \frac{dc}{dt} = \frac{dd}{dt} = 0$). अतः:

k_+ \bar{a} \bar{b}^2 = k_- \bar{c} \bar{d} \implies \frac{\bar{a} \bar{b}^2}{\bar{c} \bar{d}} = \frac{k_-}{k_+} \equiv K \quad \blacksquare

Exercise B

Coupled Two-Step Consecutive Reactions

Problem Statement: Consider the coupled reactions $A + B \underset{k_{-1}}{\stackrel{k_{+1}}{\rightleftharpoons}} C$ and $C + D \underset{k_{-2}}{\stackrel{k_{+2}}{\rightleftharpoons}} E$. Write the system of ODEs for all species and derive the overall equilibrium constant $K$ relating initial substrates $A, B, D$ to final product $E$.

Step 1: Species Dynamic Conservation ODEs

\frac{da}{dt} = -k_{+1} a b + k_{-1} c$$ $$\frac{db}{dt} = -k_{+1} a b + k_{-1} c$$ $$\frac{dc}{dt} = k_{+1} a b - k_{-1} c - k_{+2} c d + k_{-2} e$$ $$\frac{dd}{dt} = -k_{+2} c d + k_{-2} e$$ $$\frac{de}{dt} = k_{+2} c d - k_{-2} e

Step 2: Equilibrium Coupling

संतुलन पर प्रत्येक प्राथमिक चरण अलग-अलग संतुलित होता है (Principle of Microscopic Reversibility):

\text{Step 1: } k_{+1}\bar{a}\bar{b} = k_{-1}\bar{c} \implies \bar{c} = \frac{k_{+1}}{k_{-1}}\bar{a}\bar{b}$$ $$\text{Step 2: } k_{+2}\bar{c}\bar{d} = k_{-2}\bar{e} \implies \bar{e} = \frac{k_{+2}}{k_{-2}}\bar{c}\bar{d}

Step 3: Overall Equilibrium Constant

$\bar{c}$ का मान प्रतिस्थापित करने पर:

\bar{e} = \frac{k_{+2}}{k_{-2}}\left(\frac{k_{+1}}{k_{-1}}\bar{a}\bar{b}\right)\bar{d} = \frac{k_{+1}k_{+2}}{k_{-1}k_{-2}}\bar{a}\bar{b}\bar{d}$$ $$\implies K = \frac{\bar{a}\bar{b}\bar{d}}{\bar{e}} = \frac{k_{-1}k_{-2}}{k_{+1}k_{+2}} = K_1 K_2 \quad \blacksquare

Exercise C

Rigorous Derivation of Both Michaelis-Menten Models & Proof that $V_{\text{QSSA}} < V_{\text{Equil}}$

Problem Statement: Derive the Henri-Michaelis-Menten reaction velocity $V$ using both (i) the rapid equilibrium assumption ($K_s = k_{-1}/k_{+1}$), and (ii) the Briggs-Haldane quasi-steady-state approximation ($K_m = (k_{-1}+k_{+2})/k_{+1}$). Prove mathematically that for all substrate concentrations $s > 0$, $V_{\text{QSSA}} < V_{\text{Equil}}$.

Part 1: Equilibrium Assumption

मान लें $E + S \rightleftharpoons ES$ निरंतर संतुलन में है और कुल एंजाइम $e_0 = e + c$ संरक्षित है:

K_s = \frac{e \cdot s}{c} = \frac{(e_0 - c)s}{c} \implies c\left(K_s + s\right) = e_0 s \implies c = \frac{e_0 s}{K_s + s}$$ $$V = k_{+2} c = \frac{k_{+2} e_0 s}{K_s + s} = \frac{V_{\max} s}{K_s + s} \quad (\text{where } V_{\max} = k_{+2} e_0)

Part 2: Briggs-Haldane Quasi-Steady-State Approximation (QSSA)

यहाँ माना जाता है कि कॉम्प्लेक्स निर्माण के बाद उसकी सांद्रता में बदलाव बहुत धीमा होता है ($\frac{dc}{dt} \approx 0$):

\frac{dc}{dt} = k_{+1} e \cdot s - (k_{-1} + k_{+2}) c = 0$$ $$k_{+1} (e_0 - c) s - (k_{-1} + k_{+2}) c = 0 \implies c\left[k_{-1} + k_{+2} + k_{+1} s\right] = k_{+1} e_0 s$$ $$c = \frac{e_0 s}{\frac{k_{-1}+k_{+2}}{k_{+1}} + s} = \frac{e_0 s}{K_m + s}, \quad \text{where } K_m = \frac{k_{-1}+k_{+2}}{k_{+1}}$$ $$V = k_{+2} c = \frac{V_{\max} s}{K_m + s} \quad \blacksquare

Part 3: Mathematical Proof that $V_{\text{QSSA}} < V_{\text{Equil}}$

K_m = \frac{k_{-1} + k_{+2}}{k_{+1}} = \frac{k_{-1}}{k_{+1}} + \frac{k_{+2}}{k_{+1}} = K_s + \frac{k_{+2}}{k_{+1}}

चूँकि उत्प्रेरक दर $k_{+2} > 0$ तथा $k_{+1} > 0$, इसलिए $\frac{k_{+2}}{k_{+1}} > 0 \implies K_m > K_s$।

K_m > K_s \implies K_m + s > K_s + s \quad (\forall s > 0)$$ $$\implies \frac{1}{K_m + s} < \frac{1}{K_s + s} \implies \frac{V_{\max} s}{K_m + s} < \frac{V_{\max} s}{K_s + s}$$ $$\therefore V_{\text{QSSA}} < V_{\text{Equil}} \quad \blacksquare

Exercise D

Lineweaver-Burk Double-Reciprocal Linearization & Parameter Extraction

Problem Statement: Show how the Michaelis-Menten equation can be rearranged into a linear relationship between $1/V$ and $1/s$. State clearly how the maximum velocity $V_{\max}$ and Michaelis constant $K_m$ are determined from the intercepts and slope of this plot.

Step 1: Reciprocal Algebraic Inversion

माइकलिस-मेन्टन समीकरण के दोनों पक्षों का व्युत्क्रम (Reciprocal) लेने पर:

\frac{1}{V} = \frac{K_m + s}{V_{\max} s} = \frac{K_m}{V_{\max} s} + \frac{s}{V_{\max} s}$$ $$\frac{1}{V} = \left(\frac{K_m}{V_{\max}}\right) \cdot \left(\frac{1}{s}\right) + \frac{1}{V_{\max}}

Step 2: Mapping to Linear Form $y = mx + c$

y = \frac{1}{V}, \quad x = \frac{1}{s}, \quad \text{Slope } m = \frac{K_m}{V_{\max}}, \quad \text{y-intercept } c = \frac{1}{V_{\max}}

Step 3: Axis Intercept Parameter Extraction

\text{y-intercept: Set } x = 0 \implies y = \frac{1}{V_{\max}} \implies V_{\max} = \frac{1}{c}$$ $$\text{x-intercept: Set } y = 0 \implies m x_0 + c = 0 \implies x_0 = -\frac{c}{m} = -\frac{1/V_{\max}}{K_m/V_{\max}} = -\frac{1}{K_m}$$ $$\implies K_m = -\frac{1}{x_0} \quad \blacksquare

Exercise E

Hill Equation Linearization & Experimental Dataset Validation

Problem Statement: Linearize the Hill equation $V = \frac{V_{\max} s^n}{K^n + s^n}$ and describe the resulting plot. Given the experimental measurements $(s, V)$ in Chapter 1 Table 1.1, estimate the parameters $n$ and $K$ assuming $V_{\max} = 1.0\text{ mM/s}$.

Step 1: Algebraic Inversion & Ratio Separation

V(K^n + s^n) = V_{\max} s^n \implies V K^n = (V_{\max} - V) s^n$$ $$\frac{V}{V_{\max} - V} = \frac{s^n}{K^n} = \left(\frac{s}{K}\right)^n

Step 2: Logarithmic Transformation

\ln\left(\frac{V}{V_{\max} - V}\right) = n \ln s - n \ln K

Step 3: Numerical Regression Validation

पाठ्यपुस्तक के टेबल डेटा ($s \in [0.2, 4.5]\text{ mM}$, $V \in [0.01, 0.94]\text{ mM/s}$) पर लीनियर रिग्रेशन करने पर:

\text{Slope } m = n \approx 2.394$$ $$\text{Intercept } c = -n \ln K \approx -0.971 \implies \ln K = \frac{0.971}{2.394} \approx 0.4056$$ $$K = e^{0.4056} \approx 1.500\text{ mM} \quad \blacksquare

Exercise F

Competitive Inhibition Quasi-Steady-State Derivation

Problem Statement: An enzyme reaction is inhibited competitively: $E + S \underset{k_{-1}}{\stackrel{k_{+1}}{\rightleftharpoons}} ES \xrightarrow{k_{+2}} E + P$ alongside $E + I \underset{k_{-i}}{\stackrel{k_{+i}}{\rightleftharpoons}} EI$. Under the quasi-steady-state assumption for $ES$ and equilibrium for $EI$, derive the reaction velocity equation. Show that in a Lineweaver-Burk plot, the line rotates around the y-intercept with increasing inhibitor concentration $[I]$.

Step 1: Mass Balances & Conservation

e_0 = e + c + c_i \quad (\text{where } c = [ES], c_i = [EI])$$ $$K_i = \frac{e \cdot i}{c_i} \implies c_i = \frac{e \cdot i}{K_i}$$ $$e_0 = e\left(1 + \frac{i}{K_i}\right) + c \implies e = \frac{e_0 - c}{1 + \frac{i}{K_i}}

Step 2: Quasi-Steady-State Condition for $ES$

\frac{dc}{dt} = k_{+1} e \cdot s - (k_{-1} + k_{+2}) c = 0$$ $$k_{+1} \left(\frac{e_0 - c}{1 + \frac{i}{K_i}}\right) s = (k_{-1} + k_{+2}) c$$ $$(e_0 - c) s = \frac{k_{-1}+k_{+2}}{k_{+1}} \left(1 + \frac{i}{K_i}\right) c = K_m \left(1 + \frac{i}{K_i}\right) c$$ $$c = \frac{e_0 s}{K_m\left(1 + \frac{i}{K_i}\right) + s}$$ $$V = k_{+2} c = \frac{V_{\max} s}{K_m\left(1 + \frac{i}{K_i}\right) + s} \quad \blacksquare

Step 3: Lineweaver-Burk Intercept Analysis

\frac{1}{V} = \frac{K_m\left(1 + \frac{i}{K_i}\right)}{V_{\max}}\left(\frac{1}{s}\right) + \frac{1}{V_{\max}}$$ $$\text{y-intercept } = \frac{1}{V_{\max}} \quad (\text{Independent of } [I])$$ $$\text{Slope } = \frac{K_m}{V_{\max}}\left(1 + \frac{i}{K_i}\right) \quad (\text{Increases linearly with } [I])$$ $$\text{x-intercept } = -\frac{1}{K_m\left(1 + \frac{i}{K_i}\right)} \quad (\text{Shifts closer to origin as } [I] \uparrow) \quad \blacksquare