BM
Model ReferenceEvery parameter and where it came from

Provenance

What this platform knows, and how well it knows it.

Every constant in every model appears below with its value, its units, and an honest account of where it came from. The values are imported from the model code rather than retyped, and a test fails the build if a model gains a parameter that is not documented here — so this page cannot quietly drift away from what the simulations actually compute with.

Just over half of these parameters are educational estimates. That is the honest position for a teaching model of this scope, and it is the reason the platform reports directions and magnitudes rather than predictions. Select a class above to filter the tables.

Closed-loop circulation

app/hemodynamics/model.ts

Three compliant compartments exchanging a conserved blood volume, with Sunagawa ventricular-arterial coupling, a Guyton venous-return law, and a proportional baroreflex.

SymbolQuantityValueUnitsProvenanceSource
Ees0End-systolic elastance2.3mmHg/mLmeasuredSuga & Sagawa; human resting values cluster at 2–3 mmHg/mL
V0ESPVR volume intercept5mLestimatedChosen inside the reported 0–20 mL range for the unstressed volume intercept
alphaEDPVR scale1.1mmHgfittedExponential end-diastolic relation fitted to a normal human filling curve
betaEDPVR exponent0.0211/mLfittedSame fit; sets the stiffness of the passive chamber
V0dEDPVR volume intercept30mLestimatedVolume at which passive filling pressure reaches zero
CartArterial compliance1.6mL/mmHgmeasuredTotal systemic arterial compliance in healthy adults, 1.3–2.0 mL/mmHg
CvenSystemic venous compliance90mL/mmHgmeasuredVenous compliance is roughly 30–50× arterial; reported 80–200 mL/mmHg
CcpCardiopulmonary compliance35mL/mmHgestimatedLumped heart-plus-pulmonary compartment; set so resting filling pressure lands near 5.5 mmHg
Ra0Arteriolar resistance17.3Wood unitsmeasuredTotal systemic resistance ≈ 18 Wood units (about 1450 dyn·s·cm⁻⁵), split here into arteriolar and venous components
RvrVenous return resistance1.05Wood unitsfittedGuyton's resistance to venous return; set so venous return equals 5 L/min at the resting gradient
ITP0Resting intrathoracic pressure-4mmHgmeasuredEnd-expiratory intrathoracic pressure, −3 to −5 mmHg
HR0Intrinsic heart rate70bpmestimatedResting adult rate before reflex adjustment
MAPsetBaroreflex set point93mmHgestimatedNormal mean arterial pressure; the reflex operating point
gHRChronotropic reflex gain0.85bpm per mmHgestimatedOpen-loop baroreflex gains of 0.5–1.5 bpm/mmHg are reported in humans
gRVasomotor reflex gain0.6fraction at full errorestimatedChosen to give a physiological pressure response; not separately measured
gEesInotropic reflex gain0.35fraction at full errorestimatedChosen to give a physiological pressure response
gVenoVenoconstriction reflex gain260mL recruited at full errorestimatedReflex mobilisation of unstressed venous volume; order of magnitude only
tauHRChronotropic time constant2smeasuredVagal heart-rate responses settle within a few seconds
tauRVasomotor time constant6smeasuredSympathetic vascular responses develop over 5–15 s
tauEesInotropic time constant3.5sestimatedIntermediate between the vagal and vasomotor arms
tauVenoVenoconstriction time constant8sestimatedSlowest reflex effector in the model
tauFillDiastolic filling time constant90msestimatedSets how incompletely the ventricle fills when diastole is short

Molecule-to-chamber bridge

app/molecular/multiscale-model.ts

Maps myofilament activation onto ventricular mechanics using the same constitutive relations as the circulation model, so the two cannot drift apart.

SymbolQuantityValueUnitsProvenanceSource
Ees0Reference end-systolic elastance2.3mmHg/mLmeasuredShared with the circulation model
V0ESPVR volume intercept5mLestimatedShared with the circulation model
V0dEDPVR volume intercept30mLestimatedShared with the circulation model
alphaEDPVR scale1.1mmHgfittedShared with the circulation model
betaEDPVR exponent0.0211/mLfittedShared with the circulation model
EDV0Reference end-diastolic volume115mLmeasuredNormal adult left ventricular end-diastolic volume, 100–140 mL
SVR0Reference systemic resistance18.35Wood unitsmeasuredMatches Ra0 + Rvr in the circulation model
CartArterial compliance1.6mL/mmHgmeasuredShared with the circulation model
PvenDownstream venous pressure2mmHgestimatedOffset that makes the bridge's reference pressure agree with the circulation model
tauFillDiastolic filling time constant90msestimatedShared with the circulation model
elastanceExponentActivation-to-elastance exponent0.9dimensionlessillustrativeSublinear because elastance saturates as cross-bridges are already recruited; the exponent itself is not measured

Reference myofilament

app/disease-lab/disease-model.ts

The calcium transient and force-calcium relation of the reference cardiomyocyte. Every variant is expressed as a change to these.

SymbolQuantityValueUnitsProvenanceSource
calciumAmplitudePeak systolic calcium increment0.95µMmeasuredIntact ventricular myocytes reach roughly 1 µM at the peak of the transient
diastolicCalciumDiastolic calcium0.1µMmeasuredResting cytosolic calcium, 80–150 nM
calciumKdHalf-activating calcium1.65µMmeasuredpCa50 5.78; measured cardiac force-pCa relations cluster at 5.6–5.9
hillHill coefficient3.2dimensionlessmeasuredCardiac force-pCa Hill coefficients of 3–7 are reported; a lumped thin-filament cooperativity, not a site count
calciumRiseCalcium release time constant0.022smeasuredThe systolic calcium transient rises within tens of milliseconds
calciumDecayCalcium removal time constant0.185smeasuredHuman ventricular calcium decay time constants of 150–250 ms
forceActivationForce redevelopment time constant0.022sestimatedChosen so force follows calcium with a physiological lag
forceRelaxationForce relaxation time constant0.062sestimatedSet so the reference relaxation half-time lands near 120 ms
myosinAvailabilityRecruitable myosin fraction1relativeillustrativeReference value of 1 by definition; variants shift it
stiffnessPassive chamber stiffness1relativeillustrativeReference value of 1 by definition
preloadFilling relative to reference1relativeillustrativeReference value of 1 by definition; variants shift it to represent chamber remodelling
leakAmplitudeDiastolic SR calcium leak0µMillustrativeZero in the reference myofilament; the RyR2 variant raises it. Sized so the diastolic excursion stays in the 0.15–0.30 µM range
energeticFactorEnergetic cost multiplier1relativeillustrativeReference value of 1 by definition; scales the tension-time integral into a demand index
stressTriggerAdrenergic leak sensitivity0.03dimensionlessillustrativeA heuristic index, not a measured quantity
arrhythmiaBiasBaseline arrhythmic substrate0.040–1 indexillustrativeA heuristic index with no electrophysiology behind it

Vascular K_ATP channel

app/pathways/katp-model.ts

Nucleotide gating of the Kir6.1/SUR2B channel and its consequences for smooth-muscle membrane potential and arteriolar calibre.

SymbolQuantityValueUnitsProvenanceSource
atpKiATP half-inhibition90µMestimatedKir6.1-containing channels are markedly less ATP-sensitive than Kir6.2; chosen inside the reported range
atpHillATP Hill coefficient1.2dimensionlessestimatedShallow inhibition consistent with reported ATP dose-responses
adpKaMgADP half-activation180µMestimatedPlaced inside the free MgADP range spanned between perfusion and ischaemia
adpHillMgADP Hill coefficient1.5dimensionlessestimatedCooperative nucleotide activation at the SUR nucleotide-binding domains
basalOpenNucleotide-independent opening0.02fractionmeasuredNear zero: the vascular channel requires Mg-nucleotides to open (Sung et al., PNAS 2021)
maxOpenMaximum open probability0.85fractionestimatedCeiling on single-channel open probability
eKPotassium reversal potential-85mVmeasuredNernst potential for potassium at physiological gradients
vDepolarisedPotential without K_ATP-32mVestimatedWhere the myocyte would sit on its other background currents alone
gBackgroundBackground conductance1relativeillustrativeNormalised to 1; every other conductance in the module is expressed against it
gKatpMaxK_ATP conductance reserve120relative to backgroundestimatedLarge because channel density is high; opening a few percent clamps the membrane near E_K
vHalfL-type activation midpoint-34mVmeasuredSmooth-muscle L-type window current activates over −40 to −20 mV
vSlopeL-type activation slope7mVmeasuredBoltzmann slope factors of 6–9 mV are reported
toneFloorVoltage-independent tone0.2fractionestimatedStretch- and agonist-driven calcium entry that persists when the cell is hyperpolarised
toneCeilingVoltage-dependent tone ceiling1.25relativeillustrativeBounds the Boltzmann so a few millivolts of depolarisation cannot triple tone
referenceConstrictionResting arteriolar constriction0.3fraction of maximum radiusestimatedResistance arterioles hold substantial resting tone
fixedResistanceFractionNon-arteriolar resistance0.35fractionestimatedConduit vessels, capillaries and venules in series, which do not dilate

JAK-STAT signalling

app/pathways/jakstat-model.ts

STAT nucleocytoplasmic cycling on the topology of Swameye et al. (PNAS 2003), with a SOCS feedback arm added as an educational extension.

SymbolQuantityValueUnitsProvenanceSource
k1STAT phosphorylation0.42/min per unit receptorestimatedChosen to reproduce the published time to peak phosphorylation; NOT the value fitted in Swameye et al.
k2pSTAT dimerisation1.6/min per unitestimatedSecond-order dimerisation; chosen for the published timescale
k3Nuclear import0.28/minestimatedChosen so the nuclear signal peaks in the published window
k4Nuclear export0.16/minestimatedThe published fit cannot identify the export rate independently; this is an educational value
tauNuclear residence delay6minestimatedThe delay in the published cycling model; its effect on cytoplasmic pSTAT is small
kOnReceptor activation0.55/min per unit ligandestimatedEducational estimate
kOffReceptor deactivation0.09/minestimatedEducational estimate
basalActivationLigand-independent activation0/minillustrativeZero in the reference pathway; the JAK2 V617F variant raises it
socsSynthesisSOCS synthesis0.075/minestimatedPart of the educational feedback extension, not the fitted model
socsThresholdSOCS transcription threshold0.22nuclear STAT unitsestimatedEducational extension
socsDegradationSOCS degradation0.021/minestimatedEducational extension; sets how long adaptation persists
socsPotencySOCS receptor inhibition0.35SOCS unitsestimatedEducational extension
statTotalTotal STAT pool1relativeillustrativeNormalised to 1
receptorTotalTotal receptor pool1relativeillustrativeNormalised to 1

Complement alternative pathway

app/pathways/complement-model.ts

A self-amplifying convertase loop on one unit of surface, opposed by factor H, CD55, CD46 with factor I, and CD59. On-surface complement kinetics vary by orders of magnitude with surface chemistry, so every rate here is an educational estimate chosen to reproduce the published qualitative behaviour.

SymbolQuantityValueUnitsProvenanceSource
tickoverSpontaneous C3 hydrolysis3.50e-4/minestimatedSized so an untouched host surface turns over only a percent or two of plasma C3 across the window
formationConvertase assembly4.3/min per unit C3bestimatedSet so the loop sits just below instability on a fully regulated cell — the defining property of the real system
catalysisC3 cleavage by convertase2.1/minestimatedThe amplification step; chosen with formation to place the instability threshold
depositionNascent C3b reaching the surface0.3fractionestimatedMost nascent C3b hydrolyses in the fluid phase before it can bind; the surviving fraction is surface-dependent
decayIntrinsicSpontaneous convertase decay0.35/minmeasuredThe C3bBb convertase has an intrinsic half-life of roughly 90 seconds
decayCd55CD55 decay acceleration2.2/min at full densityestimatedErythrocytes lean on CD55, which is why PNH is a red-cell disease
decayFactorHFactor H decay acceleration2/min at full dockingestimatedGlomerular endothelium leans on factor H, which is why CFH variants present as renal thrombotic microangiopathy
inactivationFactor I cleavage of C3b0.62/min per unit cofactorestimatedCompetes directly against the loop for the same C3b; this rate decides the module's whole behaviour
c5DensityC3b density for half-maximal C5 convertase0.035fraction coverageestimatedLow, because the C5 convertase needs one more C3b placed locally rather than global coverage
c5CleavageC5 cleavage rate0.9/minestimatedEducational estimate
macAssemblyMAC assembly from C5b0.8/minestimatedEducational estimate
cd59PotencyCD59 pore inhibition0.94fraction blockedestimatedCD59 is a highly effective terminal inhibitor; near-complete at full density
lysisThresholdMAC density for half-maximal lysis0.045fraction coverageillustrativeMaps pore density onto a cell-population outcome; not a measured quantity
opsonicThresholdOpsonin density for half-maximal clearance0.12fraction coverageillustrativeMacrophage recognition is cooperative; the threshold itself is chosen, not measured
sheddingLoss of surface-bound fragments0.012/minillustrativeMembrane turnover and macrophage nibbling; without it a healthy cell slowly acquires a coat it would never carry

Cited variant evidence

app/disease-lab/disease-model.ts

The disease laboratory’s variants and the evidence behind each parameter shift. Directions come from these sources; the magnitudes are educational estimates chosen to be visible.

GeneVariantPhenotypeEvidenceSource
TNNT2p.Arg92Gln · R92QHypertrophic cardiomyopathydirect functional evidenceHuman R92Q cardiomyocyte study
MYBPC3loss-of-function classHypertrophic cardiomyopathyclass-level evidenceMYBPC3 functional review
TNNT2p.Lys210del · ΔK210Dilated cardiomyopathydirect functional evidenceΔK210 calcium-desensitization study
RYR2p.Arg4496Cys · R4496CCatecholaminergic polymorphic VTdirect functional evidenceR4496C calcium-wave study

Standing limitations