Explain mathematical modeling of biological systems

6.5.7 · Biology › Systems Biology & Frontiers

KIYA hai ek mathematical model?

Core object usually ek differential equation hoti hai — kyunki biology zyaadatar rates ke baare mein hai (koi cheez kitni tezi se grow, decay, produce, ya consume hoti hai).

State variable: woh cheez jo hum track karte hain, jaise $N$ = rabbits ki sankhya, $[S]$ = substrate concentration.
Parameter: ek fixed number jo system ko describe karta hai, jaise growth rate $r$ , carrying capacity $K$ .
Rate: $\frac{d N}{d t}$ = state variable kitni tezi se change ho raha hai abhi.

HUM model kaise banate hain (recipe)

State variables list karo — kya change hota hai?
Processes list karo — birth, death, binding, decay, diffusion.
Balance equation likho — har variable ke liye: $\frac{d ( stuff )}{d t} = (rate in) - (rate out)$
Rate laws chunno — jaise mass-action ( $\propto$ reactants ka product), ya saturating (Michaelis–Menten).
Solve / simulate karo, phir data se compare karo (Forecast-then-Verify).

Logistic growth model scratch se derive karna

Step 1 — unlimited growth. Har individual per-capita rate $r$ se reproduce karta hai. $N$ individuals ke saath: $\frac{d N}{d t} = r N$ Yeh step kyun? Naye individuals add hone ki rate $=$ (rate per individual) $\times$ (individuals ki sankhya). Yeh exponential growth deta hai $N (t) = N_{0} e^{r t}$ — hamesha ke liye unrealistic.

Step 2 — ek limit add karo. Carrying capacity $K$ introduce karo. Hum chahte hain ki per-capita rate $0$ tak shrink ho jaise $N \to K$ . Sabse simple factor jo yeh karta hai woh hai $(1 - \frac{N}{K})$ : $\frac{d N}{d t} = r N (1 - \frac{N}{K})$ $(1 - N / K)$ kyun? Jab $N$ bahut chota hai, factor $\approx 1$ (exponential growth). Jab $N = K$ , factor $= 0$ (growth ruk jaati hai). Yeh linearly interpolate karta hai — sabse kam assumption wala choice.

Step 3 — isse solve karo. Variables separate karo: $\int \frac{d N}{N ( 1 - N / K )} = \int r d t$ Partial fractions $\frac{1}{N ( 1 - N / K )} = \frac{1}{N} + \frac{1/ K}{1 - N / K}$ use karke aur integrate karke: $ln \frac{N}{K - N} = r t + C$ $N$ ke liye $N (0) = N_{0}$ ke saath solve karo: $N (t) = \frac{K}{1 + ( \frac{K - N _{0}}{N _{0}} ) e ^{- r t}}$ Yahi famous S-shaped (sigmoid) curve hai.

Figure — Explain mathematical modeling of biological systems

Ek doosra worked example: enzyme kinetics (Michaelis–Menten)

Worked example Saturating rate law derive karna

Kahani: enzyme $E$ substrate $S$ se bind karta hai complex $C$ banane ke liye, jo product $P$ mein badal jaata hai. $E + S k_{- 1} ⇌ k_{1} C k_{2} E + P$

Step 1 — mass action for complex. $\frac{d C}{d t} = k_{1} [E] [S] - (k_{- 1} + k_{2}) C$ Kyun? Complex rate par banta hai jo $E$ aur $S$ ke milne ke proportional hai (mass action), aur do tareekon se toot jaata hai.

Step 2 — quasi-steady-state assumption. Complex level slowly change hota hai: $\frac{d C}{d t} \approx 0$ . Kyun? Enzyme scarce hai; $C$ jaldi se balance reach kar leta hai. Yeh ek modeling simplification hai jo exactness ko solvability ke badle mein trade karta hai.

Step 3 — total enzyme $E_{T} = [E] + C$ ke saath, $C$ ke liye solve karo aur $K_{M} = \frac{k _{- 1} + k _{2}}{k _{1}}$ define karo: $v = k_{2} C = \frac{V _{ma x} [ S ]}{K _{M} + [ S ]}, V_{ma x} = k_{2} E_{T}$ Yeh kyun matter karta hai: rate saturate hoti hai — substrate ko double karna speed ko double nahi karta jab $[S] ≫ K_{M}$ . Pure mass-action endless increase predict karta; model usse correct karta hai.

Teesra example: predator–prey (Lotka–Volterra)

Common mistakes (Steel-manned)

Active recall

Recall Kya tum yeh reproduce kar sakte ho?

Biologists differential equations kyun use karte hain? → woh rates describe karte hain, aur nature tendencies specify karti hai, values nahi.
Factor $(1 - N / K)$ kya karta hai? → growth shut off kar deta hai jaise $N \to K$ .
Kaunsa assumption Michaelis–Menten ko solvable banata hai? → complex par quasi-steady-state.
Ek equilibrium ki stability kya determine karta hai? → $N^{*}$ par rate ke derivative ka sign.

Feynman: 12-saal ke bachche ko samjhao

Recall Mujhe 12 saal ke bachche ki tarah samjhao

Ek fish tank imagine karo. Kuch maachliyan daalo aur woh jaldi bachche deti hain — zyada machhli, zyada bachche, isliye sankhya tezi se badhti hai. Lekin tank mein itni hi machhliyan fit ho sakti hain; khaana aur jagah kam ho jaati hai, isliye bachche slow ho jaate hain jab tak sankhya "tank full" par ruk nahi jaati. Ek math model sirf ek rule hai jo kehta hai "abhi itni machhliyan → agle time mein itni zyada (ya kam)" — rule likho, aur ek computer usse aage play karta hai video game ki tarah bina asli machhliyan kharide future guess karne ke liye.

Connections

Differential Equations — mathematical engine.
Enzyme Kinetics — Michaelis–Menten detail mein.
Population Ecology — logistic & Lotka–Volterra applications.
Systems Biology & Frontiers — parent chapter, networks & emergence.
Stability Analysis — equilibria aur perturbations.
Feedback Loops in Gene Regulation — jahan yeh tools power up hote hain.

Biology mein mathematical model kya hai?

State variables time ke saath kaise change hote hain yeh describe karne wali (usually differential) equations ka ek set, interaction rules ke zariye.

Biology mein differential equations natural kyun hain?

Kyunki biology rates/tendencies specify karti hai (cheezein kitni tezi se change hoti hain), direct values nahi.

Logistic growth ODE likho.

dN/dt = rN(1 − N/K).

(1 − N/K) term kya accomplish karta hai?

Yeh per-capita growth ko vanish kara deta hai jaise N carrying capacity K ke paas aata hai, S-shaped curve deta hai.

Logistic equation ka solution?

N(t) = K / (1 + A e^(−rt)), jahan A = (K − N₀)/N₀.

Logistic growth ke do equilibria kya hain aur unki stability?

N* = 0 (unstable) aur N* = K (stable).

Michaelis–Menten rate law state karo.

v = Vmax[S] / (K_M + [S]).

Kaunsa assumption Michaelis–Menten yield karta hai?

Quasi-steady-state: enzyme–substrate complex ke liye dC/dt ≈ 0.

Rate constants ke terms mein K_M define karo.

K_M = (k₋₁ + k₂)/k₁.

Lotka–Volterra predator–prey equations likho.

dx/dt = αx − βxy ; dy/dt = δxy − γy.

Real populations ke liye exponential growth kyun fail karta hai?

Yeh finite resources ko ignore karta hai; valid sirf jab N ≪ K.

Tum kaise test karte ho ki equilibrium N* stable hai?

N* par d/dN(dN/dt) evaluate karo; negative ⇒ stable, positive ⇒ unstable.