TY  - CONF
AU  - Maksimović, Jelena
AU  - Ivanović-Šašić, Ana
AU  - Maćešić, Stevan
AU  - Čupić, Željko
AU  - Kolar-Anić, Ljiljana
PY  - 2023
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/7431
AB  - The impact of L-tyrosine amino acid on the kinetics of the BL oscillatory reaction was investigated
under closed reactor conditions. The study was focused on examining the sensitivity of the BL
reaction matrix to tyrosine perturbations. A high sensitivity of the BL matrix to very low tyrosine
concentrations was observed.
PB  - Serbian Society of Mechanics
C3  - 9th International Congress of the Serbian Society of Mechanics, Book of Proceedings, July 5-7, 2023, Vrnjačka Banja, Serbia
T1  - L-tyrosine influence on the reaction kinetics of iodatehydrogen peroxide oscillatory reaction
SP  - 457
EP  - 458
UR  - https://hdl.handle.net/21.15107/rcub_cer_7431
ER  - 

@conference{
author = "Maksimović, Jelena and Ivanović-Šašić, Ana and Maćešić, Stevan and Čupić, Željko and Kolar-Anić, Ljiljana",
year = "2023",
abstract = "The impact of L-tyrosine amino acid on the kinetics of the BL oscillatory reaction was investigated
under closed reactor conditions. The study was focused on examining the sensitivity of the BL
reaction matrix to tyrosine perturbations. A high sensitivity of the BL matrix to very low tyrosine
concentrations was observed.",
publisher = "Serbian Society of Mechanics",
journal = "9th International Congress of the Serbian Society of Mechanics, Book of Proceedings, July 5-7, 2023, Vrnjačka Banja, Serbia",
title = "L-tyrosine influence on the reaction kinetics of iodatehydrogen peroxide oscillatory reaction",
pages = "457-458",
url = "https://hdl.handle.net/21.15107/rcub_cer_7431"
}

Maksimović, J., Ivanović-Šašić, A., Maćešić, S., Čupić, Ž.,& Kolar-Anić, L.. (2023). L-tyrosine influence on the reaction kinetics of iodatehydrogen peroxide oscillatory reaction. in 9th International Congress of the Serbian Society of Mechanics, Book of Proceedings, July 5-7, 2023, Vrnjačka Banja, Serbia
Serbian Society of Mechanics., 457-458.
https://hdl.handle.net/21.15107/rcub_cer_7431

Maksimović J, Ivanović-Šašić A, Maćešić S, Čupić Ž, Kolar-Anić L. L-tyrosine influence on the reaction kinetics of iodatehydrogen peroxide oscillatory reaction. in 9th International Congress of the Serbian Society of Mechanics, Book of Proceedings, July 5-7, 2023, Vrnjačka Banja, Serbia. 2023;:457-458.
https://hdl.handle.net/21.15107/rcub_cer_7431 .

Maksimović, Jelena, Ivanović-Šašić, Ana, Maćešić, Stevan, Čupić, Željko, Kolar-Anić, Ljiljana, "L-tyrosine influence on the reaction kinetics of iodatehydrogen peroxide oscillatory reaction" in 9th International Congress of the Serbian Society of Mechanics, Book of Proceedings, July 5-7, 2023, Vrnjačka Banja, Serbia (2023):457-458,
https://hdl.handle.net/21.15107/rcub_cer_7431 .

TY  - CONF
AU  - Ivanović-Šašić, Ana
AU  - Maćešić, Stevan
AU  - Maksimović, Jelena
AU  - Čupić, Željko
AU  - Kolar-Anić, Ljiljana
PY  - 2023
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/7430
AB  - It is well known that almost all living or biological systems are naturally in the oscillatory dynamic states and can be considered as biochemical reaction systems. These oscillatory dynamic states can be caused by internal self-organized phenomena, but also by external periodic variations of temperature, light, food, or seasonal changes. The hypothalamic-pituitary-thyroid (HPT) axis is one such nonlinear system with feedback that is always in an oscillatory dynamic state and L-tyrosine is its main representative.
The biological importance of L-tyrosine interactions with iodine species was the motivation for modelling of the oscillatory dynamics in the Bray-Liebhafsky (BL) reaction perturbed by L- tyrosine. Also, direct experimental investigation of metabolic processes in the human body is extremely complex to be done, and therefore any alternative approach is of great importance. Therefore, the BL reaction has the potential to be used as a model of the biological system due to certain characteristics shared with the considered processes; it is characterized by its oscillatory dynamics and based on the chemistry of hydrogen peroxide and iodine compounds commonly present in the thyroid gland, where L-tyrosine is iodinated.
The impact of L-tyrosine on the dynamics of the Bray-Lienbhafsky oscillatory reaction was investigated numerically using the proposed model. The study was focused on the examination of the sensitivity of the BL reaction to L-tyrosine perturbation. The obtained results indicated possible pathways of influence
PB  - Kragujevac, Serbia : Institute for Information Technologies, University of Kragujevac
C3  - 2nd International Conference on Chemo and Bioinformatics ICCBIKG 2023, Book of proceedings, September 28-29, 2023, Kragujevac, Serbia
T1  - Numerical simulations of the oscillatory dynamics in the Bray-Liebhafsky reaction perturbed by L-tyrosine
VL  - 612
SP  - 609
DO  - 10.46793/ICCBI23.609IS
ER  - 

@conference{
author = "Ivanović-Šašić, Ana and Maćešić, Stevan and Maksimović, Jelena and Čupić, Željko and Kolar-Anić, Ljiljana",
year = "2023",
abstract = "It is well known that almost all living or biological systems are naturally in the oscillatory dynamic states and can be considered as biochemical reaction systems. These oscillatory dynamic states can be caused by internal self-organized phenomena, but also by external periodic variations of temperature, light, food, or seasonal changes. The hypothalamic-pituitary-thyroid (HPT) axis is one such nonlinear system with feedback that is always in an oscillatory dynamic state and L-tyrosine is its main representative.
The biological importance of L-tyrosine interactions with iodine species was the motivation for modelling of the oscillatory dynamics in the Bray-Liebhafsky (BL) reaction perturbed by L- tyrosine. Also, direct experimental investigation of metabolic processes in the human body is extremely complex to be done, and therefore any alternative approach is of great importance. Therefore, the BL reaction has the potential to be used as a model of the biological system due to certain characteristics shared with the considered processes; it is characterized by its oscillatory dynamics and based on the chemistry of hydrogen peroxide and iodine compounds commonly present in the thyroid gland, where L-tyrosine is iodinated.
The impact of L-tyrosine on the dynamics of the Bray-Lienbhafsky oscillatory reaction was investigated numerically using the proposed model. The study was focused on the examination of the sensitivity of the BL reaction to L-tyrosine perturbation. The obtained results indicated possible pathways of influence",
publisher = "Kragujevac, Serbia : Institute for Information Technologies, University of Kragujevac",
journal = "2nd International Conference on Chemo and Bioinformatics ICCBIKG 2023, Book of proceedings, September 28-29, 2023, Kragujevac, Serbia",
title = "Numerical simulations of the oscillatory dynamics in the Bray-Liebhafsky reaction perturbed by L-tyrosine",
volume = "612",
pages = "609",
doi = "10.46793/ICCBI23.609IS"
}

Ivanović-Šašić, A., Maćešić, S., Maksimović, J., Čupić, Ž.,& Kolar-Anić, L.. (2023). Numerical simulations of the oscillatory dynamics in the Bray-Liebhafsky reaction perturbed by L-tyrosine. in 2nd International Conference on Chemo and Bioinformatics ICCBIKG 2023, Book of proceedings, September 28-29, 2023, Kragujevac, Serbia
Kragujevac, Serbia : Institute for Information Technologies, University of Kragujevac., 612, 609.
https://doi.org/10.46793/ICCBI23.609IS

Ivanović-Šašić A, Maćešić S, Maksimović J, Čupić Ž, Kolar-Anić L. Numerical simulations of the oscillatory dynamics in the Bray-Liebhafsky reaction perturbed by L-tyrosine. in 2nd International Conference on Chemo and Bioinformatics ICCBIKG 2023, Book of proceedings, September 28-29, 2023, Kragujevac, Serbia. 2023;612:609.
doi:10.46793/ICCBI23.609IS .

Ivanović-Šašić, Ana, Maćešić, Stevan, Maksimović, Jelena, Čupić, Željko, Kolar-Anić, Ljiljana, "Numerical simulations of the oscillatory dynamics in the Bray-Liebhafsky reaction perturbed by L-tyrosine" in 2nd International Conference on Chemo and Bioinformatics ICCBIKG 2023, Book of proceedings, September 28-29, 2023, Kragujevac, Serbia, 612 (2023):609,
https://doi.org/10.46793/ICCBI23.609IS . .

TY  - JOUR
AU  - Kolar-Anić, Ljiljana
AU  - Čupić, Željko
AU  - Maćešić, Stevan
AU  - Ivanović-Šašić, Ana
AU  - Dietrich, Johannes
PY  - 2023
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/6416
AB  - 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.
PB  - Elsevier
T2  - Computers in Biology and Medicine
T1  - Modelling of the thyroid hormone synthesis as a part of nonlinear reaction mechanism with feedback
VL  - 160
SP  - 106980
DO  - 10.1016/j.compbiomed.2023.106980
ER  - 

@article{
author = "Kolar-Anić, Ljiljana and Čupić, Željko and Maćešić, Stevan and Ivanović-Šašić, Ana and Dietrich, Johannes",
year = "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.",
publisher = "Elsevier",
journal = "Computers in Biology and Medicine",
title = "Modelling of the thyroid hormone synthesis as a part of nonlinear reaction mechanism with feedback",
volume = "160",
pages = "106980",
doi = "10.1016/j.compbiomed.2023.106980"
}

Kolar-Anić, L., Čupić, Ž., Maćešić, S., Ivanović-Šašić, A.,& Dietrich, J.. (2023). Modelling of the thyroid hormone synthesis as a part of nonlinear reaction mechanism with feedback. in Computers in Biology and Medicine
Elsevier., 160, 106980.
https://doi.org/10.1016/j.compbiomed.2023.106980

Kolar-Anić L, Čupić Ž, Maćešić S, Ivanović-Šašić A, Dietrich J. Modelling of the thyroid hormone synthesis as a part of nonlinear reaction mechanism with feedback. in Computers in Biology and Medicine. 2023;160:106980.
doi:10.1016/j.compbiomed.2023.106980 .

Kolar-Anić, Ljiljana, Čupić, Željko, Maćešić, Stevan, Ivanović-Šašić, Ana, Dietrich, Johannes, "Modelling of the thyroid hormone synthesis as a part of nonlinear reaction mechanism with feedback" in Computers in Biology and Medicine, 160 (2023):106980,
https://doi.org/10.1016/j.compbiomed.2023.106980 . .

TY  - JOUR
AU  - Maćešić, Stevan
AU  - Čupić, Željko
AU  - Kolar-Anić, Ljiljana
PY  - 2023
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/7185
AB  - The effect of diffusion on the steady-state stability of an oscillatory chemical reaction model was investigated using stoichiometric network analysis and numerical simulations. Under both spatially uniform and nonuniform conditions, steady-state stability was investigated. Under spatially uniform conditions, the model can simulate oscillatory dynamics by passing through the Andronov-Hopf bifurcation. When diffusion is introduced into the system, the results have shown that two scenarios through which instabilities can occur are possible. Either, oscillations may be caused by the same instability as it was in homogeneous case, or, diffusion may cause new type of instability. Using the exponent polytope method, we derived a system of inequalities that describes the conditions for the emergence of both, oscillations, and diffusion-driven instabilities.
PB  - Elsevier
T2  - Chaos, Solitons & Fractals
T1  - Effect of diffusion on steady state stability of an oscillatory reaction model
VL  - 174
SP  - 113783
DO  - 10.1016/j.chaos.2023.113783
ER  - 

@article{
author = "Maćešić, Stevan and Čupić, Željko and Kolar-Anić, Ljiljana",
year = "2023",
abstract = "The effect of diffusion on the steady-state stability of an oscillatory chemical reaction model was investigated using stoichiometric network analysis and numerical simulations. Under both spatially uniform and nonuniform conditions, steady-state stability was investigated. Under spatially uniform conditions, the model can simulate oscillatory dynamics by passing through the Andronov-Hopf bifurcation. When diffusion is introduced into the system, the results have shown that two scenarios through which instabilities can occur are possible. Either, oscillations may be caused by the same instability as it was in homogeneous case, or, diffusion may cause new type of instability. Using the exponent polytope method, we derived a system of inequalities that describes the conditions for the emergence of both, oscillations, and diffusion-driven instabilities.",
publisher = "Elsevier",
journal = "Chaos, Solitons & Fractals",
title = "Effect of diffusion on steady state stability of an oscillatory reaction model",
volume = "174",
pages = "113783",
doi = "10.1016/j.chaos.2023.113783"
}

Maćešić, S., Čupić, Ž.,& Kolar-Anić, L.. (2023). Effect of diffusion on steady state stability of an oscillatory reaction model. in Chaos, Solitons & Fractals
Elsevier., 174, 113783.
https://doi.org/10.1016/j.chaos.2023.113783

Maćešić S, Čupić Ž, Kolar-Anić L. Effect of diffusion on steady state stability of an oscillatory reaction model. in Chaos, Solitons & Fractals. 2023;174:113783.
doi:10.1016/j.chaos.2023.113783 .

Maćešić, Stevan, Čupić, Željko, Kolar-Anić, Ljiljana, "Effect of diffusion on steady state stability of an oscillatory reaction model" in Chaos, Solitons & Fractals, 174 (2023):113783,
https://doi.org/10.1016/j.chaos.2023.113783 . .

TY  - JOUR
AU  - Čupić, Željko
AU  - Maćešić, Stevan
AU  - Anić, Slobodan
AU  - Kolar-Anić, Ljiljana
AU  - Ivanović-Šašić, Ana
AU  - Novaković, Katarina
PY  - 2022
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/5293
AB  - The study presents a mathematical model of Bruk Temkin-Gorodsky Novakovic (BT-GN) reaction system employing palladium-catalyzed oscillatory carbonylation reaction of mono alkyne-terminated poly(ethylene glycol) methyl ether. The proposed model advances simulation capabilities of this particular chemical oscillator but also BT-GN reactions in general. The model was verified using experimental data where good agreement is achieved and existence of complex pH oscillations like burst oscillations was matched. Furthermore, the model was used to predict reaction conditions capable of producing further complexities and initial conditions that would lead to desired period of pH oscillations. Having such a powerful tool, enhances our capabilities to study and further develop BT-GN oscillators with a reduced experimental effort.
PB  - Springer
T2  - Reaction Kinetics, Mechanisms and Catalysis
T1  - Oscillatory carbonylation of poly(ethylene glycol)methyl ether acetylene. Improved model of reaction mechanism
VL  - 135
IS  - 1
SP  - 3
EP  - 14
DO  - 10.1007/s11144-021-02148-9
ER  - 

@article{
author = "Čupić, Željko and Maćešić, Stevan and Anić, Slobodan and Kolar-Anić, Ljiljana and Ivanović-Šašić, Ana and Novaković, Katarina",
year = "2022",
abstract = "The study presents a mathematical model of Bruk Temkin-Gorodsky Novakovic (BT-GN) reaction system employing palladium-catalyzed oscillatory carbonylation reaction of mono alkyne-terminated poly(ethylene glycol) methyl ether. The proposed model advances simulation capabilities of this particular chemical oscillator but also BT-GN reactions in general. The model was verified using experimental data where good agreement is achieved and existence of complex pH oscillations like burst oscillations was matched. Furthermore, the model was used to predict reaction conditions capable of producing further complexities and initial conditions that would lead to desired period of pH oscillations. Having such a powerful tool, enhances our capabilities to study and further develop BT-GN oscillators with a reduced experimental effort.",
publisher = "Springer",
journal = "Reaction Kinetics, Mechanisms and Catalysis",
title = "Oscillatory carbonylation of poly(ethylene glycol)methyl ether acetylene. Improved model of reaction mechanism",
volume = "135",
number = "1",
pages = "3-14",
doi = "10.1007/s11144-021-02148-9"
}

Čupić, Ž., Maćešić, S., Anić, S., Kolar-Anić, L., Ivanović-Šašić, A.,& Novaković, K.. (2022). Oscillatory carbonylation of poly(ethylene glycol)methyl ether acetylene. Improved model of reaction mechanism. in Reaction Kinetics, Mechanisms and Catalysis
Springer., 135(1), 3-14.
https://doi.org/10.1007/s11144-021-02148-9

Čupić Ž, Maćešić S, Anić S, Kolar-Anić L, Ivanović-Šašić A, Novaković K. Oscillatory carbonylation of poly(ethylene glycol)methyl ether acetylene. Improved model of reaction mechanism. in Reaction Kinetics, Mechanisms and Catalysis. 2022;135(1):3-14.
doi:10.1007/s11144-021-02148-9 .

Čupić, Željko, Maćešić, Stevan, Anić, Slobodan, Kolar-Anić, Ljiljana, Ivanović-Šašić, Ana, Novaković, Katarina, "Oscillatory carbonylation of poly(ethylene glycol)methyl ether acetylene. Improved model of reaction mechanism" in Reaction Kinetics, Mechanisms and Catalysis, 135, no. 1 (2022):3-14,
https://doi.org/10.1007/s11144-021-02148-9 . .

TY  - CONF
AU  - Stojiljković, Aleksandra
AU  - Čupić, Željko
AU  - Maćešić, Stevan
AU  - Ivanović-Šašić, Ana
AU  - Kolar-Anić, Ljiljana
PY  - 2022
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/6413
AB  - The previously proposed stoichiometric model of the Hypothalamic-Pituitary-Adrenal (HPA)
axis activity that took into account arginine vasopressin (AVP), has been further developed to
emulate ultradian oscillations of corticotropin-releasing hormone (CRH) and AVP. With this
aim, additional coupling of HPA consisting hormones was introduced into this model by
reaction between CRH and cortisol (CORT). How additional coupling of hormones affects HPA
axis ultradian dynamics and reflects on ultradian oscillations of AVP and CRH concentrations
was examined by using numerical simulations and bifurcation analysis. Results show that the
rate constant of newly incorporated reaction alone is sufficient to be adjusted only for CRH to
exhibit oscillations with optimally prominent amplitudes. Also, oscillation frequencies of CRH
were found to be in accordance with findings in the literature under all investigated conditions
PB  - Society of Physical Chemists of Serbia
C3  - Proceedings - 16th International Conference on Fundamental and Applied Aspects of Physical Chemistry, Physical Chemistry 2022, September, 26-30, 2022, Belgrade, Serbia
T1  - Ultradian oscillations of corticotropin-releasing hormone (crh) and arginine vasopressin (avp) in modelling of hypothalamic-pituitary-adrenal axis: influence of feedback loop between crh and cortisol
VL  - I
SP  - 207
EP  - 210
UR  - https://hdl.handle.net/21.15107/rcub_cer_6413
ER  - 

@conference{
author = "Stojiljković, Aleksandra and Čupić, Željko and Maćešić, Stevan and Ivanović-Šašić, Ana and Kolar-Anić, Ljiljana",
year = "2022",
abstract = "The previously proposed stoichiometric model of the Hypothalamic-Pituitary-Adrenal (HPA)
axis activity that took into account arginine vasopressin (AVP), has been further developed to
emulate ultradian oscillations of corticotropin-releasing hormone (CRH) and AVP. With this
aim, additional coupling of HPA consisting hormones was introduced into this model by
reaction between CRH and cortisol (CORT). How additional coupling of hormones affects HPA
axis ultradian dynamics and reflects on ultradian oscillations of AVP and CRH concentrations
was examined by using numerical simulations and bifurcation analysis. Results show that the
rate constant of newly incorporated reaction alone is sufficient to be adjusted only for CRH to
exhibit oscillations with optimally prominent amplitudes. Also, oscillation frequencies of CRH
were found to be in accordance with findings in the literature under all investigated conditions",
publisher = "Society of Physical Chemists of Serbia",
journal = "Proceedings - 16th International Conference on Fundamental and Applied Aspects of Physical Chemistry, Physical Chemistry 2022, September, 26-30, 2022, Belgrade, Serbia",
title = "Ultradian oscillations of corticotropin-releasing hormone (crh) and arginine vasopressin (avp) in modelling of hypothalamic-pituitary-adrenal axis: influence of feedback loop between crh and cortisol",
volume = "I",
pages = "207-210",
url = "https://hdl.handle.net/21.15107/rcub_cer_6413"
}

Stojiljković, A., Čupić, Ž., Maćešić, S., Ivanović-Šašić, A.,& Kolar-Anić, L.. (2022). Ultradian oscillations of corticotropin-releasing hormone (crh) and arginine vasopressin (avp) in modelling of hypothalamic-pituitary-adrenal axis: influence of feedback loop between crh and cortisol. in Proceedings - 16th International Conference on Fundamental and Applied Aspects of Physical Chemistry, Physical Chemistry 2022, September, 26-30, 2022, Belgrade, Serbia
Society of Physical Chemists of Serbia., I, 207-210.
https://hdl.handle.net/21.15107/rcub_cer_6413

Stojiljković A, Čupić Ž, Maćešić S, Ivanović-Šašić A, Kolar-Anić L. Ultradian oscillations of corticotropin-releasing hormone (crh) and arginine vasopressin (avp) in modelling of hypothalamic-pituitary-adrenal axis: influence of feedback loop between crh and cortisol. in Proceedings - 16th International Conference on Fundamental and Applied Aspects of Physical Chemistry, Physical Chemistry 2022, September, 26-30, 2022, Belgrade, Serbia. 2022;I:207-210.
https://hdl.handle.net/21.15107/rcub_cer_6413 .

Stojiljković, Aleksandra, Čupić, Željko, Maćešić, Stevan, Ivanović-Šašić, Ana, Kolar-Anić, Ljiljana, "Ultradian oscillations of corticotropin-releasing hormone (crh) and arginine vasopressin (avp) in modelling of hypothalamic-pituitary-adrenal axis: influence of feedback loop between crh and cortisol" in Proceedings - 16th International Conference on Fundamental and Applied Aspects of Physical Chemistry, Physical Chemistry 2022, September, 26-30, 2022, Belgrade, Serbia, I (2022):207-210,
https://hdl.handle.net/21.15107/rcub_cer_6413 .

TY  - JOUR
AU  - Negrojević, Luka
AU  - Lončar, Aleksandar
AU  - Maksimović, Jelena
AU  - Anić, Slobodan
AU  - Čupić, Željko
AU  - Kolar-Anić, Ljiljana
AU  - Pejić, Nataša
PY  - 2022
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/5253
AB  - A kinetic method for the determination of l-tyrosine (Tyr) by analyte pulse perturbation, caused by different amounts of Tyr, in the Bray–Liebhafsky (BL) oscillatory reaction realized in a continuous-flow stirred tank reactor, is proposed. For such purposes, the BL oscillating reaction was kept in a stable non-equilibrium stationary state close to the bifurcation point and was used as the matrix reaction system. Under optimum reaction conditions, the linear relationship between the kinetic parameters, Δτ2 (the period between first two oscillations that appear after applied perturbation) as well as τend (the time elapsed between the perturbation of the BL reaction by Tyr and the termination of the oscillatory phase) and the Tyr concentrations is obtained over the ranges 1.1 × 10− 6 − 9.2 × 10− 6 mol L− 1 and 1.1 × 10− 6 − 1.4 × 10− 5 mol L− 1, with the limit of detections of 6.6 × 10− 7 mol L− 1 and 6.4 × 10− 7 mol L− 1. The described method that relies on a simple instrumental set-up, has been successfully applied to the determination of Tyr in a dietary supplement. Some aspects of the possible mechanism of Tyr action on the BL oscillator are discussed.
PB  - Springer
T2  - Reaction Kinetics, Mechanisms and Catalysis
T1  - Bray–Liebhafsky oscillatory reaction in a continuous-flow stirred tank reactor as the matrix system for determination of tyrosine
VL  - 135
SP  - 1147
EP  - 1162
DO  - 10.1007/s11144-021-02130-5
ER  - 

@article{
author = "Negrojević, Luka and Lončar, Aleksandar and Maksimović, Jelena and Anić, Slobodan and Čupić, Željko and Kolar-Anić, Ljiljana and Pejić, Nataša",
year = "2022",
abstract = "A kinetic method for the determination of l-tyrosine (Tyr) by analyte pulse perturbation, caused by different amounts of Tyr, in the Bray–Liebhafsky (BL) oscillatory reaction realized in a continuous-flow stirred tank reactor, is proposed. For such purposes, the BL oscillating reaction was kept in a stable non-equilibrium stationary state close to the bifurcation point and was used as the matrix reaction system. Under optimum reaction conditions, the linear relationship between the kinetic parameters, Δτ2 (the period between first two oscillations that appear after applied perturbation) as well as τend (the time elapsed between the perturbation of the BL reaction by Tyr and the termination of the oscillatory phase) and the Tyr concentrations is obtained over the ranges 1.1 × 10− 6 − 9.2 × 10− 6 mol L− 1 and 1.1 × 10− 6 − 1.4 × 10− 5 mol L− 1, with the limit of detections of 6.6 × 10− 7 mol L− 1 and 6.4 × 10− 7 mol L− 1. The described method that relies on a simple instrumental set-up, has been successfully applied to the determination of Tyr in a dietary supplement. Some aspects of the possible mechanism of Tyr action on the BL oscillator are discussed.",
publisher = "Springer",
journal = "Reaction Kinetics, Mechanisms and Catalysis",
title = "Bray–Liebhafsky oscillatory reaction in a continuous-flow stirred tank reactor as the matrix system for determination of tyrosine",
volume = "135",
pages = "1147-1162",
doi = "10.1007/s11144-021-02130-5"
}

Negrojević, L., Lončar, A., Maksimović, J., Anić, S., Čupić, Ž., Kolar-Anić, L.,& Pejić, N.. (2022). Bray–Liebhafsky oscillatory reaction in a continuous-flow stirred tank reactor as the matrix system for determination of tyrosine. in Reaction Kinetics, Mechanisms and Catalysis
Springer., 135, 1147-1162.
https://doi.org/10.1007/s11144-021-02130-5

Negrojević L, Lončar A, Maksimović J, Anić S, Čupić Ž, Kolar-Anić L, Pejić N. Bray–Liebhafsky oscillatory reaction in a continuous-flow stirred tank reactor as the matrix system for determination of tyrosine. in Reaction Kinetics, Mechanisms and Catalysis. 2022;135:1147-1162.
doi:10.1007/s11144-021-02130-5 .

Negrojević, Luka, Lončar, Aleksandar, Maksimović, Jelena, Anić, Slobodan, Čupić, Željko, Kolar-Anić, Ljiljana, Pejić, Nataša, "Bray–Liebhafsky oscillatory reaction in a continuous-flow stirred tank reactor as the matrix system for determination of tyrosine" in Reaction Kinetics, Mechanisms and Catalysis, 135 (2022):1147-1162,
https://doi.org/10.1007/s11144-021-02130-5 . .

TY  - JOUR
AU  - Stojiljković, Aleksandra S.
AU  - Čupić, Željko
AU  - Maćešić, Stevan
AU  - Ivanović - Šašić, Ana
AU  - Kolar-Anić, Ljiljana
PY  - 2022
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/5390
AB  - Numerous studies on humans and animals have indicated that the corticotrophin-releasing hormone (CRH) and arginine vasopressin (AVP) stimulate both individually and synergistically secretion of adrenocorticotropic hormone (ACTH) by corticotropic cells in anterior pituitary. With aim to characterize and better comprehend the mechanisms underlying the effects of AVP on Hypothalamic-Pituitary-Adrenal (HPA) axis ultradian dynamics, AVP is here incorporated into our previously proposed stoichiometric model of HPA axis in humans. This extended nonlinear network reaction model took into account AVP by: reaction steps associated with two separate inflows of AVP into pituitary portal system, that is synthesized and released from hypothalamic parvocellular and magnocellular neuronal populations, as well as summarized reaction steps related to its individual and synergistic action with CRH on corticotropic cells. To explore the properties of extended model and its capacity to emulate the effects of AVP, nonlinear dynamical systems theory and bifurcation analyses based on numerical simulations were utilized to determine the dependence of ultradian oscillations on rate constants of the inflows of CRH and AVP from parvocellular neuronal populations, the conditions under which dynamical transitions occur due to their synergistic action and, moreover, the types of these transitions. The results show that under certain conditions, HPA system could enter into oscillatory dynamic states from stable steady state and vice versa under the influence of synergy reaction rate constant. Transitions between these dynamical states were always through supercritical Andronov-Hopf bifurcation point. Also, results revealed the conditions under which amplitudes of ultradian oscillations could increase several-fold due to CRH and AVP synergistic stimulation of ACTH secretion in accordance with results reported in the literature. Moreover, results showed experimentally observed superiority of CRH as a stimulator of ACTH secretion compared to AVP in humans. The proposed model can be very useful in studies related to the role of AVP and its synergistic action with CRH in life-threatening circumstances such as acute homeostasis dynamic crisis, autoimmune inflammations or severe hypovolemia requiring instant or several-days sustained corticosteroid excess levels. Moreover, the model can be helpful for investigations of indirect AVP-induced HPA activity by exogenously administered AVP used in therapeutic treatment.
PB  - Frontiers Media S.A.
T2  - Frontiers in Endokrinology
T1  - Influence of arginine vasopressin on the ultradian dynamics of Hypothalamic-Pituitary-Adrenal axis
VL  - 13
SP  - 976323
DO  - 10.3389/fendo.2022.976323
ER  - 

@article{
author = "Stojiljković, Aleksandra S. and Čupić, Željko and Maćešić, Stevan and Ivanović - Šašić, Ana and Kolar-Anić, Ljiljana",
year = "2022",
abstract = "Numerous studies on humans and animals have indicated that the corticotrophin-releasing hormone (CRH) and arginine vasopressin (AVP) stimulate both individually and synergistically secretion of adrenocorticotropic hormone (ACTH) by corticotropic cells in anterior pituitary. With aim to characterize and better comprehend the mechanisms underlying the effects of AVP on Hypothalamic-Pituitary-Adrenal (HPA) axis ultradian dynamics, AVP is here incorporated into our previously proposed stoichiometric model of HPA axis in humans. This extended nonlinear network reaction model took into account AVP by: reaction steps associated with two separate inflows of AVP into pituitary portal system, that is synthesized and released from hypothalamic parvocellular and magnocellular neuronal populations, as well as summarized reaction steps related to its individual and synergistic action with CRH on corticotropic cells. To explore the properties of extended model and its capacity to emulate the effects of AVP, nonlinear dynamical systems theory and bifurcation analyses based on numerical simulations were utilized to determine the dependence of ultradian oscillations on rate constants of the inflows of CRH and AVP from parvocellular neuronal populations, the conditions under which dynamical transitions occur due to their synergistic action and, moreover, the types of these transitions. The results show that under certain conditions, HPA system could enter into oscillatory dynamic states from stable steady state and vice versa under the influence of synergy reaction rate constant. Transitions between these dynamical states were always through supercritical Andronov-Hopf bifurcation point. Also, results revealed the conditions under which amplitudes of ultradian oscillations could increase several-fold due to CRH and AVP synergistic stimulation of ACTH secretion in accordance with results reported in the literature. Moreover, results showed experimentally observed superiority of CRH as a stimulator of ACTH secretion compared to AVP in humans. The proposed model can be very useful in studies related to the role of AVP and its synergistic action with CRH in life-threatening circumstances such as acute homeostasis dynamic crisis, autoimmune inflammations or severe hypovolemia requiring instant or several-days sustained corticosteroid excess levels. Moreover, the model can be helpful for investigations of indirect AVP-induced HPA activity by exogenously administered AVP used in therapeutic treatment.",
publisher = "Frontiers Media S.A.",
journal = "Frontiers in Endokrinology",
title = "Influence of arginine vasopressin on the ultradian dynamics of Hypothalamic-Pituitary-Adrenal axis",
volume = "13",
pages = "976323",
doi = "10.3389/fendo.2022.976323"
}

Stojiljković, A. S., Čupić, Ž., Maćešić, S., Ivanović - Šašić, A.,& Kolar-Anić, L.. (2022). Influence of arginine vasopressin on the ultradian dynamics of Hypothalamic-Pituitary-Adrenal axis. in Frontiers in Endokrinology
Frontiers Media S.A.., 13, 976323.
https://doi.org/10.3389/fendo.2022.976323

Stojiljković AS, Čupić Ž, Maćešić S, Ivanović - Šašić A, Kolar-Anić L. Influence of arginine vasopressin on the ultradian dynamics of Hypothalamic-Pituitary-Adrenal axis. in Frontiers in Endokrinology. 2022;13:976323.
doi:10.3389/fendo.2022.976323 .

Stojiljković, Aleksandra S., Čupić, Željko, Maćešić, Stevan, Ivanović - Šašić, Ana, Kolar-Anić, Ljiljana, "Influence of arginine vasopressin on the ultradian dynamics of Hypothalamic-Pituitary-Adrenal axis" in Frontiers in Endokrinology, 13 (2022):976323,
https://doi.org/10.3389/fendo.2022.976323 . .

TY  - CONF
AU  - Ivanović-Šašić, Ana
AU  - Maćešić, Stevan
AU  - Čupić, Željko
AU  - Kolar-Anić, Ljiljana
PY  - 2021
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/6403
AB  - The Bray-Liebhafsky (BL) is one of the most analyzed oscillatory reaction both experimentally and 
numerically. Most of the experimentally obtained dynamical states of this reaction realized in a 
continuously fed well stirred tank reactor (CSTR) are successfully simulated. Beside others, 
numerous structured chaotic dynamical states were obtained between each two periodic states in the 
period doubling rout to chaos with respect to specific flow rate as the control parameter. It was an 
universal scenario throughout the whole mixed-mode region, as well as throughout other mixed-mode 
regions obtained under different initial conditions. However, the intermittent oscillations consistent 
of chaotic mixture of large-amplitude relaxation oscillations, grouped in bursts and small-amplitude 
sinusoidal ones or even quiescent parts between them known as gaps were also generated 
experimentally in the Bray–Liebhafsky reaction by varying different parameters such as temperature, 
flow rate or reactant concentrations under CSTR conditions. Nevertheless, it will be shown here that 
intermittent oscillations can be simulated by already published model of the BL reaction network.
PB  - Society of Physical Chemists of Serbia
C3  - Proceedings - 15th International Conference on Fundamental and Applied Aspects of Physical Chemistry, Physical Chemistry 2021, September 20-24 2021, Belgrade, Serbia
T1  - Various dinamical states in the Bray-Liebhafsky oscillatory reaction- from periodicity to intermittent chaos
VL  - I
SP  - 199
EP  - 199
UR  - https://hdl.handle.net/21.15107/rcub_cer_6403
ER  - 

@conference{
author = "Ivanović-Šašić, Ana and Maćešić, Stevan and Čupić, Željko and Kolar-Anić, Ljiljana",
year = "2021",
abstract = "The Bray-Liebhafsky (BL) is one of the most analyzed oscillatory reaction both experimentally and 
numerically. Most of the experimentally obtained dynamical states of this reaction realized in a 
continuously fed well stirred tank reactor (CSTR) are successfully simulated. Beside others, 
numerous structured chaotic dynamical states were obtained between each two periodic states in the 
period doubling rout to chaos with respect to specific flow rate as the control parameter. It was an 
universal scenario throughout the whole mixed-mode region, as well as throughout other mixed-mode 
regions obtained under different initial conditions. However, the intermittent oscillations consistent 
of chaotic mixture of large-amplitude relaxation oscillations, grouped in bursts and small-amplitude 
sinusoidal ones or even quiescent parts between them known as gaps were also generated 
experimentally in the Bray–Liebhafsky reaction by varying different parameters such as temperature, 
flow rate or reactant concentrations under CSTR conditions. Nevertheless, it will be shown here that 
intermittent oscillations can be simulated by already published model of the BL reaction network.",
publisher = "Society of Physical Chemists of Serbia",
journal = "Proceedings - 15th International Conference on Fundamental and Applied Aspects of Physical Chemistry, Physical Chemistry 2021, September 20-24 2021, Belgrade, Serbia",
title = "Various dinamical states in the Bray-Liebhafsky oscillatory reaction- from periodicity to intermittent chaos",
volume = "I",
pages = "199-199",
url = "https://hdl.handle.net/21.15107/rcub_cer_6403"
}

Ivanović-Šašić, A., Maćešić, S., Čupić, Ž.,& Kolar-Anić, L.. (2021). Various dinamical states in the Bray-Liebhafsky oscillatory reaction- from periodicity to intermittent chaos. in Proceedings - 15th International Conference on Fundamental and Applied Aspects of Physical Chemistry, Physical Chemistry 2021, September 20-24 2021, Belgrade, Serbia
Society of Physical Chemists of Serbia., I, 199-199.
https://hdl.handle.net/21.15107/rcub_cer_6403

Ivanović-Šašić A, Maćešić S, Čupić Ž, Kolar-Anić L. Various dinamical states in the Bray-Liebhafsky oscillatory reaction- from periodicity to intermittent chaos. in Proceedings - 15th International Conference on Fundamental and Applied Aspects of Physical Chemistry, Physical Chemistry 2021, September 20-24 2021, Belgrade, Serbia. 2021;I:199-199.
https://hdl.handle.net/21.15107/rcub_cer_6403 .

Ivanović-Šašić, Ana, Maćešić, Stevan, Čupić, Željko, Kolar-Anić, Ljiljana, "Various dinamical states in the Bray-Liebhafsky oscillatory reaction- from periodicity to intermittent chaos" in Proceedings - 15th International Conference on Fundamental and Applied Aspects of Physical Chemistry, Physical Chemistry 2021, September 20-24 2021, Belgrade, Serbia, I (2021):199-199,
https://hdl.handle.net/21.15107/rcub_cer_6403 .

TY  - CONF
AU  - Čupić, Željko
AU  - Ivanović-Šašić, Ana
AU  - Maćešić, Stevan
AU  - Anić, Slobodan
AU  - Kolar-Anić, Ljiljana
PY  - 2021
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/6404
AB  - The first discovered homogeneous oscillatory reaction was the Bray-Liebhafsky (BL) one, 
described in a paper published exactly 100 years ago. However, the applicability of oscillatory 
reactions in chemical computing was recently discovered. Here we intend to expose the native 
computing concept applied to intermittent states of the BL reaction, because we believe that 
this particular state may have some advantages. For this purpose, numerical simulations will be
used based on the known model. Sequences of perturbations will be introduced by adding iodate 
(IO3
-
) and hydrogen peroxide (H2O2), separately, as well as in various combinations with one 
another. It will be shown that dynamic states obtained after perturbations with same species 
depend very much on the sequence in which these species were used in perturbations. 
Additionally, it will be shown that obtained dynamic states shift the system from chaotic 
intermittent dynamic state to different complex periodic states. Hence, the applicability of the 
BL reaction system in chemical computing was demonstrated.
PB  - Kragujevac, Serbia : Institute for Information Technologies, University of Kragujevac
C3  - Book of proceedings - 1st International Conference on Chemo and BioInformatics  ICCBIKG 2021, October 26-27, 2021, Kragujevac, Serbia
T1  - Applicability of Bray-Liebhafsky reaction for chemical computing
SP  - 431
EP  - 434
DO  - 10.46793/ICCBI21.431C
ER  - 

@conference{
author = "Čupić, Željko and Ivanović-Šašić, Ana and Maćešić, Stevan and Anić, Slobodan and Kolar-Anić, Ljiljana",
year = "2021",
abstract = "The first discovered homogeneous oscillatory reaction was the Bray-Liebhafsky (BL) one, 
described in a paper published exactly 100 years ago. However, the applicability of oscillatory 
reactions in chemical computing was recently discovered. Here we intend to expose the native 
computing concept applied to intermittent states of the BL reaction, because we believe that 
this particular state may have some advantages. For this purpose, numerical simulations will be
used based on the known model. Sequences of perturbations will be introduced by adding iodate 
(IO3
-
) and hydrogen peroxide (H2O2), separately, as well as in various combinations with one 
another. It will be shown that dynamic states obtained after perturbations with same species 
depend very much on the sequence in which these species were used in perturbations. 
Additionally, it will be shown that obtained dynamic states shift the system from chaotic 
intermittent dynamic state to different complex periodic states. Hence, the applicability of the 
BL reaction system in chemical computing was demonstrated.",
publisher = "Kragujevac, Serbia : Institute for Information Technologies, University of Kragujevac",
journal = "Book of proceedings - 1st International Conference on Chemo and BioInformatics  ICCBIKG 2021, October 26-27, 2021, Kragujevac, Serbia",
title = "Applicability of Bray-Liebhafsky reaction for chemical computing",
pages = "431-434",
doi = "10.46793/ICCBI21.431C"
}

Čupić, Ž., Ivanović-Šašić, A., Maćešić, S., Anić, S.,& Kolar-Anić, L.. (2021). Applicability of Bray-Liebhafsky reaction for chemical computing. in Book of proceedings - 1st International Conference on Chemo and BioInformatics  ICCBIKG 2021, October 26-27, 2021, Kragujevac, Serbia
Kragujevac, Serbia : Institute for Information Technologies, University of Kragujevac., 431-434.
https://doi.org/10.46793/ICCBI21.431C

Čupić Ž, Ivanović-Šašić A, Maćešić S, Anić S, Kolar-Anić L. Applicability of Bray-Liebhafsky reaction for chemical computing. in Book of proceedings - 1st International Conference on Chemo and BioInformatics  ICCBIKG 2021, October 26-27, 2021, Kragujevac, Serbia. 2021;:431-434.
doi:10.46793/ICCBI21.431C .

Čupić, Željko, Ivanović-Šašić, Ana, Maćešić, Stevan, Anić, Slobodan, Kolar-Anić, Ljiljana, "Applicability of Bray-Liebhafsky reaction for chemical computing" in Book of proceedings - 1st International Conference on Chemo and BioInformatics  ICCBIKG 2021, October 26-27, 2021, Kragujevac, Serbia (2021):431-434,
https://doi.org/10.46793/ICCBI21.431C . .

TY  - CONF
AU  - Stojiljković, Aleksandra
AU  - Stanojević, Ana
AU  - Ivanović-Šašić, Ana
AU  - Čupić, Željko
AU  - Kolar-Anić, Ljiljana
PY  - 2021
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/6410
AB  - An upgraded model of the Hypothalamic-Pituitary-Adrenal (HPA) axis has been developed that 
advances our previously proposed low-dimensional HPA model by including the effects of arginine 
vasopressin (AVP) that is a key modulator of HPA axis function. The upgraded model allows us to 
emulate AVP effects on HPA axis dynamics individually and in synergy with the corticotropin-releasing hormone (CRH). In this work, we examine how coupling of the circadian function through 
summarised reaction steps describing CRH and AVP biosynthesis in the same neuronal cell group of the hypothalamic paraventricular nucleus (PVN) affects HPA axis dynamics. Results of numerical simulations show that coupling of the circadian function through both, CRH and AVP summarised biosynthesis reaction steps simultaneously, emulates best the HPA axis dynamics, in line with literature findings.
PB  - Society of Physical Chemists of Serbia
C3  - Proceedings - 15th International Conference on Fundamental and Applied Aspects of Physical Chemistry, Physical Chemistry 2021, September 20-24 2021, Belgrade, Serbia
T1  - Circadian rhythm function coupling to the upgraded hypothalamic-pituitary-adrenal (hpa) axis with incorporated arginine vasopressin
VL  - I
SP  - 246
EP  - 249
UR  - https://hdl.handle.net/21.15107/rcub_cer_6410
ER  - 

@conference{
author = "Stojiljković, Aleksandra and Stanojević, Ana and Ivanović-Šašić, Ana and Čupić, Željko and Kolar-Anić, Ljiljana",
year = "2021",
abstract = "An upgraded model of the Hypothalamic-Pituitary-Adrenal (HPA) axis has been developed that 
advances our previously proposed low-dimensional HPA model by including the effects of arginine 
vasopressin (AVP) that is a key modulator of HPA axis function. The upgraded model allows us to 
emulate AVP effects on HPA axis dynamics individually and in synergy with the corticotropin-releasing hormone (CRH). In this work, we examine how coupling of the circadian function through 
summarised reaction steps describing CRH and AVP biosynthesis in the same neuronal cell group of the hypothalamic paraventricular nucleus (PVN) affects HPA axis dynamics. Results of numerical simulations show that coupling of the circadian function through both, CRH and AVP summarised biosynthesis reaction steps simultaneously, emulates best the HPA axis dynamics, in line with literature findings.",
publisher = "Society of Physical Chemists of Serbia",
journal = "Proceedings - 15th International Conference on Fundamental and Applied Aspects of Physical Chemistry, Physical Chemistry 2021, September 20-24 2021, Belgrade, Serbia",
title = "Circadian rhythm function coupling to the upgraded hypothalamic-pituitary-adrenal (hpa) axis with incorporated arginine vasopressin",
volume = "I",
pages = "246-249",
url = "https://hdl.handle.net/21.15107/rcub_cer_6410"
}

Stojiljković, A., Stanojević, A., Ivanović-Šašić, A., Čupić, Ž.,& Kolar-Anić, L.. (2021). Circadian rhythm function coupling to the upgraded hypothalamic-pituitary-adrenal (hpa) axis with incorporated arginine vasopressin. in Proceedings - 15th International Conference on Fundamental and Applied Aspects of Physical Chemistry, Physical Chemistry 2021, September 20-24 2021, Belgrade, Serbia
Society of Physical Chemists of Serbia., I, 246-249.
https://hdl.handle.net/21.15107/rcub_cer_6410

Stojiljković A, Stanojević A, Ivanović-Šašić A, Čupić Ž, Kolar-Anić L. Circadian rhythm function coupling to the upgraded hypothalamic-pituitary-adrenal (hpa) axis with incorporated arginine vasopressin. in Proceedings - 15th International Conference on Fundamental and Applied Aspects of Physical Chemistry, Physical Chemistry 2021, September 20-24 2021, Belgrade, Serbia. 2021;I:246-249.
https://hdl.handle.net/21.15107/rcub_cer_6410 .

Stojiljković, Aleksandra, Stanojević, Ana, Ivanović-Šašić, Ana, Čupić, Željko, Kolar-Anić, Ljiljana, "Circadian rhythm function coupling to the upgraded hypothalamic-pituitary-adrenal (hpa) axis with incorporated arginine vasopressin" in Proceedings - 15th International Conference on Fundamental and Applied Aspects of Physical Chemistry, Physical Chemistry 2021, September 20-24 2021, Belgrade, Serbia, I (2021):246-249,
https://hdl.handle.net/21.15107/rcub_cer_6410 .

TY  - JOUR
AU  - Bubanja, Itana Nuša
AU  - Ivanović-Šašić, Ana
AU  - Čupić, Željko
AU  - Anić, Slobodan
AU  - Kolar-Anić, Ljiljana
PY  - 2020
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/3701
AB  - Dynamic states with intermittent oscillations consist of a chaotic mixture of large
amplitude relaxation oscillations grouped in bursts, and between them, small-amplitude
sinusoidal oscillations, or even the quiescent parts, known as gaps. In this study,
intermittent dynamic states were generated in Bray–Liebhafsky (BL) oscillatory reaction in
an isothermal continuously-fed, well-stirred tank reactor (CSTR) controled by changes of
specific flow rate. The intermittent states were found between two regular periodic states
and obtained for specific flow rate values from 0.020 to 0.082 min−1 . Phenomenological
analysis based on the quantitative characteristics of intermittent oscillations, as well
as, the largest Lyapunov exponents calculated from experimentally obtained time
series, both indicated the same type of behavior. Namely, fully developed chaos
arises when approaching to the vertical asymptote which is somewhere between two
bifurcations. Hence, this study proposes described route to fully developed chaos in
the Bray-Liebhafsky oscillatory reaction as an explanation for experimentally observed
intermittent dynamics. This is in correlation with our previously obtained results where the
most chaotic intermittent chaos was achieved between the periodic oscillatory dynamic
state and stable steady state, generated in BL under CSTR conditions by varying
temperature and inflow potassium iodate concentration. Moreover, it was shown that,
besides the largest Lyapunov exponent, analysis of chaos in experimentally obtained
intermittent states can be achieved by a simpler approach which involves using the
quantitative characteristics of the BL reaction evolution, that is, the number and length
of gaps and bursts obtained for the various values of specific flow rates.
PB  - Frontiers Media S.A.
T2  - Frontiers in chemistry
T1  - Intermittent Chaos in the CSTR Bray–Liebhafsky Oscillator-Specific Flow Rate Dependence
VL  - 8
SP  - 560274
DO  - 10.3389/fchem.2020.560274
ER  - 

@article{
author = "Bubanja, Itana Nuša and Ivanović-Šašić, Ana and Čupić, Željko and Anić, Slobodan and Kolar-Anić, Ljiljana",
year = "2020",
abstract = "Dynamic states with intermittent oscillations consist of a chaotic mixture of large
amplitude relaxation oscillations grouped in bursts, and between them, small-amplitude
sinusoidal oscillations, or even the quiescent parts, known as gaps. In this study,
intermittent dynamic states were generated in Bray–Liebhafsky (BL) oscillatory reaction in
an isothermal continuously-fed, well-stirred tank reactor (CSTR) controled by changes of
specific flow rate. The intermittent states were found between two regular periodic states
and obtained for specific flow rate values from 0.020 to 0.082 min−1 . Phenomenological
analysis based on the quantitative characteristics of intermittent oscillations, as well
as, the largest Lyapunov exponents calculated from experimentally obtained time
series, both indicated the same type of behavior. Namely, fully developed chaos
arises when approaching to the vertical asymptote which is somewhere between two
bifurcations. Hence, this study proposes described route to fully developed chaos in
the Bray-Liebhafsky oscillatory reaction as an explanation for experimentally observed
intermittent dynamics. This is in correlation with our previously obtained results where the
most chaotic intermittent chaos was achieved between the periodic oscillatory dynamic
state and stable steady state, generated in BL under CSTR conditions by varying
temperature and inflow potassium iodate concentration. Moreover, it was shown that,
besides the largest Lyapunov exponent, analysis of chaos in experimentally obtained
intermittent states can be achieved by a simpler approach which involves using the
quantitative characteristics of the BL reaction evolution, that is, the number and length
of gaps and bursts obtained for the various values of specific flow rates.",
publisher = "Frontiers Media S.A.",
journal = "Frontiers in chemistry",
title = "Intermittent Chaos in the CSTR Bray–Liebhafsky Oscillator-Specific Flow Rate Dependence",
volume = "8",
pages = "560274",
doi = "10.3389/fchem.2020.560274"
}

Bubanja, I. N., Ivanović-Šašić, A., Čupić, Ž., Anić, S.,& Kolar-Anić, L.. (2020). Intermittent Chaos in the CSTR Bray–Liebhafsky Oscillator-Specific Flow Rate Dependence. in Frontiers in chemistry
Frontiers Media S.A.., 8, 560274.
https://doi.org/10.3389/fchem.2020.560274

Bubanja IN, Ivanović-Šašić A, Čupić Ž, Anić S, Kolar-Anić L. Intermittent Chaos in the CSTR Bray–Liebhafsky Oscillator-Specific Flow Rate Dependence. in Frontiers in chemistry. 2020;8:560274.
doi:10.3389/fchem.2020.560274 .

Bubanja, Itana Nuša, Ivanović-Šašić, Ana, Čupić, Željko, Anić, Slobodan, Kolar-Anić, Ljiljana, "Intermittent Chaos in the CSTR Bray–Liebhafsky Oscillator-Specific Flow Rate Dependence" in Frontiers in chemistry, 8 (2020):560274,
https://doi.org/10.3389/fchem.2020.560274 . .

TY  - JOUR
AU  - Maksimović, Jelena
AU  - Čupić, Željko
AU  - Manojlović, Nedeljko
AU  - Đerić, Aleksandra
AU  - Anić, Slobodan
AU  - Kolar-Anić, Ljiljana
PY  - 2020
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/3603
AB  - Two hydroxyanthraquinones, alizarin and purpurin, that have been used throughout
history as a natural pigment, were extracted from roots of Rubia tinctorum from
Serbia. As substances with important chemical activities and, therefore, with wide
applications (for example in dyeing textile fabrics as well as in pharmacy because of
their anti-inflammatory, anti-cancer, antiviral, antimicrobial and antioxidant activities,
etc.) they were analyzed in the kinetically very sensitive Bray–Liebhafsky (BL)
oscillatory reaction. However, although they are both, hydroxyanthraquinones it is
shown that their interactions with BL nonlinear reaction system differ significantly.
Consequently, two different reactions were used to explain the mechanism of their
chemical activities. The numerical simulations based on a standard model of the BL
oscillatory reaction together with proposed reactions due to alizarin/purpurin interactions
with a matrix are correlated with experimental investigations. Moreover, it
is shown that very small amounts of alizarin and purpurin (from about 1 × 10–7 M)
produce the response of the BL matrix such that micro-quantitative analysis based
on the BL oscillatory reaction can be successfully performed in this reaction system.
The linear response of the BL matrix on the presence of alizarin and purpurin
(necessary for microquantitative determination) is analyzed as a function of two
concentration sensitive parameters: pre-oscillatory period τ1 and potential shift after
perturbation ΔE.
PB  - Springer Nature Switzerland AG 2020
T2  - Reaction Kinetics, Mechanisms and Catalysis
T1  - Bray–Liebhafsky oscillatory reaction as the matrix system for the kinetic determination of microquantities of alizarin and purpurin
VL  - 130
SP  - 655
EP  - 688
DO  - 10.1007/s11144-020-01798-5
ER  - 

@article{
author = "Maksimović, Jelena and Čupić, Željko and Manojlović, Nedeljko and Đerić, Aleksandra and Anić, Slobodan and Kolar-Anić, Ljiljana",
year = "2020",
abstract = "Two hydroxyanthraquinones, alizarin and purpurin, that have been used throughout
history as a natural pigment, were extracted from roots of Rubia tinctorum from
Serbia. As substances with important chemical activities and, therefore, with wide
applications (for example in dyeing textile fabrics as well as in pharmacy because of
their anti-inflammatory, anti-cancer, antiviral, antimicrobial and antioxidant activities,
etc.) they were analyzed in the kinetically very sensitive Bray–Liebhafsky (BL)
oscillatory reaction. However, although they are both, hydroxyanthraquinones it is
shown that their interactions with BL nonlinear reaction system differ significantly.
Consequently, two different reactions were used to explain the mechanism of their
chemical activities. The numerical simulations based on a standard model of the BL
oscillatory reaction together with proposed reactions due to alizarin/purpurin interactions
with a matrix are correlated with experimental investigations. Moreover, it
is shown that very small amounts of alizarin and purpurin (from about 1 × 10–7 M)
produce the response of the BL matrix such that micro-quantitative analysis based
on the BL oscillatory reaction can be successfully performed in this reaction system.
The linear response of the BL matrix on the presence of alizarin and purpurin
(necessary for microquantitative determination) is analyzed as a function of two
concentration sensitive parameters: pre-oscillatory period τ1 and potential shift after
perturbation ΔE.",
publisher = "Springer Nature Switzerland AG 2020",
journal = "Reaction Kinetics, Mechanisms and Catalysis",
title = "Bray–Liebhafsky oscillatory reaction as the matrix system for the kinetic determination of microquantities of alizarin and purpurin",
volume = "130",
pages = "655-688",
doi = "10.1007/s11144-020-01798-5"
}

Maksimović, J., Čupić, Ž., Manojlović, N., Đerić, A., Anić, S.,& Kolar-Anić, L.. (2020). Bray–Liebhafsky oscillatory reaction as the matrix system for the kinetic determination of microquantities of alizarin and purpurin. in Reaction Kinetics, Mechanisms and Catalysis
Springer Nature Switzerland AG 2020., 130, 655-688.
https://doi.org/10.1007/s11144-020-01798-5

Maksimović J, Čupić Ž, Manojlović N, Đerić A, Anić S, Kolar-Anić L. Bray–Liebhafsky oscillatory reaction as the matrix system for the kinetic determination of microquantities of alizarin and purpurin. in Reaction Kinetics, Mechanisms and Catalysis. 2020;130:655-688.
doi:10.1007/s11144-020-01798-5 .

Maksimović, Jelena, Čupić, Željko, Manojlović, Nedeljko, Đerić, Aleksandra, Anić, Slobodan, Kolar-Anić, Ljiljana, "Bray–Liebhafsky oscillatory reaction as the matrix system for the kinetic determination of microquantities of alizarin and purpurin" in Reaction Kinetics, Mechanisms and Catalysis, 130 (2020):655-688,
https://doi.org/10.1007/s11144-020-01798-5 . .

TY  - CONF
AU  - Čupić, Željko
AU  - Maćešić, Stevan
AU  - Kolar-Anić, Ljiljana
PY  - 2020
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/3746
UR  - http://imqf.ekof.bg.ac.rs/ConfNonlin/index.html
AB  - Stability analysis of reaction systems is described by the application of the Stoichiometric Network Analysis to the three-variable-autocatalator. Although simple, this model is complex enough to describe complex forms of nonlinear dynamics phenomena, like mixed-mode oscillations and chaos. Therefore, stability analysis of such model is not a trivial task. Using the Stoichiometric Network Analysis for this purpose makes the process clear and leads to the reliable result.
The method is described briefly in few general steps and all of them are further clarified through the application to the chosen example. First, the reaction rates in steady state are decomposed to contributions of independent pathways, called extreme currents. Then, linearized operator is constructed. Finally, through the analysis of the principal minors of the essential part of this operator, simple stability criterion is identified.
PB  - Serbian academy of nonlinear sciences
C3  - Proceedings - 1st Conference on nonlinearity
T1  - Instability region in models of nonlinear reaction systems. The Stoichiometric Network Analysis
SP  - 145
EP  - 159
UR  - https://hdl.handle.net/21.15107/rcub_cer_3746
ER  - 

@conference{
author = "Čupić, Željko and Maćešić, Stevan and Kolar-Anić, Ljiljana",
year = "2020",
abstract = "Stability analysis of reaction systems is described by the application of the Stoichiometric Network Analysis to the three-variable-autocatalator. Although simple, this model is complex enough to describe complex forms of nonlinear dynamics phenomena, like mixed-mode oscillations and chaos. Therefore, stability analysis of such model is not a trivial task. Using the Stoichiometric Network Analysis for this purpose makes the process clear and leads to the reliable result.
The method is described briefly in few general steps and all of them are further clarified through the application to the chosen example. First, the reaction rates in steady state are decomposed to contributions of independent pathways, called extreme currents. Then, linearized operator is constructed. Finally, through the analysis of the principal minors of the essential part of this operator, simple stability criterion is identified.",
publisher = "Serbian academy of nonlinear sciences",
journal = "Proceedings - 1st Conference on nonlinearity",
title = "Instability region in models of nonlinear reaction systems. The Stoichiometric Network Analysis",
pages = "145-159",
url = "https://hdl.handle.net/21.15107/rcub_cer_3746"
}

Čupić, Ž., Maćešić, S.,& Kolar-Anić, L.. (2020). Instability region in models of nonlinear reaction systems. The Stoichiometric Network Analysis. in Proceedings - 1st Conference on nonlinearity
Serbian academy of nonlinear sciences., 145-159.
https://hdl.handle.net/21.15107/rcub_cer_3746

Čupić Ž, Maćešić S, Kolar-Anić L. Instability region in models of nonlinear reaction systems. The Stoichiometric Network Analysis. in Proceedings - 1st Conference on nonlinearity. 2020;:145-159.
https://hdl.handle.net/21.15107/rcub_cer_3746 .

Čupić, Željko, Maćešić, Stevan, Kolar-Anić, Ljiljana, "Instability region in models of nonlinear reaction systems. The Stoichiometric Network Analysis" in Proceedings - 1st Conference on nonlinearity (2020):145-159,
https://hdl.handle.net/21.15107/rcub_cer_3746 .

TY  - JOUR
AU  - Jakšić, Olga
AU  - Jakšić, Zoran
AU  - Rašljić, Milena
AU  - Kolar-Anić, Ljiljana
PY  - 2019
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/2487
AB  - Starting from the fact that monocomponent adsorption, whether modeled by Lagergren or nonlinear Riccati equation, does not sustain oscillations, we speculate about the nature of multiple steady state states in multicomponent adsorption with second-order kinetics and about the possibility that multicomponent adsorption might exhibit oscillating behavior, in order to provide a tool for better discerning possible oscillations from inevitable fluctuations in experimental results or a tool for a better control of adsorption process far from equilibrium. We perform an analysis of stability of binary adsorption with second-order kinetics in multiple ways. We address perturbations around the steady state analytically, first in a classical way, then by introducing Langevin forces and analyzing the reaction flux and cross-correlations, then by applying the stochastic chemical master equation approach, and finally, numerically, by using stochastic simulation algorithms. Our results show that stationary states in this model are stable nodes. Hence, experimental results with purported oscillations in response should be addressed from the point of view of fluctuations and noise analysis.
PB  - Hindawi
T2  - Advances in Mathematical Physics
T1  - On Oscillations and Noise in Multicomponent Adsorption: The Nature of Multiple Stationary States
VL  - 2019
SP  - 7687643
DO  - 10.1155/2019/7687643
ER  - 

@article{
author = "Jakšić, Olga and Jakšić, Zoran and Rašljić, Milena and Kolar-Anić, Ljiljana",
year = "2019",
abstract = "Starting from the fact that monocomponent adsorption, whether modeled by Lagergren or nonlinear Riccati equation, does not sustain oscillations, we speculate about the nature of multiple steady state states in multicomponent adsorption with second-order kinetics and about the possibility that multicomponent adsorption might exhibit oscillating behavior, in order to provide a tool for better discerning possible oscillations from inevitable fluctuations in experimental results or a tool for a better control of adsorption process far from equilibrium. We perform an analysis of stability of binary adsorption with second-order kinetics in multiple ways. We address perturbations around the steady state analytically, first in a classical way, then by introducing Langevin forces and analyzing the reaction flux and cross-correlations, then by applying the stochastic chemical master equation approach, and finally, numerically, by using stochastic simulation algorithms. Our results show that stationary states in this model are stable nodes. Hence, experimental results with purported oscillations in response should be addressed from the point of view of fluctuations and noise analysis.",
publisher = "Hindawi",
journal = "Advances in Mathematical Physics",
title = "On Oscillations and Noise in Multicomponent Adsorption: The Nature of Multiple Stationary States",
volume = "2019",
pages = "7687643",
doi = "10.1155/2019/7687643"
}

Jakšić, O., Jakšić, Z., Rašljić, M.,& Kolar-Anić, L.. (2019). On Oscillations and Noise in Multicomponent Adsorption: The Nature of Multiple Stationary States. in Advances in Mathematical Physics
Hindawi., 2019, 7687643.
https://doi.org/10.1155/2019/7687643

Jakšić O, Jakšić Z, Rašljić M, Kolar-Anić L. On Oscillations and Noise in Multicomponent Adsorption: The Nature of Multiple Stationary States. in Advances in Mathematical Physics. 2019;2019:7687643.
doi:10.1155/2019/7687643 .

Jakšić, Olga, Jakšić, Zoran, Rašljić, Milena, Kolar-Anić, Ljiljana, "On Oscillations and Noise in Multicomponent Adsorption: The Nature of Multiple Stationary States" in Advances in Mathematical Physics, 2019 (2019):7687643,
https://doi.org/10.1155/2019/7687643 . .

TY  - JOUR
AU  - Maćešić, Stevan
AU  - Čupić, Željko
AU  - Novaković, Katarina
AU  - Parker, Julie
AU  - Anić, Slobodan
AU  - Kolar-Anić, Ljiljana
PY  - 2019
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/2999
AB  - The reaction mechanism of palladium-catalysed oscillatory carbonylation of
poly(ethylene glycol)methyl ether acetylene (PEGA) was examined by stoichiometric
network analysis (SNA) of the stability of the dynamic states. The previously published
simplified reaction network proposed to account for the experimentally observed results was
modified in order to obtain a more realistic reaction network for the investigated process. In the modified reaction network, the direct autocatalytic steps were replaced with autocatalytic
loops. In this procedure some complex net reactions had to be resolved, while some new
reactions had to be taken into account, altering the reaction mechanism responsible for the
appearance of oscillation. Moreover, in the modified reaction network, the expressions for
reaction rates correspond formally to their stoichiometry in accordance to mass action
kinetics. The SNA identified the instability condition which was further tested and confirmed
by numerical simulation.
PB  - University of Kragujevac
T2  - Match: Communications in Mathematical and in Computer Chemistry
T1  - Oscillatory Carbonylation of Poly(Ethylene Glycol)Methyl Ether Acetylene. Modelling of Reaction Mechanism and Stoichiometric Network Stability Analysis
VL  - 81
SP  - 5
EP  - 34
UR  - https://hdl.handle.net/21.15107/rcub_cer_2999
ER  - 

@article{
author = "Maćešić, Stevan and Čupić, Željko and Novaković, Katarina and Parker, Julie and Anić, Slobodan and Kolar-Anić, Ljiljana",
year = "2019",
abstract = "The reaction mechanism of palladium-catalysed oscillatory carbonylation of
poly(ethylene glycol)methyl ether acetylene (PEGA) was examined by stoichiometric
network analysis (SNA) of the stability of the dynamic states. The previously published
simplified reaction network proposed to account for the experimentally observed results was
modified in order to obtain a more realistic reaction network for the investigated process. In the modified reaction network, the direct autocatalytic steps were replaced with autocatalytic
loops. In this procedure some complex net reactions had to be resolved, while some new
reactions had to be taken into account, altering the reaction mechanism responsible for the
appearance of oscillation. Moreover, in the modified reaction network, the expressions for
reaction rates correspond formally to their stoichiometry in accordance to mass action
kinetics. The SNA identified the instability condition which was further tested and confirmed
by numerical simulation.",
publisher = "University of Kragujevac",
journal = "Match: Communications in Mathematical and in Computer Chemistry",
title = "Oscillatory Carbonylation of Poly(Ethylene Glycol)Methyl Ether Acetylene. Modelling of Reaction Mechanism and Stoichiometric Network Stability Analysis",
volume = "81",
pages = "5-34",
url = "https://hdl.handle.net/21.15107/rcub_cer_2999"
}

Maćešić, S., Čupić, Ž., Novaković, K., Parker, J., Anić, S.,& Kolar-Anić, L.. (2019). Oscillatory Carbonylation of Poly(Ethylene Glycol)Methyl Ether Acetylene. Modelling of Reaction Mechanism and Stoichiometric Network Stability Analysis. in Match: Communications in Mathematical and in Computer Chemistry
University of Kragujevac., 81, 5-34.
https://hdl.handle.net/21.15107/rcub_cer_2999

Maćešić S, Čupić Ž, Novaković K, Parker J, Anić S, Kolar-Anić L. Oscillatory Carbonylation of Poly(Ethylene Glycol)Methyl Ether Acetylene. Modelling of Reaction Mechanism and Stoichiometric Network Stability Analysis. in Match: Communications in Mathematical and in Computer Chemistry. 2019;81:5-34.
https://hdl.handle.net/21.15107/rcub_cer_2999 .

Maćešić, Stevan, Čupić, Željko, Novaković, Katarina, Parker, Julie, Anić, Slobodan, Kolar-Anić, Ljiljana, "Oscillatory Carbonylation of Poly(Ethylene Glycol)Methyl Ether Acetylene. Modelling of Reaction Mechanism and Stoichiometric Network Stability Analysis" in Match: Communications in Mathematical and in Computer Chemistry, 81 (2019):5-34,
https://hdl.handle.net/21.15107/rcub_cer_2999 .

TY  - JOUR
AU  - Macesic, Stevan
AU  - Čupić, Željko
AU  - Ivanović-Šašić, Ana
AU  - Anić, Slobodan
AU  - Radenkovic, Mirjana
AU  - Pejić, Nataša
AU  - Kolar-Anić, Ljiljana
PY  - 2018
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/2461
AB  - In this paper, we intend to show the importance of the bifurcation analysis in understanding of an oscillatory process. Hence, we use the bifurcation diagram of the Bray-Liebhafsky reaction performed in continuous well-stirred tank reactor under controlled temperature variations for the determination of the activation energies as well as rate constants of particular steps appearing in the kinetic model of oscillatory reaction mechanism. This approach has led us to the development of general procedure for treatment of experimentally obtained data and extracting kinetic parameters from them, which was very important considering that some rate constants of the already proposed model could not be determined experimentally and have to be fitted (or guessed). Also, the proposed approach has the potential to inspire the refinement of already proposed models and the development of a new one that will be able to reproduce experimentally obtained system's dynamical features more successfully. In particular, the dynamic states of the Bray-Liebhafsky oscillatory reaction have been analyzed experimentally and numerically using already proposed model together with qualitative and quantitative analysis of bifurcation diagrams in both cases.
PB  - Springer Netherlands
T2  - Reaction Kinetics, Mechanisms and Catalysis
T1  - Bifurcation analysis: a tool for determining model parameters of the considered process
VL  - 123
IS  - 1
SP  - 31
EP  - 45
DO  - 10.1007/s11144-017-1324-6
ER  - 

@article{
author = "Macesic, Stevan and Čupić, Željko and Ivanović-Šašić, Ana and Anić, Slobodan and Radenkovic, Mirjana and Pejić, Nataša and Kolar-Anić, Ljiljana",
year = "2018",
abstract = "In this paper, we intend to show the importance of the bifurcation analysis in understanding of an oscillatory process. Hence, we use the bifurcation diagram of the Bray-Liebhafsky reaction performed in continuous well-stirred tank reactor under controlled temperature variations for the determination of the activation energies as well as rate constants of particular steps appearing in the kinetic model of oscillatory reaction mechanism. This approach has led us to the development of general procedure for treatment of experimentally obtained data and extracting kinetic parameters from them, which was very important considering that some rate constants of the already proposed model could not be determined experimentally and have to be fitted (or guessed). Also, the proposed approach has the potential to inspire the refinement of already proposed models and the development of a new one that will be able to reproduce experimentally obtained system's dynamical features more successfully. In particular, the dynamic states of the Bray-Liebhafsky oscillatory reaction have been analyzed experimentally and numerically using already proposed model together with qualitative and quantitative analysis of bifurcation diagrams in both cases.",
publisher = "Springer Netherlands",
journal = "Reaction Kinetics, Mechanisms and Catalysis",
title = "Bifurcation analysis: a tool for determining model parameters of the considered process",
volume = "123",
number = "1",
pages = "31-45",
doi = "10.1007/s11144-017-1324-6"
}

Macesic, S., Čupić, Ž., Ivanović-Šašić, A., Anić, S., Radenkovic, M., Pejić, N.,& Kolar-Anić, L.. (2018). Bifurcation analysis: a tool for determining model parameters of the considered process. in Reaction Kinetics, Mechanisms and Catalysis
Springer Netherlands., 123(1), 31-45.
https://doi.org/10.1007/s11144-017-1324-6

Macesic S, Čupić Ž, Ivanović-Šašić A, Anić S, Radenkovic M, Pejić N, Kolar-Anić L. Bifurcation analysis: a tool for determining model parameters of the considered process. in Reaction Kinetics, Mechanisms and Catalysis. 2018;123(1):31-45.
doi:10.1007/s11144-017-1324-6 .

Macesic, Stevan, Čupić, Željko, Ivanović-Šašić, Ana, Anić, Slobodan, Radenkovic, Mirjana, Pejić, Nataša, Kolar-Anić, Ljiljana, "Bifurcation analysis: a tool for determining model parameters of the considered process" in Reaction Kinetics, Mechanisms and Catalysis, 123, no. 1 (2018):31-45,
https://doi.org/10.1007/s11144-017-1324-6 . .

TY  - CONF
AU  - Maksimović, Jelena
AU  - Čupić, Željko
AU  - Anić, Slobodan
AU  - Kolar-Anić, Ljiljana
AU  - Greco, Emanuela
AU  - Cervellati, Rinaldo
AU  - Pagnacco, Maja
PY  - 2018
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/6965
AB  - The effect of different concentration of gallic acid (GA) on the oscillatory behavior of the Briggs-Rauscher (BR) reaction is reported. In this work, the sample was injected before the start of oscillations. The main effect is the appearance of an induction time that increases with the concentration of GA added. Other parameters as a number of oscillations and oscillatory period are 
affected. These latter linearly depend on the concentration and can be used for quantitative GA analytical assay in a GA concentration range 3.00×10-5 M − 4.00×10-4 M. The perturbation is very different from that previously reported, in which samples of phenolics were added to an active BR mixture, indicating that the mechanism of perturbation largely depends on the injection point.
PB  - Society of Physical Chemists of Serbia
C3  - Proceedings - 14th International Conference on Fundamental and Applied Aspects of Physical Chemistry, Physical Chemistry 2018, September 24-28, 2018 Belgrade, Serbia
T1  - Gallic acid effect on the Briggs-Rauscher reaction dynamics
VL  - 1
SP  - 337
EP  - 340
UR  - https://hdl.handle.net/21.15107/rcub_cer_6965
ER  - 

@conference{
author = "Maksimović, Jelena and Čupić, Željko and Anić, Slobodan and Kolar-Anić, Ljiljana and Greco, Emanuela and Cervellati, Rinaldo and Pagnacco, Maja",
year = "2018",
abstract = "The effect of different concentration of gallic acid (GA) on the oscillatory behavior of the Briggs-Rauscher (BR) reaction is reported. In this work, the sample was injected before the start of oscillations. The main effect is the appearance of an induction time that increases with the concentration of GA added. Other parameters as a number of oscillations and oscillatory period are 
affected. These latter linearly depend on the concentration and can be used for quantitative GA analytical assay in a GA concentration range 3.00×10-5 M − 4.00×10-4 M. The perturbation is very different from that previously reported, in which samples of phenolics were added to an active BR mixture, indicating that the mechanism of perturbation largely depends on the injection point.",
publisher = "Society of Physical Chemists of Serbia",
journal = "Proceedings - 14th International Conference on Fundamental and Applied Aspects of Physical Chemistry, Physical Chemistry 2018, September 24-28, 2018 Belgrade, Serbia",
title = "Gallic acid effect on the Briggs-Rauscher reaction dynamics",
volume = "1",
pages = "337-340",
url = "https://hdl.handle.net/21.15107/rcub_cer_6965"
}

Maksimović, J., Čupić, Ž., Anić, S., Kolar-Anić, L., Greco, E., Cervellati, R.,& Pagnacco, M.. (2018). Gallic acid effect on the Briggs-Rauscher reaction dynamics. in Proceedings - 14th International Conference on Fundamental and Applied Aspects of Physical Chemistry, Physical Chemistry 2018, September 24-28, 2018 Belgrade, Serbia
Society of Physical Chemists of Serbia., 1, 337-340.
https://hdl.handle.net/21.15107/rcub_cer_6965

Maksimović J, Čupić Ž, Anić S, Kolar-Anić L, Greco E, Cervellati R, Pagnacco M. Gallic acid effect on the Briggs-Rauscher reaction dynamics. in Proceedings - 14th International Conference on Fundamental and Applied Aspects of Physical Chemistry, Physical Chemistry 2018, September 24-28, 2018 Belgrade, Serbia. 2018;1:337-340.
https://hdl.handle.net/21.15107/rcub_cer_6965 .

Maksimović, Jelena, Čupić, Željko, Anić, Slobodan, Kolar-Anić, Ljiljana, Greco, Emanuela, Cervellati, Rinaldo, Pagnacco, Maja, "Gallic acid effect on the Briggs-Rauscher reaction dynamics" in Proceedings - 14th International Conference on Fundamental and Applied Aspects of Physical Chemistry, Physical Chemistry 2018, September 24-28, 2018 Belgrade, Serbia, 1 (2018):337-340,
https://hdl.handle.net/21.15107/rcub_cer_6965 .

TY  - JOUR
AU  - Čupić, Željko
AU  - Macesic, Stevan
AU  - Novaković, Katarina
AU  - Anić, Slobodan
AU  - Kolar-Anić, Ljiljana
PY  - 2018
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/2431
AB  - Stoichiometric Network Analysis (SNA) is a powerful method that can be used to examine instabilities in modelling a broad range of reaction systems without knowing the explicit values of reaction rate constants. Due to a lack of understanding, SNA is rarely used and its full potential is not yet fulfilled. Using the oscillatory carbonylation of a polymeric substrate [poly(ethylene glycol) methyl ether acetylene] as a case study, in this work, we consider two mathematical methods for the application of SNA to the reaction models when conservation constraints between species have an important role. The first method takes conservation constraints into account and uses only independent intermediate species, while the second method applies to the full set of intermediate species, without the separation of independent and dependent variables. Both methods are used for examination of steady state stability by means of a characteristic polynomial and related Jacobian matrix. It was shown that both methods give the same results. Therefore, as the second method is simpler, we suggest it as a more straightforward method for the applications. Published by AIP Publishing.
PB  - Amer Inst Physics, Melville
T2  - Chaos
T1  - Stoichiometric network analysis of a reaction system with conservation constraints
VL  - 28
IS  - 8
SP  - 083114
DO  - 10.1063/1.5026791
ER  - 

@article{
author = "Čupić, Željko and Macesic, Stevan and Novaković, Katarina and Anić, Slobodan and Kolar-Anić, Ljiljana",
year = "2018",
abstract = "Stoichiometric Network Analysis (SNA) is a powerful method that can be used to examine instabilities in modelling a broad range of reaction systems without knowing the explicit values of reaction rate constants. Due to a lack of understanding, SNA is rarely used and its full potential is not yet fulfilled. Using the oscillatory carbonylation of a polymeric substrate [poly(ethylene glycol) methyl ether acetylene] as a case study, in this work, we consider two mathematical methods for the application of SNA to the reaction models when conservation constraints between species have an important role. The first method takes conservation constraints into account and uses only independent intermediate species, while the second method applies to the full set of intermediate species, without the separation of independent and dependent variables. Both methods are used for examination of steady state stability by means of a characteristic polynomial and related Jacobian matrix. It was shown that both methods give the same results. Therefore, as the second method is simpler, we suggest it as a more straightforward method for the applications. Published by AIP Publishing.",
publisher = "Amer Inst Physics, Melville",
journal = "Chaos",
title = "Stoichiometric network analysis of a reaction system with conservation constraints",
volume = "28",
number = "8",
pages = "083114",
doi = "10.1063/1.5026791"
}

Čupić, Ž., Macesic, S., Novaković, K., Anić, S.,& Kolar-Anić, L.. (2018). Stoichiometric network analysis of a reaction system with conservation constraints. in Chaos
Amer Inst Physics, Melville., 28(8), 083114.
https://doi.org/10.1063/1.5026791

Čupić Ž, Macesic S, Novaković K, Anić S, Kolar-Anić L. Stoichiometric network analysis of a reaction system with conservation constraints. in Chaos. 2018;28(8):083114.
doi:10.1063/1.5026791 .

Čupić, Željko, Macesic, Stevan, Novaković, Katarina, Anić, Slobodan, Kolar-Anić, Ljiljana, "Stoichiometric network analysis of a reaction system with conservation constraints" in Chaos, 28, no. 8 (2018):083114,
https://doi.org/10.1063/1.5026791 . .

TY  - JOUR
AU  - Stanojević, A
AU  - Marković, V M
AU  - Maćešić, Stevan
AU  - Kolar-Anić, Ljiljana
AU  - Vukojević, V
PY  - 2018
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/2920
AB  - The sex hormone testosterone (TTS) and the hypothalamic–pituitary–adrenal (HPA) axis mutually control one another’s activity, wherein TTS suppresses corticotrophin releasing hormone (CRH) stimulated HPA axis activity, whereas the activation of HPA axis has an inhibitory effect on TTS secretion. With an intention to explain these phenomena, a network reaction model is developed from the previously postulated stoichiometric models for HPA activity where main dynamic behaviors are controlled by two catalytic steps (one autocatalytic and one autoinhibitory) with respect to cortisol, both found experimentally. The capacity of the model to emulate TTS effects on HPA axis dynamics and its response to acute CRH-induced stress is examined using numerical simulations. Model predictions are compared with empirically obtained results reported in the literature. Thus, the reaction kinetic examinations of nonlinear biochemical transformations that constitute the HPA axis, including the negative feedback effect of TTS on HPA axis activity, recapitulates the well-established fact that TTS dampens HPA axis basal activity, decreasing both cortisol level and the amplitude of ultradian cortisol oscillations. The model also replicates TTS inhibitory action on the HPA axis response to acute environmental challenges, particularly CRH-induced stress. In addition, kinetic modelling revealed that TTS induced reduction in ultradian cortisol amplitude arises because the system moves towards a supercritical Hopf bifurcation as TTS is being increased. © 2017, The Author(s).
PB  - Springer Netherlands
T2  - Reaction Kinetics, Mechanisms and Catalysis
T1  - Kinetic modelling of testosterone-related differences in the hypothalamic–pituitary–adrenal axis response to stress
VL  - 123
IS  - 1
SP  - 17
EP  - 30
DO  - 10.1007/s11144-017-1315-7
ER  - 

@article{
author = "Stanojević, A and Marković, V M and Maćešić, Stevan and Kolar-Anić, Ljiljana and Vukojević, V",
year = "2018",
abstract = "The sex hormone testosterone (TTS) and the hypothalamic–pituitary–adrenal (HPA) axis mutually control one another’s activity, wherein TTS suppresses corticotrophin releasing hormone (CRH) stimulated HPA axis activity, whereas the activation of HPA axis has an inhibitory effect on TTS secretion. With an intention to explain these phenomena, a network reaction model is developed from the previously postulated stoichiometric models for HPA activity where main dynamic behaviors are controlled by two catalytic steps (one autocatalytic and one autoinhibitory) with respect to cortisol, both found experimentally. The capacity of the model to emulate TTS effects on HPA axis dynamics and its response to acute CRH-induced stress is examined using numerical simulations. Model predictions are compared with empirically obtained results reported in the literature. Thus, the reaction kinetic examinations of nonlinear biochemical transformations that constitute the HPA axis, including the negative feedback effect of TTS on HPA axis activity, recapitulates the well-established fact that TTS dampens HPA axis basal activity, decreasing both cortisol level and the amplitude of ultradian cortisol oscillations. The model also replicates TTS inhibitory action on the HPA axis response to acute environmental challenges, particularly CRH-induced stress. In addition, kinetic modelling revealed that TTS induced reduction in ultradian cortisol amplitude arises because the system moves towards a supercritical Hopf bifurcation as TTS is being increased. © 2017, The Author(s).",
publisher = "Springer Netherlands",
journal = "Reaction Kinetics, Mechanisms and Catalysis",
title = "Kinetic modelling of testosterone-related differences in the hypothalamic–pituitary–adrenal axis response to stress",
volume = "123",
number = "1",
pages = "17-30",
doi = "10.1007/s11144-017-1315-7"
}

Stanojević, A., Marković, V. M., Maćešić, S., Kolar-Anić, L.,& Vukojević, V.. (2018). Kinetic modelling of testosterone-related differences in the hypothalamic–pituitary–adrenal axis response to stress. in Reaction Kinetics, Mechanisms and Catalysis
Springer Netherlands., 123(1), 17-30.
https://doi.org/10.1007/s11144-017-1315-7

Stanojević A, Marković VM, Maćešić S, Kolar-Anić L, Vukojević V. Kinetic modelling of testosterone-related differences in the hypothalamic–pituitary–adrenal axis response to stress. in Reaction Kinetics, Mechanisms and Catalysis. 2018;123(1):17-30.
doi:10.1007/s11144-017-1315-7 .

Stanojević, A, Marković, V M, Maćešić, Stevan, Kolar-Anić, Ljiljana, Vukojević, V, "Kinetic modelling of testosterone-related differences in the hypothalamic–pituitary–adrenal axis response to stress" in Reaction Kinetics, Mechanisms and Catalysis, 123, no. 1 (2018):17-30,
https://doi.org/10.1007/s11144-017-1315-7 . .

TY  - JOUR
AU  - Stanojevic, Ana
AU  - Marković, Vladimir M.
AU  - Čupić, Željko
AU  - Kolar-Anić, Ljiljana
AU  - Vukojevic, Vladana
PY  - 2018
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/2365
AB  - Stress is a physiological reaction of an organism to a demand for change that is imposed by external factors or is coming from within by way of physiological strains or self-perceived mental and/or emotional threats (internal factors). It manifests itself through the sudden release of a flood of hormones, including corticosteroids, into the blood, which rouse the body for action. Normally, stress is beneficial, but when lasting or being very strong, it causes major damage to our mind and body. Despite intense research, we still do not understand fully how the stress response axis, whose main function is to respond to challenges while maintaining the normal physiological balance, loses under prolonged exposure to stressors its capacity to maintain homeostasis. Recent applications of mathematical modelling and dynamical systems theory have enabled us to emulate complex neurochemical transformations that underlie the stress response, and help us to acquire deeper understanding of this dynamical regulatory network.
PB  - Elsevier Sci Ltd, Oxford
T2  - Current Opinion in Chemical Engineering
T1  - Advances in mathematical modelling of the hypothalamic-pituitary-adrenal (HPA) axis dynamics and the neuroendocrine response to stress
VL  - 21
SP  - 84
EP  - 95
DO  - 10.1016/j.coche.2018.04.003
ER  - 

@article{
author = "Stanojevic, Ana and Marković, Vladimir M. and Čupić, Željko and Kolar-Anić, Ljiljana and Vukojevic, Vladana",
year = "2018",
abstract = "Stress is a physiological reaction of an organism to a demand for change that is imposed by external factors or is coming from within by way of physiological strains or self-perceived mental and/or emotional threats (internal factors). It manifests itself through the sudden release of a flood of hormones, including corticosteroids, into the blood, which rouse the body for action. Normally, stress is beneficial, but when lasting or being very strong, it causes major damage to our mind and body. Despite intense research, we still do not understand fully how the stress response axis, whose main function is to respond to challenges while maintaining the normal physiological balance, loses under prolonged exposure to stressors its capacity to maintain homeostasis. Recent applications of mathematical modelling and dynamical systems theory have enabled us to emulate complex neurochemical transformations that underlie the stress response, and help us to acquire deeper understanding of this dynamical regulatory network.",
publisher = "Elsevier Sci Ltd, Oxford",
journal = "Current Opinion in Chemical Engineering",
title = "Advances in mathematical modelling of the hypothalamic-pituitary-adrenal (HPA) axis dynamics and the neuroendocrine response to stress",
volume = "21",
pages = "84-95",
doi = "10.1016/j.coche.2018.04.003"
}

Stanojevic, A., Marković, V. M., Čupić, Ž., Kolar-Anić, L.,& Vukojevic, V.. (2018). Advances in mathematical modelling of the hypothalamic-pituitary-adrenal (HPA) axis dynamics and the neuroendocrine response to stress. in Current Opinion in Chemical Engineering
Elsevier Sci Ltd, Oxford., 21, 84-95.
https://doi.org/10.1016/j.coche.2018.04.003

Stanojevic A, Marković VM, Čupić Ž, Kolar-Anić L, Vukojevic V. Advances in mathematical modelling of the hypothalamic-pituitary-adrenal (HPA) axis dynamics and the neuroendocrine response to stress. in Current Opinion in Chemical Engineering. 2018;21:84-95.
doi:10.1016/j.coche.2018.04.003 .

Stanojevic, Ana, Marković, Vladimir M., Čupić, Željko, Kolar-Anić, Ljiljana, Vukojevic, Vladana, "Advances in mathematical modelling of the hypothalamic-pituitary-adrenal (HPA) axis dynamics and the neuroendocrine response to stress" in Current Opinion in Chemical Engineering, 21 (2018):84-95,
https://doi.org/10.1016/j.coche.2018.04.003 . .

TY  - JOUR
AU  - Bubanja, Itana Nuša
AU  - Ivanović-Šašić, Ana
AU  - Čupić, Željko
AU  - Anić, Slobodan
AU  - Kolar-Anić, Ljiljana
PY  - 2017
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/2159
AB  - Chaotic dynamic states with intermittent oscillations were generated in a Bray-Liebhafsky (BL) oscillatory reaction in an isothermal open reactor i.e., in the continuously-fed well-stirred tank reactor (CSTR) when the inflow concentration of potassium iodate was the control parameter. They are found between periodic oscillations obtained when [KIO3](0)  LT  3.00 x 10(-2) M and stable steady states when [KIO3](0) > 4.10 x 10(-2) M. It was shown that the most chaotic states obtained experimentally somewhere in the middle of this region are in high correlation with results obtained by means of largest Lyapunov exponents and phenomenological analysis based on the quantitative characteristics of intermittent oscillations.
PB  - Maik Nauka/Interperiodica/Springer, New York
T2  - Russian Journal of Physical Chemistry A
T1  - Intermittent Chaos in the Bray-Liebhafsky Oscillator. Dependence of Dynamic States on the Iodate Concentration
VL  - 91
IS  - 13
SP  - 2525
EP  - 2529
DO  - 10.1134/S0036024417130076
ER  - 

@article{
author = "Bubanja, Itana Nuša and Ivanović-Šašić, Ana and Čupić, Željko and Anić, Slobodan and Kolar-Anić, Ljiljana",
year = "2017",
abstract = "Chaotic dynamic states with intermittent oscillations were generated in a Bray-Liebhafsky (BL) oscillatory reaction in an isothermal open reactor i.e., in the continuously-fed well-stirred tank reactor (CSTR) when the inflow concentration of potassium iodate was the control parameter. They are found between periodic oscillations obtained when [KIO3](0)  LT  3.00 x 10(-2) M and stable steady states when [KIO3](0) > 4.10 x 10(-2) M. It was shown that the most chaotic states obtained experimentally somewhere in the middle of this region are in high correlation with results obtained by means of largest Lyapunov exponents and phenomenological analysis based on the quantitative characteristics of intermittent oscillations.",
publisher = "Maik Nauka/Interperiodica/Springer, New York",
journal = "Russian Journal of Physical Chemistry A",
title = "Intermittent Chaos in the Bray-Liebhafsky Oscillator. Dependence of Dynamic States on the Iodate Concentration",
volume = "91",
number = "13",
pages = "2525-2529",
doi = "10.1134/S0036024417130076"
}

Bubanja, I. N., Ivanović-Šašić, A., Čupić, Ž., Anić, S.,& Kolar-Anić, L.. (2017). Intermittent Chaos in the Bray-Liebhafsky Oscillator. Dependence of Dynamic States on the Iodate Concentration. in Russian Journal of Physical Chemistry A
Maik Nauka/Interperiodica/Springer, New York., 91(13), 2525-2529.
https://doi.org/10.1134/S0036024417130076

Bubanja IN, Ivanović-Šašić A, Čupić Ž, Anić S, Kolar-Anić L. Intermittent Chaos in the Bray-Liebhafsky Oscillator. Dependence of Dynamic States on the Iodate Concentration. in Russian Journal of Physical Chemistry A. 2017;91(13):2525-2529.
doi:10.1134/S0036024417130076 .

Bubanja, Itana Nuša, Ivanović-Šašić, Ana, Čupić, Željko, Anić, Slobodan, Kolar-Anić, Ljiljana, "Intermittent Chaos in the Bray-Liebhafsky Oscillator. Dependence of Dynamic States on the Iodate Concentration" in Russian Journal of Physical Chemistry A, 91, no. 13 (2017):2525-2529,
https://doi.org/10.1134/S0036024417130076 . .

TY  - JOUR
AU  - Čupić, Željko
AU  - Stanojevic, Ana
AU  - Marković, Vladimir M.
AU  - Kolar-Anić, Ljiljana
AU  - Terenius, Lars
AU  - Vukojevic, Vladana
PY  - 2017
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/2087
AB  - Stress and alcohol use are interrelatedstress contributes to the initiation and upholding of alcohol use and alcohol use alters the way we perceive and respond to stress. Intricate mechanisms through which ethanol alters the organism's response to stress remain elusive. We have developed a stoichiometric network model to succinctly describe neurochemical transformations underlying the stress response axis and use numerical simulations to model ethanol effects on complex daily changes of blood levels of cholesterol, 6 peptide and 8 steroid hormones. Modelling suggests that ethanol alters the dynamical regulation of hypothalamic-pituitary-adrenal (HPA) axis activity by affecting the amplitude of ultradian oscillations of HPA axis hormones, which defines the threshold with respect to which the response to stress is being set. These effects are complexlow/moderate acute ethanol challenge ( LT  8mM) may reduce, leave unaltered or increase the amplitude of ultradian cortisol (CORT) oscillations, giving rise to an intricate response at the organism level, offering also a potential explanation as to why apparently discordant results were observed in experimental studies. In contrast, high-dose acute ethanol challenge (> 8mM) increases instantaneous CORT levels and the amplitude of ultradian CORT oscillations in a dose-dependent manner, affecting the HPA axis activity also during the following day(s). Chronic exposure to ethanol qualitatively changes the HPA axis dynamics, whereas ethanol at intoxicating levels shuts down this dynamic regulation mechanism. Mathematical modelling gives a quantitative biology-based framework that can be used for predicting how the integral HPA axis response is perturbed by alcohol.
PB  - Wiley, Hoboken
T2  - Addiction Biology
T1  - The HPA axis and ethanol: a synthesis of mathematical modelling and experimental observations
VL  - 22
IS  - 6
SP  - 1486
EP  - 1500
DO  - 10.1111/adb.12409
ER  - 

@article{
author = "Čupić, Željko and Stanojevic, Ana and Marković, Vladimir M. and Kolar-Anić, Ljiljana and Terenius, Lars and Vukojevic, Vladana",
year = "2017",
abstract = "Stress and alcohol use are interrelatedstress contributes to the initiation and upholding of alcohol use and alcohol use alters the way we perceive and respond to stress. Intricate mechanisms through which ethanol alters the organism's response to stress remain elusive. We have developed a stoichiometric network model to succinctly describe neurochemical transformations underlying the stress response axis and use numerical simulations to model ethanol effects on complex daily changes of blood levels of cholesterol, 6 peptide and 8 steroid hormones. Modelling suggests that ethanol alters the dynamical regulation of hypothalamic-pituitary-adrenal (HPA) axis activity by affecting the amplitude of ultradian oscillations of HPA axis hormones, which defines the threshold with respect to which the response to stress is being set. These effects are complexlow/moderate acute ethanol challenge ( LT  8mM) may reduce, leave unaltered or increase the amplitude of ultradian cortisol (CORT) oscillations, giving rise to an intricate response at the organism level, offering also a potential explanation as to why apparently discordant results were observed in experimental studies. In contrast, high-dose acute ethanol challenge (> 8mM) increases instantaneous CORT levels and the amplitude of ultradian CORT oscillations in a dose-dependent manner, affecting the HPA axis activity also during the following day(s). Chronic exposure to ethanol qualitatively changes the HPA axis dynamics, whereas ethanol at intoxicating levels shuts down this dynamic regulation mechanism. Mathematical modelling gives a quantitative biology-based framework that can be used for predicting how the integral HPA axis response is perturbed by alcohol.",
publisher = "Wiley, Hoboken",
journal = "Addiction Biology",
title = "The HPA axis and ethanol: a synthesis of mathematical modelling and experimental observations",
volume = "22",
number = "6",
pages = "1486-1500",
doi = "10.1111/adb.12409"
}

Čupić, Ž., Stanojevic, A., Marković, V. M., Kolar-Anić, L., Terenius, L.,& Vukojevic, V.. (2017). The HPA axis and ethanol: a synthesis of mathematical modelling and experimental observations. in Addiction Biology
Wiley, Hoboken., 22(6), 1486-1500.
https://doi.org/10.1111/adb.12409

Čupić Ž, Stanojevic A, Marković VM, Kolar-Anić L, Terenius L, Vukojevic V. The HPA axis and ethanol: a synthesis of mathematical modelling and experimental observations. in Addiction Biology. 2017;22(6):1486-1500.
doi:10.1111/adb.12409 .

Čupić, Željko, Stanojevic, Ana, Marković, Vladimir M., Kolar-Anić, Ljiljana, Terenius, Lars, Vukojevic, Vladana, "The HPA axis and ethanol: a synthesis of mathematical modelling and experimental observations" in Addiction Biology, 22, no. 6 (2017):1486-1500,
https://doi.org/10.1111/adb.12409 . .

TY  - JOUR
AU  - Stanojevic, A.
AU  - Markovic, V. M.
AU  - Čupić, Željko
AU  - Vukojevic, V.
AU  - Kolar-Anić, Ljiljana
PY  - 2017
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/2254
AB  - A model was developed that can be used to study the effect of gradual cholesterol intake by food on the HPA axis dynamics. Namely, well defined oscillatory dynamics of vital neuroendocrine hypothalamic-pituitary-adrenal (HPA) axis has proven to be necessary for maintaining regular basal physiology and formulating appropriate stress response to various types of perturbations. Cholesterol, as a precursor of all steroid HPA axis hormones, can alter the dynamics of HPA axis. To analyse its particular influence on the HPA axis dynamics we used stoichiometric model of HPA axis activity, and simulate cholesterol perturbations in the form of finite duration pulses, with asymmetrically distributed concentration profile. Our numerical simulations showed that there is a complex, nonlinear dependence between the HPA axis responsiveness and different forms of applied cholesterol concentration pulses, indicating the significance of kinetic modelling, and dynamical systems theory for the understanding of large-scale self-regulatory, and homeostatic processes within this neuroendocrine system.
PB  - Maik Nauka/Interperiodica/Springer, New York
T2  - Russian Journal of Physical Chemistry A
T1  - Modelling of the Hypothalamic-Pituitary-Adrenal Axis Perturbations by Externally Induced Cholesterol Pulses of Finite Duration and with Asymmetrically Distributed Concentration Profile
VL  - 91
IS  - 13
SP  - 2600
EP  - 2607
DO  - 10.1134/S0036024417130027
ER  - 

@article{
author = "Stanojevic, A. and Markovic, V. M. and Čupić, Željko and Vukojevic, V. and Kolar-Anić, Ljiljana",
year = "2017",
abstract = "A model was developed that can be used to study the effect of gradual cholesterol intake by food on the HPA axis dynamics. Namely, well defined oscillatory dynamics of vital neuroendocrine hypothalamic-pituitary-adrenal (HPA) axis has proven to be necessary for maintaining regular basal physiology and formulating appropriate stress response to various types of perturbations. Cholesterol, as a precursor of all steroid HPA axis hormones, can alter the dynamics of HPA axis. To analyse its particular influence on the HPA axis dynamics we used stoichiometric model of HPA axis activity, and simulate cholesterol perturbations in the form of finite duration pulses, with asymmetrically distributed concentration profile. Our numerical simulations showed that there is a complex, nonlinear dependence between the HPA axis responsiveness and different forms of applied cholesterol concentration pulses, indicating the significance of kinetic modelling, and dynamical systems theory for the understanding of large-scale self-regulatory, and homeostatic processes within this neuroendocrine system.",
publisher = "Maik Nauka/Interperiodica/Springer, New York",
journal = "Russian Journal of Physical Chemistry A",
title = "Modelling of the Hypothalamic-Pituitary-Adrenal Axis Perturbations by Externally Induced Cholesterol Pulses of Finite Duration and with Asymmetrically Distributed Concentration Profile",
volume = "91",
number = "13",
pages = "2600-2607",
doi = "10.1134/S0036024417130027"
}

Stanojevic, A., Markovic, V. M., Čupić, Ž., Vukojevic, V.,& Kolar-Anić, L.. (2017). Modelling of the Hypothalamic-Pituitary-Adrenal Axis Perturbations by Externally Induced Cholesterol Pulses of Finite Duration and with Asymmetrically Distributed Concentration Profile. in Russian Journal of Physical Chemistry A
Maik Nauka/Interperiodica/Springer, New York., 91(13), 2600-2607.
https://doi.org/10.1134/S0036024417130027

Stanojevic A, Markovic VM, Čupić Ž, Vukojevic V, Kolar-Anić L. Modelling of the Hypothalamic-Pituitary-Adrenal Axis Perturbations by Externally Induced Cholesterol Pulses of Finite Duration and with Asymmetrically Distributed Concentration Profile. in Russian Journal of Physical Chemistry A. 2017;91(13):2600-2607.
doi:10.1134/S0036024417130027 .

Stanojevic, A., Markovic, V. M., Čupić, Željko, Vukojevic, V., Kolar-Anić, Ljiljana, "Modelling of the Hypothalamic-Pituitary-Adrenal Axis Perturbations by Externally Induced Cholesterol Pulses of Finite Duration and with Asymmetrically Distributed Concentration Profile" in Russian Journal of Physical Chemistry A, 91, no. 13 (2017):2600-2607,
https://doi.org/10.1134/S0036024417130027 . .

TY  - JOUR
AU  - Abulseoud, Osama A.
AU  - Ho, Man Choi
AU  - Choi, Doo-Sup
AU  - Stanojevic, Ana
AU  - Čupić, Željko
AU  - Kolar-Anić, Ljiljana
AU  - Vukojevic, Vladana
PY  - 2017
UR  - https://cer.ihtm.bg.ac.rs/handle/123456789/2269
AB  - Changes in the hypothalamic-pituitary-adrenal (HPA) axis activity constitute a key component of bipolar mania, but the extent and nature of these alterations are not fully understood. We use here the lateral hypothalamic-kindled (LHK) rat model to deliberately induce an acute manic-like episode and measure serum corticosterone concentrations to assess changes in HPA axis activity. A mathematical model is developed to succinctly describe the entwined biochemical transformations that underlay the HPA axis and emulate by numerical simulations the considerable increase in serum corticosterone concentration induced by LHK. Synergistic combination of the LHK rat model and dynamical systems theory allows us to quantitatively characterize changes in HPA axis activity under controlled induction of acute manic-like states and provides a framework to study in silico how the dynamic integration of neurochemical transformations underlying the HPA axis is disrupted in these states.
PB  - Public Library of Science (PLoS)
T2  - Plos One
T1  - Corticosterone oscillations during mania induction in the lateral hypothalamic kindled rat-Experimental observations and mathematical modeling
VL  - 12
IS  - 5
DO  - 10.1371/journal.pone.0177551
ER  - 

@article{
author = "Abulseoud, Osama A. and Ho, Man Choi and Choi, Doo-Sup and Stanojevic, Ana and Čupić, Željko and Kolar-Anić, Ljiljana and Vukojevic, Vladana",
year = "2017",
abstract = "Changes in the hypothalamic-pituitary-adrenal (HPA) axis activity constitute a key component of bipolar mania, but the extent and nature of these alterations are not fully understood. We use here the lateral hypothalamic-kindled (LHK) rat model to deliberately induce an acute manic-like episode and measure serum corticosterone concentrations to assess changes in HPA axis activity. A mathematical model is developed to succinctly describe the entwined biochemical transformations that underlay the HPA axis and emulate by numerical simulations the considerable increase in serum corticosterone concentration induced by LHK. Synergistic combination of the LHK rat model and dynamical systems theory allows us to quantitatively characterize changes in HPA axis activity under controlled induction of acute manic-like states and provides a framework to study in silico how the dynamic integration of neurochemical transformations underlying the HPA axis is disrupted in these states.",
publisher = "Public Library of Science (PLoS)",
journal = "Plos One",
title = "Corticosterone oscillations during mania induction in the lateral hypothalamic kindled rat-Experimental observations and mathematical modeling",
volume = "12",
number = "5",
doi = "10.1371/journal.pone.0177551"
}

Abulseoud, O. A., Ho, M. C., Choi, D., Stanojevic, A., Čupić, Ž., Kolar-Anić, L.,& Vukojevic, V.. (2017). Corticosterone oscillations during mania induction in the lateral hypothalamic kindled rat-Experimental observations and mathematical modeling. in Plos One
Public Library of Science (PLoS)., 12(5).
https://doi.org/10.1371/journal.pone.0177551

Abulseoud OA, Ho MC, Choi D, Stanojevic A, Čupić Ž, Kolar-Anić L, Vukojevic V. Corticosterone oscillations during mania induction in the lateral hypothalamic kindled rat-Experimental observations and mathematical modeling. in Plos One. 2017;12(5).
doi:10.1371/journal.pone.0177551 .

Abulseoud, Osama A., Ho, Man Choi, Choi, Doo-Sup, Stanojevic, Ana, Čupić, Željko, Kolar-Anić, Ljiljana, Vukojevic, Vladana, "Corticosterone oscillations during mania induction in the lateral hypothalamic kindled rat-Experimental observations and mathematical modeling" in Plos One, 12, no. 5 (2017),
https://doi.org/10.1371/journal.pone.0177551 . .