Mathematical Modeling of the Pituitary-Thyroid Feedback Loop: Role of a TSH-T3-Shunt and Sensitivity Analysis

Impaired autonomic function and somatosensory disturbance in patients with treated autoimmune thyroiditis

Despite treatment with levothyroxine, hypothyroidism and autoimmune thyroiditis (AIT) may be associated with reduced quality of life (QoL), an enigmatic condition referred to as "syndrome T". Peripheral neuropathy, described in untreated thyroid disease, could be a contributing mechanism. We analysed autonomic and somatosensory function in 29 patients with AIT and treated hypothyroidism and 27 healthy volunteers. They underwent heart rate variability (HRV) analysis and quantitative sensory testing (n = 28), comprising 13 parameters of small and large nerve fibre function and pain thresholds. Autonomic cardiovascular function was assessed in rest, deep respiration and orthostasis. Additionally, biomarkers for autoimmunity and thyroid function were measured. Anxiety, depression and QoL were assessed using validated questionnaires. 36% of the patients showed at least one sign of somatosensory small or large fibre dysfunction. 57% presented with mild hyperalgesia to at least one stimulus. Several markers of autonomic function and some detection thresholds were related to the antibody titres. Anxiety, depression scores and QoL correlated to antibody titres and HRV measures. Autonomic and somatosensory dysfunction indicate that in treated hypothyroidism and AIT a subgroup of patients suffers from neuropathic symptoms leading to impaired QoL. Additionally, mild hyperalgesia as a possible sensitisation phenomenon should be considered a target for symptomatic treatment.

Spatially Dependent Tissue Distribution of Thyroid Hormones by Plasma Thyroid Hormone Binding Proteins

Plasma thyroid hormone (TH) binding proteins (THBPs), including thyroxine-binding globulin (TBG), transthyretin (TTR), and albumin (ALB), carry THs to extrathyroidal sites, where THs are unloaded locally and then taken up via membrane transporters into the tissue proper. The respective roles of THBPs in supplying THs for tissue uptake are not completely understood. To investigate this, we developed a spatial human physiologically based kinetic (PBK) model of THs, which produces several novel findings. (1) Contrary to postulations that TTR and/or ALB are the major local T4 contributors, the three THBPs may unload comparable amounts of T4 in Liver, a rapidly perfused organ; however, their contributions in slowly perfused tissues follow the order of abundances of T4TBG, T4TTR, and T4ALB. The T3 amounts unloaded from or loaded onto THBPs in a tissue acting as a T3 sink or source respectively follow the order of abundance of T3TBG, T3ALB, and T3TTR regardless of perfusion rate. (2) Any THBP alone is sufficient to maintain spatially uniform TH tissue distributions. (3) The TH amounts unloaded by each THBP species are spatially dependent and nonlinear in a tissue, with ALB being the dominant contributor near the arterial end but conceding to TBG near the venous end. (4) Spatial gradients of TH transporters and metabolic enzymes may modulate these contributions, producing spatially invariant or heterogeneous TH tissue concentrations depending on whether the blood-tissue TH exchange operates in near-equilibrium mode. In summary, our modeling provides novel insights into the differential roles of THBPs in local TH tissue distribution.

Automatic Levothyroxine Dosing Algorithm for Patients Suffering from Hashimoto's Thyroiditis

Hypothyroidism is a condition where the patient's thyroid gland cannot produce sufficient thyroid hormones (mainly triiodothyronine and thyroxine). The primary cause of hypothyroidism is autoimmune-mediated destruction of the thyroid gland, referred to as Hashimoto's thyroiditis. A patient's desired thyroid hormone concentration is achieved by oral administration of thyroid hormone, usually levothyroxine. Establishing individual levothyroxine doses to achieve desired thyroid hormone concentrations requires several patient visits. Additionally, clear guidance for the dosing regimen is lacking, and significant inter-individual differences exist. This study aims to design a digital automatic dosing algorithm for patients suffering from Hashimoto's thyroiditis. The dynamic behaviour of the relevant thyroid function is mathematically modelled. Methods of automatic control are exploited for the design of the proposed robust model-based levothyroxine dosing algorithm. Numerical simulations are performed to evaluate the mathematical model and the dosing algorithm. With the help of the developed controller thyroid hormone concentrations of patients, emulated using Thyrosim, have been regulated under the euthyroid state. The proposed concept demonstrates reliable responses amidst varying patient parameters. Our developed model provides a useful basis for the design of automatic levothyroxine dosing algorithms. The proposed robust feedback loop contributes to the first results for computer-assisted thyroid dosing algorithms.

A minimal human physiologically based kinetic model of thyroid hormones and chemical disruption of plasma thyroid hormone binding proteins

The thyroid hormones (THs), thyroxine (T4) and triiodothyronine (T3), are under homeostatic control by the hypothalamic-pituitary-thyroid axis and plasma TH binding proteins (THBPs), including thyroxine-binding globulin (TBG), transthyretin (TTR), and albumin (ALB). THBPs buffer free THs against transient perturbations and distribute THs to tissues. TH binding to THBPs can be perturbed by structurally similar endocrine-disrupting chemicals (EDCs), yet their impact on circulating THs and health risks are unclear. In the present study, we constructed a human physiologically based kinetic (PBK) model of THs and explored the potential effects of THBP-binding EDCs. The model describes the production, distribution, and metabolism of T4 and T3 in the Body Blood, Thyroid, Liver, and Rest-of-Body (RB) compartments, with explicit consideration of the reversible binding between plasma THs and THBPs. Rigorously parameterized based on literature data, the model recapitulates key quantitative TH kinetic characteristics, including free, THBP-bound, and total T4 and T3 concentrations, TH productions, distributions, metabolisms, clearance, and half-lives. Moreover, the model produces several novel findings. (1) The blood-tissue TH exchanges are fast and nearly at equilibrium especially for T4, providing intrinsic robustness against local metabolic perturbations. (2) Tissue influx is limiting for transient tissue uptake of THs when THBPs are present. (3) Continuous exposure to THBP-binding EDCs does not alter the steady-state levels of THs, while intermittent daily exposure to rapidly metabolized TBG-binding EDCs can cause much greater disruptions to plasma and tissue THs. In summary, the PBK model provides novel insights into TH kinetics and the homeostatic roles of THBPs against thyroid disrupting chemicals.

Modelling of the thyroid hormone synthesis as a part of nonlinear reaction mechanism with feedback

The synthesis of thyroid hormones in the hypothalamic-pituitary-thyroid (HPT) axis was studied. For this purpose, a reaction model for HPT axis with stoichiometric relations between the main reaction species was postulated. Using the law of mass action, this model has been transformed into a set of nonlinear ordinary differential equations. This new model has been examined by stoichiometric network analysis (SNA) with the aim to see if it possesses the ability to reproduce oscillatory ultradian dynamics founded on the internal feedback mechanism. In particular, a feedback regulation of TSH production based on the interplay between TRH, TSH, somatostatin and thyroid hormones was proposed. Besides, the ten times larger amount of produced T₄ with respect to T₃ in the thyroid gland was successfully simulated. The properties of SNA in combination with experimental results, were used to determine the unknown parameters (19 rate constants of particular reaction steps) necessary for numerical investigations. The steady-state concentrations of 15 reactive species were tuned to be consistent with the experimental data. The predictive potential of the proposed model was illustrated on numerical simulations of somatostatin influence on TSH dynamics investigated experimentally by Weeke et al. in 1975. In addition, all programs for SNA analysis were adapted for this kind of a large model. The procedure of calculating rate constants from steady-state reaction rates and very limited available experimental data was developed. For this purpose, a unique numerical method was developed to fine-tune model parameters while preserving the fixed rate ratios and using the magnitude of the experimentally known oscillation period as the only target value. The postulated model was numerically validated by perturbation simulations with somatostatin infusion and the results were compared with experiments available in literature. Finally, as far as we know, this reaction model with 15 variables is the most dimensional one that have been analysed mathematically to obtain instability region and oscillatory dynamic states. Among the existing models of thyroid homeostasis this theory represents a new class that may improve our understanding of basic physiological processes and helps to develop new therapeutic approaches. Additionally, it may pave the way to improved diagnostic methods for pituitary and thyroid disorders.

Autoimmune thyroid diseases as a cost of physiological autoimmune surveillance

Graves' disease (GD) and Hashimoto's thyroiditis (HT) are common autoimmune diseases of the thyroid gland, causing hyperthyroidism and hypothyroidism, respectively. Despite their opposing clinical manifestation, they have several enigmatic links. Here, we propose that GD and HT have the same fundamental origin: both diseases are the cost of a beneficial physiological process called autoimmune surveillance of hypersecreting mutants. Autoreactive T cells selectively eliminate mutant cells that hypersecrete the hormones and threaten to become toxic nodules. These T cells can trigger a humoral response in susceptible individuals, leading to the production of antibodies against thyroid antigens. This shared origin can explain similarities in incidence and risk factors between HT and GD, despite their opposite clinical phenotypes.

Clinically practical pharmacometrics computer model to evaluate and personalize pharmacotherapy in pediatric rare diseases: application to Graves' disease

Objectives

Graves' disease (GD) with onset in childhood or adolescence is a rare disease (ORPHA:525731). Current pharmacotherapeutic approaches use antithyroid drugs, such as carbimazole, as monotherapy or in combination with thyroxine hormone substitutes, such as levothyroxine, as block-and-replace therapy to normalize thyroid function and improve patients' quality of life. However, in the context of fluctuating disease activity, especially during puberty, a considerable proportion of pediatric patients with GD is suffering from thyroid hormone concentrations outside the therapeutic reference ranges. Our main goal was to develop a clinically practical pharmacometrics computer model that characterizes and predicts individual disease activity in children with various severity of GD under pharmacotherapy.

Methods

Retrospectively collected clinical data from children and adolescents with GD under up to two years of treatment at four different pediatric hospitals in Switzerland were analyzed. Development of the pharmacometrics computer model is based on the non-linear mixed effects approach accounting for inter-individual variability and incorporating individual patient characteristics. Disease severity groups were defined based on free thyroxine (FT4) measurements at diagnosis.

Results

Data from 44 children with GD (75% female, median age 11 years, 62% receiving monotherapy) were analyzed. FT4 measurements were collected in 13, 15, and 16 pediatric patients with mild, moderate, or severe GD, with a median FT4 at diagnosis of 59.9 pmol/l (IQR 48.4, 76.8), and a total of 494 FT4 measurements during a median follow-up of 1.89 years (IQR 1.69, 1.97). We observed no notable difference between severity groups in terms of patient characteristics, daily carbimazole starting doses, and patient years. The final pharmacometrics computer model was developed based on FT4 measurements and on carbimazole or on carbimazole and levothyroxine doses involving two clinically relevant covariate effects: age at diagnosis and disease severity.

Discussion

We present a tailored pharmacometrics computer model that is able to describe individual FT4 dynamics under both, carbimazole monotherapy and carbimazole/levothyroxine block-and-replace therapy accounting for inter-individual disease progression and treatment response in children and adolescents with GD. Such clinically practical and predictive computer model has the potential to facilitate and enhance personalized pharmacotherapy in pediatric GD, reducing over- and underdosing and avoiding negative short- and long-term consequences. Prospective randomized validation trials are warranted to further validate and fine-tune computer-supported personalized dosing in pediatric GD and other rare pediatric diseases.

The role of supporting and disruptive mechanisms of FT3 homeostasis in regulating the hypothalamic-pituitary-thyroid axis

BACKGROUND

Thyroid hormones are controlled by the hypothalamic-pituitary-thyroid (HPT) axis through a complex network of regulatory loops, involving the hormones TRH, TSH, FT4, and FT3. The relationship between TSH and FT4 is widely used for diagnosing thyroid diseases. However, mechanisms of FT3 homeostasis are not well understood.

OBJECTIVE

We used mathematical modelling to further examine mechanisms that exist in the HPT axis regulation for protecting circulating FT3 levels.

METHODS

A mathematical model consisting of a system of four coupled first-order parameterized non-linear ordinary differential equations (ODEs) was developed, accounting for the interdependencies between the hormones in the HPT axis regulation. While TRH and TSH feed forward to the pituitary and thyroid, respectively, FT4 and FT3 feed backward to both the pituitary and hypothalamus. Stable equilibrium solutions of the ODE system express homeostasis for a particular variable, such as FT3, if this variable stays in a narrow range while certain other parameter(s) and system variable(s) may vary substantially.

RESULTS

The model predicts that (1) TSH-feedforward protects FT3 levels if the FT4 production rate declines and (2) combined negative feedback by FT4 and FT3 on both TSH and TRH production rates keeps FT3 levels insensitive to moderate changes in FT4 production rates and FT4 levels. The optimum FT4 and FT3 feedback and TRH and TSH-feedforward ranges that preserve FT3 homeostasis were found by numerical continuation analysis. Model predictions were in close agreement with clinical studies and individual patient examples of hypothyroidism and hyperthyroidism.

CONCLUSIONS

These findings further extend the concept of HPT axis regulation beyond TSH and FT4 to integrate the more active sister hormone FT3 and mechanisms of FT3 homeostasis. Disruption of homeostatic mechanisms leads to disease. This provides a perspective for novel testable concepts in clinical studies to therapeutically target the disruptive mechanisms.

The progression of secondary diabetes: A review of modeling studies

Mathematical modeling has provided quantitative information consistent with experimental data, greatly improving our understanding of the progression of type 1 and type 2 diabetes. However, diabetes is a complex metabolic disease and has been found to be involved in crosstalk interactions with diverse endocrine diseases. Mathematical models have also been developed to investigate the quantitative impact of various hormonal disorders on glucose imbalance, advancing the precision treatment for secondary diabetes. Here we review the models established for the study of dysglycemia induced by hormonal disorders, such as excessive glucocorticoids, epinephrine, and growth hormone. To investigate the influence of hyperthyroidism on the glucose regulatory system, we also propose a hyperthyroid-diabetes progression model. Model simulations indicate that timely thyroid treatment can halt the progression of hyperglycemia and prevent beta-cell failure. This highlights the diagnosis of hormonal disorders, together withblood sugar tests, as significant measures for the early diagnosis and treatment of diabetes. The work recapitulates updated biological research on the interactions between the glucose regulatory system and other endocrine axes. Further mathematical modeling of secondary diabetes is desired to promote the quantitative study of the disease and the development of individualized diabetic therapies.

Intrathyroidal feedforward and feedback network regulating thyroid hormone synthesis and secretion

Thyroid hormones (THs), including T4 and T3, are produced and released by the thyroid gland under the stimulation of thyroid-stimulating hormone (TSH). The homeostasis of THs is regulated via the coordination of the hypothalamic-pituitary-thyroid axis, plasma binding proteins, and local metabolism in tissues. TH synthesis and secretion in the thyrocytes-containing thyroid follicles are exquisitely regulated by an elaborate molecular network comprising enzymes, transporters, signal transduction machineries, and transcription factors. In this article, we synthesized the relevant literature, organized and dissected the complex intrathyroidal regulatory network into structures amenable to functional interpretation and systems-level modeling. Multiple intertwined feedforward and feedback motifs were identified and described, centering around the transcriptional and posttranslational regulations involved in TH synthesis and secretion, including those underpinning the Wolff-Chaikoff and Plummer effects and thyroglobulin-mediated feedback regulation. A more thorough characterization of the intrathyroidal network from a systems biology perspective, including its topology, constituent network motifs, and nonlinear quantitative properties, can help us to better understand and predict the thyroidal dynamics in response to physiological signals, therapeutic interventions, and environmental disruptions.

Dynamics of thyroid diseases and thyroid-axis gland masses

Thyroid disorders are common and often require lifelong hormone replacement. Treating thyroid disorders involves a fascinating and troublesome delay, in which it takes many weeks for serum thyroid-stimulating hormone (TSH) concentration to normalize after thyroid hormones return to normal. This delay challenges attempts to stabilize thyroid hormones in millions of patients. Despite its importance, the physiological mechanism for the delay is unclear. Here, we present data on hormone delays from Israeli medical records spanning 46 million life-years and develop a mathematical model for dynamic compensation in the thyroid axis, which explains the delays. The delays are due to a feedback mechanism in which peripheral thyroid hormones and TSH control the growth of the thyroid and pituitary glands; enlarged or atrophied glands take many weeks to recover upon treatment due to the slow turnover of the tissues. The model explains why thyroid disorders such as Hashimoto's thyroiditis and Graves' disease have both subclinical and clinical states and explains the complex inverse relation between TSH and thyroid hormones. The present model may guide approaches to dynamically adjust the treatment of thyroid disorders.

Optimal Hormone Replacement Therapy in Hypothyroidism - A Model Predictive Control Approach

In this paper, we address the problem of optimal thyroid hormone replacement strategy development for hypothyroid patients. This is challenging for the following reasons. First, it is difficult to determine the correct dosage leading to normalized serum thyroid hormone concentrations of a patient. Second, it remains unclear whether a levothyroxine L-T4) monotherapy or a liothyronine/levothyroxine (L-T3/L-T4) combined therapy is more suitable to treat hypothyroidism. Third, the optimal intake frequency of L-T3/L-T4 is unclear. We address these issues by extending a mathematical model of the pituitary-thyroid feedback loop to be able to consider an oral intake of L-T3/L-T4. A model predictive controller (MPC) is employed to determine optimal dosages with respect to the thyroid hormone concentrations for each type of therapy. The results indicate that the L-T3/L-T4 combined therapy is slightly better (in terms of the achieved hormone concentrations) to treat hypothyroidism than the L-T4 monotherapy. In case of a specific genetic variant, namely genotype CC in polymorphism rs2235544 of gene DIO1, the simulation results suggest that the L-T4 monotherapy is better to treat hypothyroidism. In turn, when genotype AA is considered, the L-T3/L-T4 combined therapy is better to treat hypothyroidism. Furthermore, when genotype CC of polymorphism rs225014 (also referred to as c.274A>G or p.Thr92Ala) in the DIO2 gene is considered, the outcome of the L-T3/L-T4 combined therapy is better in terms of the steady-state hormone concentrations (for a triiodothyronine setpoint at the upper limit of the reference range of healthy individuals). Finally, the results suggest that two daily intakes of L-T3 could be the best trade-off between stable hormone concentrations and inconveniences for the patient.

Principles of Endocrine Regulation: Reconciling Tensions Between Robustness in Performance and Adaptation to Change

Endocrine regulation in the hypothalamic-pituitary-thyroid (HPT) axis is orchestrated by physiological circuits which integrate multiple internal and external influences. Essentially, it provides either of the two responses to overt biological challenges: to defend the homeostatic range of a target hormone or adapt it to changing environmental conditions. Under certain conditions, such flexibility may exceed the capability of a simple feedback control loop, rather requiring more intricate networks of communication between the system's components. A new minimal mathematical model, in the form of a parametrized nonlinear dynamical system, is here formulated as a proof-of-concept to elucidate the principles of the HPT axis regulation. In particular, it allows uncovering mechanisms for the homeostasis of the key biologically active hormone free triiodothyronine (FT3). One mechanism supports the preservation of FT3 homeostasis, whilst the other is responsible for the adaptation of the homeostatic state to a new level. Together these allow optimum resilience in stressful situations. Preservation of FT3 homeostasis, despite changes in FT4 and TSH levels, is found to be an achievable system goal by joining elements of top-down and bottom-up regulation in a cascade of targeted feedforward and feedback loops. Simultaneously, the model accounts for the combination of properties regarded as essential to endocrine regulation, namely sensitivity, the anticipation of an adverse event, robustness, and adaptation. The model therefore offers fundamental theoretical insights into the effective system control of the HPT axis.

Clustering Based Parameter Estimation of Thyroid Hormone Pathway

An ordinary differential equation (ODE) model of the working of the thyroid system for euthyroidism has been presented. As clinical data for thyroid hormones is relatively scarce, such modelling offers potential benefits over wet lab procedures. Genetic algorithms developed for determining of parameters of the ODE system using the available data have been presented and evaluated. This approach enables subject specific parameter estimation towards characterisation of individual thyroid operation. Initially, a simple steady state model was used. Later a cosinor model for the circadian variation of thyroid hormones was used to obtain more reliable results, as indicated through sensitivity analysis in conjunction with other statistical methods. Our parameter determination method has been tested on groups of patients with similar observed values of thyroid stimulating hormone (TSH), free T ₃ and free T ₄ (identified through clustering) to determine their parameter values jointly. This approach appears to produce parameter sets with lower variation than parameters determined independently, thus leading to better parameter determination.

Mathematical modeling and simulation of thyroid homeostasis: Implications for the Allan-Herndon-Dudley syndrome

INTRODUCTION

A mathematical model of the pituitary-thyroid feedback loop is extended to deepen the understanding of the Allan-Herndon-Dudley syndrome (AHDS). The AHDS is characterized by unusual thyroid hormone concentrations and a mutation in the SLC16A2 gene encoding for the monocarboxylate transporter 8 (MCT8). This mutation leads to a loss of thyroid hormone transport activity. One hypothesis to explain the unusual hormone concentrations of AHDS patients is that due to the loss of thyroid hormone transport activity, thyroxine (T ₄) is partially retained in thyroid cells.

METHODS

This hypothesis is investigated by extending a mathematical model of the pituitary-thyroid feedback loop to include a model of the net effects of membrane transporters such that the thyroid hormone transport activity can be considered. A nonlinear modeling approach based on the Michaelis-Menten kinetics and its linear approximation are employed to consider the membrane transporters. The unknown parameters are estimated through a constrained parameter optimization.

RESULTS

In dynamic simulations, damaged membrane transporters result in a retention of T ₄ in thyroid cells and ultimately in the unusual hormone concentrations of AHDS patients. The Michaelis-Menten modeling approach and its linear approximation lead to similar results.

DISCUSSION

The results support the hypothesis that a partial retention of T ₄ in thyroid cells represents one mechanism responsible for the unusual hormone concentrations of AHDS patients. Moreover, our results suggest that the retention of T ₄ in thyroid cells could be the main reason for the unusual hormone concentrations of AHDS patients.

Enhancing the Natural Biological Control in the Thyroid Hormone Homeostasis As a First-Order Control System

This study explores the natural control system that exists within the pituitary gland. More specifically, this study investigates the regulation of the thyroid stimulating hormone (TSH), released by the anterior pituitary, with regards to the thyroid releasing hormone (TRH), which is released by the hypothalamus. Using appropriate assumptions on the behavior of the hormones, along with relevant boundary conditions, we modeled an output of TSH using constant TRH input over the course of a six-hour period. Other relevant hormones such as thyroxine (T4), triiodothyronine (T3), and their relevant intermediaries were also modeled as a means to complete the natural feedback found physiologically. Due to our boundary conditions, we do not consider the consumption or final function of these hormones since they leave the pituitary gland, our control system; instead, we consider a constant TRH since it is produced by the hypothalamus. Finally, we explore the results of reducing the TRH input while observing the TSH response. We append a short loop controller feedback that uses the TSH output to regulate a TRH input to remedy the reduction of TRH. The open-loop transfer function derived presented three poles at the clearance exponents for T4, TSH, and central T3, with a phase margin of 74.1°, characterizing a stable but slow system that can be improved with a simple proportional control.

Modeling of levothyroxine in newborns and infants with congenital hypothyroidism: challenges and opportunities of a rare disease multi-center study

Modeling of retrospectively collected multi-center data of a rare disease in pediatrics is challenging because laboratory data can stem from several decades measured with different assays. Here we present a retrospective pharmacometrics (PMX) based data analysis of the rare disease congenital hypothyroidism (CH) in newborns and infants. Our overall aim is to develop a model that can be applied to optimize dosing in this pediatric patient population since suboptimal treatment of CH during the first 2 years of life is associated with a reduced intelligence quotient between 10 and 14 years. The first goal is to describe a retrospectively collected dataset consisting of 61 newborns and infants with CH up to 2 years of age. Overall, 505 measurements of free thyroxine (FT4) and 510 measurements of thyrotropin or thyroid-stimulating hormone were available from patients receiving substitution treatment with levothyroxine (LT4). The second goal is to introduce a scale/location-scale normalization method to merge available FT4 measurements since 34 different postnatal age- and assay-specific laboratory reference ranges were applied. This method takes into account the change of the distribution of FT4 values over time, i.e. a transformation from right-skewed towards normality during LT4 treatment. The third goal is to develop a practical and useful PMX model for LT4 treatment to characterize FT4 measurements, which is applicable within a clinical setting. In summary, a time-dependent normalization method and a practical PMX model are presented. Since there is no on-going or planned development of new pharmacological approaches for CH, PMX based modeling and simulation can be leveraged to personalize dosing with the goal to enhance longer-term neurological outcome in children with the rare disease CH.

A unified mathematical model of thyroid hormone regulation and implication for personalized treatment of thyroid disorders

Current clinician practice for thyroid hormone regulation of patients is based upon guesswork and experience rather than quantified analysis, which exposes patients under longer risk and discomfort. To quantitatively analyze the thyroid regulation for patients of different thyroid states, we develop a two-dimensional mathematical model that can be applied to analyze the dynamic behaviors of thyroid hormones with or without drug intervention. The unified model can be employed to study the regulation of TSH (thyroid-stimulating hormone) and FT4 (free thyroxine) for euthyroid (normal thyroid) subjects, Hashimoto's thyroiditis, and Graves' disease patients, respectively. The results suggest that the level of TPOAb (thyroid peroxidase antibody) may be a factor determining whether the patient would progress from euthyroid state to subclinical or clinical hypothyroidism, and that increased TRAb (TSH receptor antibody) may lead Graves' disease to deteriorate from the early stage to overt hyperthyroidism. Given the early blood-test data, we demonstrate the feasibility for healthcare professionals to apply our model in choosing an appropriate dosage regimen for patients to achieve the desired TSH and FT4 levels within a specified time frame. This proposed model has the potential to optimize personalized treatment and shorten the therapeutic time for patients suffering from Hashimoto's thyroiditis and Graves' disease.

Mechanistic Computational Model for Extrapolating In vitro Thyroid Peroxidase (TPO) Inhibition Data to Predict Serum Thyroid Hormone Levels in Rats

High throughput (HTP) in vitro assays are developed to screen chemicals for their potential to inhibit thyroid hormones (THs) synthesis. Some of these experiments, such as the thyroid peroxidase (TPO) inhibition assay, are based on thyroid microsomal extracts. However, the regulation of thyroid disruption chemicals (TDCs) is based on THs in vivo serum levels. This necessitates the estimation of TDCs in vivo tissue levels in the thyroid where THs synthesis inhibition by TPO takes place. The in vivo tissue levels of chemicals are controlled by pharmacokinetic determinants such as absorption, distribution, metabolism and excretion (ADME), and can be described quantitatively in physiologically based pharmacokinetic (PBPK) models. An integrative computational model including chemical specific PBPK and TH kinetics models provides a mechanistic quantitative approach to translate thyroidal HTP in vitro assays to in vivo measures of circulating THs serum levels. This computational framework is developed to quantitatively establish the linkage between applied dose, chemical thyroid tissue levels, thyroid TPO inhibition potential, and in vivo TH serum levels. Once this link is established quantitively, the overall model is used to calibrate the TH kinetics parameters using experimental data for THs levels in thyroid tissue and serum for the two drugs Propylthiouracil (PTU) and Methimazole (MMI). The calibrated quantitative framework is then evaluated against literature data for the environmental chemical ethylenethiourea (ETU). The linkage of PBPK and TH kinetics models illustrates a computational framework that can be extrapolated to humans to screen chemicals based on their exposure levels and potential to disrupt serum THs levels in vivo.

The Quest for System-Theoretical Medicine in the COVID-19 Era

Precision medicine and molecular systems medicine (MSM) are highly utilized and successful approaches to improve understanding, diagnosis, and treatment of many diseases from bench-to-bedside. Especially in the COVID-19 pandemic, molecular techniques and biotechnological innovation have proven to be of utmost importance for rapid developments in disease diagnostics and treatment, including DNA and RNA sequencing technology, treatment with drugs and natural products and vaccine development. The COVID-19 crisis, however, has also demonstrated the need for systemic thinking and transdisciplinarity and the limits of MSM: the neglect of the bio-psycho-social systemic nature of humans and their context as the object of individual therapeutic and population-oriented interventions. COVID-19 illustrates how a medical problem requires a transdisciplinary approach in epidemiology, pathology, internal medicine, public health, environmental medicine, and socio-economic modeling. Regarding the need for conceptual integration of these different kinds of knowledge we suggest the application of general system theory (GST). This approach endorses an organism-centered view on health and disease, which according to Ludwig von Bertalanffy who was the founder of GST, we call Organismal Systems Medicine (OSM). We argue that systems science offers wider applications in the field of pathology and can contribute to an integrative systems medicine by (i) integration of evidence across functional and structural differentially scaled subsystems, (ii) conceptualization of complex multilevel systems, and (iii) suggesting mechanisms and non-linear relationships underlying the observed phenomena. We underline these points with a proposal on multi-level systems pathology including neurophysiology, endocrinology, immune system, genetics, and general metabolism. An integration of these areas is necessary to understand excess mortality rates and polypharmacological treatments. In the pandemic era this multi-level systems pathology is most important to assess potential vaccines, their effectiveness, short-, and long-time adverse effects. We further argue that these conceptual frameworks are not only valid in the COVID-19 era but also important to be integrated in a medicinal curriculum.

Could Glyphosate and Glyphosate-Based Herbicides Be Associated With Increased Thyroid Diseases Worldwide?

The increased incidence of thyroid diseases raises a series of questions about what the main predisposing factors are nowadays. If dietary restriction of iodine was once a major global health concern, today, the processes of industrialization of food and high exposure to a wide variety of environmental chemicals may be affecting, directly or indirectly, thyroid function. The homeostasis of hypothalamus-pituitary-thyroid (HPT) axis is finely regulated through the negative feedback mechanism exerted by thyroid hormones. Allostatic mechanisms are triggered to adjust the physiology of HPT axis in chronic conditions. Glyphosate and glyphosate-based herbicides are pesticides with controversial endocrine disrupting activities and only few studies have approached their effects on HPT axis and thyroid function. However, glyphosate has an electrophilic and nucleophilic zwitterion chemical structure that may affect the mechanisms involved in iodide oxidation and organification, as well as the oxidative phosphorylation in the ATP synthesis. Thus, in this review, we aimed to: (1) discuss the critical points in the regulation of HPT axis and thyroid hormones levels balance, which may be susceptible to the toxic action of glyphosate and glyphosate-based herbicides, correlating the molecular mechanisms involved in glyphosate toxicity described in the literature that may, directly or indirectly, be associated to the higher incidence of thyroid diseases; and (2) present the literature regarding glyphosate toxicity in HPT axis.

Lipidomic Perturbations in Cynomolgus Monkeys are Regulated by Thyroid Stimulating Hormone

Thyroid disease affects an estimated 200 million people worldwide, and is commonly associated with increased blood lipid levels. However, the mechanism by which thyroid-stimulating hormone (TSH) affects lipid profiles is not clear. Twenty-four cynomolgus monkeys were treated with a novel exogenous recombinant human TSH (rhTSH) (SNA001) at 9 μg kg⁻¹, 22 μg kg⁻¹, or 54 μg kg⁻¹, and reference rhTSH (Thyrogen^®) at 22 μg kg⁻¹. The primary TSH (SNA001) pharmacokinetic (PK) parameters increased in a dose-dependent manner across the dose range of 9 μg kg⁻¹, 22 μg kg⁻¹, or 54 μg kg⁻¹. Peak triiodothyronine (T3) and thyroxine (T4) levels were reached within 24 h after rhTSH administration, which was delayed by approximately 20 h. In total, 420 lipid species were detected and quantified by ultra-performance liquid chromatography high resolution spectrometry (UPLC-HR-MS)-based lipidomics. Notably, peak levels of lipid accumulation, particularly sphingomyelin (SM) and triglycerides (TG), appeared at 4 and 24 h, which was consistent with the pattern of TSH and T3/T4 levels, respectively. According to weighted correlation network analysis (WGCNA), perturbations of many lipid species were strongly correlated with TSH and T3/T4 levels. TSH and the stimulated T3/T4 levels and lipid profiles following SNA001 administration were comparable to those after administration of the reference rhTSH (Thyrogen^®). The plasma lipidome and changes in lipid levels after rhTSH stimulation were associated with TSH and T3/T4 concentrations. T3/T4 and lipid profiles were delayed after TSH stimulation. Such phenomena require further exploration.

Endocrine Autoimmune Disease as a Fragility of Immune Surveillance against Hypersecreting Mutants

Some endocrine organs are frequent targets of autoimmune attack. Here, we addressed the origin of autoimmune disease from the viewpoint of feedback control. Endocrine tissues maintain mass through feedback loops that balance cell proliferation and removal according to hormone-driven regulatory signals. We hypothesized the existence of a dedicated mechanism that detects and removes mutant cells that missense the signal and therefore hyperproliferate and hypersecrete with potential to disrupt organismal homeostasis. In this mechanism, hypersecreting cells are preferentially eliminated by autoreactive T cells at the cost of a fragility to autoimmune disease. The "autoimmune surveillance of hypersecreting mutants" (ASHM) hypothesis predicts the presence of autoreactive T cells in healthy individuals and the nature of self-antigens as peptides from hormone secretion pathway. It explains why some tissues get prevalent autoimmune disease, whereas others do not and instead show prevalent mutant-expansion disease (e.g., hyperparathyroidism). The ASHM hypothesis is testable, and we discuss experimental follow-up.

Heterogenous biochemical expression of hormone activity in subclinical/overt hyperthyroidism and exogenous thyrotoxicosis

Background

Subclinical hyperthyroidism/thyrotoxicosis originates from different causes and clinical conditions, sharing the laboratory constellation of a suppressed TSH in the presence of thyroid hormone concentrations within the reference range.

Aim

Presentation of hyperthyroidism can manifest itself in several ways. We questioned whether there is either a consistent biochemical equivalence of thyroid hormone response to these diagnostic categories, or a high degree of heterogeneity may exist both within and between the different clinical manifestations.

Methods

This secondary analysis of a former prospective cross-sectional trial involved 461 patients with untreated thyroid autonomy, Graves’ disease or on levothyroxine (LT4) after thyroidectomy for thyroid carcinoma. TSH response and biochemical equilibria between TSH and thyroid hormones were contrasted between endogenous hyperthyroidism and thyrotoxicosis (LT4 overdose).

Results

Concentrations of FT4, FT3, TSH, deiodinase activity and BMI differed by diagnostic category. Over various TSH strata, FT4 concentrations were significantly higher in LT4-treated thyroid carcinoma patients, compared to the untreated diseases, though FT3 levels remained comparable. They were concentrated in the upper FT4- but low deiodinase range, distinguishing them from patients with thyroid autonomy and Graves’ disease. In exogenous thyrotoxicosis, TSH and FT3 were less responsive to FT4 concentrations approaching its upper normal/hyperthyroid range.

Conclusions

The presence or lack of TSH feedforward activity determines the system response in the thyroid-active (hyperthyroidism) and no-thyroid response to treatment (thyrotoxicosis). This rules out a consistent thread of thyroid hormone response running through the different diagnostic categories. TSH measurements should therefore be interpreted conditionally and differently in subclinical hyperthyroidism and thyrotoxicosis.

Triiodothyronine secretion in early thyroid failure: The adaptive response of central feedforward control

BACKGROUND

Defined by thyroid-pituitary feedback control, clinical diagnosis of hypothyroidism and hyperthyroidism has become synonymous with TSH measurement. We combined in silico analysis and in vivo data to explore the central influences on thyroidal T3 production.

MATERIALS & METHODS

A system of five coupled first-order nonlinear parameterised ordinary differential equations (ODEs) is used to model the feedback control of TSH and TRH by thyroid hormones together with the feedforward control of thyroidal T3 secretion and enzymatic T4-T3 conversion. Dependencies of the stable equilibrium solutions of this ODE system, that is the homeostasis of the underlying physiological process, on the system parameters were investigated whether they accounted for clinical observations.

RESULTS

During the modelled transition to hypothyroidism, central control imposed an increasing influence in maintaining serum FT3 levels, compared to peripheral conversion efficiency. Numerical continuation analysis revealed dependencies of T3 production on different elements of TSH feedforward control. While T4-T3 conversion provided the main T3 source in euthyroidism, this was overtaken by increasing glandular T3 secretion when thyroid reserve declined. The computational results were in good agreement with data from untreated patients with autoimmune thyroiditis.

CONCLUSIONS

Dependencies revealed in the expression of control differ in thyroid health and disease, using a physiologically based mathematical model of combined feedback-feedforward control of the hypothalamic-pituitary-thyroid regulation. Strong T3-protective mechanisms of the control system emerge with declining thyroid function, when glandular T3 secretion becomes increasingly influential over conversion efficiency. This has wide-ranging implications for the utility of TSH in clinical decision-making.

Individualised requirements for optimum treatment of hypothyroidism: complex needs, limited options

Levothyroxine (LT4) therapy has a long history, a well-defined pharmacological profile and a favourable safety record in the alleviation of hypothyroidism. However, questions remain in defining the threshold for the requirement of treatment in patients with subclinical hypothyroidism, assessing the dose adequacy of the drug, and selecting the best treatment mode (LT4 monotherapy versus liothyronine [LT3]/LT4 combinations) for subpopulations with persisting complaints. Supplied as a prodrug, LT4 is enzymatically converted into the biologically more active thyroid hormone, triiodothyronine (T3). Importantly, tetraiodothyronine (T4) to T3 conversion efficiency may be impaired in patients receiving LT4, resulting in a loss of thyroid-stimulating hormone (TSH)-mediated feed-forward control of T3, alteration of the interlocking equilibria between serum concentrations of TSH, free thyroxine (FT4), and free triiodothyonine (FT3), and a decrease in FT3 to FT4 ratios. This downgrades the value of the TSH reference system derived in thyroid health for guiding the replacement dose in the treatment situation. Individualised conditionally defined setpoints may therefore provide appropriate biochemical targets to be clinically tested, together with a stronger focus on clinical presentation and future endpoint markers of tissue thyroid state. This cautionary note encompasses the use of aggregated statistical data from clinical trials which are not safely applicable to the individual level of patient care under these circumstances.

Functional and Symptomatic Individuality in the Response to Levothyroxine Treatment

Background: For significant numbers of patients dissatisfied on standard levothyroxine (LT4) treatment for hypothyroidism, patient-specific responses to T4 could play a significant role. Aim: To assess response heterogeneity to LT4 treatment, identifying confounders and hidden clusters within a patient panel, we performed a secondary analysis using data from a prospective cross-sectional and retrospective longitudinal study. Methods: Multivariate and multivariable linear models adjusted for covariates (gender, age, and BMI) were stratified by disease-specific treatment indication. During follow-up, pooled observations were compared from the same patient presenting either with or without self-reported symptoms. Statistical analysis was extended to multilevel models to derive intra-class correlation coefficients and reliability measures during follow-up. Results: Equilibria between TSH, FT4, and FT3 serum concentrations in 342 patients were examined by treatment indication (benign goiter, autoimmune thyroiditis, thyroid carcinoma), consequently displaying complex interactive response patterns. Seventy-seven patients treated with LT4 and monitored for thyroid carcinoma presented, in association with changes in LT4 dose, either with hypothyroid symptoms or symptom-free. Significant biochemical differences appeared between the different presentations. Leveled trajectories by subject to relief from hypothyroid symptoms differed significantly, indicating distinct responses, and denying a single shared outcome. These were formally defined by a high coefficient of the intraclass correlation (ICC1, exceeding 0.60 in all thyroid parameters) during follow-up on multiple visits at the same LT4 dose, when lacking symptoms. The intra-personal clusters were clearly differentiated from random variability by random group resampling. Symptomatic change in these patients was strongly associated with serum FT3, but not with FT4 or TSH concentrations. In 25 patients transitioning from asymptomatic to symptomatically hyperthyroid, FT3 concentrations remained within the reference limits, whilst at the same time marked biochemical differences were apparent between the presentations. Conclusions: Considerable intra-individual clustering occurred in the biochemical and symptomatic responses to LT4 treatment, implying statistically multileveled response groups. Unmasking individual differences in the averaged treatment response hereby highlights clinically distinguishable subgroups within an indiscriminate patient panel. This, through well-designed larger clinical trials will better target the different therapeutic needs of individual patients.

The role of functional thyroid capacity in pituitary thyroid feedback regulation

BACKGROUND

Thyroid feedback regulation and equilibria between thyroid hormones differ in the presence or absence of a functioning thyroid remnant.

MATERIALS AND METHODS

This study examines the relationship between the sensitivity of TSH feedback and thyroid capacity in untreated patients with thyroid autoimmune disease (n = 86) and healthy controls (n = 271). Functional capacity was estimated at maximum TSH stimulation, and pituitary TSH response was FT4-standardised with two established indices, the TSH index and the thyrotroph thyroid hormone resistance index.

RESULTS

The two indices correlated inversely with thyroid volume and functional thyroid capacity. Relationships were shifted upwards in patients with thyroid autoimmune disease. This positioned patients with thyroid autoimmune disease predominantly at the lower capacity range and upper part of TSH index. The relationship was modulated by serum FT3 concentrations, shifting 0.19 [95%CI: 0.12, 0.26] mIU/L per pmol FT3 increase. FT3 correlated with TSH index in total group ( τ  = 0.09, P = 0.009) and both subgroups. FT3 levels were maintained despite a substantial capacity loss by progressively increasing conversion rates of T3 from T4, only collapsing at capacities below <1.5 pmol/s.

CONCLUSION

Functional thyroid capacity and preferential T3 generation are essential elements in adjusting the sensitivity of hypothalamic-pituitary-thyroid feedback control and balancing system equilibria. This suggests that the indirect regulatory role of glandular T3 co-secretion exceeds its quantitative contribution to the thyroid hormone pool. Implications for clinical practice extend to the diagnostic use of TSH in patients with impaired thyroid reserve.

Lessons from Randomised Clinical Trials for Triiodothyronine Treatment of Hypothyroidism: Have They Achieved Their Objectives?

Randomised controlled trials are deemed to be the strongest class of evidence in evidence-based medicine. Failure of trials to prove superiority of T3/T4 combination therapy over standard LT4 monotherapy has greatly influenced guidelines, while not resolving the ongoing debate. Novel studies have recently produced more evidence from the examination of homeostatic equilibria in humans and experimental treatment protocols in animals. This has exacerbated a serious disagreement with evidence from the clinical trials. We contrasted the weight of statistical evidence against strong physiological counterarguments. Revisiting this controversy, we identify areas of improvement for trial design related to validation and sensitivity of QoL instruments, patient selection, statistical power, collider stratification bias, and response heterogeneity to treatment. Given the high individuality expressed by thyroid hormones, their interrelationships, and shifted comfort zones, the response to LT4 treatment produces a statistical amalgamation bias (Simpson's paradox), which has a key influence on interpretation. In addition to drug efficacy, as tested by RCTs, efficiency in clinical practice and safety profiles requires reevaluation. Accordingly, results from RCTs remain ambiguous and should therefore not prevail over physiologically based counterarguments. In giving more weight to other forms of valid evidence which contradict key assumptions of historic trials, current treatment options should remain open and rely on personalised biochemical treatment targets. Optimal treatment choices should be guided by strict requirements of organizations such as the FDA, demanding treatment effects to be estimated under actual conditions of use. Various improvements in design and analysis are recommended for future randomised controlled T3/T4 combination trials.