Thomsen-Friedenreich antigens in bladder cancer: evaluation of their prognostic value

RESUMO: Introdução. O cancro de bexiga é uma patologia comum que representa o 6° e o 5° cancro mais incidente em Portugal e na Itália, respetivamente. Em mais de metade dos casos ocorre reincidência durante o primeiro ano, requerendo acompanhamento clínico ao longo da vida. A instilação intravesical de Bacillus Calmette-Guérin (BCG) (uma estirpe atenuada do Mycobacterium bovis) representa uma imunoterapia eficaz no combate ao cancro de bexiga, no entanto, muitos aspetos da interação de BCG com as células tumorais bem como com as células do sistema imunitário permanecem por desvendar. As células tumorais de bexiga expressam frequentemente as formas sialiladas dos antigénios de Thomsen-Friedenreich (TF), i.e., sialil-T (sT) e sialil-Tn (sTn). Contudo ainda se desconhece o significado da sua expressão na malignidade tumoral e se afeta a eficácia da terapêutica BCG. Objetivo do estudo. Investigar o papel dos antigénios sT e sTn no fenótipo maligno de células de cancro de bexiga bem como na resposta mediada pelo sistema imunitário à terapia com BCG. Metodologia. Para tal, foram utilizadas as linhas celulares de cancro da bexiga HT1376 e MCR, geneticamente modificadas por transdução com vetores codificantes para as sialiltransferases ST3GAL1 ou ST6GALNAC1, de forma a expressar homogeneamente os antigénios sT ou sTn respetivamente. Estes modelos celulares foram estudados após confronto com BCG. O nível de BCG internalizado foi avaliado por citometria de fluxo. O perfil global de expressão genética dos modelos celulares antes e após incubação com BCG foi analisado pela tecnologia de microarray. O perfil de citocinas secretadas pelos modelos celulares após incubação com BCG, bem como de macrófagos estimulados pelo secretoma de células de cancro de bexiga que por sua vez foram estimuladas previamente por BCG, foi estudado pelo sistema multiplex de “imuno-esferas”. Resultados. A análise do transcritoma dos modelos celulares revelou que grupos de genes envolvidos em funções específicas foram modulados em paralelo nos dois modelos celulares, após transdução, independentemente da sialiltransferase expressa. Ou seja, em células que expressavam a sialiltransferase ST3GAL1 ou ST6GALNAC1, os genes envolvidos na regulação da segregação cromossómica e na reparação do DNA foram consistentemente regulados negativamente. Genes descritos na literatura como marcadores para o cancro de bexiga foram também modulados. A incubação com BCG resultou numa tendência ao aumento da expressão de genes relevantes na preservação e estabilidade genómica e menor malignidade, no entanto, apenas em células que expressavam sT ou sTn. Entre as dez citocinas testadas, apenas a IL-6 e IL-8 foram expressas pelas linhas celulares de cancro da bexiga, com indução destas após estimulação com BCG, e principalmente em células que expressavam ST3GAL1 ou ST6GALNAC1. Em macrófagos, citocinas inflamatórias, tais como IL-1β, IL-6 e TNFα, e a citocina anti-inflamatória IL-10, foram induzidas apenas pelo secretoma de células de cancro da bexiga confrontadas com BCG, com maior relevância quando estas expressavam ST3GAL1 ou ST6GALNAC1, prevendo a estimulação de macrófagos semelhantes aos de tipo M1 e uma melhor resposta à terapia com BCG. Conclusões. O efeito geral da expressão destas sialiltransferases e dos produtos enzimáticos sT ou sTn nas células de cancro de bexiga conduz a um fenótipo de maior malignidade. Contudo, a maior avidez de estas na produção de citocinas inflamatórias após confronto com BCG, bem como a maior capacidade de estimulação de macrófagos, predirá uma resposta à terapia com BCG mais eficaz em tumores que expressem os antigénios de TF sialilados. Tais conclusões são totalmente concordantes com os nossos mais recentes dados clínicos obtidos em colaboração, que mostram que em doentes com cancro de bexiga que expressam sTn respondem melhor a terapia BCG. ----------ABSTRACT: Background. Bladder cancer is a common malignancy representing the 6th and the 5th most incident cancer in Portugal and in Italy, respectively. More than half of the cases relapse within one year, requiring though a lifelong follow-up. Intravesical instillation of Bacillus Calmette-Guérin (BCG) (an attenuated strain of Mycobacterium bovis) represents an effective immunotherapy of bladder cancer, although many aspects of the interaction of BCG with cancer cells and host immune cells remain obscure. Bladder cancer cells often express the sialylated forms of the Thomsen-Friedenreich (TF), i.e., sialil-T (sT) e sialil-Tn (sTn). However, it’s still unknown the sense of such expression in tumour malignancy and in the BCG therapy efficacy. Aim of the study. To investigate the role of the sT and sTn antigens on the malignant phenotype of bladder cancer cells and the immune mediated response to BCG therapy. Experimental. We have utilized populations of the bladder cancer cell lines HT1376 and MCR, genetically modified by transduction with the sialyltransferases ST3GAL1 or ST6GALNAC1 to express homogeneously sT or sTn antigens. The level of BCG internalized was assessed by flow cytometry. The whole gene expression profile of BCG-challenged or unchallenged bladder cancer cell lines was studied by microarray technology. The profile of cytokines secreted by BCG-challenged bladder cancer cells and that of macrophages challenged by the secretome of BCG-challenged bladder cancer cells was studied by multiplex immune-beads assay. Results. Transcriptome analysis of the sialyltransferase-transduced cells revealed that groups of genes involved in specific functions were regulated in parallel in the two cell lines, regardless the sialyltransferase expressed. Namely, in sialyltransferase-expressing cells, genes involved in the proper chromosomal segregation and in the DNA repair were consistently down-regulated, while genes reported in literature as markers for bladder cancer were modulated. BCG-challenging induced a tendency to up-regulation of the genes preserving genomic stability and reducing malignancy, but only in cells expressing either sT or sTn. Among the ten cytokines tested, only IL-6 and IL-8 were expressed by bladder cancer cell lines and up-regulated by BCG-challenging, mainly in sialyltransferases-expressing cells. In macrophages, inflammatory cytokines, such as IL-1β, IL-6 and TNFα, and the antinflammatory IL-10 were induced only by the secretome of BCG-challenged bladder cancer cells, particularly when expressing either sialyltransferase, predicting the stimulation of M1-like macrophages and a better response to BCG therapy. Conclusions. The general effect of the expression of the two sialyltransferases and their products in the bladder cancer cells is toward a more malignant phenotype. However, the stronger ability of sialyltransferase expressing cells to produce inflammatory cytokines upon BCG-challenging and to stimulate macrophages predicts a more effective response to BCG in tumours expressing the sialylated TF antigens. This is fully consistent with our recent clinical data obtained in collaboration, showing that patients with bladder cancer expressing sTn respond better to BCG therapy.

Structural stability of adenylate kinase from the sulfate-reducing bacteria Desulfovibrio gigas

Biophysical Chemistry 110 (2004) 83–92

Automatic learning for the classification of chemical reactions and in statistical thermodynamics

This Thesis describes the application of automatic learning methods for a) the classification of organic and metabolic reactions, and b) the mapping of Potential Energy Surfaces(PES). The classification of reactions was approached with two distinct methodologies: a representation of chemical reactions based on NMR data, and a representation of chemical reactions from the reaction equation based on the physico-chemical and topological features of chemical bonds. NMR-based classification of photochemical and enzymatic reactions. Photochemical and metabolic reactions were classified by Kohonen Self-Organizing Maps (Kohonen SOMs) and Random Forests (RFs) taking as input the difference between the 1H NMR spectra of the products and the reactants. The development of such a representation can be applied in automatic analysis of changes in the 1H NMR spectrum of a mixture and their interpretation in terms of the chemical reactions taking place. Examples of possible applications are the monitoring of reaction processes, evaluation of the stability of chemicals, or even the interpretation of metabonomic data. A Kohonen SOM trained with a data set of metabolic reactions catalysed by transferases was able to correctly classify 75% of an independent test set in terms of the EC number subclass. Random Forests improved the correct predictions to 79%. With photochemical reactions classified into 7 groups, an independent test set was classified with 86-93% accuracy. The data set of photochemical reactions was also used to simulate mixtures with two reactions occurring simultaneously. Kohonen SOMs and Feed-Forward Neural Networks (FFNNs) were trained to classify the reactions occurring in a mixture based on the 1H NMR spectra of the products and reactants. Kohonen SOMs allowed the correct assignment of 53-63% of the mixtures (in a test set). Counter-Propagation Neural Networks (CPNNs) gave origin to similar results. The use of supervised learning techniques allowed an improvement in the results. They were improved to 77% of correct assignments when an ensemble of ten FFNNs were used and to 80% when Random Forests were used. This study was performed with NMR data simulated from the molecular structure by the SPINUS program. In the design of one test set, simulated data was combined with experimental data. The results support the proposal of linking databases of chemical reactions to experimental or simulated NMR data for automatic classification of reactions and mixtures of reactions. Genome-scale classification of enzymatic reactions from their reaction equation. The MOLMAP descriptor relies on a Kohonen SOM that defines types of bonds on the basis of their physico-chemical and topological properties. The MOLMAP descriptor of a molecule represents the types of bonds available in that molecule. The MOLMAP descriptor of a reaction is defined as the difference between the MOLMAPs of the products and the reactants, and numerically encodes the pattern of bonds that are broken, changed, and made during a chemical reaction. The automatic perception of chemical similarities between metabolic reactions is required for a variety of applications ranging from the computer validation of classification systems, genome-scale reconstruction (or comparison) of metabolic pathways, to the classification of enzymatic mechanisms. Catalytic functions of proteins are generally described by the EC numbers that are simultaneously employed as identifiers of reactions, enzymes, and enzyme genes, thus linking metabolic and genomic information. Different methods should be available to automatically compare metabolic reactions and for the automatic assignment of EC numbers to reactions still not officially classified. In this study, the genome-scale data set of enzymatic reactions available in the KEGG database was encoded by the MOLMAP descriptors, and was submitted to Kohonen SOMs to compare the resulting map with the official EC number classification, to explore the possibility of predicting EC numbers from the reaction equation, and to assess the internal consistency of the EC classification at the class level. A general agreement with the EC classification was observed, i.e. a relationship between the similarity of MOLMAPs and the similarity of EC numbers. At the same time, MOLMAPs were able to discriminate between EC sub-subclasses. EC numbers could be assigned at the class, subclass, and sub-subclass levels with accuracies up to 92%, 80%, and 70% for independent test sets. The correspondence between chemical similarity of metabolic reactions and their MOLMAP descriptors was applied to the identification of a number of reactions mapped into the same neuron but belonging to different EC classes, which demonstrated the ability of the MOLMAP/SOM approach to verify the internal consistency of classifications in databases of metabolic reactions. RFs were also used to assign the four levels of the EC hierarchy from the reaction equation. EC numbers were correctly assigned in 95%, 90%, 85% and 86% of the cases (for independent test sets) at the class, subclass, sub-subclass and full EC number level,respectively. Experiments for the classification of reactions from the main reactants and products were performed with RFs - EC numbers were assigned at the class, subclass and sub-subclass level with accuracies of 78%, 74% and 63%, respectively. In the course of the experiments with metabolic reactions we suggested that the MOLMAP / SOM concept could be extended to the representation of other levels of metabolic information such as metabolic pathways. Following the MOLMAP idea, the pattern of neurons activated by the reactions of a metabolic pathway is a representation of the reactions involved in that pathway - a descriptor of the metabolic pathway. This reasoning enabled the comparison of different pathways, the automatic classification of pathways, and a classification of organisms based on their biochemical machinery. The three levels of classification (from bonds to metabolic pathways) allowed to map and perceive chemical similarities between metabolic pathways even for pathways of different types of metabolism and pathways that do not share similarities in terms of EC numbers. Mapping of PES by neural networks (NNs). In a first series of experiments, ensembles of Feed-Forward NNs (EnsFFNNs) and Associative Neural Networks (ASNNs) were trained to reproduce PES represented by the Lennard-Jones (LJ) analytical potential function. The accuracy of the method was assessed by comparing the results of molecular dynamics simulations (thermal, structural, and dynamic properties) obtained from the NNs-PES and from the LJ function. The results indicated that for LJ-type potentials, NNs can be trained to generate accurate PES to be used in molecular simulations. EnsFFNNs and ASNNs gave better results than single FFNNs. A remarkable ability of the NNs models to interpolate between distant curves and accurately reproduce potentials to be used in molecular simulations is shown. The purpose of the first study was to systematically analyse the accuracy of different NNs. Our main motivation, however, is reflected in the next study: the mapping of multidimensional PES by NNs to simulate, by Molecular Dynamics or Monte Carlo, the adsorption and self-assembly of solvated organic molecules on noble-metal electrodes. Indeed, for such complex and heterogeneous systems the development of suitable analytical functions that fit quantum mechanical interaction energies is a non-trivial or even impossible task. The data consisted of energy values, from Density Functional Theory (DFT) calculations, at different distances, for several molecular orientations and three electrode adsorption sites. The results indicate that NNs require a data set large enough to cover well the diversity of possible interaction sites, distances, and orientations. NNs trained with such data sets can perform equally well or even better than analytical functions. Therefore, they can be used in molecular simulations, particularly for the ethanol/Au (111) interface which is the case studied in the present Thesis. Once properly trained, the networks are able to produce, as output, any required number of energy points for accurate interpolations.

Gammaretroviral and lentiviral vectors for gene therapy: stability and inactivation mechanisms

Dissertation presented to obtain a Ph.D degree in Engineering and Technology Sciences, Gene Therapy at the Instituto de Tecnologia Quimica e Biológica, Universidade Nova de Lisboa

Genomic and transcriptomic features of Chlamydia trachomatis. Tracking the basis for the ecological success

Dissertation presented to obtain a Ph.D degree in Biology, speciality in Molecular Genetics, by Universidade Nova de Lisboa, Faculdade de Ciências e Tecnologia

A proteomic approach toward the selection of proteins with enhanced intrinsic conformational stability

Journal of Proteome Research (2006)5: 2720-2726

Metal ions modulate the folding and stability of the tumor suppressor protein S100A2

Febs Journal (2009)276:1776-1786

Penicillium glabrum cork colonising isolates – preliminary analysis of their genomic similari

Cork stopper manufacturing process includes an operation, known as stabilisation, by which humid cork slabs are extensively colonised by fungi. The effects of fungal growth on cork are yet to be completely understood and are considered to be involved in the so called “cork taint” of bottled wine. It is essential to identify environmental constraints which define the appearance of the colonising fungal species and to trace their origin to the forest and/or as residents in the manufacturing space. The present article correlates two sets of data, from consecutive years and the same season, of systematic biologic sampling of two manufacturing units, located in the North and South of Portugal. Chrysonilia sitophila dominance was identified, followed by a high diversity of Penicillium species. Penicillium glabrum, found in all samples, was the most frequent isolated species. P. glabrum intra-species variability was investigated using DNA fingerprinting techniques revealing highly discriminative polymorphic markers in the genome. Cluster analysis of P. glabrum data was discussed in relation to the geographical location of strains, and results suggest that P. glabrum arise from predominantly the manufacturing space, although cork resident fungi can also contrib

Ankylosing spondylitis: genomic and functional characterization of candidate genes and their repercussion in clinical practice

RESUMO: Introdução: A espondilite anquilosante (EA) é uma doença inflamatória crónica caracterizada pela inflamação das articulações sacroilíacas e da coluna. A anquilose progressiva motiva uma deterioração gradual da função física e da qualidade de vida. O diagnóstico e o tratamento precoces podem contribuir para um melhor prognóstico. Neste contexto, a identificação de biomarcadores, assume-se como sendo muito útil para a prática clínica e representa hoje um grande desafio para a comunidade científica. Objetivos: Este estudo teve como objetivos: 1 - caracterizar a EA em Portugal; 2 - investigar possíveis associações entre genes, MHC e não-MHC, com a suscetibilidade e as características fenotípicas da EA; 3 - identificar genes candidatos associados a EA através da tecnologia de microarray. Material e Métodos: Foram recrutados doentes com EA, de acordo com os critérios modificados de Nova Iorque, nas consultas de Reumatologia dos diferentes hospitais participantes. Colecionaram-se dados demográficos, clínicos e radiológicos e colhidas amostras de sangue periférico. Selecionaram-se de forma aleatória, doentes HLA-B27 positivos, os quais foram tipados em termos de HLA classe I e II por PCR-rSSOP. Os haplótipos HLA estendidos foram estimados pelo algoritmo Expectation Maximization com recurso ao software Arlequin v3.11. As variantes alélicas dos genes IL23R, ERAP1 e ANKH foram estudadas através de ensaios de discriminação alélica TaqMan. A análise de associação foi realizada utilizando testes da Cochrane-Armitage e de regressão linear, tal como implementado pelo PLINK, para variáveis qualitativas e quantitativas, respetivamente. O estudo de expressão génica foi realizado por Illumina HT-12 Whole-Genome Expression BeadChips. Os genes candidatos foram validados usando qPCR-based TaqMan Low Density Arrays (TLDAs). Resultados: Foram incluídos 369 doentes (62,3% do sexo masculino, com idade média de 45,4 ± 13,2 anos, duração média da doença de 11,4 ± 10,5 anos). No momento da avaliação, 49,9% tinham doença axial, 2,4% periférica, 40,9% mista e 7,1% entesopática. A uveíte anterior aguda (33,6%) foi a manifestação extra-articular mais comum. Foram positivos para o HLA-B27, 80,3% dos doentes. Os haplótipo A*02/B*27/Cw*02/DRB1*01/DQB1*05 parece conferir suscetibilidade para a EA, e o A*02/B*27/Cw*01/DRB1*08/DQB1*04 parece conferir proteção em termos de atividade, repercussão funcional e radiológica da doença. Três variantes (2 para IL23R e 1 para ERAP1) mostraram significativa associação com a doença, confirmando a associação destes genes com a EA na população Portuguesa. O mesmo não se verificou com as variantes estudadas do ANKH. Não se verificou associação entre as variantes génicas não-MHC e as manifestações clínicas da EA. Foi identificado um perfil de expressão génica para a EA, tendo sido validados catorze genes - alguns têm um papel bem documentado em termos de inflamação, outros no metabolismo da cartilagem e do osso. Conclusões: Foi estabelecido um perfil demográfico e clínico dos doentes com EA em Portugal. A identificação de variantes génicas e de um perfil de expressão contribuem para uma melhor compreensão da sua fisiopatologia e podem ser úteis para estabelecer modelos com relevância em termos de diagnóstico, prognóstico e orientação terapêutica dos doentes. -----------ABSTRACT: Background: Ankylosing Spondylitis (AS) is a chronic inflammatory disorder characterized by inflammation in the spine and sacroiliac joints leading to progressive joint ankylosis and in progressive deterioration of physical function and quality of life. An early diagnosis and early therapy may contribute to a better prognosis. The identification of biomarkers would be helpful and represents a great challenge for the scientific community. Objectives: The present study had the following aims: 1- to characterize the pattern of AS in Portuguese patients; 2- to investigate MHC and non-MHC gene associations with susceptibility and phenotypic features of AS and; 3- to identify candidate genes associated with AS by means of whole-genome microarray. Material and Methods: AS was defined in accordance to the modified New York criteria and AS cases were recruited from hospital outcares patient clinics. Demographic and clinical data were recorded and blood samples collected. A random group of HLA-B27 positive patients and controls were selected and typed for HLA class I and II by PCR-rSSOP. The extended HLA haplotypes were estimated by Expectation Maximization Algorithm using Arlequin v3.11 software. Genotyping of IL23R, ERAP1 and ANKH allelic variants was carried out with TaqMan allelic discrimination assays. Association analysis was performed using the Cochrane-Armitage and linear regression tests as implemented in PLINK, for dichotomous and quantitative variables, respectively. Gene expression profile was carried out using Illumina HT-12 Whole-Genome Expression BeadChips and candidate genes were validated using qPCR-based TaqMan Low Density Arrays (TLDAs). Results: A total of 369 patients (62.3% male; mean age 45.4±13.2 years; mean disease duration 11.4±10.5 years), were included. Regarding clinical disease pattern, at the time of assessment, 49.9% had axial disease, 2.4% peripheral disease, 40.9% mixed disease and 7.1% isolated enthesopathic disease. Acute anterior uveitis (33.6%) was the most common extra-articular manifestation. 80.3% of AS patients were HLA-B27 positive. The haplotype A*02/B*27/Cw*02/DRB1*01/DQB1*05 seems to confer susceptibility to AS, whereas A*02/B*27/Cw*01/DRB1*08/DQB1*04 seems to provide protection in terms of disease activity, functional and radiological repercussion. Three markers (two for IL23R and one for ERAP1) showed significant single-locus disease associations. Association of these genes with AS in the Portuguese population was confirmed, whereas ANKH markers studied did not show an association with AS. No association was seen between non-MHC genes and clinical manifestations of AS. A gene expression signature for AS was established; among the fourteen validated genes, a number of them have a well-documented inflammatory role or in modulation of cartilage and bone metabolism. Conclusions: A demographic and clinical profile of patients with AS in Portugal was established. Identification of genetic variants of target genes as well as gene expression signatures could provide a better understanding of AS pathophysiology and could be useful to establish models with relevance in terms of susceptibility, prognosis, and potential therapeutic guidance.

Insights into human carboxylesterase 2 stability and activity in vivo and in vitro

Dissertation to obtain a Master Degree in Biotechnology

Insight into the multicopper oxidases stability

Dissertation presented to obtain the PhD degree in Biochemistry

Weekly tracking of stability of the flow of conversations into the subprocesses of the Last Planner System

21st Annual Conference of the International Group for Lean Construction – IGLC 21 – Fortaleza, Brazil

Protein stability in a proteomic perspective

Dissertation presented to obtain the Ph.D degree in Biochemistry

Functional genomic analysis of heat stress in Vitis vinifera

Grapevine (Vitis vinifera) is one of most agro-economically important fruit crops worldwide, with a special relevance in Portugal where over 300 varieties are used for wine production. Due to global warming, temperature stress is currently a serious issue affecting crop production especially in temperate climates. Mobile genetic elements such as retrotransposons have been shown to be involved in environmental stress induced genetic and epigenetic modifications. In this study, sequences related to Grapevine Retrotransposon 1 (Gret1) were utilized to determine heat induced genomic and transcriptomic modifications in Touriga Nacional, a traditional Portuguese grapevine variety. For this purpose, growing canes were treated to 42 oC for four hours and leaf genomic DNA and RNA was utilized for various techniques to observe possible genomic alterations and variation in transcription levels of coding and non-coding sequences between non-treated plants and treated plants immediately after heat stress (HS-0 h) or after a 24 hour recovery period (HS-24 h). Heat stress was found to induce a significant decrease in Gret1 related sequences in HS-24 h leaves, indicating an effect of heat stress on genomic structure. In order to identify putative heat induced DNA modifications, genome wide approaches such as Amplified Fragment Length Polymorphism were utilized. This resulted in the identification of a polymorphic DNA fragment in HS-0 h and HS-24 h leaves whose sequence mapped to a genomic region flanking a house keeping gene (NADH) that is represented in multiple copies in the Vitis vinifera genome. Heat stress was also found to affect the transcript levels of various non-coding and gene coding sequences. Accordingly, quantitative real time PCR results established that Gret1 related sequences are up regulated immediately after heat stress whereas the level of transcript of genes involved in identification and repair of double strand breaks are significantly down regulated in HS-0 h plants. Taken together, the results of this work demonstrated heat stress affects both genomic integrity and transcription levels.

Genomic and environmental factors influence Wolbachia-Drosophila symbiosis

Intracellular, vertically transmitted bacteria form complex and intimate relationships with their hosts. Wolbachia, maternally transmitted α- proteobacteria, live within the cells of numerous arthropod species. Wolbachia are famous master manipulators of insect reproduction: to favour their own spread they can induce male killing, parthenogenesis or cytoplasmic incompatibility. Wolbachia can also protect various insects from pathogens, which makes them a promising tool for the control of vector-borne diseases. Mosquitoes with Wolbachia have already been released in the wild to eliminate dengue. Yet, how Wolbachia manipulate their hosts remains largely unknown.(...)