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.

The calcium-modulated proteins, S100A1 and S100B, as potential regulators of the dynamics of type III intermediate filaments

The Ca2+-modulated, dimeric proteins of the EF-hand (helix-loop-helix) type, S100A1 and S100B, that have been shown to inhibit microtubule (MT) protein assembly and to promote MT disassembly, interact with the type III intermediate filament (IF) subunits, desmin and glial fibrillary acidic protein (GFAP), with a stoichiometry of 2 mol of IF subunit/mol of S100A1 or S100B dimer and an affinity of 0.5-1.0 µM in the presence of a few micromolar concentrations of Ca2+. Binding of S100A1 and S100B results in inhibition of desmin and GFAP assemblies into IFs and stimulation of the disassembly of preformed desmin and GFAP IFs. S100A1 and S100B interact with a stretch of residues in the N-terminal (head) domain of desmin and GFAP, thereby blocking the head-to-tail process of IF elongation. The C-terminal extension of S100A1 (and, likely, S100B) represents a critical part of the site that recognizes desmin and GFAP. S100B is localized to IFs within cells, suggesting that it might have a role in remodeling IFs upon elevation of cytosolic Ca2+ concentration by avoiding excess IF assembly and/or promoting IF disassembly in vivo. S100A1, that is not localized to IFs, might also play a role in the regulation of IF dynamics by binding to and sequestering unassembled IF subunits. Together, these observations suggest that S100A1 and S100B may be regarded as Ca2+-dependent regulators of the state of assembly of two important elements of the cytoskeleton, IFs and MTs, and, potentially, of MT- and IF-based activities.

The role of keratin intermediate filaments in the colon epithelial cells

Keratins (K) are cytoskeletal proteins mainly expressed in the epithelium and constitute the largest subgroup of intermediate filaments (IFs). Simple epithelial keratins (SEKs) K7-K8 and K18-K20 are the major IF elements in the colon. SEK mutations are known to cause around 30 human diseases, mainly affecting liver and skin. However, so far no strong associations between K8 mutations and the development of human colitis have been found. The keratin contribution to colonic health comes from the K8 knock-out (K8-/-) mouse model, which develops an early chronic inflammation and hyperproliferation in the colon. The aim of this thesis was to investigate how keratins contribute to intestinal health and disease mainly by the experimental analysis using the K8-/- mouse colon and cell culture models. The work described here is divided into three studies. The first study revealed involvement of keratins in Notch1 signaling, which is the master regulator of cell fate in the colon. Immunoprecipitation and immunostaining, both in vitro and in vivo showed that K8 binds and co-localizes with Notch1. Interestingly, overexpression of keratins enhanced Notch1 levels and stabilized Notch intracellular domain (NICD), leading to higher activity of Notch signaling. The dramatic decrease in Notch activity in the K8-/- colon resulted in a differentiation shift towards goblet and enteroendocrine cells. The second study focused on the involvement of keratins in colitis-associated cancer (CAC). Although, the K8-/- inflamed colon did not develop colorectal cancer (CRC) spontaneously, it was dramatically more susceptible to induced CRC in two CRC models: azoxymethane (AOM) and multiple intestinal neoplasia (ApcMin/+). To understand how the loss of K8 contributes to CAC, the epithelial inflammasome signaling pathway was analyzed. The released component of active inflammasome, cleaved caspase-1 and its downstream protein, interleukin (IL)-18, were significantly increased in K8-/- and K8-/-ApcMin/+ colons. The inflammasome pathway has recently been suggested to control the levels of IL-22 binding protein (IL-22BP), which is a negative regulator of IL-22 activity. Interestingly, the activated inflammasome correlated with an upregulation of IL-22 and a complete loss of IL-22BP in the K8-null colons. The activation of IL-22 was confirmed by increased levels of downstream signaling, which is phosphorylated signal transducer and activator of transcription 3 (P-STAT3), a transcription factor promoting proliferation and tissue regeneration in the colon. The objective of the third study, was to examine the role of keratins in colon energy metabolism. A proteomic analysis identified mitochondrial 3-hydroxy-3-methylglutaryl-CoA synthase 2 (HMGCS2) as the major ownregulated protein in the K8-/- colonocytes. HMGCS2 is the rate-limiting enzyme in ketogenesis, where energy from bacterially produced short chain fatty acids (SCFAs), mainly butyrate, is converted into ketone bodies in colonic epithelium. Lower levels and activity of HMGCS2 in the K8-/- colon resulted in a blunted ketogenesis. The studies upstream from HMGCS2, identified decreased levels of the SCFA-transporter monocarboxylate transporter 1 (MCT1), which led to increased SCFA content in the stool suggesting impaired butyrate transport through the colonic epithelium. Taken together, the results of the herein thesis indicate that keratins are essential regulators of colon homeostasis, in particular epithelial differentiation, tumorigenesis and energy metabolism.

Instability of surface-temperature filaments in strain and shear

The effects of uniform straining and shearing on the stability of a surface quasi-geostrophic temperature filament are investigated. Straining is shown to stabilize perturbations for wide filaments but only for a finite time until the filament thins to a critical width, after which some perturbations can grow. No filament can be stabilized in practice, since there are perturbations that can grow large for any strain rate. The optimally growing perturbations, defined as solutions that reach a certain threshold amplitude first, are found numerically for a wide range of parameter values. The radii of the vortices formed through nonlinear roll-up are found to be proportional to θ/s, where θ is the temperature anomaly of the filament and s the strain rate, and are not dependent on the initial size of the filament. Shearing is shown to reduce the normal-mode growth rates, but it cannot stabilize them completely when there are temperature discontinuities in the basic state; smooth filaments can be stabilized completely by shearing and a simple scaling argument provides the shear rate required. Copyright © 2010 Royal Meteorological Society

Supramolecular helix and beta-sheet through self-assembly of two isomeric tetrapeptides in crystals and formation of filaments and ribbons in the solid state

Single crystal X-ray diffraction studies show that the beta-turn structure of tetrapeptide I, Boc-Gly-Phe-Aib-Leu-OMe (Aib: alpha-amino isobutyric acid) self-assembles to a supramolecular helix through intermolecular hydrogen bonding along the crystallographic a axis. By contrast the beta-turn structure of an isomeric tetrapeptide II, Boc-Gly-Leu-Aib-Phe-OMe self-assembles to a supramolecular beta-sheet-like structure via a two-dimensional (a, b axis) intermolecular hydrogen bonding network and pi-pi interactions. FT-IR studies of the peptides revealed that both of them form intermolecularly hydrogen bonded supramolecular structures in the solid state. Field emission scanning electron micrographs (FE-SEM) of the dried fibrous materials of the peptides show different morphologies, non-twisted filaments in case of peptide I and non-twisted filaments and ribbon-like structures in case of peptide II.

Conductance of amyloid β based peptide filaments: structure–function relations.

Controlling the morphology of self-assembled peptide nanostructures, particularly those based on amyloid peptides, has been the focus of intense research. In order to exploit these structures in electronic applications, further understanding of their electronic behavior is required. In this work, the role of peptide morphology in determining electronic conduction along self-assembled peptide nanofilament networks is demonstrated. The peptides used in this work were based on the sequence AAKLVFF, which is an extension of a core sequence from the amyloid b peptide. We show that the incorporation of a non-natural amino acid, 2-thienylalanine, instead of phenylalanine improves the obtained conductance with respect to that obtained for a similar structure based on the native sequence, which was not the case for the incorporation of 3-thienylalanine. Furthermore, we demonstrate that the morphology of the self-assembled structures, which can be controlled by the solvent used in the assembly process, strongly affects the conductance, with larger conduction obtained for a morphology of long, straight filaments. Our results demonstrate that, similar to natural systems, the assembly and folding of peptides could be of great importance for optimizing their function as components of electronic devices. Hence, sequence design and assembly conditions can be used to control the performance of peptide based structures in such electronic applications.

Heterogeneous surface expression of EspA translocon filaments by Escherichia coli O157:H7 is controlled at the posttranscriptional level

Type III secretion systems of enteric bacteria enable translocation of effector proteins into host cells. Secreted proteins of verotoxigenic Escherichia coli O157 strains include components of a translocation apparatus, EspA, -B, and -D, as well as "effectors" such as the translocated intimin receptor (Tir) and the mitochondrion-associated protein (Map). This research has investigated the regulation of LEE4 translocon proteins, in particular EspA. EspA filaments could not be detected on the bacterial cell surface when E. coli O157:H7 was cultured in M9 minimal medium but were expressed from only a proportion of the bacterial population when cultured in minimal essential medium modified with 25 mM HEPES. The highest proportions of EspA-filamented bacteria were detected in late exponential phase, after which filaments were lost rapidly from the bacterial cell surface. Our previous research had shown that human and bovine E. coli O157:H7 strains exhibit marked differences in EspD secretion levels. Here it is demonstrated that the proportion of the bacterial population expressing EspA filaments was associated with the level of EspD secretion. The ability of individual bacteria to express EspA filaments was not controlled at the level of LEE1-4 operon transcription, as demonstrated by using both beta-galactosidase and green fluorescent protein (GFP) promoter fusions. All bacteria, whether expressing EspA filaments or not, showed equivalent levels of GFP expression when LEEI-4 translational fusions were used. Despite this, the LEE4-espADB mRNA was more abundant from populations with a high proportion of nonsecreting bacteria (low secretors) than from populations with a high proportion of secreting and therefore filamented bacteria (high secretors). This research demonstrates that while specific environmental conditions are required to induce LEEI-4 expression, a further checkpoint exists before EspA filaments are produced on the bacterial surface and secretion of effector proteins occurs. This checkpoint in E. coli O157:H7 translocon expression is controlled by a posttranscriptional mechanism acting on LEE4-espADB mRNA. The heterogeneity in EspA filamentation could arise from phase-variable expression of regulators that control this posttranscriptional mechanism.

The significance of changes in Mytella falcata (Orbigny, 1842) gill filaments chronically exposed to polluted environments

The impact of pollutants in an organism can be observed by changes in functional complexity at different levels. Bivalve gills are suitable for histopathological analysis because of their structure and function. This study aimed at examining the morphology of Mytella falcata gill filaments from three sites in the Santos estuary (São Paulo, Brazil) with different levels of environmental degradation to identify possible changes in gill structure and discuss the significance of these alterations. For this purpose, histological, histochemical and ultrastructural techniques were used. The filaments of animals from site A (less impacted site) were intact, while in sites B and C, pathological changes were observed, such as: detachment of the epithelium in the intermediate zone, morphological changes of this epithelium, inflammatory process, increase in the number of mucous cells and cell turnover processes. These results suggest that the related changes are an attempt to prevent the entrance of pollutants through gill filaments into the entire organism and that cell turnover is the final way to compensate cell injury. (C) 2008 Elsevier Ltd. All rights reserved

Surface morphology of Diplodon expansus (Kuster, 1856; Mollusca, Bivalvia, Hyriidae) gill filaments after exposure to environmentally relevant concentrations of atrazine herbicide

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Fine structure of Mytella falcata (Bivalvia) gill filaments

Bivalve filter feeders are sessile animals that live in constant contact with water and its pollutants. Their gill is an organ highly exposed to these conditions due to its large surface and its involvement in gas exchanges and feeding. The bivalve Mytella falcata is found in estuaries of Latin America, on the Atlantic as well as the Pacific Coast. It is commonly consumed, and sometimes is the only source of protein of low-income communities. In this study, gill filaments of M. falcata were characterized using histology, histochemistry and transmission electron microscopy for future comparative studies among animals exposed to environmental pollutants. Gill filaments may be divided into abfrontal, intermediate and frontal zones. Filaments are interconnected by ciliary discs. In the center of filaments, haemocytes circulate through a haemolymph vessel internally lined by an endothelium and supported by an acellular connective tissue rich in polysaccharides and collagen. The abfrontal zone contains cuboidal cells, while the intermediate zone consists of a simple squamous epithelium. The frontal zone is composed of five columnar cell types: one absorptive, mainly characterized by the presence of pinocytic vesicles in the apical region of the cell; one secretory, rarely observed and three ciliated with abundant mitochondria. All cells lining the filament exhibit numerous microvilli and seem to absorb substances from the environment. PAS staining was observed in mucous cells in the frontal and abfrontal zones. Bromophenol blue allowed the distinction of haemocytes and detection of a glycoprotein secretion in the secretory cells of the frontal region. The characteristics of M. falcata gill filaments observed in this study were very similar to those of other bivalves, especially other Mytilidae, and are suitable for histopathological studies on the effect of water-soluble pollutants. (C) 2007 Elsevier Ltd. All rights reserved.

The pre-emptive effect of epidural ketamine on wound sensitivity in horses tested by using von Frey filaments

Objective To evaluate the pre-emptive analgesic effect of pre-incisional epidural ketamine.Study Design A blinded, randomized experimental study.Animals Sixteen mixed breed mares, 17.6 +/- 2.8 years old, weighing 352 +/- 32 kg.Methods In a pilot study, an incision was made on one lateral thigh using a lidocaine block and no further analgesics, and it was verified that the nociceptive threshold was lower on the incised side than nonincised side (p < 0.05), and that von Frey filaments evoked a pain response. The 16 animals were divided into group A (ketamine, n = 9) and B (saline, n = 7). An epidural catheter was inserted 24 hours before the trials, the thigh was shaved bilaterally, and the right side was blocked (incised side) using lidocaine. Twenty-five minutes later, ketamine (A) or saline (B) was administered epidurally. Five minutes later, a 10-cm. skin incision was made on the right side, and then sutured. Nociceptive threshold was determined with von Frey filaments at 1, 3, and 5 cm. around the incision at 15-minute intervals for 2 hours, then at 4, 6, and 8 hours. Behavioral alterations, heart and respiratory rates were recorded. Nociceptive thresholds from these points were averaged to obtain mean values at each time, converted to a logarithmic scale, and submitted to a nonparametric analysis (Mann-Whitney and one-way repeated measures anova test,p less than or equal to 0.05).Results After 8 hours, the global range score revealed reduced hyperalgesia (p < 0.01) around the incision in 92% (4.65-4.27) of evaluated intervals in group A (ketamine). There were no significant changes in behavior, heart and respiratory rates,Conclusions It was concluded that pre-emptive epidural ketamine reduced post-incisional pain in the horse, and that von Frey filaments were able to quantify cutaneous sensitivity after tissue damage.Clinical relevance Epidural ketamine injection can reduce post-incisional sensitivity in the horse.

Model for electrodes' filaments of hot cathode fluorescent lamps, during preheating with constant rms current

This paper presents a new model for the representation of electrodes' filaments of hot-cathode fluorescent lamps, during preheating processes based on the injection of currents with constant root mean square (rms) values. The main improvement obtained with this model is the prediction of the R-h/R-c ratio during the preheating process, as a function of the preheating time and of the rms current injected in the electrodes. Using the proposed model, it is possible to obtain an estimate of the time interval and the current that should be provided by the electronic ballast, in order to ensure a suitable preheating process. is estimate of time and current can be used as input data in the design of electronic ballasts with programmed lamp start, permitting the prediction of the R-h/R-c ratio during the initial steps of the design (theoretical analysis and digital simulation). Therefore, the use of the proposed model permits to reduce the necessity of several empirical adjustments in the prototype, in order to set the operation of electronic ballasts during the preheating process. This fact reduces time and costs associated to the global design procedure of electronic ballasts.