Quantitative refinement of reaction-based biomodels

In the field of molecular biology, scientists adopted for decades a reductionist perspective in their inquiries, being predominantly concerned with the intricate mechanistic details of subcellular regulatory systems. However, integrative thinking was still applied at a smaller scale in molecular biology to understand the underlying processes of cellular behaviour for at least half a century. It was not until the genomic revolution at the end of the previous century that we required model building to account for systemic properties of cellular activity. Our system-level understanding of cellular function is to this day hindered by drastic limitations in our capability of predicting cellular behaviour to reflect system dynamics and system structures. To this end, systems biology aims for a system-level understanding of functional intraand inter-cellular activity. Modern biology brings about a high volume of data, whose comprehension we cannot even aim for in the absence of computational support. Computational modelling, hence, bridges modern biology to computer science, enabling a number of assets, which prove to be invaluable in the analysis of complex biological systems, such as: a rigorous characterization of the system structure, simulation techniques, perturbations analysis, etc. Computational biomodels augmented in size considerably in the past years, major contributions being made towards the simulation and analysis of large-scale models, starting with signalling pathways and culminating with whole-cell models, tissue-level models, organ models and full-scale patient models. The simulation and analysis of models of such complexity very often requires, in fact, the integration of various sub-models, entwined at different levels of resolution and whose organization spans over several levels of hierarchy. This thesis revolves around the concept of quantitative model refinement in relation to the process of model building in computational systems biology. The thesis proposes a sound computational framework for the stepwise augmentation of a biomodel. One starts with an abstract, high-level representation of a biological phenomenon, which is materialised into an initial model that is validated against a set of existing data. Consequently, the model is refined to include more details regarding its species and/or reactions. The framework is employed in the development of two models, one for the heat shock response in eukaryotes and the second for the ErbB signalling pathway. The thesis spans over several formalisms used in computational systems biology, inherently quantitative: reaction-network models, rule-based models and Petri net models, as well as a recent formalism intrinsically qualitative: reaction systems. The choice of modelling formalism is, however, determined by the nature of the question the modeler aims to answer. Quantitative model refinement turns out to be not only essential in the model development cycle, but also beneficial for the compilation of large-scale models, whose development requires the integration of several sub-models across various levels of resolution and underlying formal representations.

Proteomic analysis of cytosolic proteins associated with petite mutations in Candida glabrata

The incidence of superficial or deep-seated infections due to Candida glabrata has increased markedly, probably because of the low intrinsic susceptibility of this microorganism to azole antifungals and its relatively high propensity to acquire azole resistance. To determine changes in the C. glabrata proteome associated with petite mutations, cytosolic extracts from an azole-resistant petite mutant of C. glabrata induced by exposure to ethidium bromide, and from its azole-susceptible parent isolate were compared by two-dimensional polyacrylamide gel electrophoresis. Proteins of interest were identified by peptide mass fingerprinting or sequence tagging using a matrix-assisted laser desorption/ionization tandem time-of-flight mass spectrometer. Tryptic peptides from a total of 160 Coomassie-positive spots were analyzed for each strain. Sixty-five different proteins were identified in the cytosolic extracts of the parent strain and 58 in the petite mutant. Among the proteins identified, 10 were higher in the mutant strain, whereas 23 were lower compared to the parent strain. The results revealed a significant decrease in the enzymes associated with the metabolic rate of mutant cells such as aconitase, transaldolase, and pyruvate kinase, and changes in the levels of specific heat shock proteins. Moreover, transketolase, aconitase and catalase activity measurements decreased significantly in the ethidium bromide-induced petite mutant. These data may be useful for designing experiments to obtain a better understanding of the nuclear response to impairment of mitochondrial function associated with this mutation in C. glabrata.

Reduction of carcass weight loss in swines

The carcass fast freezing is one of the aspects of great prominence to the final quality of pork. In order to reduce weight loss, two experiments were performed, in which the carcasses were monitored during 20 hours to evaluate the main variables involved during two different freezing processes (standard and proposed) as follows: microbiological quality, storage temperature, relative humidity (RH) and air velocity. In experiment I, the carcasses were submitted to a system using heat shock (2 hours in static tunnel at - 25 °C) and subsequently sent to the equalization chamber. In experiment II, the carcasses were submitted to the heat shock and stored in a chamber with RH between 80-85%. The chambers used in both experiments showed no change in the variables studied (internal temperature of 5 °C and air velocity of approximately 0.3 m/s). However, the relative humidity in the three chambers was evaluated and significant differences were found; as a consequence, high levels of weight loss were observed in both chambers In experiment II there was an increase of RH, which reduced the weight loss of the carcasses.

Computational models for and from biology : simple gene assembly and reaction systems

The advancement of science and technology makes it clear that no single perspective is any longer sufficient to describe the true nature of any phenomenon. That is why the interdisciplinary research is gaining more attention overtime. An excellent example of this type of research is natural computing which stands on the borderline between biology and computer science. The contribution of research done in natural computing is twofold: on one hand, it sheds light into how nature works and how it processes information and, on the other hand, it provides some guidelines on how to design bio-inspired technologies. The first direction in this thesis focuses on a nature-inspired process called gene assembly in ciliates. The second one studies reaction systems, as a modeling framework with its rationale built upon the biochemical interactions happening within a cell. The process of gene assembly in ciliates has attracted a lot of attention as a research topic in the past 15 years. Two main modelling frameworks have been initially proposed in the end of 1990s to capture ciliates’ gene assembly process, namely the intermolecular model and the intramolecular model. They were followed by other model proposals such as templatebased assembly and DNA rearrangement pathways recombination models. In this thesis we are interested in a variation of the intramolecular model called simple gene assembly model, which focuses on the simplest possible folds in the assembly process. We propose a new framework called directed overlap-inclusion (DOI) graphs to overcome the limitations that previously introduced models faced in capturing all the combinatorial details of the simple gene assembly process. We investigate a number of combinatorial properties of these graphs, including a necessary property in terms of forbidden induced subgraphs. We also introduce DOI graph-based rewriting rules that capture all the operations of the simple gene assembly model and prove that they are equivalent to the string-based formalization of the model. Reaction systems (RS) is another nature-inspired modeling framework that is studied in this thesis. Reaction systems’ rationale is based upon two main regulation mechanisms, facilitation and inhibition, which control the interactions between biochemical reactions. Reaction systems is a complementary modeling framework to traditional quantitative frameworks, focusing on explicit cause-effect relationships between reactions. The explicit formulation of facilitation and inhibition mechanisms behind reactions, as well as the focus on interactions between reactions (rather than dynamics of concentrations) makes their applicability potentially wide and useful beyond biological case studies. In this thesis, we construct a reaction system model corresponding to the heat shock response mechanism based on a novel concept of dominance graph that captures the competition on resources in the ODE model. We also introduce for RS various concepts inspired by biology, e.g., mass conservation, steady state, periodicity, etc., to do model checking of the reaction systems based models. We prove that the complexity of the decision problems related to these properties varies from P to NP- and coNP-complete to PSPACE-complete. We further focus on the mass conservation relation in an RS and introduce the conservation dependency graph to capture the relation between the species and also propose an algorithm to list the conserved sets of a given reaction system.

Selecting Relevant Genes with a Spectral Approach

Array technologies have made it possible to record simultaneously the expression pattern of thousands of genes. A fundamental problem in the analysis of gene expression data is the identification of highly relevant genes that either discriminate between phenotypic labels or are important with respect to the cellular process studied in the experiment: for example cell cycle or heat shock in yeast experiments, chemical or genetic perturbations of mammalian cell lines, and genes involved in class discovery for human tumors. In this paper we focus on the task of unsupervised gene selection. The problem of selecting a small subset of genes is particularly challenging as the datasets involved are typically characterized by a very small sample size ?? the order of few tens of tissue samples ??d by a very large feature space as the number of genes tend to be in the high thousands. We propose a model independent approach which scores candidate gene selections using spectral properties of the candidate affinity matrix. The algorithm is very straightforward to implement yet contains a number of remarkable properties which guarantee consistent sparse selections. To illustrate the value of our approach we applied our algorithm on five different datasets. The first consists of time course data from four well studied Hematopoietic cell lines (HL-60, Jurkat, NB4, and U937). The other four datasets include three well studied treatment outcomes (large cell lymphoma, childhood medulloblastomas, breast tumors) and one unpublished dataset (lymph status). We compared our approach both with other unsupervised methods (SOM,PCA,GS) and with supervised methods (SNR,RMB,RFE). The results clearly show that our approach considerably outperforms all the other unsupervised approaches in our study, is competitive with supervised methods and in some case even outperforms supervised approaches.

A Comparative study of germination strategies of two species of genus Allium sect. Allium

Germination experiments were performed with seeds of two species of genus Allium section Allium, a rare and endangered species A. pyrenaicum and a common A. sphaerocephalon. Different pre-treatments and a photoperiod of 24 h darkness were applied in order to simulate different germination conditions. Both species showed a high percentage of viable seeds a part of which were dormant. An elevate percentage of dormant seeds could be caused by a later collection time. Low altitude populations had more mortality than the others, possibly caused by the hard summer conditions during flowering and fruiting time. Comparisons between dates of species coexistence localities only show inter-population variability and it could be caused by the detected dormancy. Darkness accelerates germination, possibly for elongation radicle stimulation. Heat-shock pre-treatments decreased germination time in seeds from localities where fire is a probable event. The rarity of A. Pyrenaicum not seems to be caused by restricted germination requirements but is attributable to distinct habitat preferences, related to his altitudinal range of distribution

The adsorbed conformation of globular proteins at the air/water interface

External reflection FTIR spectroscopy and surface pressure measurements were used to compare conformational changes in the adsorbed structures of three globular proteins at the air/water interface. Of the three proteins studied, lysozyme, bovine serum albumin and P-lactoglobulin, lysozyme was unique in its behaviour. Lysozyme adsorption was slow, taking approximately 2.5 h to reach a surface pressure plateau (from a 0.07 mM solution), and led to significant structural change. The FTIR spectra revealed that lysozyme formed a highly networked adsorbed layer of unfolded protein with high antiparallel beta-sheet content and that these changes occurred rapidly (within 10 min). This non-native secondary structure is analogous to that of a 3D heat-set protein gel, suggesting that the adsorbed protein formed a highly networked interfacial layer. Albumin and P-lactoglobulin adsorbed rapidly (reaching a plateau within 10 min) and with little chance to their native secondary structure.

Contingency locus in ctsR of Listeria monocytogenes Scott A: a strategy for occurrence of abundant piezotolerant isolates within clonal populations

In a recent study we demonstrated that a high-hydrostatic-pressure-tolerant isolate of Listeria monocytogenes lacks a codon in the class 3 heat shock regulator gene ctsR. This mutation in the region that encodes four consecutive glycines was directly responsible for the observed piezotolerance, increased stress resistance, and reduced virulence. The aim of the present study was to determine whether mutations in ctsR are frequently associated with piezotolerance in L. monocytogenes. Wild-type cultures of L. monocytogenes were therefore exposed to 350 MPa for 20 min, and the piezotolerance of individual surviving isolates was assessed. This rendered 33 isolates with a stable piezotolerant phenotype from a total of 84 survivors. Stable piezotolerant mutants were estimated to be present in the initial wild-type population at frequencies of >10�5. Subsequent sequencing of the ctsR gene of all stable piezotolerant isolates revealed that two-thirds of the strains (i.e., n � 21) had mutations in this gene. The majority of the mutations (16 of 21 strains) consisted of a triplet deletion in the glycine-encoding region of ctsR, identical to what was found in our previous study. Interestingly, 2 of 21 mutants contained a codon insertion in this repeat region. The remaining three stable piezotolerant strains showed a 19-bp insertion in the glycine repeat region, a 16-bp insertion downstream of the glycine repeat area (both leading to frameshifts and a truncated ctsR), and an in-frame 114-bp deletion encoding a drastically shortened carboxy terminus of CtsR. In four instances it was not possible to generate a PCR product. A piezotolerant phenotype could not be linked to mutations in ctsR in 8 of 33 isolates, indicating that other thus-far-unknown mechanisms also lead to stable piezotolerance. The present study highlights the importance of ctsR in piezotolerance and stress tolerance of L. monocytogenes, and it demonstrates that short-sequence repeat regions contribute significantly to the occurrence of a piezotolerant and stress-tolerant subpopulation within L. monocytogenes cultures, thus playing an important role in survival.

Optimisation of recombinant production of active human cardiac SERCA2a ATPase

Methods for recombinant production of eukaryotic membrane proteins, yielding sufficient quantity and quality of protein for structural biology, remain a challenge. We describe here, expression and purification optimisation of the human SERCA2a cardiac isoform of Ca2+ translocating ATPase, using Saccharomyces cerevisiae as the heterologous expression system of choice. Two different expression vectors were utilised, allowing expression of C-terminal fusion proteins with a biotinylation domain or a GFP- His8 tag. Solubilised membrane fractions containing the protein of interest were purified onto Streptavidin-Sepharose, Ni-NTA or Talon resin, depending on the fusion tag present. Biotinylated protein was detected using specific antibody directed against SERCA2 and, advantageously, GFP-His8 fusion protein was easily traced during the purification steps using in-gel fluorescence. Importantly, talon resin affinity purification proved more specific than Ni-NTA resin for the GFP-His8 tagged protein, providing better separation of oligomers present, during size exclusion chromatography. The optimised method for expression and purification of human cardiac SERCA2a reported herein, yields purified protein (> 90%) that displays a calcium-dependent thapsigargin-sensitive activity and is suitable for further biophysical, structural and physiological studies. This work provides support for the use of Saccharomyces cerevisiae as a suitable expression system for recombinant production of multi-domain eukaryotic membrane proteins.

Changes in gene expression induced by H(2)O(2) in cardiac myocytes.

Oxidative stress induces cardiac myocyte apoptosis. At least some effects are probably mediated through changes in gene expression. Using Affymetrix arrays, we examined the changes in gene expression induced by H(2)O(2) (0.04, 0.1, and 0.2mM; 2 and 4h) in rat neonatal ventricular myocytes. Changes in selected upregulated genes were confirmed by ratiometric RT-PCR. p21(Cip1/Waf1) was one of the only two genes upregulated in all conditions studied. Of the heat shock proteins, only Hsp70/70.1 was induced by H(2)O(2) with no change in the expression of Hsp25, Hsp60 or Hsp90. Heme oxygenase 1 was also potently upregulated, but not heme oxygenases 2 or 3. Of the intercellular adhesion proteins, syndecan-1 was significantly upregulated in response to H(2)O(2), with little change in the expression of other syndecans and no change in expression of any of the integrins studied. Thus, oxidative stress, exemplified by H(2)O(2), selectively promotes the expression of specific gene family members.

Characterization of Caulobacter crescentus response to low temperature and identification of genes involved in freezing resistance

Free-living bacteria must respond to a wide range of temperature changes, and have developed specific mechanisms to survive in extreme environments. In this work we describe a remarkable resistance of mesophilic bacterium Caulobacter crescentus to several cycles of freezing at -80 degrees C, which was able to grow at low temperatures. Exponentially growing cells and late stationary-phase cells presented higher freezing resistance at both -20 and -80 degrees C than early stationary-phase cells. Cryotolerance was observed when log-phase cultures grown at 30 degrees C were preincubated at 5, 15 or 20 degrees C before freezing at -20 degrees C. A transposon library was screened to identify mutants sensitive to freezing at -80 degrees C and three strains presenting < 10% survival were isolated. Identification of genes disrupted in each mutant showed that they encoded an AddA family DNA helicase, a DEAD/DEAH box RNA helicase and a putative RND (resistance, nodulation, cell division) efflux system component. These strains showed longer generation times than wild-type cells when growing at 15 degrees C, with the RNA helicase mutant presenting a severe growth defect. These analyses suggest that the singular intrinsic resistance to freezing of C. crescentus is in fact a consequence of several independent traits, especially the maintenance of a proper degree of supercoiling of nucleic acids.

CspC and CspD are essential for Caulobacter crescentus stationary phase survival

The cold shock response in bacteria involves the expression of low-molecular weight cold shock proteins (CSPs) containing a nucleic acid-binding cold shock domain (CSD), which are known to destabilize secondary structures on mRNAs, facilitating translation at low temperatures. Caulobacter crescentus cspA and cspB are induced upon cold shock, while cspC and cspD are induced during stationary phase. In this work, we determined a new coding sequence for the cspC gene, revealing that it encodes a protein containing two CSDs. The phenotypes of C. crescentus csp mutants were analyzed, and we found that cspC is important for cells to maintain viability during extended periods in stationary phase. Also, cspC and cspCD strains presented altered morphology, with frequent non-viable filamentous cells, and cspCD also showed a pronounced cell death at late stationary phase. In contrast, the cspAB mutant presented increased viability in this phase, which is accompanied by an altered expression of both cspC and cspD, but the triple cspABD mutant loses this characteristic. Taken together, our results suggest that there is a hierarchy of importance among the csp genes regarding stationary phase viability, which is probably achieved by a fine tune balance of the levels of these proteins.

SpdR, a Response Regulator Required for Stationary-Phase Induction of Caulobacter crescentus cspD

The cold shock protein (CSP) family includes small polypeptides that are induced upon temperature downshift and stationary phase. The genome of the alphaproteobacterium Caulobacter crescentus encodes four CSPs, with two being induced by cold shock and two at the onset of stationary phase. In order to identify the environmental signals and cell factors that are involved in cspD expression at stationary phase, we have analyzed cspD transcription during growth under several nutrient conditions. The results showed that expression of cspD was affected by the medium composition and was inversely proportional to the growth rate. The maximum levels of expression were decreased in a spoT mutant, indicating that ppGpp may be involved in the signalization for carbon starvation induction of cspD. A Tn5 mutant library was screened for mutants with reduced cspD expression, and 10 clones that showed at least a 50% reduction in expression were identified. Among these, a strain with a transposon insertion into a response regulator of a two-component system showed no induction of cspD at stationary phase. This protein (SpdR) was able to acquire a phosphate group from its cognate histidine kinase, and gel mobility shift assay and DNase I footprinting experiments showed that it binds to an inverted repeat sequence of the cspD regulatory region. A mutated SpdR with a substitution of the conserved aspartyl residue that is the probable phosphorylation site is unable to bind to the cspD regulatory region and to complement the spdR mutant phenotype.

The iron stimulon of Xylella fastidiosa includes genes for type IV pilus and colicin V-like bacteriocins

Xylella fastidiosa is the etiologic agent of a wide range of plant diseases, including citrus variegated chlorosis (CVC), a major threat to citrus industry. The genomes of several strains of this phytopathogen were completely sequenced, enabling large-scale functional studies. DNA microarrays representing 2,608 (91.6%) coding sequences (CDS) of X. fastidiosa CVC strain 9a5c were used to investigate transcript levels during growth with different iron availabilities. When treated with the iron chelator 2,2`-dipyridyl, 193 CDS were considered up-regulated and 216 were considered down-regulated. Upon incubation with 100 mu M ferric pyrophosphate, 218 and 256 CDS were considered up- and down-regulated, respectively. Differential expression for a subset of 44 CDS was further evaluated by reverse transcription-quantitative PCR. Several CDS involved with regulatory functions, pathogenicity, and cell structure were modulated under both conditions assayed, suggesting that major changes in cell architecture and metabolism occur when X. fastidiosa cells are exposed to extreme variations in iron concentration. Interestingly, the modulated CDS include those related to colicin V-like bacteriocin synthesis and secretion and to functions of pili/fimbriae. We also investigated the contribution of the ferric uptake regulator Fur to the iron stimulon of X. fastidiosa. The promoter regions of the strain 9a5c genome were screened for putative Fur boxes, and candidates were analyzed by electrophoretic mobility shift assays. Taken together, our data support the hypothesis that Fur is not solely responsible for the modulation of the iron stimulon of X fastidiosa, and they present novel evidence for iron regulation of pathogenicity determinants.

Xylella fastidiosa: An in vivo system to study possible survival strategies within citrus xylem vessels based on global gene expression analysis

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)