1
|
Maćešić S, Novakovic K. Diffusion-driven instabilities in the BT-GN oscillatory carbonylation reaction network. CHAOS (WOODBURY, N.Y.) 2024; 34:073152. [PMID: 39047161 DOI: 10.1063/5.0211536] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/31/2024] [Accepted: 07/01/2024] [Indexed: 07/27/2024]
Abstract
This study explores the role of diffusion in creating instabilities in the Bruk Temkin-Gorodsky Novakovic (BT-GN) oscillatory carbonylation reaction network. Stoichiometric network analysis and numerical methods revealed the presence of two destabilizing feedback cycles responsible for these instabilities. Analysis of a spatially uniform system showed that the saddle-node bifurcation can be simulated within the reaction network. The introduction of diffusion results in two types of instabilities: one occurs when a spatially uniform system is already unstable, leading to a reaction-diffusion front; and another involves diffusion-driven instabilities where introducing diffusion destabilizes a stable spatially uniform system. Slower PdI2 diffusion plays a key role in inducing these instabilities. Equations describing conditions for the emergence of the instabilities in both cases were derived.
Collapse
Affiliation(s)
- Stevan Maćešić
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, Belgrade, Serbia
| | - Katarina Novakovic
- School of Engineering, Newcastle University, Newcastle upon Tyne NE1 7RU, United Kingdom
| |
Collapse
|
2
|
Kolar-Anić L, Čupić Ž, Maćešić S, Ivanović-Šašić A, Dietrich JW. Modelling of the thyroid hormone synthesis as a part of nonlinear reaction mechanism with feedback. Comput Biol Med 2023; 160:106980. [PMID: 37141650 DOI: 10.1016/j.compbiomed.2023.106980] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2022] [Revised: 04/04/2023] [Accepted: 04/23/2023] [Indexed: 05/06/2023]
Abstract
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 T4 with respect to T3 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.
Collapse
Affiliation(s)
| | - Željko Čupić
- University of Belgrade, Institute of Chemistry, Technology and Metallurgy, Department of Catalysis and Chemical Engineering, Serbia.
| | - Stevan Maćešić
- University of Belgrade, Faculty of Physical Chemistry, Serbia
| | - Ana Ivanović-Šašić
- University of Belgrade, Institute of Chemistry, Technology and Metallurgy, Department of Catalysis and Chemical Engineering, Serbia
| | - Johannes W Dietrich
- Diabetes, Endocrinology and Metabolism Section, Department of Internal Medicine I, St. Josef Hospital, Ruhr University Bochum, Bochum, NRW, Germany; Diabetes Centre Bochum/Hattingen, St. Elisabeth Hospital Blankenstein, Hattingen, NRW, Germany; Centre for Rare Endocrine Diseases, Ruhr Centre for Rare Diseaeses (CeSER), Ruhr University of Bochum and Witten/Herdecke University, Bochum, NRW, Germany; Centre for Diabetes Technology, Catholic Hospitals Bochum, Hattingen, NRW, Germany
| |
Collapse
|
3
|
|
4
|
Improvement of the stoichiometric network analysis for determination of instability conditions of complex nonlinear reaction systems. Chem Eng Sci 2010. [DOI: 10.1016/j.ces.2010.03.008] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
5
|
Schmitz G, Kolar-Anić LZ, Anić SR, Čupić ŽD. Stoichiometric Network Analysis and Associated Dimensionless Kinetic Equations. Application to a Model of the Bray−Liebhafsky Reaction. J Phys Chem A 2008; 112:13452-7. [DOI: 10.1021/jp8056674] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Guy Schmitz
- Faculté des Sciences Appliquées, Universté Libre de Bruxelles, CP165/63, Av. F. Roosevelt 50, 1050 Bruxelles, Belgium, Faculty of Physical Chemistry, University of Belgrade, P.O.Box 47, 11000 Belgrade, Serbia, and Department of Catalysis and Chemical Engineering, IChTM, University of Belgrade, Njegoševa 12, 11000 Belgrade, Serbia
| | - Ljiljana Z. Kolar-Anić
- Faculté des Sciences Appliquées, Universté Libre de Bruxelles, CP165/63, Av. F. Roosevelt 50, 1050 Bruxelles, Belgium, Faculty of Physical Chemistry, University of Belgrade, P.O.Box 47, 11000 Belgrade, Serbia, and Department of Catalysis and Chemical Engineering, IChTM, University of Belgrade, Njegoševa 12, 11000 Belgrade, Serbia
| | - Slobodan R. Anić
- Faculté des Sciences Appliquées, Universté Libre de Bruxelles, CP165/63, Av. F. Roosevelt 50, 1050 Bruxelles, Belgium, Faculty of Physical Chemistry, University of Belgrade, P.O.Box 47, 11000 Belgrade, Serbia, and Department of Catalysis and Chemical Engineering, IChTM, University of Belgrade, Njegoševa 12, 11000 Belgrade, Serbia
| | - Željko D. Čupić
- Faculté des Sciences Appliquées, Universté Libre de Bruxelles, CP165/63, Av. F. Roosevelt 50, 1050 Bruxelles, Belgium, Faculty of Physical Chemistry, University of Belgrade, P.O.Box 47, 11000 Belgrade, Serbia, and Department of Catalysis and Chemical Engineering, IChTM, University of Belgrade, Njegoševa 12, 11000 Belgrade, Serbia
| |
Collapse
|
6
|
Samadi-Maybodi A, Akhoondi R. Studies of visible oscillating chemiluminescence in luminol–H2O2–KSCN–CuSO4 system using (2-hydroxyethyl) trimethylammonium hydroxide. LUMINESCENCE 2008; 23:42-8. [DOI: 10.1002/bio.1015] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
7
|
Samadi-Maybodi A, Naseri N, Chaichi MJ. Fluorimetry Studies of Oscillating Chemilumnescence in the Luminol= H2O2= KSCN= CuSO4= TMAOH System. J Fluoresc 2005; 15:117-22. [PMID: 15883765 DOI: 10.1007/s10895-005-2518-2] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2004] [Accepted: 12/02/2004] [Indexed: 10/25/2022]
Abstract
Oscillating chemical reactions are complex systems, involving a large number of chemical species. In oscillating chemical reaction, some species, usually a reaction intermediate, exhibit fluctuation in its concentration. In this report, visible oscillating chemiluminescence produced by the addition of luminol (3-aminophthalhydrazide) to the oscillating system of H2O2-KSCN-CuSO4-TMAOH was investigated using spectrofluorimetry method. The effects of ingredient concentration of the oscillating system and complexing agents like citric acid and cysteine on the behavior of the oscillating system were investigated. Moreover, the influence of nonaqueous solvents such as ethanol and ethylene glycol has been studied.
Collapse
|
8
|
Katare S, Caruthers JM, Delgass WN, Venkatasubramanian V. An Intelligent System for Reaction Kinetic Modeling and Catalyst Design. Ind Eng Chem Res 2004. [DOI: 10.1021/ie034067h] [Citation(s) in RCA: 58] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Santhoji Katare
- Center for Catalyst Design, School of Chemical Engineering, Purdue University, West Lafayette, Indiana 47907
| | - James M. Caruthers
- Center for Catalyst Design, School of Chemical Engineering, Purdue University, West Lafayette, Indiana 47907
| | - W. Nicholas Delgass
- Center for Catalyst Design, School of Chemical Engineering, Purdue University, West Lafayette, Indiana 47907
| | - Venkat Venkatasubramanian
- Center for Catalyst Design, School of Chemical Engineering, Purdue University, West Lafayette, Indiana 47907
| |
Collapse
|
9
|
Samadi-Maybodi A, Ourad SM. Studies of visible oscillating chemiluminescence with a luminol-H2O2-KSCN-CuSO4-NaOH system in batch reactor. LUMINESCENCE 2003; 18:42-8. [PMID: 12536379 DOI: 10.1002/bio.705] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
Oscillating chemical reactions are complex systems involving a large number of chemical species. In oscillating chemical reactions some species, usually reaction intermediates, exhibit fluctuation in concentration. Visible oscillating chemiluminescence, produced by the addition of luminol (3-aminophthalhydrazide) to the oscillating system H(2)O(2)-KSCN-CuSO(4)-NaOH, was investigated. In this study the effect of varying the concentration of H(2)O(2), KSCN, CuSO(4), NaOH and luminol was investigated in a batch reactor. We showed that the concentration of all components involved in the oscillating chemilumenscent reaction influenced the light intensity and the oscillation period.
Collapse
|
10
|
Clarke BL, Jiang W. Method for deriving Hopf and saddle‐node bifurcation hypersurfaces and application to a model of the Belousov–Zhabotinskii system. J Chem Phys 1993. [DOI: 10.1063/1.466073] [Citation(s) in RCA: 22] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
11
|
|
12
|
Clarke BL. General method for simplifying chemical networks while preserving overall stoichiometry in reduced mechanisms. J Chem Phys 1992. [DOI: 10.1063/1.463911] [Citation(s) in RCA: 23] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|