Stevan Maćešić

Modelling cholesterol effects on the dynamics of the hypothalamic-pituitary-adrenal (HPA) axis.

A mathematical model of the hypothalamic-pituitary-adrenal (HPA) axis with cholesterol as a dynamical variable was derived to investigate the effects of cholesterol, the primary precursor of all steroid hormones, on the ultradian and circadian HPA axis activity. To develop the model, the parameter space was systematically examined by stoichiometric network analysis to identify conditions for ultradian oscillations, determine conditions under which dynamic transitions, i.e. bifurcations occur and identify bifurcation types. The bifurcations were further characterized using numerical simulations. Model predictions agree well with empirical findings reported in the literature, indicating that cholesterol levels may critically affect the global dynamics of the HPA axis. The proposed model provides a base for better understanding of experimental observations, it may be used as a tool for designing experiments and offers useful insights into the characteristics of basic dynamic regulatory mechanisms that, when impaired, may lead to the development of some modern-lifestyle-associated diseases.

Dynamic transitions in a model of the hypothalamic-pituitary-adrenal axis

Dynamic properties of a nonlinear five-dimensional stoichiometric model of the hypothalamic-pituitary-adrenal (HPA) axis were systematically investigated. Conditions under which qualitative transitions between dynamic states occur are determined by independently varying the rate constants of all reactions that constitute the model. Bifurcation types were further characterized using continuation algorithms and scale factor methods. Regions of bistability and transitions through supercritical Andronov-Hopf and saddle loop bifurcations were identified. Dynamic state analysis predicts that the HPA axis operates under basal (healthy) physiological conditions close to an Andronov-Hopf bifurcation. Dynamic properties of the stress-control axis have not been characterized experimentally, but modelling suggests that the proximity to a supercritical Andronov-Hopf bifurcation can give the HPA axis both, flexibility to respond to external stimuli and adjust to new conditions and stability, i.e., the capacity to return to the original dynamic state afterwards, which is essential for maintaining homeostasis. The analysis presented here reflects the properties of a low-dimensional model that succinctly describes neurochemical transformations underlying the HPA axis. However, the model accounts correctly for a number of experimentally observed properties of the stress-response axis. We therefore regard that the presented analysis is meaningful, showing how in silico investigations can be used to guide the experimentalists in understanding how the HPA axis activity changes under chronic disease and/or specific pharmacological manipulations.

Current rates and reaction rates in the Stoichiometric Network Analysis (SNA)

Abstract In stoichiometric network analysis (SNA) the instability condition is calculated by the current rates. Recently, we have shown that in the final result the current rates can be substituted by reaction rates, which is a more appropriate value for the examination of instability from experimental point of view. Here, we elaborate the problem of whether the current rates are necessary parameters in the calculation, with the aim of obtaining the region of instability. All calculations are performed on a model for Belousov-Zhabotinsky (BZ) reaction, which has not been examined by SNA.

Bifurcation analysis: a tool for determining model parameters of the considered process

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.

Stoichiometric network analysis of a reaction system with conservation constraints

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.