Different approaches to modelling the cost-effectiveness of schistosomiasis control

This paper reviews three different approaches to modelling the cost-effectiveness of schistosomiasis control. Although these approaches vary in their assessment of costs, the major focus of the paper is on the evaluation of effectiveness. The first model presented is a static economic model which assesses effectiveness in terms of the proportion of cases cured. This model is important in highlighting that the optimal choice of chemotherapy regime depends critically on the level of budget constraint, the unit costs of screening and treatment, the rates of compliance with screening and chemotherapy and the prevalence of infection. The limitations of this approach is that it models the cost-effectiveness of only one cycle of treatment, and effectiveness reflects only the immediate impact of treatment. The second model presented is a prevalence-based dynamic model which links prevalence rates from one year to the next, and assesses effectiveness as the proportion of cases prevented. This model was important as it introduced the concept of measuring the long-term impact of control by using a transmission model which can assess reduction in infection through time, but is limited to assessing the impact only on the prevalence of infection. The third approach presented is a theoretical framework which describes the dynamic relationships between infection and morbidity, and which assesses effectiveness in terms of case-years prevented of infection and morbidity. The use of this model in assessing the cost-effectiveness of age-targeted treatment in controlling Schistosoma mansoni is explored in detail, with respect to varying frequencies of treatment and the interaction between drug price and drug efficacy.

Squamous cell carcinoma in 6 patients with chronic discharging osteomyelitis: Is it really rare?

Introduction: Development of a squamous cell carcinoma (Marjolin's ulcer) is a rare but well-known complication of chronic discharging osteomyelitis. A high index of suspicion and highquality of histopathological examination are paramount in order to make the correct diagnosis. Methods: During a 15-year period (1993 and 2008), patients with long-standing chronic osteomyelitis with clinical symptoms of >1 year of duration, were retrospectively reviewed. Included were patients with histologically confirmed squamous-cell carcinoma associated with chronic wound overlying the site of chronic osteomyelitis. Clinical features and treatment approaches of these patients were analyzed. Results: During the study period, 6 patients were identified (2 women and 4 men) aged 52 to 67 years (mean 59 years). All patients had a long history of chronic discharging osteomyelitis (12, 19, 21, 30, 39 and 40 years), localized in the lower (5 patients) or upper extremitiy (1 patient). All tumors were histologically highly-differentiated squamous-cell carcinomas involving the deep soft tissues and the bone. 5 of 6 patients were initially misdiagnosed as chronic bone infection since bacteria were isolated in wound swabs, including Staphylococcus aureus (n = 3), Escherichia coli (n = 1) and Staphylococcus epidermidis (n = 1). Treatment consisted of major amputation in 4 patients and radical surgical excision in 2 patients refusing amputation. 4 patients were lost to follow-up due to return to their country of origin, the remaining 2 patients were without signs of tumor recurrence (both after major amputation). Conclusion: Malignant transformation of is a rare, but serious complication of long-standing discharging chronic osteomyelitis. All 6 patients were diagnosed after >10 years of persistent or recurrent wounds. It should be particularly suspected in case of a pathological fracture, development of exophytic mass and changes in local pattern of the ulcers itself. Multiple biopsies, including deep soft tissues and bone, are recommended to distinguish between chronic osteomyelitis, pseudoepitheliomatous hyperplasia and highly-differentiated squamous cell carcinoma. Early diagnosis and a team approach (orthopaedic and plastic surgeon) are crucial for optimal management of the patient

Oxygen exchange in galls of Eurosta solidaginis (Diptera: Tephritidae)

Hypoxia in plant tissue should affect animals living within. Gallmakers stimulate their plant hosts to produce the gall they inhabit and feed on, and also influence the gall phenotype for other adaptations, such as defense against predators. The potential for hypoxia in galls of Eurosta solidaginis was studied in the context of potential adaptations to gall oxygen level, using a combination of direct measurement, mathematical modelling, and respirometry on both gallmakers and hosts. Modelling results suggested mild hypoxia tolerable to the larva persists for most of the growth season, whereas more severe hypoxia may occur earlier in fully-grown young galls. Field data from one of the two years studied showed hypoxia more severe than expected, and coincided with adverse weather conditions and high larval mortality. The hypoxia may be related to host response to adverse weather. Whether hypoxia directly caused larval mortality requires further study.

A "Trojan horse" strategy to reverse drug-resistance in brain tumors

Los gliomas malignos representan una de las formas más agresivas de los tumores del sistema nervioso central (SNC). De acuerdo con la clasificación de los tumores cerebrales de la Organización Mundial de la Salud (OMS), los astrocitomas han sido categorizados en cuatro grados, determinados por la patología subyacente. Es así como los gliomas malignos (o de alto grado) incluyen el glioma anaplásico (grado III) así como el glioblastoma multiforme (GBM, grado IV),estos últimos los más agresivos con el peor pronóstico (1). El manejo terapéutico de los tumores del SNC se basa en la cirugía, la radioterapia y la quimioterapia, dependiendo de las características del tumor, el estadio clínico y la edad (2),(3), sin embargo ninguno de los tratamientos estándar es completamente seguro y compatible con una calidad de vida aceptable (3), (4). En general, la quimioterapia es la primera opción en los tumores diseminados, como el glioblastoma invasivo y el meduloblastoma de alto riesgo o con metástasis múltiple, pero el pronóstico en estos pacientes es muy pobre (2),(3). Solamente nuevas terapias dirigidas (2) como las terapias anti-angiogénicas (4); o terapias génicas muestran un beneficio real en grupos limitados de pacientes con defectos moleculares específicos conocidos (4). De este modo, se hace necesario el desarrollo de nuevas terapias farmacológicas para atacar los tumores cerebrales. Frente a las terapias los gliomas malignos son con frecuencia quimioresistentes, y esta resistencia parece depender de al menos dos mecanismos: en primer lugar, la pobre penetración de muchas drogas anticáncer a través de la barrera hematoencefálica (BBB: Blood Brain Barrier), la barrera del fluido sangre-cerebroespinal (BCSFB: Blood-cerebrospinal fluid barrier) y la barrera sangre-tumor (BTB: blood-tumor barrier). Dicha resistencia se debe a la interacción de la droga con varios transportadores o bombas de eflujo de droga ABC (ABC: ATP-binding cassette) que se sobre expresan en las células endoteliales o epiteliales de estas barreras. En segundo lugar, estos transportadores de eflujo de drogas ABC propios de las células tumorales confieren un fenotipo conocido como resistencia a multidrogas (MDR: multidrug resistance), el cual es característico de varios tumores sólidos. Este fenotipo también está presente en los tumores del SNC y su papel en gliomas es objeto de investigación (5). Por consiguiente el suministro de medicamentos a través de la BBB es uno de los problemas vitales en los tratamientos de terapia dirigida. Estudios recientes han demostrado que algunas moléculas pequeñas utilizadas en estas terapias son sustratos de la glicoproteína P (Pgp: P-gycoprotein), así como también de otras bombas de eflujo como las proteínas relacionadas con la resistencia a multidrogas (MRPs: multidrug resistance-related proteins (MRPs) o la proteína relacionada con cáncer de seno (BCRP: breast-cancer resistance related protein)) que no permiten que las drogas de este tipo alcancen el tumor (1). Un sustrato de Pgp y BCRP es la DOXOrubicina (DOXO), un fármaco utilizado en la terapia anti cáncer, el cual es muy eficaz para atacar las células del tumor cerebral in vitro, pero con un uso clínico limitado por la poca entrega a través de la barrera hematoencefálica (BBB) y por la resistencia propia de los tumores. Por otra parte las células de BBB y las células del tumor cerebral tienen también proteínas superficiales, como el receptor de la lipoproteína de baja densidad (LDLR), que podría utilizarse como blanco terapéutico en BBB y tumores cerebrales. Es asi como la importancia de este estudio se basa en la generación de estrategias terapéuticas que promuevan el paso de las drogas a través de la barrera hematoencefalica y tumoral, y a su vez, se reconozcan mecanismos celulares que induzcan el incremento en la expresión de los transportadores ABC, de manera que puedan ser utilizados como blancos terapéuticos.Este estudio demostró que el uso de una nueva estrategia basada en el “Caballo de Troya”, donde se combina la droga DOXOrubicina, la cual es introducida dentro de un liposoma, salvaguarda la droga de manera que se evita su reconocimiento por parte de los transportadores ABC tanto de la BBB como de las células del tumor. La construcción del liposoma permitió utilizar el receptor LDLR de las células asegurando la entrada a través de la BBB y hacia las células tumorales a través de un proceso de endocitosis. Este mecanismo fue asociado al uso de estatinas o drogas anticolesterol las cuales favorecieron la expresión de LDLR y disminuyeron la actividad de los transportadores ABC por nitración de los mismos, incrementando la eficiencia de nuestro Caballo de Troya. Por consiguiente demostramos que el uso de una nueva estrategia o formulación denominada ApolipoDOXO más el uso de estatinas favorece la administración de fármacos a través de la BBB, venciendo la resistencia del tumor y reduciendo los efectos colaterales dosis dependiente de la DOXOrubicina. Además esta estrategia del "Caballo de Troya", es un nuevo enfoque terapéutico que puede ser considerado como una nueva estrategia para aumentar la eficacia de diferentes fármacos en varios tumores cerebrales y garantiza una alta eficiencia incluso en un medio hipóxico,característico de las células cancerosas, donde la expresión del transportador Pgp se vió aumentada. Teniendo en cuenta la relación entre algunas vías de señalización reconocidas como moduladores de la actividad de Pgp, este estudio presenta no solo la estrategia del Caballo de Troya, sino también otra propuesta terapéutica relacionada con el uso de Temozolomide más DOXOrubicina. Esta estrategia demostró que el temozolomide logra penetrar la BBB por que interviene en la via de señalización de la Wnt/GSK3/β-catenina, la cual modula la expresión del transportador Pgp. Se demostró que el TMZ disminuye la proteína y el mRNA de Wnt3 permitiendo plantear la hipótesis de que la droga al disminuir la transcripción del gen Wnt3 en células de BBB, incrementa la activación de la vía fosforilando la β-catenina y conduciendo a disminuir la β-catenina nuclear y por tanto su unión al promotor del gen mdr1. Con base en los resultados este estudio permitió el reconocimiento de tres mecanismos básicos relacionados con la expresión de los transportadores ABC y asociados a las estrategias empleadas: el primero fue el uso de las estatinas, el cual condujo a la nitración de los transportadores disminuyendo su actividad por la via del factor de transcripción NFκB; el segundo a partir del uso del temozolomide, el cual metila el gen de Wnt3 reduciendo la actividad de la via de señalización de la la β-catenina, disminuyendo la expresión del transportador Pgp. El tercero consistió en la determinación de la relación entre el eje RhoA/RhoA quinasa como un modulador de la via (no canónica) GSK3/β-catenina. Se demostró que la proteína quinasa RhoA promovió la activación de la proteína PTB1, la cual al fosforilar a GSK3 indujo la fosforilación de la β-catenina, lo cual dio lugar a su destrucción por el proteosoma, evitando su unión al promotor del gen mdr1 y por tanto reduciendo su expresión. En conclusión las estrategias propuestas en este trabajo incrementaron la citotoxicidad de las células tumorales al aumentar la permeabilidad no solo de la barrera hematoencefálica, sino también de la propia barrera tumoral. Igualmente, la estrategia del “Caballo de Troya” podría ser útil para la terapia de otras enfermedades asociadas al sistema nervioso central. Por otra parte estos estudios indican que el reconocimiento de mecanismos asociados a la expresión de los transportadores ABC podría constituir una herramienta clave en el desarrollo de nuevas terapias anticáncer.

Mathematical Optimization Applied to Thermal and Electrical Energy Systems

This Thesis aims at building and discussing mathematical models applications focused on Energy problems, both on the thermal and electrical side. The objective is to show how mathematical programming techniques developed within Operational Research can give useful answers in the Energy Sector, how they can provide tools to support decision making processes of Companies operating in the Energy production and distribution and how they can be successfully used to make simulations and sensitivity analyses to better understand the state of the art and convenience of a particular technology by comparing it with the available alternatives. The first part discusses the fundamental mathematical background followed by a comprehensive literature review about mathematical modelling in the Energy Sector. The second part presents mathematical models for the District Heating strategic network design and incremental network design. The objective is the selection of an optimal set of new users to be connected to an existing thermal network, maximizing revenues, minimizing infrastructure and operational costs and taking into account the main technical requirements of the real world application. Results on real and randomly generated benchmark networks are discussed with particular attention to instances characterized by big networks dimensions. The third part is devoted to the development of linear programming models for optimal battery operation in off-grid solar power schemes, with consideration of battery degradation. The key contribution of this work is the inclusion of battery degradation costs in the optimisation models. As available data on relating degradation costs to the nature of charge/discharge cycles are limited, we concentrate on investigating the sensitivity of operational patterns to the degradation cost structure. The objective is to investigate the combination of battery costs and performance at which such systems become economic. We also investigate how the system design should change when battery degradation is taken into account.

Modelling of transient heat transfer in annealing furnaces

This is a study of heat transfer in a lift-off furnace which is employed in the batch annealing of a stack of coils of steel strip. The objective of the project is to investigate the various factors which govern the furnace design and the heat transfer resistances, so as to reduce the time of the annealing cycle, and hence minimize the operating costs. The work involved mathematical modelling of patterns of gas flow and modes of heat transfer. These models are: Heat conduction and its conjectures in the steel coils;Convective heat transfer in the plates separating the coils in the stack and in other parts of the furnace; and Radiative and convective heat transfer in the furnace by using the long furnace model. An important part of the project is the development of numerical methods and computations to solve the transient models. A limited number of temperature measurements was available from experiments on a test coil in an industrial furnace. The mathematical model agreed well with these data. The model has been used to show the following characteristics of annealing furnaces, and to suggest further developments which would lead to significant savings: - The location of the limiting temperature in a coil is nearer to the hollow core than to the outer periphery. - Thermal expansion of the steel tends to open the coils, reduces their thermal conductivity in the radial direction, and hence prolongs the annealing cycle. Increasing the tension in the coils and/or heating from the core would overcome this heat transfer resistance. - The shape and dimensions of the convective channels in the plates have significant effect on heat convection in the stack. An optimal design of a channel is shown to be of a width-to-height ratio equal to 9. - Increasing the cooling rate, by using a fluidized bed instead of the normal shell and tube exchanger, would shorten the cooling time by about 15%, but increase the temperature differential in the stack. - For a specific charge weight, a stack of different-sized coils will have a shorter annealing cycle than one of equally-sized coils, provided that production constraints allow the stacking order to be optimal. - Recycle of hot flue gases to the firing zone of the furnace would produce a. decrease in the thermal efficiency up to 30% but decreases the heating time by about 26%.

Hypoxia-inducible factor (HIF) network:insights from mathematical models

Oxygen is a crucial molecule for cellular function. When oxygen demand exceeds supply, the oxygen sensing pathway centred on the hypoxia inducible factor (HIF) is switched on and promotes adaptation to hypoxia by up-regulating genes involved in angiogenesis, erythropoiesis and glycolysis. The regulation of HIF is tightly modulated through intricate regulatory mechanisms. Notably, its protein stability is controlled by the oxygen sensing prolyl hydroxylase domain (PHD) enzymes and its transcriptional activity is controlled by the asparaginyl hydroxylase FIH (factor inhibiting HIF-1).To probe the complexity of hypoxia-induced HIF signalling, efforts in mathematical modelling of the pathway have been underway for around a decade. In this paper, we review the existing mathematical models developed to describe and explain specific behaviours of the HIF pathway and how they have contributed new insights into our understanding of the network. Topics for modelling included the switch-like response to decreased oxygen gradient, the role of micro environmental factors, the regulation by FIH and the temporal dynamics of the HIF response. We will also discuss the technical aspects, extent and limitations of these models. Recently, HIF pathway has been implicated in other disease contexts such as hypoxic inflammation and cancer through crosstalking with pathways like NF?B and mTOR. We will examine how future mathematical modelling and simulation of interlinked networks can aid in understanding HIF behaviour in complex pathophysiological situations. Ultimately this would allow the identification of new pharmacological targets in different disease settings.

The Basic Reproduction Number Obtained From Jacobian And Next Generation Matrices - A Case Study Of Dengue Transmission Modelling.

The basic reproduction number is a key parameter in mathematical modelling of transmissible diseases. From the stability analysis of the disease free equilibrium, by applying Routh-Hurwitz criteria, a threshold is obtained, which is called the basic reproduction number. However, the application of spectral radius theory on the next generation matrix provides a different expression for the basic reproduction number, that is, the square root of the previously found formula. If the spectral radius of the next generation matrix is defined as the geometric mean of partial reproduction numbers, however the product of these partial numbers is the basic reproduction number, then both methods provide the same expression. In order to show this statement, dengue transmission modelling incorporating or not the transovarian transmission is considered as a case study. Also tuberculosis transmission and sexually transmitted infection modellings are taken as further examples.

Structural alterations in adult rat carotid bodies exposed to hyperbaric oxygenation

Inhibition of carotid body (CB) function is the main mechanism involved in the attenuation of respiratory drive observed during hyperoxia However, only a few studies at 5 0 atmospheres absolutes (ATA) have analyzed carotid body structure or function in hyperbaric oxygenation (HBO(2)) situations We hypothesized that rats will present CB structural alterations when exposed to different lower hyperbaric oxygen doses enough to alter their chemosensory response to hypoxia Methods - Twenty-one adult male Wistar rats, divided into three groups, were maintained in room air or exposed to O(2) at 2 4 or 3 0 ATA for six hours Histological, ultrastructural and immunohistochemical analyses for neuronal nitric oxide synthase (nNOS) and F2-isoprostane were performed in the excised CBs Results - Histological analyses revealed signs of intracellular edema in animals exposed to both conditions, but this was more marked in the 3 0 ATA group, which showed ultrastructural alterations at the mitochondrial level There was a significant increase in the volume density of intraglomic-congested capillaries in the 3 0 ATA group associated with an arteriolar vasoconstriction In the 2 4 ATA group, there was a relative increase of glomic light cells and a decrease of glomic progenitor cells Additionally, there was a stronger immunoreactivity for F2-isoprostane in the 3 0 ATA O(2)-exposed carotid bodies The glomic cells stained positive for nNOS, but no difference was observed between the groups Our results show that high 02 exposures may induce structural alterations in glomic cells with signs of lipid peroxidation We further suggest that deviation of blood flow toward intraglomic capillaries occurs in hyperbaric hyperoxia

Use of laser phototherapy on a delayed wound healing of oral mucosa previously submitted to radiotherapy: case report

Radiotherapy produces both acute and delayed effects on mucosal tissues, disturbing their healing. This report shows a successful treatment with laser phototherapy (LPT) on a delayed wound healing in oral mucosa previously submitted to radiotherapy with a follow up of 3 years. A 47-year-old patient treated 6 months earlier for tongue squamous cell carcinoma by surgery and radiotherapy presented with a mass in the operated area. Biopsy showed chronic inflammatory infiltrate around a residual polyglactin suture. After 2 months there was a painful mucosal dehiscence on the biopsy site. LPT was performed using a semiconductor laser with 660-nm wavelength (InGaAlP) and spot size of 0.04 cm(2). The parameters applied were 40 mW, 4 Jcm(2)/point, 0.16 J/point, 2.4 J/session. The irradiation was performed punctually, through contact mode in 15 points (4 seconds/point), on top of and around the lesion, during ten sessions. The wound healed completely after ten sessions. This treatment proved to be conservative and effective, inducing healing of a chronic wound in a tissue previously submitted to radiotherapy.

A mathematical model of concentrate solids content for the wet drum magnetic separator

Low concentrate density from wet drum magnetic separators in dense medium circuits can cause operating difficulties due to inability to obtain the required circulating medium density and, indirectly, high medium solids losses. The literature is almost silent on the processes controlling concentrate density. However, the common name for the region through which concentrate is discharged-the squeeze pan gap-implies that some extrusion process is thought to be at work. There is no model of magnetics recovery in a wet drum magnetic separator, which includes as inputs all significant machine and operating variables. A series of trials, in both factorial experiments and in single variable experiments, was done using a purpose built rig which featured a small industrial scale (700 mm lip length, 900 turn diameter) wet drum magnetic separator. A substantial data set of 191 trials was generated in this work. The results of the factorial experiments were used to identify the variables having a significant effect on magnetics recovery. It is proposed, based both on the experimental observations of the present work and on observations reported in the literature, that the process controlling magnetic separator concentrate density is one of drainage. Such a process should be able to be defined by an initial moisture, a drainage rate and a drainage time, the latter being defined by the volumetric flowrate and the volume within the drainage zone. The magnetics can be characterised by an experimentally derived ultimate drainage moisture. A model based on these concepts and containing adjustable parameters was developed. This model was then fitted to a randomly chosen 80% of the data, and validated by application to the remaining 20%. The model is shown to be a good fit to data over concentrate solids content values from 40% solids to 80% solids and for both magnetite and ferrosilicon feeds. (C) 2003 Elsevier Science B.V. All rights reserved.

Peroxisome proliferator-activated receptor-beta as a target for wound healing drugs.

Healing of cutaneous wounds, which is crucial for survival after an injury, proceeds via a well-tuned pattern of events including inflammation, re-epithelialisation, and matrix and tissue remodelling. These events are regulated spatio-temporally by a variety of growth factors and cytokines. The inflammation that immediately follows injury increases the expression of peroxisome proliferator-activated receptor (PPAR)-beta in the wound edge keratinocytes and triggers the production of endogenous PPARbeta ligands that activate the newly produced receptor. This elevated PPARbeta activity results in increased resistance of the keratinocytes to the apoptotic signals released during wounding, allowing faster re-epithelialisation. The authors speculate that, in parallel, ligand activation of PPARbeta in infiltrated macrophages attenuates the inflammatory response, which also promotes repair. Thus, current understanding of the roles of PPARbeta in different cell types implicated in tissue repair has revealed an intriguing intercellular cross-talk that coordinates, spatially and temporally, inflammation, keratinocyte survival, proliferation and migration, which are all essential for efficient wound repair. These novel insights into the orchestrating roles of PPARbeta during wound healing may be helpful in the development of drugs for acute and chronic wound disorders.

Peroxisome proliferator-activated receptors (PPARs): from metabolic control to epidermal wound healing.

Peroxisome proliferator-activated receptors control many cellular and metabolic processes. They are transcription factors belonging to the family of ligand-inducible nuclear receptors. Three isotypes called PPARalpha, PPARbeta/delta and PPARgamma have been identified in lower vertebrates and mammals. They display differential tissue distribution and each of the three isotypes fulfills specific functions. PPARalpha and PPARgamma control energy homoeostasis and inflammatory responses. Their activity can be modulated by drugs such as the hypolipidaemic fibrates and the insulin sensitising thiazolidinediones (pioglitazone and rosiglitazone). Thus, these receptors are involved in the control of chronic diseases such as diabetes, obesity, and atherosclerosis. Little is known about the main function of PPARbeta, but it has been implicated in embryo implantation, tumorigenesis in the colon, reverse cholesterol transport, and recently in skin wound healing. Here, we present recent developments in the PPAR field with particular emphasis on both the function of PPARs in lipid metabolism and energy homoeostasis (PPARalpha and PPARgamma), and their role in epidermal maturation and skin wound repair (PPARalpha and PPARbeta).

Methods for construction and analysis of computational models in systems biology : applications to the modelling of the heat shock response and the self-assembly of intermediate filaments

Systems biology is a new, emerging and rapidly developing, multidisciplinary research field that aims to study biochemical and biological systems from a holistic perspective, with the goal of providing a comprehensive, system- level understanding of cellular behaviour. In this way, it addresses one of the greatest challenges faced by contemporary biology, which is to compre- hend the function of complex biological systems. Systems biology combines various methods that originate from scientific disciplines such as molecu- lar biology, chemistry, engineering sciences, mathematics, computer science and systems theory. Systems biology, unlike “traditional” biology, focuses on high-level concepts such as: network, component, robustness, efficiency, control, regulation, hierarchical design, synchronization, concurrency, and many others. The very terminology of systems biology is “foreign” to “tra- ditional” biology, marks its drastic shift in the research paradigm and it indicates close linkage of systems biology to computer science. One of the basic tools utilized in systems biology is the mathematical modelling of life processes tightly linked to experimental practice. The stud- ies contained in this thesis revolve around a number of challenges commonly encountered in the computational modelling in systems biology. The re- search comprises of the development and application of a broad range of methods originating in the fields of computer science and mathematics for construction and analysis of computational models in systems biology. In particular, the performed research is setup in the context of two biolog- ical phenomena chosen as modelling case studies: 1) the eukaryotic heat shock response and 2) the in vitro self-assembly of intermediate filaments, one of the main constituents of the cytoskeleton. The range of presented approaches spans from heuristic, through numerical and statistical to ana- lytical methods applied in the effort to formally describe and analyse the two biological processes. We notice however, that although applied to cer- tain case studies, the presented methods are not limited to them and can be utilized in the analysis of other biological mechanisms as well as com- plex systems in general. The full range of developed and applied modelling techniques as well as model analysis methodologies constitutes a rich mod- elling framework. Moreover, the presentation of the developed methods, their application to the two case studies and the discussions concerning their potentials and limitations point to the difficulties and challenges one encounters in computational modelling of biological systems. The problems of model identifiability, model comparison, model refinement, model inte- gration and extension, choice of the proper modelling framework and level of abstraction, or the choice of the proper scope of the model run through this thesis.