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 use of protein structure/activity relationships in the rational design of stable particulate delivery systems

The recombinant heat shock protein (18 kDa-hsp) from Mycobacterium leprae was studied as a T-epitope model for vaccine development. We present a structural analysis of the stability of recombinant 18 kDa-hsp during different processing steps. Circular dichroism and ELISA were used to monitor protein structure after thermal stress, lyophilization and chemical modification. We observed that the 18 kDa-hsp is extremely resistant to a wide range of temperatures (60% of activity is retained at 80ºC for 20 min). N-Acylation increased its ordered structure by 4% and decreased its ß-T1 structure by 2%. ELISA demonstrated that the native conformation of the 18 kDa-hsp was preserved after hydrophobic modification by acylation. The recombinant 18 kDa-hsp resists to a wide range of temperatures and chemical modifications without loss of its main characteristic, which is to be a source of T epitopes. This resistance is probably directly related to its lack of organization at the level of tertiary and secondary structures.

Anabolic steroid- and exercise-induced cardiac stress protein (HSP72) in the rat

The present study investigated the effects of exercise and anabolic-androgenic steroids on cardiac HSP72 expression. Male Wistar rats were divided into experimental groups: nandrolone exercise (NE, N = 6), control exercise (CE, N = 6), nandrolone sedentary (NS, N = 6), and control sedentary (CS, N = 6). Animals in the NE and NS groups received a weekly intramuscular injection (6.5 mg/kg of body weight) of nandrolone decanoate, while those in the CS and CE groups received mineral oil as vehicle. Animals in the NE and CE groups were submitted to a progressive running program on a treadmill, for 8 weeks. Fragments of the left ventricle were collected at sacrifice and the relative immunoblot contents of HSP72 were determined. Heart weight to body weight ratio was higher in exercised than in sedentary animals (P < 0.05, 4.65 ± 0.38 vs 4.20 ± 0.47 mg/g, respectively), independently of nandrolone, and in nandrolone-treated than untreated animals (P < 0.05, 4.68 ± 0.47 vs 4.18 ± 0.32 mg/g, respectively), independently of exercise. Cardiac HSP72 accumulation was higher in exercised than in sedentary animals (P < 0.05, 677.16 ± 129.14 vs 246.24 ± 46.30 relative unit, respectively), independently of nandrolone, but not different between nandrolone-treated and untreated animals (P > 0.05, 560.88 ± 127.53 vs 362.52 ± 95.97 relative unit, respectively) independently of exercise. Exercise-induced HSP72 expression was not affected by nandrolone. These levels of HSP72 expression in response to nandrolone administration suggest either a low intracellular stress or a possible less protection to the myocardium.

Identification of the divergent calmodulin binding motif in yeast Ssb1/Hsp75 protein and in other HSP70 family members

Yeast soluble proteins were fractionated by calmodulin-agarose affinity chromatography and the Ca2+/calmodulin-binding proteins were analyzed by SDS-PAGE. One prominent protein of 66 kDa was excised from the gel, digested with trypsin and the masses of the resultant fragments were determined by MALDI/MS. Twenty-one of 38 monoisotopic peptide masses obtained after tryptic digestion were matched to the heat shock protein Ssb1/Hsp75, covering 37% of its sequence. Computational analysis of the primary structure of Ssb1/Hsp75 identified a unique potential amphipathic alpha-helix in its N-terminal ATPase domain with features of target regions for Ca2+/calmodulin binding. This region, which shares 89% similarity to the experimentally determined calmodulin-binding domain from mouse, Hsc70, is conserved in near half of the 113 members of the HSP70 family investigated, from yeast to plant and animals. Based on the sequence of this region, phylogenetic analysis grouped the HSP70s in three distinct branches. Two of them comprise the non-calmodulin binding Hsp70s BIP/GR78, a subfamily of eukaryotic HSP70 localized in the endoplasmic reticulum, and DnaK, a subfamily of prokaryotic HSP70. A third heterogeneous group is formed by eukaryotic cytosolic HSP70s containing the new calmodulin-binding motif and other cytosolic HSP70s whose sequences do not conform to those conserved motif, indicating that not all eukaryotic cytosolic Hsp70s are target for calmodulin regulation. Furthermore, the calmodulin-binding domain found in eukaryotic HSP70s is also the target for binding of Bag-1 - an enhancer of ADP/ATP exchange activity of Hsp70s. A model in which calmodulin displaces Bag-1 and modulates Ssb1/Hsp75 chaperone activity is discussed.

Heat-shock response in Arabidopsis thaliana explored by multiplexed quantitative proteomics using differential metabolic labeling

We have developed a general method for multiplexed quantitative proteomics using differential metabolic stable isotope labeling and mass spectrometry. The method was successfully used to study the dynamics of heat-shock response in Arabidopsis thaliana. A number of known heat-shock proteins were confirmed, and some proteins not previously associated with heat shock were discovered. The method is applicable in stable isotope labeling and allows for high degrees of multiplexing.

Physiological and proteomic approaches to address heat tolerance during anthesis in rice (Oryza sativa L.)

Episodes of high temperature at anthesis, which in rice is the most sensitive stage to temperature, are expected to occur more frequently in future climates. The morphology of the reproductive organs and pollen number, and changes in anther protein expression, were studied in response to high temperature at anthesis in three rice (Oryza sativa L.) genotypes. Plants were exposed to 6 h of high (38 °C) and control (29 °C) temperature at anthesis and spikelets collected for morphological and proteomic analysis. Moroberekan was the most heat-sensitive genotype (18% spikelet fertility at 38 °C), while IR64 (48%) and N22 (71%) were moderately and highly heat tolerant, respectively. There were significant differences among the genotypes in anther length and width, apical and basal pore lengths, apical pore area, and stigma and pistil length. Temperature also affected some of these traits, increasing anther pore size and reducing stigma length. Nonetheless, variation in the number of pollen on the stigma could not be related to measured morphological traits. Variation in spikelet fertility was highly correlated (r=0.97, n=6) with the proportion of spikelets with ≥20 germinated pollen grains on the stigma. A 2D-gel electrophoresis showed 46 protein spots changing in abundance, of which 13 differentially expressed protein spots were analysed by MS/MALDI-TOF. A cold and a heat shock protein were found significantly up-regulated in N22, and this may have contributed to the greater heat tolerance of N22. The role of differentially expressed proteins and morphology during anther dehiscence and pollination in shaping heat tolerance and susceptibility is discussed.

GTP binding controls complex formation by the human ROCO protein MASL1

The human ROCO proteins are a family of multi-domain proteins sharing a conserved ROC-COR supra-domain. The family has four members: leu- cine-rich repeat kinase 1 (LRRK1), leucine-rich repeat kinase 2 (LRRK2), death-associated protein kinase 1 (DAPK1) and malignant fibrous histiocy- toma amplified sequences with leucine-rich tandem repeats 1 (MASL1). Previous studies of LRRK1/2 and DAPK1 have shown that the ROC (Ras of complex proteins) domain can bind and hydrolyse GTP, but the cellular consequences of this activity are still unclear. Here, the first biochemical characterization of MASL1 and the impact of GTP binding on MASL1 complex formation are reported. The results demonstrate that MASL1, similar to other ROCO proteins, can bind guanosine nucleotides via its ROC domain. Furthermore, MASL1 exists in two distinct cellular com- plexes associated with heat shock protein 60, and the formation of a low molecular weight pool of MASL1 is modulated by GTP binding. Finally, loss of GTP enhances MASL1 toxicity in cells. Taken together, these data point to a central role for the ROC/GTPase domain of MASL1 in the reg- ulation of its cellular function.

Stimulation of the p38 mitogen-activated protein kinase pathway in neonatal rat ventricular myocytes by the G protein-coupled receptor agonists, endothelin-1 and phenylephrine: a role in cardiac myocyte hypertrophy?

We examined the activation of the p38 mitogen-activated protein kinase (p38-MAPK) pathway by the G protein-coupled receptor agonists, endothelin-1 and phenylephrine in primary cultures of cardiac myocytes from neonatal rat hearts. Both agonists increased the phosphorylation (activation) of p38-MAPK by approximately 12-fold. A p38-MAPK substrate, MAPK-activated protein kinase 2 (MAPKAPK2), was activated approximately fourfold and 10 microM SB203580, a p38-MAPK inhibitor, abolished this activation. Phosphorylation of the MAPKAPK2 substrate, heat shock protein 25/27, was also increased. Using selective inhibitors, activation of the p38-MAPK pathway by endothelin-1 was shown to involve protein kinase C but not Gi/Go nor the extracellularly responsive kinase (ERK) pathway. SB203580 failed to inhibit the morphological changes associated with cardiac myocyte hypertrophy induced by endothelin-1 or phenylephrine between 4 and 24 h. However, it decreased the myofibrillar organization and cell profile at 48 h. In contrast, inhibition of the ERK cascade with PD98059 prevented the increase in myofibrillar organization but not cell profile. These data are not consistent with a role for the p38-MAPK pathway in the immediate induction of the morphological changes of hypertrophy but suggest that it may be necessary over a longer period to maintain the response.

Effect of heat exposure on the thermoregulatory responses of selected naked neck chickens

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

Heat shock proteins (HSP): dermatological implications and perspectives

In recent years, several studies have demonstrated the protective effect of Heat Shock Proteins (HSP) on different organs and tissues under stressful conditions. However, most research explores the performance of those molecular chaperones during immune responses or pathological conditions like cancer, whereas the number of studies related to the performance of HSPs in the skin during diverse natural or physiopathological conditions is very low. Therefore, the aim of this article was to summarize the main concepts concerning the expression and performance of HSPs, from analysis of current medicine and cosmetics publications, as well as exploring the importance of these proteins in the dermatological area in physiological events such as cutaneous aging, skin cancer and wound healing and to present final considerations related to biotechnology performance in this area.

Proline content and protein patterns in Eucalyptus grandis shoots submitted to high and low temperature shocks

As respostas às mudanças de temperatura de plantas aclimatadas e não aclimatadas de E. grandis cultivadas in vitro foram avaliadas considerando alterações dos níveis de prolina e proteínas solúveis totais. Análises de proteínas solúveis através de SDS-PAGE e prolina foram realizadas após 12h a 12ºC (aclimatação ao frio) ou a 33ºC (aclimatação ao calor), e imediatamente depois dos choques térmicos a 41ºC e 0ºC. Análises também foram realizadas após um período de 24h depois dos choques térmicos (período de recuperação). O tratamento de temperatura a 0ºC não alterou o padrão de proteínas nas plantas aclimatadas e não aclimatadas, entretanto a temperatura baixa induziu altos níveis de prolina, que se mantiveram relativamente altos após o período de recuperação. Três novas proteínas (90,5, 75 e 39 kDa), provavelmente HSPs, foram observadas nas plantas aclimatadas e não aclimatadas submetidas às temperaturas altas. As plantas expostas a 41ºC foram capazes de recuperar-se dos choques após o período de recuperação, entretanto não houve recuperação completa das plantas expostas às baixas temperaturas. O efeito da aclimatação sobre a recuperação (homeostasis) pode variar dependendo do parâmetro avaliado, tipo e duração do choque térmico.

Genomic organization of the Neurospora crassa gsn gene: possible involvement of the STRE and HSE elements in the modulation of transcription during heat shock

Glycogen synthase, an enzyme involved in glycogen biosynthesis, is regulated by phosphorylation and by the allosteric ligand glucose-6-phosphate (G6P). In addition, enzyme levels can be regulated by changes in gene expression. We recently cloned a cDNA for glycogen synthase (gsn) from Neurospora crassa, and showed that gsn transcription decreased when cells were exposed to heat shock (shifted from 30degreesC to 45degreesC). In order to understand the mechanisms that control gsn expression, we isolated the gene, including its 5' and 3' flanking regions, from the genome of N. crassa. An ORF of approximately 2.4 kb was identified, which is interrupted by four small introns (II-V). Intron I (482 bp) is located in the 5'UTR region. Three putative Transcription Initiation Sites (TISs) were mapped, one of which lies downstream of a canonical TATA-box sequence (5'-TGTATAAA-3'). Analysis of the 5'-flanking region revealed the presence of putative transcription factor-binding sites, including Heat Shock Elements (HSEs) and STress Responsive Elements (STREs). The possible involvement of these motifs in the negative regulation of gsn transcription was investigated using Electrophoretic Mobility Shift Assays (EMSA) with nuclear extracts of N. crassa mycelium obtained before and after heat shock, and DNA fragments encompassing HSE and STRE elements from the 5'-flanking region. While elements within the promoter region are involved in transcription under heat shock, elements in the 5'UTR intron may participate in transcription during vegetative growth. The results thus suggest that N. crassa possesses trans-acting elements that interact with the 5'-flanking region to regulate gsn transcription during heat shock and vegetative growth.

Single dose of a vaccine based on DNA encoding mycobacterial hsp65 protein plus TDM-loaded PLGA microspheres protects mice against a virulent strain of Mycobacterium tuberculosis

The high incidence of tuberculosis around the world and the inability of BCG to protect certain populations clearly indicate that an improved vaccine against tuberculosis is needed. A single antigen, the mycobacterial heat shock protein hsp65, is sufficient to protect BALB/c mice against challenge infection when administered as DNA vaccine in a three-dose-based schedule. In order to simplify the vaccination schedule, we coencapsulated hsp65-DNA and trehalose dimicolate (TDM) into biodegradable poly(DL-lactide-co-glycolide) (PLGA) microspheres. BALB/c mice immunized with a single dose of DNA-hsp65/TDM-1oaded microspheres produced high levels of IgG2a subtype antibody and high amounts of IFN-gamma in the supernatant of spleen cell cultures. DNA-hsp65/TDM-loaded microspheres were also able to induce high IFN-gamma production in bulk lung cells from challenged mice and confer protection as effective as that attained after three doses of naked DNA administration. This new formulation also allowed a ten-fold reduction in the DNA dose when compared to naked DNA. Thus, this combination of DNA vaccine and adjuvants with immunomodulatory and carrier properties holds the potential for an improved vaccine against tuberculosis.