How to Use SimThyr for Thyroid Function Modeling

How to Use SimThyr for Thyroid Function Modeling

Overview

SimThyr is a software tool for simulating the hypothalamic–pituitary–thyroid (HPT) axis to study thyroid hormone dynamics and feedback. Use it to build, run, and analyze models of TSH, T4, T3 production and regulation under physiological or pathological conditions.

Installation

1. Download: Get the latest SimThyr release from the official project page (assume Windows/Linux/Mac builds available).
2. Install: Run the installer or extract the archive; place the executable in a folder with read/write permission for model files.
3. Dependencies: Ensure required runtime (e.g., Java or Python runtime) is installed if specified by the release notes.

Getting Started (basic workflow)

1. Open SimThyr and create a new model project.
2. Define compartments/variables: Add variables for TRH, TSH, free T4 (fT4), free T3 (fT3), and peripheral conversion if not preset.
3. Set equations: Use differential equations or built-in kinetic blocks to define production, clearance, and feedback loops (e.g., negative feedback of fT4/fT3 on TSH).
4. Parameterize: Enter parameter values — secretion rates, half-lives, feedback sensitivities, conversion rates, volumes of distribution. Use literature values for baseline models.
5. Initial conditions: Set initial hormone concentrations consistent with normal physiology.
6. Simulation settings: Choose time span, solver (stiff/nonstiff), tolerances, and output intervals.
7. Run simulation and observe time-series plots for all key variables.

Common Modeling Tasks

• Baseline physiology: Reproduce steady-state hormone levels by tuning secretion and clearance to match typical TSH, fT4, fT3.
• Thyroid hormone replacement: Simulate levothyroxine dosing schedules by adding dosing events to the thyroid hormone pool and observe TSH response over weeks.
• Central hypothyroidism: Reduce pituitary TSH responsiveness or TRH drive and simulate consequent hormone profiles.
• Thyrotoxicosis: Increase thyroidal hormone production rate; include changes in peripheral conversion if needed.
• Deiodinase alterations: Modify T4→T3 conversion rate to study effects on fT3 despite normal fT4.

Parameter Estimation & Fitting

• Use any built-in parameter estimation tools or export simulation outputs to CSV for external fitting (e.g., with R or Python).
• Fit model parameters to clinical time-series (TSH, fT4, fT3) using least-squares or Bayesian methods; constrain parameters within physiological ranges.

Sensitivity Analysis

• Perform local sensitivity (one-parameter-at-a-time) or global sensitivity (e.g., Sobol) to identify influential parameters (feedback strength, clearance rates).
• Use results to prioritize parameters for fitting or targeted experiments.

Validation & Good Practices

• Compare to clinical data: Validate steady-state and dynamic responses (e.g., after dosing or perturbation).
• Check mass balance and units: Ensure consistency across equations.
• Solver checks: Try multiple solvers/tolerances to confirm numerical stability.
• Document assumptions: Record simplifications (e.g., lumped compartments, ignored delays).

Troubleshooting

• Oscillatory or unstable solutions: Reduce time-step, switch solver, or check for missing negative feedback/delays.
• No change after perturbation: Ensure event/dose is correctly linked to the variable and units match.
• Poor fit to data: Re-examine parameter identifiability; perform sensitivity analysis; add measurement noise model if fitting.

Exporting & Reporting

• Export plots and time-series CSV for publications or further analysis.
• Include model equations, parameter sets, and simulation scripts with reports for reproducibility.

Example (conceptual)

• Model equations (conceptual):
  d[TSH]/dt = secretion_TSH(TRH, fT4) – clearance_TSH[TSH]
  d[fT4]/dt = production_T4(Thyroid, TSH) – conversion_T4_to_T3 – clearance_T4 * [fT4]
  d[fT3]/dt = conversion_T4_to_T3 – clearance_T3 * [fT3]

Further resources

• Check official SimThyr documentation and user guides for version-specific instructions and example models.