I now have rich content from multiple textbooks. Let me compile a thorough, well-grounded answer.
Here are model ideas drawn directly from renal physiology, organized by the system or mechanism each model would represent:
1. Nephron Filtration Barrier Model
What to model: The three-layer glomerular filtration barrier - fenestrated capillary endothelium (pores 70-90 nm), the glomerular basement membrane (GBM), and the podocyte filtration slits (~25 nm wide).
Model concepts:
- Size selectivity: free passage of neutral solutes <4 nm, near-total exclusion >8 nm
- Charge selectivity: the GBM and slit diaphragm are negatively charged, repelling anionic molecules (albumin)
- Mesangial cell contractility modulating filtration surface area
- ~0.8 m² total filtration area in a human kidney
Modeling approaches: A layered membrane model with variable pore sizes and charge distribution; could be implemented as a computational sieve or as a physical 3D-printed layered structure for teaching.
2. GFR Autoregulation Model (Myogenic + Tubuloglomerular Feedback)
What to model: The kidney maintains GFR nearly constant over a mean arterial pressure range of roughly 80-180 mmHg through two mechanisms:
- Myogenic response: afferent arteriole stretches with increased pressure → smooth muscle contracts → raises resistance → limits flow
- Tubuloglomerular feedback (TGF): macula densa senses NaCl delivery → signals juxtaglomerular cells → adjusts afferent arteriole tone
Starling forces at the glomerulus:
- GFR = Kf × [(Pgc - Pbs) - (πgc - πbs)]
- Where Kf = filtration coefficient, Pgc = glomerular capillary hydrostatic pressure, Pbs = Bowman's space pressure, π = oncotic pressures
Modeling approaches: Differential equation or feedback control system (control theory model); the two feedback loops can be modeled as parallel proportional controllers. Good for demonstrating how NSAIDS (block prostaglandin-mediated afferent dilation) or ACE inhibitors alter GFR.
3. Countercurrent Multiplier and Exchanger Model
What to model: The loop of Henle builds a medullary osmotic gradient (300 mOsm/kg at cortex → ~1200 mOsm/kg at papilla):
- Thin descending limb: freely permeable to water, impermeable to Na⁺/Cl⁻ → tubular fluid becomes concentrated
- Thin/thick ascending limb: impermeable to water; actively transports Na⁺, K⁺, Cl⁻ out (NKCC2 in thick limb) → dilutes tubular fluid, concentrates interstitium
- Countercurrent multiplier effect: the hairpin arrangement amplifies a small single-effect into a large axial gradient
- Vasa recta: act as countercurrent exchangers - descending limb picks up solute and loses water, ascending limb returns solute and gains water, preserving the gradient
Modeling approaches: A classic compartmental or finite-difference numerical model with transport coefficients; you can model what happens when loop length changes (short loops = less concentration ability), or when NKCC2 is blocked (furosemide), or when ADH is absent (diabetes insipidus - dilute urine produced).
4. Tubular Transport Maximum (Tm) and Glucose Reabsorption Model
What to model: The proximal tubule reabsorbs glucose via SGLT2 (Na⁺-glucose cotransporter) with a transport maximum (Tm ~375 mg/min in humans). Below plasma glucose ~180 mg/dL, all filtered glucose is reabsorbed. Above the threshold, glucosuria appears.
Key features:
- Splay in the titration curve (due to nephron heterogeneity)
- SGLT2 inhibitors (gliflozins) reduce Tm and induce intentional glucosuria
Modeling approaches: A saturation kinetics (Michaelis-Menten) model; vary filtered load (GFR × plasma [glucose]) against Tm; model the splay with a distribution of Tm values across nephrons.
5. ADH-Aquaporin Urine Concentration Model
What to model: How ADH (vasopressin) acting on collecting duct V2 receptors → cAMP → PKA phosphorylation → vesicular trafficking of aquaporin-2 (AQP2) to the apical membrane → increased water permeability → concentrated urine.
Key states to model:
- No ADH: dilute urine (~50 mOsm/kg), high urine volume - central/nephrogenic diabetes insipidus
- Max ADH: urine up to ~1200 mOsm/kg in humans
- The gradient established by the countercurrent multiplier is the "battery" that ADH "switches on"
Modeling approaches: Signal transduction cascade model (kinetic equations for cAMP, PKA, AQP2 trafficking); or a simpler osmotic flow model coupling collecting duct water permeability to the medullary gradient.
6. RAAS Feedback Model (Renin-Angiotensin-Aldosterone System)
What to model: A classic negative-feedback pressure/volume control loop:
- Decreased renal perfusion pressure → JG cells release renin → angiotensinogen → Ang I → Ang II (via ACE)
- Ang II: efferent arteriole vasoconstriction (maintains GFR), aldosterone secretion, Na⁺ reabsorption in collecting duct, thirst, ADH release
- Volume/pressure restoration feeds back to suppress renin
Clinically powerful: Renal artery stenosis (one-kidney Goldblatt model) creates a decoupled loop - the stenosed kidney keeps secreting renin despite systemic hypertension because its local pressure is low. ACE inhibitors interrupt the loop at a defined node.
Modeling approaches: A nonlinear ODE feedback model with inputs (renal perfusion pressure, Na⁺ intake) and outputs (GFR, blood pressure, aldosterone); this is the foundation of Guyton's classic cardiovascular-renal model of long-term blood pressure control.
7. Acid-Base Buffering and Bicarbonate Reabsorption Model
What to model: The kidney regulates plasma HCO₃⁻ and excretes acid:
- Proximal tubule reclaims ~85% of filtered HCO₃⁻ via carbonic anhydrase (H⁺ secreted, combines with HCO₃⁻ in lumen to form H₂CO₃ → CO₂ + H₂O, then CO₂ re-enters cell)
- Collecting duct type A intercalated cells secrete H⁺ via H⁺-ATPase and H⁺/K⁺-ATPase; titratable acid and NH₄⁺ are the final urine buffers
Modeling approaches: A Henderson-Hasselbalch steady-state model linking tubular H⁺ secretion rates to plasma pH; model responses to metabolic acidosis/alkalosis; excellent for showing why carbonic anhydrase inhibitors (acetazolamide) cause proximal RTA.
8. Potassium Secretion Model in the Collecting Duct
What to model: The principal cell of the collecting duct is the main site of K⁺ secretion. It is driven by:
- High intracellular [K⁺] from Na⁺/K⁺-ATPase on basolateral side
- Apical ROMK channels allowing K⁺ to exit into lumen
- Na⁺ reabsorption (via ENaC) creating a lumen-negative potential that drives K⁺ secretion
- Aldosterone upregulates both ENaC and Na⁺/K⁺-ATPase → more K⁺ loss
Modeling approaches: An electrophysiological/electrochemical model with membrane potentials and channel conductances; model hyperkalemia in Addison's disease (no aldosterone) vs. hypokalemia in primary hyperaldosteronism.
Summary Table
| Model | Core Mechanism | Key Parameters |
|---|
| Filtration barrier | Size + charge selectivity | Pore diameter, charge density, Kf |
| GFR autoregulation | Myogenic + TGF feedback | MAP, afferent/efferent tone, oncotic pressure |
| Countercurrent multiplier | Na⁺ gradient amplification | Loop length, NKCC2 activity, ADH |
| Glucose Tm | Saturable SGLT2 transport | Plasma glucose, GFR, Tm, splay |
| ADH-AQP2 | Signal-driven water permeability | cAMP, AQP2 density, medullary gradient |
| RAAS loop | Pressure-volume feedback | Renin, Ang II, aldosterone, Na⁺ balance |
| Acid-base | Carbonic anhydrase + H⁺ secretion | HCO₃⁻ filtered load, NH₄⁺ production |
| K⁺ secretion | ENaC/ROMK electrochemical gradient | Aldosterone, luminal flow rate, plasma [K⁺] |
Sources: Ganong's Review of Medical Physiology (26th ed.), Brenner and Rector's The Kidney, Costanzo Physiology (7th ed.), Histology: A Text and Atlas (Pawlina).
Would you like me to go deeper on any of these - for example, the mathematical formulation, the specific transport proteins, or a clinical disease model built around one of these systems?