Next Article in Journal
From Green Remediation to Polymer Hybrid Fabrication with Improved Optical Band Gaps
Next Article in Special Issue
Mechanistic Insight of Na/K-ATPase Signaling and HO-1 into Models of Obesity and Nonalcoholic Steatohepatitis
Previous Article in Journal
Significantly Elevated FMR1 mRNA and Mosaicism for Methylated Premutation and Full Mutation Alleles in Two Brothers with Autism Features Referred for Fragile X Testing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Systems Approach to Study Associations between OxLDL and Abdominal Aortic Aneurysms

by
Łukasz Gutowski
1,
Kaja Gutowska
2,3,
Maria Pioruńska-Stolzmann
1,
Piotr Formanowicz
2,3 and
Dorota Formanowicz
1,*
1
Department of Clinical Biochemistry and Laboratory Medicine, Poznan University of Medical Sciences, Rokietnicka 8, 60-806 Poznan, Poland
2
Institute of Computing Science, Poznan University of Technology, Piotrowo 2, 60-965 Poznan, Poland
3
Institute of Bioorganic Chemistry, Polish Academy of Sciences, Noskowskiego 12/14, 61-704 Poznan, Poland
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(16), 3909; https://doi.org/10.3390/ijms20163909
Submission received: 4 June 2019 / Revised: 4 August 2019 / Accepted: 8 August 2019 / Published: 11 August 2019
(This article belongs to the Special Issue Redox Regulation of Cardiovascular Metabolism and Inflammation)

Abstract

:
Although abdominal aortic aneurysm (AAA) is a common vascular disease and is associated with high mortality, the full pathogenesis of AAA remains unknown to researchers. Abdominal aortic aneurysms and atherosclerosis are strongly related. Currently, it is more often suggested that development of AAA is not a result of atherosclerosis, however, individual factors can act independently or synergistically with atherosclerosis. One of such factors is low-density lipoprotein (LDL) and its oxidized form (oxLDL). It is known that oxLDL plays an important role in the pathogenesis of atherosclerosis, thus, we decided to examine oxLDL impact on the development of AAA by creating two models using Petri-nets. The first, full model, contains subprocess of LDL oxidation and all subprocesses in which it participates, while the second, reduced model, does not contain them. The analysis of such models can be based on t-invariants. They correspond to subprocesses which do not change the state of the modeled system. Moreover, the knockout analysis has been used to estimate how crucial a selected transition (representing elementary subprocess) is, based on the number of excluded subprocesses as a result of its knockout. The results of the analysis of our models show that oxLDL affects 55.84% of subprocesses related to AAA development, but the analysis of the nets based on knockouts and simulation has shown that the influence of oxLDL on enlargement and rupture of AAA is negligible.

Graphical Abstract

1. Introduction

An abdominal aortic aneurysm can occur in 5% of men aged 65–74 [1] and is associated with over 80% mortality rate after its rupture [2]. In 2013 rupture of the aortic aneurysm caused more than 150,000 deaths, making this disease the thirteenth leading cause of death in the world [1]. An abdominal aortic aneurysm is defined as an enlargement exceeding 1.5 times the normal size in a given segment [3]. The full mechanism of pathogenesis of the abdominal aortic aneurysm is still unknown.
Atherosclerosis and abdominal aortic aneurysms are strongly interrelated, and therefore, they cause controversy as to the individuality of these diseases [4,5]. Currently, it is more often suggested that the development of abdominal aortic aneurysm is not a result of atherosclerosis, and there are also separate factors that can act independently or synergistically with atherosclerosis [2]. Both disorders, on the other hand, share some risk factors and some similar pathological characteristics [5,6]. There are, therefore, questions about the factors that can be shared in the development of both disorders. One of such factors is LDL and its oxidized form—oxLDL. It is known that oxLDL plays an important role in the pathogenesis of atherosclerosis [7]. Nonetheless, studies about its role in AAA are inconclusive. Some of these studies do not indicate a relationship [2,8,9], while others report LDL participation in the development of abdominal aortic aneurysms [10,11,12]. According to our best knowledge there is no animal model describing influence of oxLDL on AAA, but the impact of oxLDL on the development of other aneurysms has been proven [13]. Therefore, there is a possibility of synergistic development of both diseases, which we decided to examine by creating a model of oxLDL influence on the development of AAA using the Petri-nets.
It is not always easy to study certain aspects of aneurysm on human tissues, that is why animal AAA models are used. There are several types of animal AAA models, however three of them are mainly used: elastase perfusion, CaCl 2 application, and infusion of angiotensin II into either ApoE / . or LDL receptor / mice. In the elastase model abdominal aorta is perfused with this elastase for a short period of time. It leads to dilatation after about 14 days and development of a chronic inflammatory response in aortic wall. Another model uses CaCl 2 . Calcium chloride is applied to mice by the peri-aortic incubation. In this model mice develop aneurysm at the end of the third week, and it is characterized by a chronic inflammatory response. In the third model, subcutaneous infusion of AngII into either LDL receptor / or ApoE / mice leads to the formation of AAAs in the suprarenal region within the 28-day infusion period. This method is the most widely used because it does not require abdominal surgery. Depending on the studies mice were fed with normal or high fat diet [6,11,12].
The model proposed by us describes interdependencies between many processes, such as the formation of oxLDL, production of reactive oxygen species by nitric oxide synthase (NOS) and NADPH oxidase enzymes (NOX) and their effect on the development of inflammation, increase in the amount of the metalloproteinases, and as a result - development and rupture of AAA.
The most important reactive form of oxygen in aneurysm development is superoxide anion radical (O 2 ) [14,15]. It is synthesized by uncoupled inducible NOS (iNOS) and endothelial NOS (eNOS), and NOX [15]. NOS with cofactor tetrahydrobiopterin (BH 4 )produces nitric oxide (NO) in physiological conditions. However, the oxidation of BH 4 to BH 2 under oxidative stress conditions leads to uncoupling of the enzyme and production of O 2 . This radical intensifies the uncoupling of the enzyme, although it comes to it much faster under the influence of peroxynitrite (ONOO ), which is produced as a result of the reaction of O 2 and NO [16,17]. OxLDL also affects NOS, leading to its dephosphorylation, [18,19,20], and consequently to the production of superoxide anion radical instead of NO. Superoxide anion radical is reduced to hydrogen peroxide (H 2 O 2 ) by enzyme dismutase (SOD) [21]. H 2 O 2 is decomposed to water and oxygen by the enzyme catalase. O 2 , H 2 O 2 and ONOO are reactive oxygen and nitrogen species (ROS/RNS) that have a significant impact on the development of oxidative stress, which plays a significant role in the pathogenesis of AAA. Polymorphonuclear leukocytes (PMNs) release myeloperoxidase (MPO) which produces another ROS—hypochlorous acid (HOCl) from H 2 O 2 and Cl . HOCl indirectly stimulates proteolytic activity. MPO products can lead to lipid peroxidation [22]. ROS regulate Akt kinase activity involved in the process of vascular smooth muscle cells (VSMC) apoptosis [23,24,25], and activate activator protein-1 (AP-1) through MAPK, which leads to the production of chemokines and proinflammatory cytokines [24,26]. The most important cytokines in this disease entity are IL-6, IL-1 β , IL-8 and TNF α . Production of some of them, like IL-6, could be increased, amongst others, by aortic wall stretch [27]. ROS damage DNA, which leads to in nuclear factor κ B (NF κ B) activity. NF κ B also leads to the production of pro-inflammatory chemokines and cytokines as well as adhesion molecules: intercellular adhesion molecule-1 (ICAM-1) and vascular cell adhesion molecule-1 (VCAM-1) [24]. Oxidative stress increases the expression of an angiotensin-converting enzyme (ACE) [23]. It participates in the conversion process of angiotensin I in II. Angiotensin II activates NF κ B and NOX by activating its subunits [24]. NOX is also stimulated by tumor necrosis factor α (TNF α ), oxLDL and hemodynamic stress [28,29]. TNF α also activates NF κ B [30]. ROS can oxidize LDL to oxLDL. OxLDL activates superoxide anion generation by NOX by activating its subunits [7,31]. OxLDL also affects on increase of O 2 concentration in other ways: it leads to dephosphorylation of NOS (as a result of which O 2 is produced, instead of NO) [18,19,20], and inhibits dismutase activity (decreases O 2 reduction to H 2 O 2 ) [7]. OxLDL increases activity of metalloproteinase 2 (MMP2) [31] and the production of adhesion particles ICAM and VCAM [20].
Increase in adhesive particles and chemokines amount foster infiltration of the aortic wall by inflammatory cells. These cells intensify the synthesis of proinflammatory cytokines, which deepens inflammation and leads to the synthesis of proteases, mostly metalloproteinases [32,33,34]. Reactive oxygen species also take part in the expression of genes and activation of latent MMPs [10]. The essential metalloproteinases in the pathogenesis of AAA are 1, 2, 7, 9, 12, 13. They degrade connective tissues. The proteolysis products are chemoattractants, intensifying infiltration by inflammatory cells [35]. MMPs 2, 7, 9, and 12 are mainly responsible for the degradation of elastin which results in aortic dilatation [35]. This process is intensified by previously mentioned VSMC apoptosis because these cells are responsible for the synthesis of the elastin [17,35]. Enlargement of the aneurysm increases the hemodynamic stress, which leads to compensatory deposition of collagen in the aortic wall. This makes it possible to maintain resistance to the pressure exerted on the vessel wall [35]. In turn, MMPs 1, 8, and 13 are mainly responsible for the degradation of collagen. If collagen predominate in the aortic wall, the long-term effects of metalloproteinases 1, 8 and 13 may result in rupture of the aneurysm.
The description above shows very complex biological phenomenon, which should also be treated as a complex biological system. To understand this phenomenon better, systematize knowledge, discover new properties, and confirm specific facts, the systems approach has been used. In this study, a systems approach based on Petri-net, and its analysis based on t-invariants has been applied.

2. Analysis

The presented analysis of the proposed Petri-net models is based on t-invariants. Two main parts of such an analysis can be distinguished, i.e., an analysis of Maximum Common Transition sets (MCT sets) and an analysis of t-clusters. In addition to classical analysis, significance and knockout analyses were performed. The former one determines percentage of all subprocesses in which selected elementary subprocesses are involved. Such an analysis allows to distinguish which subprocesses are more crucial for functioning of the modeled system. As a complement, the knockout analysis allows to estimate how crucial selected transition is, based on the number of excluded subprocesses. The knockout simulation collects more detailed data about a behavior of the net when selected transitions are excluded.
The analysis of MCT sets is associated with determination of certain functional blocks. For the presented full model, 15 MCT sets were described in Table A3 in Appendix B. A need of the analysis of t-clusters follows from a large number of t-invariants. They were grouped in such clusters and a biological meaning of each of them was determined. For the proposed Petri-net 17 t-clusters described in Table A4 in Appendix C were calculated. In this case, analysis of t-clusters does not allow for accurate analysis of all subprocesses. Therefore, a significance analysis was performed for selected transitions (elementary subprocess). This means that for each transition, its attendance frequency is determined in all supports of t-invariants.
To complement the analysis of selected t-invariants a knockout analysis was performed (cf. [36,37]). This analysis relies on disabling selected transition (elementary subprocess/selected reaction). Turning off certain active component of a Petri-net consequently leads to exclusion of others. Knockouted transition (certain reaction) may disable large fragments of the net, when it plays a key role in the functioning of the entire model. After turning off the selected transition or several transitions, t-invariants are calculated again. On this basis it is possible to estimate which subprocesses have been excluded in consequence of knockout of selected transitions. Thus, the knockout analysis allows to estimate how crucial selected transition (elementary subprocess) is, based on the number of excluded subprocesses.
These analyses focuses on subprocesses having an indirect and direct impact on the development and rupture of the abdominal aortic aneurysm. Distinguished subprocesses are: LDL oxidation, oxidative stress, production of inflammatory cytokines, influence of MMPs, and impact of NOX. The results were presented for both full model and reduced model, to assess the impact of oxLDL on the development of AAA on this basis. The results of the significance analysis of selected subprocesses for full model and reduced model are presented in Table 1. Selected subprocesses may often consist of more than one elementary subprocess (transition). On this basis, it is possible to determine the percentage of all subprocesses in which a specific elementary subprocess is involved. From a mathematical point of view it is determined in how many t-invariants a selected transition occurs. In consequence, conclusions which subprocesses are more crucial for functioning of the modeled system than others can be drawn. In Table 1, it can be seen that subprocess of oxLDL appears in 55.84% of all modeled subprocesses in full model. This result could suggest that oxLDL has high significance for whole system. Nonetheless, knockout analysis was also conducted to achieve full analysis.
Table 1 includes information about percentage contribution of particular elementary subprocess (single transition) in the whole system, while Table 2 contains percentage contribution of particular subprocess. To be precise, Table 2 contains results of knockout analysis of selected subprocesses for two models: full and reduced. As mentioned before, some subprocesses consist of several elementary subprocesses (several transitions), and all of them are turned off in this analysis. Knockout of such set of transitions allows to determine contribution of certain subprocess in the whole system. Excluded sets of transitions are indicated in the row called “knockouted transitions”. It should be noted that the reduced model does not include certain transitions, inter alia, t 0 which is associated with LDL oxidation, hence no result in this field of Table 2. Thus, comparing full model (with oxLDL) and reduced model (without oxLDL) allows to estimate the percentage contribution of oxLDL in the selected subprocesses, i.e., oxidative stress, influence of MMPs, influence of NOX, production of inflammatory cytokines, enlargement and rupture of AAA.
Table 3 summarizes the most important results of the analysis of the proposed full model. For subprocesses associated with enlargement and rupture of aneurysm, the significance of a given subprocess in the full network (736 t-invariants) and in the knockouted network (325 t-invariants) is presented. In the last column, the difference between percentage contribution in the whole system for full model, and full model with knockout was calculated. On this basis influence of oxLDL on aneurysm enlargement (1.15 p.p.) and rupture (2.71 p.p.) has been determined. It can be noticed that the knockout analysis for the full model (excluded transition t 0 is associated with LDL oxidation) gives the same results as the significance analysis (Table 1) for the reduced Petri-net model. As can be seen, results suggest that these analyzes are consistent.
To complete the above analyses and confirm the results that the influence of oxLDL on aneurysm enlargement is only 1.15 p.p. and aneurysm rupture is 2.71 p.p., a knockout simulation was performed. Its purpose was collection of more detailed data about behavior of Petri-net in situation of excluding various transitions and comparing them with each other. This simulation can be performed in the Java application called Holmes [38].
The simulation was carried out for full model and for full model with knockout (knockouted transition t 0 —LDL oxidation). Simulation properties are as follows: 10,000 steps, 1000 repetitions. On this basis, we obtained average number of firing of transition (AvgF) in all steps in all simulations for these two models. Our results focus on transitions t 52 (enlargement of AAA) and t 54 (rupture of AAA):
  • AvgF for t 52 in full model is 19.40, while AvgF for t 52 in full model with knockout of t 0 is 19.46.
  • AvgF for t 54 in full model is 4.85, while AvgF for t 54 in full model with knockout of t 0 is 4.86.
As it can be noticed, differences in the average number of firings of transitions t 52 and t 54 before and after knockout are insignificant. This confirms the observations summarized in Table 3. All key subprocesses have not changed after the knockout, which means that the subprocesses of enlargement and rupture of AAA can occur without influence of oxLDL.

3. Petri-Net-Based Models

Two Petri-net-based models of aneurysm formation and rupture have been proposed (they are extended versions of the net presented in [39]). These models have been created using Holmes. The first, full model, contains subprocess of LDL oxidation to oxLDL and all subprocesses in which they participate, i.e., stimulation of NOX, stimulation of adhesion particles, inhibition of SOD, activation of elastin proteolysis by MMP2, and synthesis of superoxide anion radical together with NOS. While the second, reduced model, does not contain these subprocesses. These two models of Petri-nets are compared to each other. The main purpose of such comparison is evaluation how large impact on the enlargement and rupture of AAA have subprocesses related to oxLDL. A Petri-net model, despite its intuitive graphical representation, is complicated due to a large number of passive and active components. Therefore, a schematic diagram is shown in Figure 1 [3,7,10,14,15,16,17,18,19,20,21,22,23,24,25,27,28,29,30,31,33,35,40,41,42,43,44,45,46]. In this scheme of AAA development green color indicates reactive oxygen/nitrogen species, NOX and NOSs are marked blue, participation of angiotensin is marked orange, processes associated with oxLDL are marked yellow, factors affecting enlargement and rupture of AAA are marked purple, and AAA enlargement and rupture are marked with pink color.
The proposed full model of Petri-net is shown in Figure A1 in Appendix A. In that figure oxLDL subprocesses, marked with red color, are included in the full model, and they were removed from reduced model.
Petri-net-based full model of AAA development contains 64 transitions and 40 places. Places correspond to passive components of the modeled system, while transitions correspond to its active components. Descriptions of all places and transitions of the net are included in Appendix A in Table A1 and Table A2, respectively. Full model is covered by 736 t-invariants (which correspond to biological subprocesses. Reduced model, which does not contains subprocesses related to oxLDL, contains 57 transitions and 38 places. Transitions and places that have been removed from this model are marked with red in Figure A1 and their names are marked in bold font in Table A1 and Table A2 in Appendix A. Reduced model is covered by 325 t-invariants.

4. Methods

Petri-nets are mathematical objects suitable for modeling a wide class of systems, especially those ones which contain concurrent processes. They were proposed in 1962 by Carl A. Petri in the context of computer science [47]. Such nets have been used for years in modeling and analysis of technical systems. In the mid of 1990s it has been realised that nets of this type can be used also for investigations of properties of biological systems (cf. [48]).
Petri-nets have a structure of a directed bipartite graph what means that they are composed of two disjoint subsets of vertices, called places and transitions. These vertices are connected by arcs in such a way that an arc connects a place with a transition or a transition with a place (i.e., no two places nor two transitions are connected). When a Petri-net is a model of a biological system places correspond to its passive components, as chemical compounds, while transitions are counterparts of active components, as chemical reactions. Arcs describe causal relations between the passive and active components and they are labeled by positive integer numbers called weights [49,50].
There is one more type of components of Petri-nets, i.e., tokens. They bring into the net dynamics (which is crucial for modeling systems) not only the biological ones. Tokens flow from one place to another through transitions. This flow corresponds to a flow of substances, information etc. through the modeled system. It is governed by a simple rule of transition firing. According to this rule transition t j is active if in every place p i directly preceding it (such a place is called a pre-place of transition t j ) the number of tokens is equal to at least the weight of arc ( p i , t j ) , i.e., the arc connecting p i with t j . An active transition can be fired, what means that tokens flow from its pre-places to its post-places, i.e, the places directly succeeding t j , and the numbers of flowing tokens are equal to weights of the appropriate arcs [49,50].
Petri-nets have an intuitive graphical representation. In this representation places are depicted as circles, transitions as rectangles, arcs as arrows and tokens as dots or numbers residing in places. When a weight of an arc is equal to one, it is not shown in the graphical representation of the net. While this representation is very helpful in understanding the structure of the modeled system and supports simulations of its behavior, it is not very well suited for a formal analysis of its properties. For this purpose another representation, called an incidence matrix, can be used. In such a matrix A = ( a i j ) n × m , where n is a number of places and m is a number of transitions, entry a i j is equal to a difference between numbers of tokens in place p i before and after firing transition t j [48].
In the analysis of Petri-net-based models of biological systems especially important are t-invariants, which are vectors x Z m being solutions to equation A · x = 0 . To every t-invariant there corresponds set of transitions s ( x ) = { t j : x j > 0 } , called its support. When a Petri-nets is a model of a biological system, usually it should be covered by t-invariants, what means that each transition should be an element of a support of some such an invariant. t-invariants are counterparts of subprocesses which do not change a state of the modeled system. More precisely, if every transition t j s ( x ) is fired x j times, a distribution of tokens in places (called a marking of the net) does not change t j [48].
From this follows that an analysis of similarities among t-invariants may lead to discoveries of previously unknown properties of the modeled system. Indeed, such similarities may be a reason of interactions between the above mentioned subprocesses, since supports of similar t-invariants have a non-empty intersection. This intersection contains some transitions corresponding to elementary processes being components of the subprocesses. Since these subprocesses contain some common elementary processes they can interact with each other through them. Hence, looking for similar t-invariants may lead to identifying unknown interactions of subprocesses which may be a source of some important properties of the analyzed biological system [51,52,53].
Searching for similarities among t-invariants can be done using standard clustering algorithms. They divide the set of all t-invariants into a disjoint subset containing t-invariants which are similar to each other according to some similarity measure. However, it is not a trivial task since a proper algorithm, a proper similarity measure and a proper number of clusters should be chosen. All of these three components of the clustering procedure should be adjusted to the modeled system and to the problem to be solved on the basis of the model. The resulting clusters are called t-clusters and correspond to some functional modules of the biological system.
Moreover, also the set of transitions can be divided into disjoint subsets called MCT sets. Such a set contains transitions being elements of supports of the same t-invariants and corresponds to some functional block of the modeled biological system. During the analysis biological meaning of MCT sets as well as t-clusters should be determined (cf. [53,54]).
Figure 2 presents a general scheme of work, where three stages can be identified. The first—creation of model, the second—analyses of model, and the third—result obtainment.

5. Conclusions

There is no agreement in the scientific community about the influence of LDL and its oxidized form on the formation of abdominal aortic aneurysms. According to popular belief, it might seem that oxLDL should be involved, not only in the pathogenesis of atherosclerosis, but also in the pathogenesis of abdominal aortic aneurysm, however knockout analysis and simulation of the proposed Petri-net-based models with and without oxLDL have shown that there is almost no influence of oxLDL on aneurysm enlargement and rupture. Nevertheless, it should be noted that in our model oxLDL is associated with 55.84% of all modeled subprocesses involved in the pathogenesis of aneurysm. Despite the involvement of oxLDL in more than half of the subprocesses, its participation does not significantly translate into enlargement or rupture of aneurysm. Our results seem to reconcile some of the contradictory reports: oxLDL is indeed involved in the formation of aneurysms, but it is not reflected significantly in its enlargement or rupture, and the whole process can take place without its participation.

Author Contributions

Conceptualization—Ł.G., K.G., M.P.-S., P.F. and D.F.; methodology—Ł.G., K.G., P.F.; validation—Ł.G., K.G, P.F. and D.F.; formal analysis—Ł.G., K.G., and P.F.; investigation—Ł.G., K.G., P.F. and D.F.; writing—original draft preparation—Ł.G., K.G, M.P.-S., P.F. and D.F.; writing—review and editing—Ł.G., K.G, P.F. and D.F.; supervision—P.F. and D.F; funding acquisition—P.F. and D.F.

Funding

This research has been partially supported by the National Science Centre (Poland) Grant No. 2012/07/B/ST6/01537 and by the statutory funds of Poznan University of Medical Sciences No. 502-01-22283690-08800.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. The Proposed Petri Net Model

Table A1. List of places for full model.
Table A1. List of places for full model.
PlaceBiological MeaningPlaceBiological Meaning
p 0 superoxide anion radical p 20 PKC
p 1 oxidized LDL p 21 NF κ B
p 2 LDL p 22 chemokines
p 3 NO p 23 inflammatory cytokines
p 4 peroxynitrite p 24 MMPs
p 5 BH 4 p 25 adhesion particles
p 6 NOSs p 26 circulating inflammatory cells
p 7 BH 2 p 27 TNF α
p 8 H 2 O 2 p 28 p38MAPK
p 9 ROS p 29 active inflammatory cells
p 10 dismutase p 30 hemodynamic stress
p 11 catalase p 31 COX2
p 12 NOX p 32 PGE2
p 13 ANGII p 33 less elastin
p 14 ANGI p 34 less collagen
p 15 ACE p 35 enlarged AAA
p 16 tyrosine kinase p 36 deposited collagen
p 17 PARP polymerase p 37 HOCl
p 18 AP1 p 38 MPO
p 19 Akt p 39 IL-6
Table A2. List of transitions for full model.
Table A2. List of transitions for full model.
TransitionBiological MeaningTransitionBiological Meaning
t 0 LDL oxidation t 32 expression of MMPs by inflammatory cytokines
t 1 blood as LDL source t 33 activation of MMPs by ROS
t 2 peroxynitrite synthesis t 34 production of adhesion particles
t 3 BH 4 synthesis t 35 diapedesis
t 4 NO synthesis t 36 source of inflammatory cells
t 5 NOSs sources t 37 source of ANGI
t 6 BH 4 oxidation to BH 2 through peroxynitrite t 38 elastine proteolysis by MMPs
t 7 BH 4 oxidation to BH 2 through O 2 t 39 source of TNF α
t 8 O 2 synthesis via NOSs and BH 2 t 40 pool of inflamatory cytokines
t 9 pool of ROS t 41 activation of p38MAPK
t 10 dismutation t 42 production of inflammatory cytokines by inflmmatory cells
t 11 source of dismutase t 43 direct stimulation of NF κ B via ANGII
t 12 H 2 O 2 reduction t 44 O 2 synthesis via NOSs and oxLDL
t 13 source of catalase t 45 stimulation of adhesion particles
t 14 O 2 synthesis through NOX t 46 stimulation of NF κ B by TNF α
t 15 indirect stimulation of NOX through oxLDL t 47 activation of NOX by hemodynamic stress
t 16 activation of NOX through TNF α t 48 activation of COX2 by hemodynamic stress
t 17 activation of NOX by ANGII t 49 production of PGE2 by COX2
t 18 conversion of ANGI to ANGII t 50 activation of NF κ B via PGE2
t 19 source of ACE t 51 collagen proteolysis by MMPs
t 20 activation of ACE by ROS t 52 enlargement of AAA
t 21 stimulation of tyrosine kinase t 53 compensating collagen deposition
t 22 DNA damage and release of PARP t 54 rupture of AAA
t 23 stimulation of AP1 t 55 degradation of chemokines
t 24 indirect stimulation of Akt t 56 additional activation of elastin proteolysis by MMP2
t 25 VSMC apoptosis t 57 increase of hemodynamic stress
t 26 stimulation of PKC by PARP t 58 activation of iNOS by inflammatory cytokines
t 27 activation of NF κ B by PKC t 59 inhibition of SOD
t 28 stimulation of inflammatory cytokines production by AP1 t 60 production of HOCl by MPO
t 29 stimulation of inflammatory cytokines production by NF κ B t 61 source of MPO
t 30 stimulation of chemokines production by AP1 t 62 activation of MMPs by HOCl
t 31 stimulation of chemokines production by NF κ B t 63 stimulation of IL-6
Figure A1. The Petri-net-based model of formation and rupture of abdominal aortic aneurysm including subprocesses associated with oxLDL (full model, with oxLDL). The model was divided into several blocks: (a) production of HOCl, (b) synthesis of NO, (c) synthesis of O 2 by NOSs and BH 2 , (d) compensating collagen deposition by hemodynamic stress, (e) rupture of AAA, (f) oxidation of BH 4 to BH 2 , (g) synthesis of ONOO by O 2 and NO, (h) enlargement of AAA, (i) proteolysis of elastin and collagen by MMPs, (j) additional activation of MMP2 by oxLDL, (k) VSMC apoptosis, (l) expression of MMPs by inflammatory cytokines and ROS, (m) oxidation of LDL, (n) synthesis of O 2 by NOX, (o) dismutation of O 2 into H 2 O 2 , (p) synthesis of O 2 by NOSs and oxLDL, (q) activation of NF κ B caused by DNA damage by ROS, (r) stimulation of AP1 by ROS, (s) inhibition of SOD by oxLDL, (t) stimulation of NOX by oxLDL, (u) activation of ACE by ROS, (v) reduction of H 2 O 2 by catalase, (w) stimulation of NF κ B by TNF α , (x) activation of NOX by TNF α , (y) activation of NOX, COX2 and IL-6 by hemodynamic stress, (z) activation of NOX by ANGII (conversion of ANGI to ANGII by ACE), (aa) stimulation of NF κ B by ANGII, (ab) production of adhesion particles by NF κ B, (ac) stimulation of chemokines production by NF κ B, (ad) stimulation of inflammatory cytokines production by NF κ B, (ae) activation of NF κ B via PGE2 activated by COX2, (af) stimulation of adhesion particles by oxLDL, (ag) diapedesis, (ah) infiltration by inflammatory cells. In addition, transitions and places related to LDL oxidation and their roles in the modeled process were marked with red squares and circles. These specific elements were omitted in the second model (reduced model, without oxLDL). In addition, some places and transitions that are common for the highlighted blocks are marked with a non-colored octagon.
Figure A1. The Petri-net-based model of formation and rupture of abdominal aortic aneurysm including subprocesses associated with oxLDL (full model, with oxLDL). The model was divided into several blocks: (a) production of HOCl, (b) synthesis of NO, (c) synthesis of O 2 by NOSs and BH 2 , (d) compensating collagen deposition by hemodynamic stress, (e) rupture of AAA, (f) oxidation of BH 4 to BH 2 , (g) synthesis of ONOO by O 2 and NO, (h) enlargement of AAA, (i) proteolysis of elastin and collagen by MMPs, (j) additional activation of MMP2 by oxLDL, (k) VSMC apoptosis, (l) expression of MMPs by inflammatory cytokines and ROS, (m) oxidation of LDL, (n) synthesis of O 2 by NOX, (o) dismutation of O 2 into H 2 O 2 , (p) synthesis of O 2 by NOSs and oxLDL, (q) activation of NF κ B caused by DNA damage by ROS, (r) stimulation of AP1 by ROS, (s) inhibition of SOD by oxLDL, (t) stimulation of NOX by oxLDL, (u) activation of ACE by ROS, (v) reduction of H 2 O 2 by catalase, (w) stimulation of NF κ B by TNF α , (x) activation of NOX by TNF α , (y) activation of NOX, COX2 and IL-6 by hemodynamic stress, (z) activation of NOX by ANGII (conversion of ANGI to ANGII by ACE), (aa) stimulation of NF κ B by ANGII, (ab) production of adhesion particles by NF κ B, (ac) stimulation of chemokines production by NF κ B, (ad) stimulation of inflammatory cytokines production by NF κ B, (ae) activation of NF κ B via PGE2 activated by COX2, (af) stimulation of adhesion particles by oxLDL, (ag) diapedesis, (ah) infiltration by inflammatory cells. In addition, transitions and places related to LDL oxidation and their roles in the modeled process were marked with red squares and circles. These specific elements were omitted in the second model (reduced model, without oxLDL). In addition, some places and transitions that are common for the highlighted blocks are marked with a non-colored octagon.
Ijms 20 03909 g0a1

Appendix B. Analysis of MCT Sets

Analysis of models based on Petri-net focuses mainly on t-invariants. On the basis of t-invariants, MCT sets are calculated. These sets divide Petri-net into blocks containing certain transitions. To be precise, these transitions are elements of supports of exactly the same t-invariants. From the biological point of view, MCT sets divide Petri-net into some functional blocks. Biological descriptions for all MCT sets are included in Table A3.
Table A3. Biological meaning of non-trivial MCT sets for full model.
Table A3. Biological meaning of non-trivial MCT sets for full model.
MCT-SetContained TransitionsBiological Meaning
m 1 t 2 , t 4 , t 6 Oxidation of BH 4 to BH 2 by ONOO .
m 2 t 10 , t 11 , t 14 Dismutation of NOX-derived O 2 .
m 3 t 21 , t 24 , t 25 Apoptosis of VSMC.
m 4 t 22 , t 26 , t 27 Activation of NF κ B caused by DNA damage by ROS.
m 5 t 35 , t 36 , t 42 Infiltration by inflammatory cells.
m 6 t 48 , t 49 , t 50 Activation of NF κ B via PGE2 activated by COX2.
m 7 t 51 , t 53 , t 54 Collagen proteolysis and rupture of aneurysm.
m 8 t 0 , t 1 Oxidation of LDL.
m 9 t 3 , t 8 Synthesis of O 2 by NOS and BH 2 .
m 10 t 12 , t 13 Reduction of H 2 O 2 by catalase.
m 11 t 18 , t 37 Conversion of ANGI to ANGII by ACE.
m 12 t 23 , t 41 Stimulation of AP1 by ROS.
m 13 t 40 , t 63 Stimulation of IL-6 by hemodynamic stress.
m 14 t 52 , t 57 Influence of aneurysm enlargement on hemodynamic stress.
m 15 t 60 , t 61 Production of HOCl by MPO and H 2 O 2 .

Appendix C. Analysis of t-Clusters

As it has been already mentioned, analysis of Petri-nets can be based on t-invariants. Such invariants are crucial from biological point of view because they correspond to modeled subprocesses. In the case of a large number of t-invariants, the biological significance is not assigned to each of them. The proposed full model contains 736 t-invariants, therefore clustering was used. Biological significance was assigned for groups of t-invariants similar to each other (called t-clusters) and their descriptions are given in Table A4. 17 t-clusters for the proposed Petri-net were distinguished.
Clustering of t-invariants was performed by using a Java application called Holmes [38]. For clustering 7 different algorithms and 8 different distance measures has been used. The results presented in this paper have been obtained using a clustering algorithm based on the average linkage method and Pearson’s measure of distance. The clustering method, the distance measure and the number of t-clusters were selected on the basis of the MSS index (Mean Split Silhouette) and the C-H coefficient (Calinski-Harabasz) [55,56].
Table A4. Biological meaning of t-clusters for full model.
Table A4. Biological meaning of t-clusters for full model.
t-ClusterBiological Meaning
c 1 Uncoupling of NOS by ONOO and O 2 .
c 2 Influence of hemodynamic stress on activation of NF κ B, which activates production of chemokines and inflammatory cytokines, and results in inflammatory cells infiltration. Inflammatory state leads to production of MMPs and enlargement of aneurysm.
c 3 This cluster includes all modeled subprocesses.
c 4 This cluster includes almost all modeled subprocesses except: activation of NOX through TNF α .
c 5 Activation of NF κ B by TNF α . NF κ B produces chemokines.
c 6 Influence of ANGII on enlargement and rupture of aneurysm.
c 7 Synthesis of HOCL by MPO leads to activation of MMPs which results in enlargement and rupture of aneurysm. Activation of NOX by ANGII, TNF α , hemodynamic stress. PGE2 stimulates NF κ B which activates production of chemokines.
c 8 Activation of inflammatory cytokines (TNF α and IL-6) resulting in enlargement of aneurysm.
c 9 This cluster includes almost all modeled subprocesses except:
  • Activation of NF κ B via PGE2 activated by COX2.
  • Activation of NOX through TNF α .
c 10 Activation of NF κ B by TNF α . NF κ B produces cytokines, which activates iNOS. iNOS and eNOS are uncoupled by ONOO and O 2 .
c 11 Activation of NF κ B by ANGII. NF κ B activates production of chemokines and inflammatory cytokines, and results in inflammatory cells infiltration. Cytokines activates iNOS, iNOS and eNOS are uncoupled by ONOO and O 2 .
c 12 Activation of NOX by ANGII and TNF α . PGE2 stimulates NF κ B, which enhances tissue infiltration by inflammatory cells. Inflammatory cells produce cytokines resulting in activation of iNOS. Part of O 2 undergoes dismutation to H 2 O 2 , and other part reacts with NO leading to ONOO synthesis. ONOO leads to uncoupling of NOS. Synthesis of HOCL by MPO leads to activation of MMPs which results in enlargement of aneurysm.
c 13 Activation of NOX by ANGII, TNF α and hemodynamic stress. ANGII and TNF α stimulates NF κ B, which enhances tissue infiltration by inflammatory cells. Inflammatory cells produce cytokines (mainly IL-6, IL-1 and TNF α ) resulting in activation of iNOS. Part of O 2 undergoes dismutation to H 2 O 2 , and other part reacts with NO leading to ONOO synthesis. ONOO leads to uncoupling of NOS. Synthesis of HOCl by MPO leads to activation of MMPs which results in enlargement and rupture of aneurysm.
c 14 Activation of NF κ B by TNF α . NF κ B activates production of chemokines and adhesion particles, which results in inflammatory cells infiltration. Cytokines activates iNOS, iNOS and eNOS are uncoupled by ONOO .
c 15 Activation of NF κ B by ANGII and TNF α . TNF α and IL-6 as inflammatory cytokines lead to expression of MMPs and enlargement of aneurysm. Adhesion particles produced by NF κ B and chemokines result in inflammatory cells infiltration. Cytokines activates iNOS, iNOS and eNOS are uncoupled by ONOO .
c 16 Activation of NF κ B by TNF α leads to inflammatory cytokines production. Inflammatory cytokines activates expression of MMPs, which results in enlargement and rupture of aneurysm.
c 17 This cluster includes almost all modeled subprocesses except:
  • Activation of NF κ B via PGE2 activated by COX2.
  • Stimulation of IL-6 by hemodynamic stress.
  • Activation of NOX through TNF α .
  • Uncoupling of NOS by ONOO and O 2 .
  • Activation of MMPs by HOCl.
  • Stimulation of chemokines production by NF κ B.
  • Production of adhesion particles.
  • Rupture of aneurysm.

References

  1. Newby, D.; Forsythe, R.; McBride, O. Aortic Wall Inflammation Predicts Abdominal Aortic Aneurysm Expansion, Rupture, and Need for Surgical Repair. Circulation 2017, 136, 787–797. [Google Scholar] [CrossRef]
  2. Toghill, B.J.; Saratzis, A.; Bown, M.J. Abdominal aortic aneurysm-an independent disease to atherosclerosis? Cardiovasc. Pathol. 2017, 27, 71–75. [Google Scholar] [CrossRef] [PubMed]
  3. Guzik, B. Rola Anionu Ponadtlenkowego w Etiopatogenezie Tętniaków Aorty Brzusznej. Ph.D. Dissertation, Institute of Cardiology, Jagiellonian University Collegium Medicum, Krakow, Poland, 2006. [Google Scholar]
  4. Golledge, J.; Norman, P. Atherosclerosis and abdominal aortic aneurysm: Cause, response or common risk factors? Arterioscler. Thromb. Vasc. Biol. 2010, 30, 1075–1077. [Google Scholar] [CrossRef] [PubMed]
  5. Trollope, A.F.; Golledge, J. Angiopoietins, abdominal aortic aneurysm and atherosclerosis. Atherosclerosis 2011, 214, 237–243. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Kuivaniemi, H.; Ryer, E.J.; Elmore, J.R.; Tromp, G. Understanding the pathogenesis of abdominal aortic aneurysms. Expert Rev. Cardiovasc. Ther. 2015, 13, 975–987. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Yansheng, Y.; Yanbing, W.; Yibo, Z.; Chang, L. Klotho ameliorates oxidized low density lipoprotein (ox-LDL)-induced oxidative stress via regulating LOX-1 and PI3K/Akt/eNOS pathways. Lipids Health Dis. 2017, 16, 77. [Google Scholar] [CrossRef]
  8. Kubota, Y.; Folsom, A.R.; Ballantyne, C.M.; Tang, W. Lipoprotein(a) and abdominal aortic aneurysm risk: The Atherosclerosis Risk in Communities study. Atherosclerosis 2018, 268, 63–67. [Google Scholar] [CrossRef]
  9. Li, W.; Luo, S.; Luo, J.; Liu, Y.; Ning, B.; Huang, W.; Xue, L.; Chen, J. Predictors Associated With Increased Prevalence of Abdominal AorticAneurysm in Chinese Patients with Atherosclerotic Risk Factors. Eur. J. Vasc. Endovasc. Surg. 2017, 54, 43–49. [Google Scholar] [CrossRef]
  10. Hobbs, S.D.; Claridge, M.W.; Quick, C.R.; Day, N.E.; Bradbury, A.W.; Wilmink, A.B. LDL cholesterol is associated with small abdominal aortic aneurysms. Eur. J. Vasc. Endovasc. Surg. 2003, 26, 618–622. [Google Scholar] [CrossRef]
  11. Liu, J.; Lu, H.; Howatt, D.A.; Balakrishnan, A.; Moorleghen, J.J.; Sorci-Thomas, M.; Cassis, L.A.; Daugherty, A. Associations of ApoAI and ApoB-containing lipoproteins with AngII-induced abdominal aortic aneurysms in mice. Arterioscler. Thromb. Vasc. Biol. 2015, 35, 1826–1834. [Google Scholar] [CrossRef]
  12. Prins, P.A.; Hill, M.F.; Airey, D.; Nwosu, S.; Perati, P.R.; Tavori, H.; Linton, M.F.; Kon, V.; Fazio, S.; Sampson, U.K. Angiotensin-induced abdominal aortic aneurysms in hypercholesterolemic mice: Role of serum cholesterol and temporal effects of exposure. PLoS ONE 2014, 9, e84517. [Google Scholar] [CrossRef] [PubMed]
  13. Frösen, J.; Tulamo, R.; Heikura, T.; Sammalkorpi, S.; Niemelä, M.; Hernesniemi, J.; Levonen, A.L.; Hörkkö, S.; Ylä-Herttuala, S. Lipid accumulation, lipid oxidation, and low plasma levels of acquired antibodies against oxidized lipids associate with degeneration and rupture of the intracranial aneurysm wall. Acta Neuropathol. Commun. 2013, 1, 71. [Google Scholar] [CrossRef] [PubMed]
  14. Gao, L.; Siu, K.L.; Chalupsky, K.; Nguyen, A.; Chen, P.; Weintraub, N.L.; Galis, Z.; Cai, H. Role of uncoupled endothelial nitric oxide synthase in abdominal aortic aneurysm formation: Treatment with folic acid. Hypertension 2012, 59, 158–166. [Google Scholar] [CrossRef] [PubMed]
  15. Muller, G.; Morawietz, H. Nitric oxide, NAD(P)H oxidase, and atherosclerosis. Antioxid. Redox Signal. 2009, 11, 1711–1731. [Google Scholar] [CrossRef] [PubMed]
  16. Schulz, E.; Gori, T.; Münzel, T. Oxidative stress and endothelial dysfunction in hypertension. Hypertens. Res. 2011, 34, 665–673. [Google Scholar] [CrossRef] [PubMed]
  17. McCormick, M.L.; Gavrila, D.; Weintraub, N.L. Role of oxidative stress in the pathogenesis of abdominal aortic aneurysms. Arterioscler. Thromb. Vasc. Biol. 2007, 27, 461–469. [Google Scholar] [CrossRef] [PubMed]
  18. Fleming, I.; Mohamed, A.; Galle, J.; Turchanowa, L.; Brandes, R.P.; Fisslthaler, B.; Busse, R. Oxidized low-density lipoprotein increases superoxide production by endothelial nitric oxide synthase by inhibiting PKCalpha. Cardiovasc. Res. 2005, 65, 897–906. [Google Scholar] [CrossRef]
  19. Honjo, T.; Otsui, K.; Shiraki, R.; Kawashima, S.; Sawamura, T.; Yokoyama, M.; Inoue, N. Essential role of NOXA1 in generation of reactive oxygen species induced by oxidized low-density lipoprotein in human vascular endothelial cells. Endothelium 2008, 15, 137–141. [Google Scholar] [CrossRef]
  20. Wang, Y.; Che, J.; Zhao, H.; Tang, J.; Shi, G. Osthole alleviates oxidized low-density lipoprotein-induced vascular endothelial injury through suppression of transforming growth factor-β1/Smad pathway. Int. Immunopharmacol. 2018, 65, 373–381. [Google Scholar] [CrossRef]
  21. Yang, H.; Shi, M.; Richardson, A.; Vijg, J.; Guo, Z. Attenuation of leukocyte-endothelium interaction by antioxidant enzymes. Free Radic. Biol. Med. 2003, 35, 266–276. [Google Scholar] [CrossRef]
  22. Kim, H.W.; Blomkalns, A.L.; Ogbi, M.; Thomas, M.; Gavrila, D.; Neltner, B.S.; Cassis, L.A.; Thompson, R.W.; Weiss, R.M.; Lindower, P.D.; et al. Role of myeloperoxidase in abdominal aortic aneurysm formation: Mitigation by taurine. Am. J. Physiol. Heart Circ. Physiol. 2017, 1, 1168–1179. [Google Scholar] [CrossRef] [PubMed]
  23. Emeto, T.I.; Moxon, J.V.; Au, M.; Golledge, J. Oxidative stress and abdominal aortic aneurysm: Potential treatment targets. Clin. Sci. 2016, 130, 301–315. [Google Scholar] [CrossRef] [PubMed]
  24. Schramm, A.; Matusik, P.; Osmenda, G.; Guzik, T.J. Targeting NADPH oxidases in vascular pharmacology. Vascul. Pharmacol. 2012, 56, 216–231. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Xiong, W.; Mactaggart, J.; Knispel, R.; Worth, J.; Zhu, Z.; Yulong, L.; Sun, Y.; Baxter, T.; Johanning, J. Inhibition of reactive oxygen species attenuates aneurysm formation in a murine model. Atherosclerosis 2009, 202, 128–134. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Fukatsu, H.; Koide, N.; Tada-Oikawa, S.; Izuoka, K.; Ikegami, A.; Ichihara, S.; Ukaji, T.; Morita, N.; Naiki, Y.; Komatsu, T.; et al. NF κB inhibitor DHMEQ inhibits titanium dioxide nanoparticle induced interleukin 1β production: Inhibition of the PM2.5 induced inflammation model. Mol. Med. Rep. 2018, 18, 5279–5285. [Google Scholar] [CrossRef]
  27. Akerman, A.W.; Stroud, R.E.; Bans, R.W.; Grespin, R.T.; McDonald, L.T.; LaRue, R.A.; Mukherjee, R.; Ikonomidis, J.S.; Jones, J.A.; Ruddy, J.M. Elevated Wall Tension Initiates Interleukin-6 Expression and Abdominal Aortic Dilation. Ann. Vasc. Surg. 2018, 46, 193–204. [Google Scholar] [CrossRef]
  28. Ray, R.; Shah, A.M. NADPH oxidase and endothelial cell function. Clin. Sci. 2005, 109, 217–226. [Google Scholar] [CrossRef] [Green Version]
  29. Siu, K.L.; Li, Q.; Zhang, Y.; Guo, J.; Youn, J.Y.; Du, J.; Cai, H. NOX isoforms in the development of abdominal aortic aneurysm. Redox Biol. 2017, 11, 118–125. [Google Scholar] [CrossRef]
  30. Ramella, M.; Boccafoschi, F.; Bellofatto, K.; Follenzi, A.; Fusaro, L.; Boldorini, R.; Casella, F.; Porta, C.; Settembrini, P.; Cannas, M. Endothelial MMP-9 drives the inflammatory response in abdominal aortic aneurysm (AAA). Am. J. Transl. Res. 2017, 9, 5485. [Google Scholar]
  31. Madamanchi, N.R.; Runge, M.S. Redox signaling in cardiovascular health and disease. Free Radic. Biol. Med. 2013, 61, 473–501. [Google Scholar] [CrossRef] [Green Version]
  32. Dawson, J.; Cockerill, G.W.; Choke, E.; Belli, A.M.; Loftus, I.; Thompson, M.M. Aortic aneurysms secrete interleukin-6 into the circulation. J. Vasc. Surg. 2007, 45, 350–356. [Google Scholar] [CrossRef] [Green Version]
  33. Sprague, A.H.; Khalil, R.A. Inflammatory cytokines in vascular dysfunction and vascular disease. Biochem. Pharmacol. 2009, 78, 539–552. [Google Scholar] [CrossRef] [Green Version]
  34. Wang, S.K.; Green, L.A.; Gutwein, A.R.; Drucker, N.A.; Motaganahalli, R.L.; Gupta, A.K.; Fajardo, A.; Murphy, M.P. Description of human AAA by cytokine and immune cell aberrations compared to risk-factor matched controls. Surgery 2018, 164, 354–358. [Google Scholar] [CrossRef]
  35. Choke, E.; Cockerill, G.; Wilson, W.R.; Sayed, S.; Dawson, J.; Loftus, I.; Thompson, M.M. A review of biological factors implicated in abdominal aortic aneurysm rupture. Eur. J. Vasc. Endovasc. Surg. 2005, 30, 227–244. [Google Scholar] [CrossRef]
  36. Formanowicz, D.; Gutowska, K.; Formanowicz, P. Theoretical Studies on the Engagement of Interleukin 18 in the Immuno-Inflammatory Processes Underlying Atherosclerosis. Int. J. Mol. Sci. 2018, 19, 3476. [Google Scholar] [CrossRef]
  37. Formanowicz, D.; Radom, M.; Zawierucha, P.; Formanowicz, P. Petri-net-based approach to modeling and analysis of selected aspects of the molecular regulation of angiogenesis. PLoS ONE 2017, 12, e0173020. [Google Scholar] [CrossRef]
  38. Radom, M.; Rybarczyk, A.; Szawulak, B.; Andrzejewski, H.; Chabelski, P.; Kozak, A.; Formanowicz, P. Holmes: A graphical tool for development, simulation and analysis of Petri-net-based models of complex biological systems. Bioinformatics 2017, 33, 3822–3823. [Google Scholar] [CrossRef]
  39. Formanowicz, P.; Nowicki, M.; Formanowicz, P. Mathematical modeling of aortic aneurysm progression. In New Approaches to Aortic Diseases from Valve to Abdominal Bifurcation; Tintoiu, I.C., Ursulescu, A., Elefteriades, J.A., Underwood, M.J., Droc, I., Eds.; Academic Press: London, UK, 2018. [Google Scholar]
  40. Zmijewski, J.W.; Landar, A.; Watanabe, N.; Dickinson, D.A.; Noguchi, N.; Darley-Usmar, V.M. Cell signalling by oxidized lipids and the role of reactive oxygen species in the endothelium. Biochem. Soc. Trans. 2005, 33, 1385–1389. [Google Scholar] [CrossRef] [Green Version]
  41. Zhang, J.; Schmidt, J.; Ryschich, E.; Mueller-Schilling, M.; Schumacher, H.; Allenberg, J.R. Inducible nitric oxide synthase is present in human abdominal aortic aneurysm and promotes oxidative vascular injury. J. Vasc. Surg. 2003, 38, 360–367. [Google Scholar] [CrossRef]
  42. Gomez, I.; Ozen, G.; Deschildre, C.; Amgoud, Y.; Boubaya, L.; Gorenne, I.; Benyahia, C.; Roger, T.; Lesèche, G.; Galardon, E.; et al. Reverse Regulatory Pathway (H2S/PGE2/MMP) in Human Aortic Aneurysm and Saphenous Vein Varicosity. PLoS ONE 2016, 11, e0158421. [Google Scholar] [CrossRef]
  43. Yokoyama, U.; Ishiwata, R.; Jin, M.H.; Kato, Y.; Suzuki, O.; Jin, H.; Ichikawa, Y.; Kumagaya, S.; Katayama, Y.; Fujita, T.; et al. Inhibition of EP4 signaling attenuates aortic aneurysm formation. PLoS ONE 2012, 7, e36724. [Google Scholar] [CrossRef]
  44. Russell-Puleri, S.; Dela Paz, N.G.; Adams, D.; Chattopadhyay, M.; Cancel, L.; Ebong, E.; Orr, A.W.; Frangos, J.A.; Tarbell, J.M. Fluid shear stress induces upregulation of COX-2 and PGI2 release in endothelial cells via a pathway involving PECAM-1, PI3K, FAK, and p38. Am. J. Physiol. Heart Circ. Physiol. 2017, 312, H485–H500. [Google Scholar] [CrossRef]
  45. Cao, R.Y.; St Amand, T.; Li, X.; Yoon, S.H.; Wang, C.P.; Song, H.; Maruyama, T.; Brown, P.M.; Zelt, D.T.; Funk, C.D. Prostaglandin receptor EP4 in abdominal aortic aneurysms. Am. J. Pathol. 2012, 181, 313–321. [Google Scholar] [CrossRef]
  46. Ozen, G.; Gomez, I.; Daci, A.; Deschildre, C.; Boubaya, L.; Teskin, O.; Uydeş-Doğan, B.S.; Jakobsson, P.J.; Longrois, D.; Topal, G.; et al. Inhibition of microsomal PGE synthase-1 reduces human vascular tone by increasing PGI2: A safer alternative to COX-2 inhibition. Br. J. Pharmacol. 2017, 174, 4087–4098. [Google Scholar] [CrossRef]
  47. Petri, C.A. Communication with Automata; Schriften des Instituts fur Instrumentelle Mathematik: Bonn, Germany, 1962. (In German) [Google Scholar]
  48. Koch, I.; Reisig, W.; Schreiber, F. (Eds.) Modeling in Systems Biology; The Petri Net Approach; Springer: London, UK, 2011; ISBN 978-1-84996-474-6. [Google Scholar]
  49. Murata, T. Petri Nets: Properties, Analysis and Applications. Proc. IEEE 1989, 77, 541–580. [Google Scholar] [CrossRef]
  50. Reising, W. Understanding Petri Nets. Modeling Techniques, Analysis Methods, Case Studies; Springer: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  51. Sackmann, A.; Heiner, M.; Koch, I. Application of Petri net-based analysis techniques to signal transduction pathways. BMC Bioinform. 2006, 7, 482. [Google Scholar] [CrossRef]
  52. Grafahrend-Belau, E.; Schreiber, F.; Heiner, M.; Sackmann, A.; Junker, B.H.; Grunwald, S.; Speer, A.; Winder, K.; Koch, I. Modularization of biochemical networks based on classification of Petri net t-invariants. BMC Bioinform. 2008, 9. [Google Scholar] [CrossRef]
  53. Formanowicz, D.; Kozak, A.; Głowacki, T.; Radom, M.; Formanowicz, P. Hemojuvelin–Hepcidin axis modeled and analyzed using Petri nets. J. Biomed. Inform. 2013, 46, 1030–1043. [Google Scholar] [CrossRef]
  54. Formanowicz, D.; Radom, M.; Rybarczyk, A.; Formanowicz, P. The role of Fenton reaction in ROS-induced toxicity underlying atherosclerosis—Modeled and analyzed using a Petri-net-based approach. BioSystems 2018, 165, 71–87. [Google Scholar] [CrossRef]
  55. Calinski, T.; Harabasz, J. A dendrite method for cluster analysis. Commun. Stat. 1974, 3, 1–27. [Google Scholar] [CrossRef]
  56. Kaufman, L.; Rousseeuw, P.J. Finding Groups in Data: An Introduction to Cluster Analysis; John Wiley and Sons: New York, NY, USA, 1990; ISBN 0-471-87876-6. [Google Scholar]
Figure 1. Scheme of the proposed full model of AAA development. References was marked above the arrows for analysis facilitation [3,7,10,14,15,16,17,18,19,20,21,22,23,24,25,27,28,29,30,31,33,35,40,41,42,43,44,45,46].
Figure 1. Scheme of the proposed full model of AAA development. References was marked above the arrows for analysis facilitation [3,7,10,14,15,16,17,18,19,20,21,22,23,24,25,27,28,29,30,31,33,35,40,41,42,43,44,45,46].
Ijms 20 03909 g001
Figure 2. The general scheme of work.
Figure 2. The general scheme of work.
Ijms 20 03909 g002
Table 1. Significance analysis of selected subprocesses for full and reduced model. The following abbreviations in the columns headings has been distinguished: trans—transition, inv— t-invariant, frequency trans/inv—frequency of occurrence of selected transition in all t-invariants.
Table 1. Significance analysis of selected subprocesses for full and reduced model. The following abbreviations in the columns headings has been distinguished: trans—transition, inv— t-invariant, frequency trans/inv—frequency of occurrence of selected transition in all t-invariants.
Full Model
with oxLDL
Reduced Model
without oxLDL
SubprocessElementary SubprocessFrequency
Trans/Inv
Percentage of Transition
in Whole System
(100% = 736 inv)
Frequency
Trans/Inv
Percentage of Transition
in Whole System
(100% = 325 inv)
production of O 2
by NOX
71396.88%30292.92%
production of O 2
by NOS and oxLDL
20527.85%--
oxidative
stress
production of O 2
by NOS and BH 2
36449.46%21867.08%
production of ONOO
by O 2 and NO
18224.73%10933.54%
production of H 2 O 2
through dismutation
71396.88%30292.92%
production of HOCl
by MPO and H 2 O 2
70796.06%29691.08%
activation by ROS19726.77%8827.08%
activation by cytokines19226.09%9930.46%
influence
of MMPs
activation of MMP2
by oxLDL
13217.93%--
activation of MPO
by HOCl
18425.00%7523.08%
LDL
oxidation
LDL oxidation41155.84%--
activation by oxLDL60.82%--
influence
of NOX
activation by TNF α 34246.47%14344.00%
activation by
hemodynamic stress
14519.70%7924.31%
activation by ANGII29139.54%11535.38%
production by
inflammatory cells
55275.00%21164.92%
production by TNF α 10814.67%5316.31%
production of
inflammatory
cytokines
production by NF κ B11515.63%5516.92%
production by AP1263.53%206.15%
stimulation of IL-610814.67%5316.31%
enlargement
of AAA
enlargement of AAA62985.46%27484.31%
rupture
of AAA
rupture of AAA26035.33%10632.62%
Table 2. Knockout analysis of selected subprocesses for full and reduced model.
Table 2. Knockout analysis of selected subprocesses for full and reduced model.
Oxidative
Stress
Influence
of MMPs
LDL
Oxidation
Influence
of NOX
Production of
Inflammatory
Cytokines
Enlargement
of AAA
Rupture
of AAA
full model
(with oxLDL)
98.78%80.84%55.84%96.88%93.75%85.46%35.33%
reduced model
(without oxLDL)
97.23%78.15%-92.92%90.15%84.31%32.62%
knockouted
transitions
t 2 , t 8 , t 10 ,
t 14 , t 44 , t 60
t 32 , t 33 ,
t 56 , t 62
t 0 t 15 , t 16 ,
t 17 , t 47
t 28 , t 29 , t 40 ,
t 42 , t 63
t 52 t 54
Table 3. The key results from analysis.
Table 3. The key results from analysis.
Full Model with oxLDL
Significance Analysis of
Selected Subprocesses
Knockout of Full Model
(Knockout Transition t 0 )
SubprocessFrequency
Trans/Inv
Percentage
Contribution in
Whole System
(100% = 736 inv)
Frequency
Trans/Inv
Percentage
Contribution in
Whole System
(100% = 325 inv)
Difference between
Percentage Contribution
(736 inv vs. 325 inv)
enlargement of AAA62985.46%27484.31%1.15 p.p.
rupture of AAA26035.33%10632.62%2.71 p.p.

Share and Cite

MDPI and ACS Style

Gutowski, Ł.; Gutowska, K.; Pioruńska-Stolzmann, M.; Formanowicz, P.; Formanowicz, D. Systems Approach to Study Associations between OxLDL and Abdominal Aortic Aneurysms. Int. J. Mol. Sci. 2019, 20, 3909. https://doi.org/10.3390/ijms20163909

AMA Style

Gutowski Ł, Gutowska K, Pioruńska-Stolzmann M, Formanowicz P, Formanowicz D. Systems Approach to Study Associations between OxLDL and Abdominal Aortic Aneurysms. International Journal of Molecular Sciences. 2019; 20(16):3909. https://doi.org/10.3390/ijms20163909

Chicago/Turabian Style

Gutowski, Łukasz, Kaja Gutowska, Maria Pioruńska-Stolzmann, Piotr Formanowicz, and Dorota Formanowicz. 2019. "Systems Approach to Study Associations between OxLDL and Abdominal Aortic Aneurysms" International Journal of Molecular Sciences 20, no. 16: 3909. https://doi.org/10.3390/ijms20163909

APA Style

Gutowski, Ł., Gutowska, K., Pioruńska-Stolzmann, M., Formanowicz, P., & Formanowicz, D. (2019). Systems Approach to Study Associations between OxLDL and Abdominal Aortic Aneurysms. International Journal of Molecular Sciences, 20(16), 3909. https://doi.org/10.3390/ijms20163909

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop