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.

Quantification de la relation séquence-activité de l’ARN par prédiction de structure tridimensionnelle

Dans un premier temps, nous avons modélisé la structure d’une famille d’ARN avec une grammaire de graphes afin d’identifier les séquences qui en font partie. Plusieurs autres méthodes de modélisation ont été développées, telles que des grammaires stochastiques hors-contexte, des modèles de covariance, des profils de structures secondaires et des réseaux de contraintes. Ces méthodes de modélisation se basent sur la structure secondaire classique comparativement à nos grammaires de graphes qui se basent sur les motifs cycliques de nucléotides. Pour exemplifier notre modèle, nous avons utilisé la boucle E du ribosome qui contient le motif Sarcin-Ricin qui a été largement étudié depuis sa découverte par cristallographie aux rayons X au début des années 90. Nous avons construit une grammaire de graphes pour la structure du motif Sarcin-Ricin et avons dérivé toutes les séquences qui peuvent s’y replier. La pertinence biologique de ces séquences a été confirmée par une comparaison des séquences d’un alignement de plus de 800 séquences ribosomiques bactériennes. Cette comparaison a soulevée des alignements alternatifs pour quelques unes des séquences que nous avons supportés par des prédictions de structures secondaires et tertiaires. Les motifs cycliques de nucléotides ont été observés par les membres de notre laboratoire dans l'ARN dont la structure tertiaire a été résolue expérimentalement. Une étude des séquences et des structures tertiaires de chaque cycle composant la structure du Sarcin-Ricin a révélé que l'espace des séquences dépend grandement des interactions entre tous les nucléotides à proximité dans l’espace tridimensionnel, c’est-à-dire pas uniquement entre deux paires de bases adjacentes. Le nombre de séquences générées par la grammaire de graphes est plus petit que ceux des méthodes basées sur la structure secondaire classique. Cela suggère l’importance du contexte pour la relation entre la séquence et la structure, d’où l’utilisation d’une grammaire de graphes contextuelle plus expressive que les grammaires hors-contexte. Les grammaires de graphes que nous avons développées ne tiennent compte que de la structure tertiaire et négligent les interactions de groupes chimiques spécifiques avec des éléments extra-moléculaires, comme d’autres macromolécules ou ligands. Dans un deuxième temps et pour tenir compte de ces interactions, nous avons développé un modèle qui tient compte de la position des groupes chimiques à la surface des structures tertiaires. L’hypothèse étant que les groupes chimiques à des positions conservées dans des séquences prédéterminées actives, qui sont déplacés dans des séquences inactives pour une fonction précise, ont de plus grandes chances d’être impliqués dans des interactions avec des facteurs. En poursuivant avec l’exemple de la boucle E, nous avons cherché les groupes de cette boucle qui pourraient être impliqués dans des interactions avec des facteurs d'élongation. Une fois les groupes identifiés, on peut prédire par modélisation tridimensionnelle les séquences qui positionnent correctement ces groupes dans leurs structures tertiaires. Il existe quelques modèles pour adresser ce problème, telles que des descripteurs de molécules, des matrices d’adjacences de nucléotides et ceux basé sur la thermodynamique. Cependant, tous ces modèles utilisent une représentation trop simplifiée de la structure d’ARN, ce qui limite leur applicabilité. Nous avons appliqué notre modèle sur les structures tertiaires d’un ensemble de variants d’une séquence d’une instance du Sarcin-Ricin d’un ribosome bactérien. L’équipe de Wool à l’université de Chicago a déjà étudié cette instance expérimentalement en testant la viabilité de 12 variants. Ils ont déterminé 4 variants viables et 8 létaux. Nous avons utilisé cet ensemble de 12 séquences pour l’entraînement de notre modèle et nous avons déterminé un ensemble de propriétés essentielles à leur fonction biologique. Pour chaque variant de l’ensemble d’entraînement nous avons construit des modèles de structures tertiaires. Nous avons ensuite mesuré les charges partielles des atomes exposés sur la surface et encodé cette information dans des vecteurs. Nous avons utilisé l’analyse des composantes principales pour transformer les vecteurs en un ensemble de variables non corrélées, qu’on appelle les composantes principales. En utilisant la distance Euclidienne pondérée et l’algorithme du plus proche voisin, nous avons appliqué la technique du « Leave-One-Out Cross-Validation » pour choisir les meilleurs paramètres pour prédire l’activité d’une nouvelle séquence en la faisant correspondre à ces composantes principales. Finalement, nous avons confirmé le pouvoir prédictif du modèle à l’aide d’un nouvel ensemble de 8 variants dont la viabilité à été vérifiée expérimentalement dans notre laboratoire. En conclusion, les grammaires de graphes permettent de modéliser la relation entre la séquence et la structure d’un élément structural d’ARN, comme la boucle E contenant le motif Sarcin-Ricin du ribosome. Les applications vont de la correction à l’aide à l'alignement de séquences jusqu’au design de séquences ayant une structure prédéterminée. Nous avons également développé un modèle pour tenir compte des interactions spécifiques liées à une fonction biologique donnée, soit avec des facteurs environnants. Notre modèle est basé sur la conservation de l'exposition des groupes chimiques qui sont impliqués dans ces interactions. Ce modèle nous a permis de prédire l’activité biologique d’un ensemble de variants de la boucle E du ribosome qui se lie à des facteurs d'élongation.

Macrophyte community structure and species occurrence in relation to environmental determinants in the ephemeral aquatic habitats of Gavdos, Greece

The aims of this study were to explore the environmental factors that determine the distribution of plant communities in temporary rock pools and provide a quantitative analysis of vegetation-environment relationships for five study sites on the island of Gavdos, southwest of Crete, Greece. Data from 99 rock pools were collected and analysed using Two-Way Indicator Species Analysis (TWINSPAN), Detrended Correspondence Analysis (DCA) and Canonical Correspondence Analysis (CCA) to identify the principal communities and environmental gradients that are linked to community distribution. A total of 46 species belonging to 21 families were recorded within the study area. The dominant families were Labiatae, Gramineae and Compositae while therophytes and chamaephytes were the most frequent life forms. The samples were classified into six community types using TWINSPAN, which were also corroborated by CCA analysis. The principal gradients for vegetation distribution, identified by CCA, were associated with water storage and water retention ability, as expressed by pool perimeter and water depth. Generalised Additive Models (GAMs) were employed to identify responses of four dominant rock pool species to water depth. The resulting species response curves showed niche differentiation in the cases of Callitriche pulchra and Tillaea vaillantii and revealed competition between Zannichellia pedunculata and Chara vulgaris. The use of classification in combination with ordination techniques resulted in a good discrimination between plant communities. Generalised Additive Models are a powerful tool in investigating species response curves to environmental gradients. The methodology adopted can be employed for improving baseline information on plant community ecology and distribution in Mediterranean ephemeral pools.

Crystal structure of a myotoxic Asp49-phospholipase A(2) with low catalytic activity: Insights into Ca2+ -independent catalytic mechanism

A myotoxic Asp49-phospholipase A(2) (Asp49-PLA(2)) with low catalytic activity (BthTX-II from Bothrops jararacussu venom) was crystallized and the molecular-replacement solution has been obtained with a dimer in the asymmetric unit. The quaternary structure of BthTX-II resembles the myotoxic Asp49-PLA2 PrTX-III (piratoxin III from B. pirajai venom) and all non-catalytic and myotoxic dimeric Lys49-PLA(2)s. Despite of this, BthTX-II is different from the highly catalytic and non-myotoxic BthA-I (acidic PLA(2) from B. jararacussu) and other Asp49-PLA(2)s. BthTX-II structure showed a severe distortion of calcium-binding loop leading to displacement of the C-terminal region. Tyr28 side chain, present in this region, is in an opposite position in relation to the same residue in the catalytic activity Asp49-PLA(2)s, making a hydrogen bond with the atom 0 delta 2 of the catalytically active Asp49, which should coordinate the calcium. This high distortion may also be confirmed by the inability of BthTX-II to bind Na+ ions at the Ca2+-binding loop, despite of the crystallization to have occurred in the presence of this ion. In contrast, other Asp49-PLA(2)s which are able to bind Ca2+ ions are also able to bind Na+ ions at this loop. The comparison with other catalytic, non-catalytic and inhibited PLA(2)s indicates that the BthTX-II is not able to bind calcium ions; consequently, we suggest that its low catalytic function is based on an alternative way compared with other PLA(2)s. (c) 2008 Elsevier B.V All rights reserved.

Activity pattern of Cuniculus paca (Rodentia: Cuniculidae) in relation to lunar illumination and other abiotic variables in the southern Brazilian Amazon

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

Optimization of molecular and crystalline forms of drugs, agrochemicals, pesticides in relation to activity, bioavailability, patentability and to the fabrication of polymorphs, solvates, co-crystals with green chemistry methods

This doctorate was funded by the Regione Emilia Romagna, within a Spinner PhD project coordinated by the University of Parma, and involving the universities of Bologna, Ferrara and Modena. The aim of the project was: - Production of polymorphs, solvates, hydrates and co-crystals of active pharmaceutical ingredients (APIs) and agrochemicals with green chemistry methods; - Optimization of molecular and crystalline forms of APIs and pesticides in relation to activity, bioavailability and patentability. In the last decades, a growing interest in the solid-state properties of drugs in addition to their solution chemistry has blossomed. The achievement of the desired and/or the more stable polymorph during the production process can be a challenge for the industry. The study of crystalline forms could be a valuable step to produce new polymorphs and/or co-crystals with better physical-chemical properties such as solubility, permeability, thermal stability, habit, bulk density, compressibility, friability, hygroscopicity and dissolution rate in order to have potential industrial applications. Selected APIs (active pharmaceutical ingredients) were studied and their relationship between crystal structure and properties investigated, both in the solid state and in solution. Polymorph screening and synthesis of solvates and molecular/ionic co-crystals were performed according to green chemistry principles. Part of this project was developed in collaboration with chemical/pharmaceutical companies such as BASF (Germany) and UCB (Belgium). We focused on on the optimization of conditions and parameters of crystallization processes (additives, concentration, temperature), and on the synthesis and characterization of ionic co-crystals. Moreover, during a four-months research period in the laboratories of Professor Nair Rodriguez-Hormedo (University of Michigan), the stability in aqueous solution at the equilibrium of ionic co-crystals (ICCs) of the API piracetam was investigated, to understand the relationship between their solid-state and solution properties, in view of future design of new crystalline drugs with predefined solid and solution properties.

Patterned library analysis: A method for the quantitative assessment of hypotheses concerning the determinants of protein structure

Site-directed mutagenesis and combinatorial libraries are powerful tools for providing information about the relationship between protein sequence and structure. Here we report two extensions that expand the utility of combinatorial mutagenesis for the quantitative assessment of hypotheses about the determinants of protein structure. First, we show that resin-splitting technology, which allows the construction of arbitrarily complex libraries of degenerate oligonucleotides, can be used to construct more complex protein libraries for hypothesis testing than can be constructed from oligonucleotides limited to degenerate codons. Second, using eglin c as a model protein, we show that regression analysis of activity scores from library data can be used to assess the relative contributions to the specific activity of the amino acids that were varied in the library. The regression parameters derived from the analysis of a 455-member sample from a library wherein four solvent-exposed sites in an α-helix can contain any of nine different amino acids are highly correlated (P < 0.0001, R2 = 0.97) to the relative helix propensities for those amino acids, as estimated by a variety of biophysical and computational techniques.

Quantitative characterization of physico-chemical properties of the active layers of reverse osmosis and nanofiltration membranes, and their relation to membrane performance

The main objectives of this dissertation were: (i) to develop experimental and analytical procedures to quantify different physico-chemical properties of the ultra-thin (~ 100 nm) active layers of reverse osmosis (RO) and nanofiltration (NF) membranes and their interactions with contaminants; (ii) to use such procedures to evaluate the similarities and differences between the active layers of different RO/NF membranes; and (iii) to relate characterization results to membrane performance. Such objectives were motivated by the current limited understanding of the physico-chemical properties of active layers as a result of traditional characterization techniques having limitations associated with the nanometer-scale spatial resolution required to study these ultra-thin films. Functional groups were chosen as the main active layer property of interest. Specific accomplishments of this study include the development of procedures to quantify in active layers as a function of pH: (1) the concentration of both negatively and positively ionized functional groups; (2) the stoichiometry of association between ions (i.e., barium) and ionized functional groups (i.e., carboxylate and sulfonate); and (3) the steric effects experienced by ions (i.e., barium). Conceptual and mathematical models were developed to describe experimental results. The depth heterogeneity of the active layer physico-chemical properties and interactions with contaminants studied in this dissertation was also characterized. Additionally, measured concentrations of ionized functional groups in the polyamide active layers of several commercial RO/NF membranes were used as input in a simplified RO/NF transport model to predict the rejection of a strong electrolyte (i.e., potassium iodide) and a weak acid (i.e., arsenious acid) at different pH values based on rejection results at one pH condition. The good agreement between predicted and experimental results showed that the characterization procedures developed in this study serve as useful tools in the advancement of the understanding of the properties and structure of the active layers of RO/NF membranes, and the mechanisms of contaminant transport through them.

Human Mitochondrial Hsp70 (mortalin): Shedding Light On Atpase Activity, Interaction With Adenosine Nucleotides, Solution Structure And Domain Organization.

The human mitochondrial Hsp70, also called mortalin, is of considerable importance for mitochondria biogenesis and the correct functioning of the cell machinery. In the mitochondrial matrix, mortalin acts in the importing and folding process of nucleus-encoded proteins. The in vivo deregulation of mortalin expression and/or function has been correlated with age-related diseases and certain cancers due to its interaction with the p53 protein. In spite of its critical biological roles, structural and functional studies on mortalin are limited by its insoluble recombinant production. This study provides the first report of the production of folded and soluble recombinant mortalin when co-expressed with the human Hsp70-escort protein 1, but it is still likely prone to self-association. The monomeric fraction of mortalin presented a slightly elongated shape and basal ATPase activity that is higher than that of its cytoplasmic counterpart Hsp70-1A, suggesting that it was obtained in the functional state. Through small angle X-ray scattering, we assessed the low-resolution structural model of monomeric mortalin that is characterized by an elongated shape. This model adequately accommodated high resolution structures of Hsp70 domains indicating its quality. We also observed that mortalin interacts with adenosine nucleotides with high affinity. Thermally induced unfolding experiments indicated that mortalin is formed by at least two domains and that the transition is sensitive to the presence of adenosine nucleotides and that this process is dependent on the presence of Mg2+ ions. Interestingly, the thermal-induced unfolding assays of mortalin suggested the presence of an aggregation/association event, which was not observed for human Hsp70-1A, and this finding may explain its natural tendency for in vivo aggregation. Our study may contribute to the structural understanding of mortalin as well as to contribute for its recombinant production for antitumor compound screenings.

Changes in the genetic structure of Bacteria and microbial activity in an agricultural soil amended with tannery sludge

The application of tannery sludge to soils is a form of recycling; however, few studies have examined the impacts of this practice on soil microbial properties. We studied effects of two applications (2006 and 2007) of tannery sludge (with a low chromium content) on the structure of the bacterial community and on the microbial activity of soils. We fertilized an agricultural area in Rolandia, Parana state, Brazil with different doses of sludge based on total N content, which ranged from 0 to 1200 kg N ha(-1). Sludge remained on the soil surface for three months before being plowed. Soils were sampled seven times during the experiment. Bacterial community structure, assessed by denaturing gradient gel electrophoresis (DGGE), was modified by the application of tannery sludge. Soon after the first application, there was clear separation between the bacterial communities in different treatments, such that each dose of sludge was associated with a specific community. These differences remained until 300 days after application and also after the second sludge application, but 666 days after the beginning of the experiment no differences were found in the bacterial communities of the lowest doses and the control. The principal response curve (PRC) analysis showed that the first sludge application strongly stimulated biological activity even 300 days after application. The second application also stimulated activity, but at a lower magnitude and for a shorter time, given that 260 days after the second application there was no difference in biological activity among treatments. PRC also showed that the properties most influenced by the application of tannery sludge were enzymatic activities related to N cycling (asparaginase and urease). The redundancy analysis (RDA) showed that tannery sludge`s influence on microbial activity is mainly related to increases in inorganic N and soil pH. Results showed that changes in the structure of the bacterial community in the studied soils were directly related to changes of their biological activity. (C) 2010 Elsevier Ltd. All rights reserved.

Candida albicans and Candida tropicalis in Oral Candidosis: Quantitative Analysis, Exoenzyme Activity, and Antifungal Drug Sensitivity

Candida albicans and C. tropicalis obtained from whole saliva of patients presenting signs of oral candidosis were assayed for quantification of colony forming units, exoenzyme activity (phospholipase and proteinase) and antifungal drug sensitivity (amphotericin B, fluconazole and itraconazole) by the reference method of the Clinical and Laboratory Standards Institute. The number of colony forming units per milliliter varied according to the Candida species involved and whether a single or mixed infection was present. Proteinase activity was observed in both C. albicans and C. tropicalis, but phospholipase activity was noted only in C. albicans. In vitro resistance to antifungals was verified in both species, but C. tropicalis appears to be more resistant to the tested antifungals than C. albicans.

Habitual physical activity and physical activity intensity: their relation to body composition in 5.0-10.5 year old children

Background: Concerns of a decrease in physical activity levels (PALs) of children and a concurrent increase in childhood obesity exist worldwide. The exact relation between these two parameters however has as yet to be fully defined in children. Objective: This study examined the relation in 47 children, aged 5–10.5 y (mean age 8.4plusminus0.9 y) between habitual physical activity, minutes spent in moderate, vigorous and hard intensity activity and body composition parameters. Design: Total energy expenditure (TEE) was calculated using the doubly labelled water technique and basal metabolic rate (BMR) was predicted from Schofield's equations. PAL was determined by PAL=TEE/BMR. Time spent in moderate, vigorous and hard intensity activity was determined by accelerometry, using the Tritrac-R3D. Body fatness and body mass index (BMI) were used as the two measures of body composition. Results: Body fat and BMI were significantly inversely correlated with PAL (r=-0.43, P=0.002 and r=-0.45, P=0.001). Times spent in vigorous activity and hard activity were significantly correlated to percentage body fat (r=-0.44, P=0.004 and r=-0.39, P=0.014), but not BMI. Children who were in the top tertiles for both vigorous activity and hard activity had significantly lower body fat percentages than those in the middle and lowest tertiles. Moderate intensity activity was not correlated with measures of body composition. Conclusions: As well as showing a significant relation between PAL and body composition, these data intimate that there may be a threshold of intensity of physical activity that is influential on body fatness. In light of world trends showing increasing childhood obesity, this study supports the need to further investigate the importance of physical activity for children.

Role of the 6-20 disulfide bridge in the structure and activity of epidermal growth factor

Two synthetic analogues of murine epidermal. growth factor, [Abu6, 20] mEGF4-48 (where Abu denotes amino-butyric acid) and [G1, M3, K21, H40] mEGF1-48, have been investigated by NMR spectroscopy. [Abu6, 20] mEGF4-48 was designed to determine the contribution of the 6-20 disulfide bridge to the structure and function of mEGF The overall structure of this analogue was similar to that of native mEGF, indicating that the loss of the 6-20 disulfide bridge did not affect the global fold of the molecule. Significant structural differences were observed near the N-terminus, however, with the direction of the polypeptide chain between residues four and nine being altered such that these residues were now located on the opposite face of the main beta-sheet from their position in native mEGF Thermal denaturation experiments also showed that the structure of [Abu6, 20] mEGF4-48 was less stable than that of mEGF. Removal of this disulfide bridge resulted in a significant loss of both mitogenic activity in Balb/c 3T3 cells and receptor binding on A431 cells compared with native mEGF and mEGF4-48, implying that the structural changes in [Abu6, 20] mEGF4-48, although limited to the N-terminus, were sufficient to interfere with receptor binding. The loss of binding affinity probably arose mainly from steric interactions of the dislocated N-terminal region with part of the receptor binding surface of EGF [G1, M3, K21, H40] mEGF1-48 was also synthesized in order to compare the synthetic polypeptide with the corresponding product of recombinant expression. Its mitogenic activity in Balb/c 3T3 cells was similar to that of native mEGF and analysis of its H-1 chemical shifts suggested that its structure was also very similar to native.

Crystal and molecular structure of 2-hydroxy-1-naphthaldehyde isonicotinoyl hydrazone (NIH) and its iron(III) complex: an iron chelator with anti-tumour activity

Previous studies have demonstrated that 2-hydroxy-1-naphthaldehyde isonicotinoyl hydrazone (NIH) and several other aroylhydrazone chelators possess anti-neoplastic activity due to their ability to bind intracellular iron. In this study we have examined the structure and properties of NIH and its Fe-III complex in order to obtain further insight into its anti-tumour activity. Two tridentate NIH ligands deprotonate upon coordination to Fe-III in a meridional fashion to form a distorted octahedral, high-spin complex. Solution electrochemistry of [Fe(NIH-H)(2)](+) shows that the trivalent oxidation state is dominant over a wide potential range and that the Fe-II analogue is not a stable form of this complex. The fact that [Fe(NIH-H)(2)](+) cannot-cycle between the Fe-II and Fe-III states suggests that the production of toxic free- radical species, e.g. OH. or O2(.-),is not part of this ligand's cytotoxic action. This suggestion is supported by cell culture experiments demonstrating that the addition of Fe-III to NIH prevents its anti-proliferative effect. The chemistry of this chelator and its Fe-III complex are discussed in the context of understanding its anti-tumour activity.

Carbon structure and porosity of carbonaceous adsorbents in relation to their adsorption properties

A number of carbonaceous adsorbents were prepared by carbonisation at 600 degrees C following acidic oxidation under various conditions. Effects of the chemical nature of the precursor, such as the ratio of aromatic to aliphatic carbons and oxygen content, on the chemical and structural characteristics of the resultant chars were investigated using C-13 NMR and Raman spectroscopy, respectively. The C-13 NMR spectral parameters of the coal samples show that as the severity of oxidation conditions increased, the ratio of aromatic to aliphatic carbons increased. Furthermore, it was also found that the amount of disorganised carbon affects both the pore structure and the adsorption properties of carbonaceous adsorbents. It is demonstrated that higher amount of the disorganised carbon indicates smaller micropore size. (C) 1999 Elsevier Science Ltd. All rights reserved.