Analysis of the organic matter associated to sea salt : definition of potential molecular markers based on the volatile composition and presence of glycosidic derivatives and polysaccharides

Sea salt is a natural product obtained from the evaporation of seawater in saltpans due to the combined effect of wind and sunlight. Nowadays, there is a growing interest for protection and re-valorisation of saltpans intrinsically associated to the quality of sea salt that can be evaluated by its physico-chemical properties. These man-made systems can be located in different geographical areas presenting different environmental surroundings. During the crystallization process, organic compounds coming from these surroundings can be incorporated into sea salt crystals, influencing their final composition. The organic matter associated to sea salt arises from three main sources: algae, surrounding bacterial community, and anthropogenic activity. Based on the hypothesis that sea salt contains associated organic compounds that can be used as markers of the product, including saltpans surrounding environment, the aim of this PhD thesis was to identify these compounds. With this purpose, this work comprised: 1) a deep characterisation of the volatile composition of sea salt by headspace solid phase microextraction combined with comprehensive two-dimensional gas chromatography time-of-flight mass spectrometry (HS-SPME/GCGC–ToFMS) methodology, in search of potential sea salt volatile markers; 2) the development of a methodology to isolate the polymeric material potentially present in sea salt, in amounts that allow its characterisation in terms of polysaccharides and protein; and 3) to explore the possible presence of triacylglycerides. The high chromatographic resolution and sensitivity of GC×GC–ToFMS enabled the separation and identification of a higher number of volatile compounds from sea salt, about three folds, compared to unidimentional chromatography (GC–qMS). The chromatographic contour plots obtained revealed the complexity of marine salt volatile composition and confirmed the relevance of GC×GC–ToFMS for this type of analysis. The structured bidimentional chromatographic profile arising from 1D volatility and 2D polarity was demonstrated, allowing more reliable identifications. Results obtained for analysis of salt from two locations in Aveiro and harvested over three years suggest the loss of volatile compounds along the time of storage of the salt. From Atlantic Ocean salts of seven different geographical origins, all produced in 2007, it was possible to identify a sub-set of ten compounds present in all salts, namely 6-methyl-5-hepten-2-one, 2,2,6-trimethylcyclohexanone, isophorone, ketoisophorone, β-ionone-5,6-epoxide, dihydroactinidiolide, 6,10,14-trimethyl-2-pentadecanone, 3-hydroxy-2,4,4-trimethylpentyl 2-methylpropanoate, 2,4,4-trimethylpentane-1,3-diyl bis(2-methylpropanoate), and 2-ethyl-1-hexanol. These ten compounds were considered potential volatile markers of sea salt. Seven of these compounds are carotenoid-derived compounds, and the other three may result from the integration of compounds from anthropogenic activity as metabolites of marine organisms. The present PhD work also allowed the isolation and characterisation, for the first time, of polymeric material from sea salt, using 16 Atlantic Ocean salts. A dialysis-based methodology was developed to isolate the polymeric material from sea salt in amounts that allowed its characterisation. The median content of polymeric material isolated from the 16 salts was 144 mg per kg of salt, e.g. 0.014% (w/w). Mid-infrared spectroscopy and thermogravimetry revealed the main occurrence of sulfated polysaccharides, as well as the presence of protein in the polymeric material from sea salt. Sea salt polysaccharides were found to be rich in uronic acid residues (21 mol%), glucose (18), galactose (16), and fucose (13). Sulfate content represented a median of 45 mol%, being the median content of sulfated polysaccharides 461 mg/g of polymeric material, which accounted for 66 mg/kg of dry salt. Glycosidic linkage composition indicates that the main sugar residues that could carry one or more sulfate groups were identified as fucose and galactose. This fact allowed to infer that the polysaccharides from sea salt arise mainly from algae, due to their abundance and composition. The amino acid profile of the polymeric material from the 16 Atlantic Ocean salts showed as main residues, as medians, alanine (25 mol%), leucine (14), and valine (14), which are hydrophobic, being the median protein content 35 mg/g, i.e. 4,9 mg per kg of dry salt. Beside the occurrence of hydrophobic volatile compounds in sea salt, hydrophobic non-volatile compounds were also detected. Triacylglycerides were obtained from sea salt by soxhlet extraction with n-hexane. Fatty acid composition revealed palmitic acid as the major residue (43 mol%), followed by stearic (13), linolenic (13), oleic (12), and linoleic (9). Sea salt triacylglycerides median content was 1.5 mg per kg of dry salt. Both protein and triacylglycerides seem to arise from macro and microalgae, phytoplankton and cyanobacteria, due to their abundance and composition. Despite the variability resulting from saltpans surrounding environment, this PhD thesis allowed the identification of a sea salt characteristic organic compounds profile based on volatile compounds, polysaccharides, protein, and triacylglycerides.

Hydraulics of stepped chutes: experimental-numerical-theoretical study

O estudo de escoamentos turbulentos em descarregadores em degraus tem sido um desafio para os investigadores. A macro-rugosidade do leito, a ondulação da superfície livre, a intermitência da localização da secção inicial de entrada de ar e o escoamento bi-fásico a jusante da secção inicial de entrada de ar fazem com que a caracterização do escoamento deslizante sobre turbilhões em descarregadores em degraus não seja simples. Actualmente, é possível combinar técnicas de medição fiáveis com simulações numéricas e análise teórica. Nesta dissertação, o estudo experimental baseia-se em resultados experimentais obtidos em duas instalações experimentais: a instalação A, do Laboratório Nacional de Engenharia Civil (LNEC), dotada de um descarregador em degraus com declive de 1V:0.75H e a instalação B, do Instituto Superior Técnico (IST), dotada de um descarregador em degraus com declive de 1V:2H. Uma sonda de ar, um tubo de Pitot modificado e vários hidrómetros permitiram o estudo do escoamento deslizante sobre turbilhões mono- e bifásico ao longo do descarregador em degraus da instalação A. Na bacia de dissipação de energia, a observação visual do escoamento foi auxiliada pelas leituras efectuadas em tomadas de pressão localizadas na soleira da bacia e pelas leituras da altura do escoamento efectuadas em réguas graduadas localizadas nas paredes da bacia. Na instalação B foram utilizados um tubo de Pitot e vários hidrómetros para estudar a região não arejada do escoamento deslizante sobre turbilhões no descarregador em degraus. O estudo numérico é baseado em simulações numéricas da região não arejada do escoamento deslizante sobre turbilhões sobre descarregadores em degraus com o código comercial de CFD FLOW-3D® de modo a reproduzir as condições ensaiadas experimentalmente. As simulações beneficiaram da técnica de blocos múltiplos (multi-block) num sistema de coordenadas cartesianas, da determinação da superfície livre pelo método TruVOF e da utilização de dois modelos de turbulência: os modelos k- e RNG k-. Por último, o estudo teórico consistiu em desenvolver um modelo simplificado 1D para determinar as características hidráulicas principais do trecho não arejado do escoamento deslizante sobre turbilhões em descarregadores em degraus. O modelo foi desenvolvido a partir das equações de Navier-Stokes, conjuntamente com resultados experimentais e numéricos. Os resultados apresentados nesta dissertação contribuem para o conhecimento do escoamento deslizante sobre turbilhões em descarregadores em degraus, nomeadamente na região não arejada, na secção inicial de entrada de ar e na região arejada. A hidráulica dos dissipadores de energia a jusante de descarregadores em degraus, em particular das bacias tipo III do USBR, é também objecto de estudo. Em relação à região não arejada do escoamento em descarregadores com declive acentuado, e com base em resultados experimentais e numéricos, são propostas expressões para estimar o desenvolvimento da altura equivalente de água, da espessura da camada limite, da concentração média de ar, do coeficiente de energia cinética, da dissipação de energia, do factor de resistência e do coeficiente n da fórmula de Manning. São ainda propostas expressões adimensionais para a energia cinética turbulenta e sua dissipação. Para declives moderados, são propostas expressões para estimar o desenvolvimento da altura equivalente de água, do coeficiente de energia cinética e da energia específica residual. São ainda apresentados valores do expoente 1/N da expressão adimensional da distribuição de velocidades, quer para descarregadores com declive acentuado quer com declive moderado. Em conformidade com outros estudos centrados em escoamentos de parede e com derivações teóricas, para a região não arejada do escoamento em descarregadores em degraus de acentuado declive, observa-se que o factor de resistência depende da macro-rugosidade criada pelos degraus e da geometria da secção transversal e que o coeficiente n da fórmula de Manning aumenta com a rugosidade. A descrição estatística da turbulência do escoamento é igualmente explorada, contribuindo para o conhecimento da estrutura do escoamento. Observou-se que para números de Reynolds rugoso não superiores a 6.8x104 a energia cinética turbulência e a sua dissipação cumprem leis de semelhança. Estas expressões adimensionais estão de acordo com os resultados obtidos por outros autores para escoamentos completamente desenvolvidos em canais abertos e no escoamento em rios com leito de gravilha. Em acréscimo, a taxa de dissipação de energia, quer para descarregadores de declive acentuado quer de moderado declive, é baixa. Por último, observa-se que os valores da média temporal da concentração de ar entre 0 e 1 medidos na região não arejada do escoamento dizem respeito não só ao ar capturado entre ondas de água, na zona de ondulação da superfície livre, mas também ao ar emulsionado no escoamento, i.e., sob a forma de bolhas de ar, quando perto da secção média inicial de entrada de ar, devido à diferença entre localizações instantânea e média temporal. Foram revistas metodologias e fórmulas para estimar a localização da secção inicial de entrada de ar e apresentadas expressões para estimar a concentração média de ar e a altura equivalente de água nessa secção. Relativamente à região de escoamento arejado em descarregadores em degraus com declive acentuado, os resultados experimentais apresentados nesta dissertação permitiram estimar a influência da definição da superfície livre nos parâmetros hidráulicos da região do escoamento arejado e estimar a máxima elevação do escoamento nesta região do escoamento. Com base nos resultados experimentais obtidos na bacia de dissipação de energia do tipo III do USBR localizada a jusante do descarregador em degraus da instalação A, observou-se que os perfis da altura piezométrica e da altura do escoamento tendem a seguir o perfil recomendado pelo USBR para bacias tipo III. A excepção ocorre à entrada da bacia, onde as alturas piezométricas apresentadas nesta dissertação excedem largamente as apresentadas pelo USBR. É ainda observado que, tal como entre as bacias tipo I e tipo III do USBR, o ressalto hidráulico estabiliza muito mais rapidamente numa bacia tipo III a jusante de um descarregador em degraus do que uma bacia tipo I a jusante do mesmo descarregador em degraus. Finalmente, observa-se que os blocos de amortecimento a colocar no descarregador não têm influência visível nos resultados da altura piezométrica nem da altura do escoamento ao longo da bacia. Relativamente às simulações numéricas do escoamento não arejado, a proximidade entre resultados experimentais e numéricos permite validar o modelo teórico e a integração numérica usados no FLOW-3D®. As simulações desenvolvidas também mostraram que o modelo de turbulência k- permite representar as características do escoamento não arejado em descarregadores em degraus, uma vez que não foram observadas diferenças significativas entre as simulações com este modelo e com o modelo RNG k-. Finalmente, observou-se que o modelo de entrada de ar usado no FLOW-3D® é válido para estimar a localização da secção inicial de entrada de ar. Por último, a proximidade entre os resultados obtidos da aplicação do modelo teórico desenvolvido no âmbito desta dissertação e os resultados experimentais indica que as hipóteses e simplificações consideradas no desenvolvimento do modelo são adequadas.

A study of singular integral operators with shift

Nesta tese, consideram-se operadores integrais singulares com a acção extra de um operador de deslocacamento de Carleman e com coeficientes em diferentes classes de funções essencialmente limitadas. Nomeadamente, funções contínuas por troços, funções quase-periódicas e funções possuíndo factorização generalizada. Nos casos dos operadores integrais singulares com deslocamento dado pelo operador de reflexão ou pelo operador de salto no círculo unitário complexo, obtêm-se critérios para a propriedade de Fredholm. Para os coeficientes contínuos, uma fórmula do índice de Fredholm é apresentada. Estes resultados são consequência das relações de equivalência explícitas entre aqueles operadores e alguns operadores adicionais, tais como o operador integral singular, operadores de Toeplitz e operadores de Toeplitz mais Hankel. Além disso, as relações de equivalência permitem-nos obter um critério de invertibilidade e fórmulas para os inversos laterais dos operadores iniciais com coeficientes factorizáveis. Adicionalmente, aplicamos técnicas de análise numérica, tais como métodos de colocação de polinómios, para o estudo da dimensão do núcleo dos dois tipos de operadores integrais singulares com coeficientes contínuos por troços. Esta abordagem permite também a computação do inverso no sentido Moore-Penrose dos operadores principais. Para operadores integrais singulares com operadores de deslocamento do tipo Carleman preservando a orientação e com funções contínuas como coeficientes, são obtidos limites superiores da dimensão do núcleo. Tal é implementado utilizando algumas estimativas e com a ajuda de relações (explícitas) de equivalência entre operadores. Focamos ainda a nossa atenção na resolução e nas soluções de uma classe de equações integrais singulares com deslocamento que não pode ser reduzida a um problema de valor de fronteira binomial. De forma a atingir os objectivos propostos, foram utilizadas projecções complementares e identidades entre operadores. Desta forma, as equações em estudo são associadas a sistemas de equações integrais singulares. Estes sistemas são depois analisados utilizando um problema de valor de fronteira de Riemann. Este procedimento tem como consequência a construção das soluções das equações iniciais a partir das soluções de problemas de valor de fronteira de Riemann. Motivados por uma grande diversidade de aplicações, estendemos a definição de operador integral de Cauchy para espaços de Lebesgue sobre grupos topológicos. Assim, são investigadas as condições de invertibilidade dos operadores integrais neste contexto.

Development of numerical methodologies for parameter identification and shape optimization in metal forming simulations

Por parte da indústria de estampagem tem-se verificado um interesse crescente em simulações numéricas de processos de conformação de chapa, incluindo também métodos de engenharia inversa. Este facto ocorre principalmente porque as técnicas de tentativa-erro, muito usadas no passado, não são mais competitivas a nível económico. O uso de códigos de simulação é, atualmente, uma prática corrente em ambiente industrial, pois os resultados tipicamente obtidos através de códigos com base no Método dos Elementos Finitos (MEF) são bem aceites pelas comunidades industriais e científicas Na tentativa de obter campos de tensão e de deformação precisos, uma análise eficiente com o MEF necessita de dados de entrada corretos, como geometrias, malhas, leis de comportamento não-lineares, carregamentos, leis de atrito, etc.. Com o objetivo de ultrapassar estas dificuldades podem ser considerados os problemas inversos. No trabalho apresentado, os seguintes problemas inversos, em Mecânica computacional, são apresentados e analisados: (i) problemas de identificação de parâmetros, que se referem à determinação de parâmetros de entrada que serão posteriormente usados em modelos constitutivos nas simulações numéricas e (ii) problemas de definição geométrica inicial de chapas e ferramentas, nos quais o objetivo é determinar a forma inicial de uma chapa ou de uma ferramenta tendo em vista a obtenção de uma determinada geometria após um processo de conformação. São introduzidas e implementadas novas estratégias de otimização, as quais conduzem a parâmetros de modelos constitutivos mais precisos. O objetivo destas estratégias é tirar vantagem das potencialidades de cada algoritmo e melhorar a eficiência geral dos métodos clássicos de otimização, os quais são baseados em processos de apenas um estágio. Algoritmos determinísticos, algoritmos inspirados em processos evolucionários ou mesmo a combinação destes dois são usados nas estratégias propostas. Estratégias de cascata, paralelas e híbridas são apresentadas em detalhe, sendo que as estratégias híbridas consistem na combinação de estratégias em cascata e paralelas. São apresentados e analisados dois métodos distintos para a avaliação da função objetivo em processos de identificação de parâmetros. Os métodos considerados são uma análise com um ponto único ou uma análise com elementos finitos. A avaliação com base num único ponto caracteriza uma quantidade infinitesimal de material sujeito a uma determinada história de deformação. Por outro lado, na análise através de elementos finitos, o modelo constitutivo é implementado e considerado para cada ponto de integração. Problemas inversos são apresentados e descritos, como por exemplo, a definição geométrica de chapas e ferramentas. Considerando o caso da otimização da forma inicial de uma chapa metálica a definição da forma inicial de uma chapa para a conformação de um elemento de cárter é considerado como problema em estudo. Ainda neste âmbito, um estudo sobre a influência da definição geométrica inicial da chapa no processo de otimização é efetuado. Este estudo é realizado considerando a formulação de NURBS na definição da face superior da chapa metálica, face cuja geometria será alterada durante o processo de conformação plástica. No caso dos processos de otimização de ferramentas, um processo de forjamento a dois estágios é apresentado. Com o objetivo de obter um cilindro perfeito após o forjamento, dois métodos distintos são considerados. No primeiro, a forma inicial do cilindro é otimizada e no outro a forma da ferramenta do primeiro estágio de conformação é otimizada. Para parametrizar a superfície livre do cilindro são utilizados diferentes métodos. Para a definição da ferramenta são também utilizados diferentes parametrizações. As estratégias de otimização propostas neste trabalho resolvem eficientemente problemas de otimização para a indústria de conformação metálica.

Volatile terpenoids and C13 norisoprenoids in wines: development of rapid methods of analysis and evaluation of the sesquiterpenoids biological potential

Vitis vinifera L., the most widely cultivated fruit crop in the world, was the starting point for the development of this PhD thesis. This subject was exploited following on two actual trends: i) the development of rapid, simple, and high sensitive methodologies with minimal sample handling; and ii) the valuation of natural products as a source of compounds with potential health benefits. The target group of compounds under study were the volatile terpenoids (mono and sesquiterpenoids) and C13 norisoprenoids, since they may present biological impact, either from the sensorial point of view, as regards to the wine aroma, or by the beneficial properties for the human health. Two novel methodologies for quantification of C13 norisoprenoids in wines were developed. The first methodology, a rapid method, was based on the headspace solid-phase microextraction combined with gas chromatography-quadrupole mass spectrometry operating at selected ion monitoring mode (HS-SPME/GC-qMS-SIM), using GC conditions that allowed obtaining a C13 norisoprenoid volatile signature. It does not require any pre-treatment of the sample, and the C13 norisoprenoid composition of the wine was evaluated based on the chromatographic profile and specific m/z fragments, without complete chromatographic separation of its components. The second methodology, used as reference method, was based on the HS-SPME/GC-qMS-SIM, allowing the GC conditions for an adequate chromatographic resolution of wine components. For quantification purposes, external calibration curves were constructed with β-ionone, with regression coefficient (r2) of 0.9968 (RSD 12.51 %) and 0.9940 (RSD of 1.08 %) for the rapid method and for the reference method, respectively. Low detection limits (1.57 and 1.10 μg L-1) were observed. These methodologies were applied to seventeen white and red table wines. Two vitispirane isomers (158-1529 L-1) and 1,1,6-trimethyl-1,2-dihydronaphthalene (TDN) (6.42-39.45 μg L-1) were quantified. The data obtained for vitispirane isomers and TDN using the two methods were highly correlated (r2 of 0.9756 and 0.9630, respectively). A rapid methodology for the establishment of the varietal volatile profile of Vitis vinifera L. cv. 'Fernão-Pires' (FP) white wines by headspace solid-phase microextraction combined with comprehensive two-dimensional gas chromatography with time-of-flight mass spectrometry (HS-SPME/GCxGC-TOFMS) was developed. Monovarietal wines from different harvests, Appellations, and producers were analysed. The study was focused on the volatiles that seem to be significant to the varietal character, such as mono and sesquiterpenic compounds, and C13 norisoprenoids. Two-dimensional chromatographic spaces containing the varietal compounds using the m/z fragments 93, 121, 161, 175 and 204 were established as follows: 1tR = 255-575 s, 2tR = 0,424-1,840 s, for monoterpenoids, 1tR = 555-685 s, 2tR = 0,528-0,856 s, for C13 norisoprenoids, and 1tR = 695-950 s, 2tR = 0,520-0,960 s, for sesquiterpenic compounds. For the three chemical groups under study, from a total of 170 compounds, 45 were determined in all wines, allowing defining the "varietal volatile profile" of FP wine. Among these compounds, 15 were detected for the first time in FP wines. This study proposes a HS-SPME/GCxGC-TOFMS based methodology combined with classification-reference sample to be used for rapid assessment of varietal volatile profile of wines. This approach is very useful to eliminate the majority of the non-terpenic and non-C13 norisoprenic compounds, allowing the definition of a two-dimensional chromatographic space containing these compounds, simplifying the data compared to the original data, and reducing the time of analysis. The presence of sesquiterpenic compounds in Vitis vinifera L. related products, to which are assigned several biological properties, prompted us to investigate the antioxidant, antiproliferative and hepatoprotective activities of some sesquiterpenic compounds. Firstly, the antiradical capacity of trans,trans-farnesol, cis-nerolidol, α-humulene and guaiazulene was evaluated using chemical (DPPH• and hydroxyl radicals) and biological (Caco-2 cells) models. Guaiazulene (IC50= 0.73 mM) was the sesquiterpene with higher scavenger capacity against DPPH•, while trans,trans-farnesol (IC50= 1.81 mM) and cis-nerolidol (IC50= 1.48 mM) were more active towards hydroxyl radicals. All compounds, with the exception of α-humulene, at non-cytotoxic levels (≤ 1 mM), were able to protect Caco-2 cells from oxidative stress induced by tert-butyl hydroperoxide. The activity of the compounds under study was also evaluated as antiproliferative agents. Guaiazulene and cis-nerolidol were able to more effectively arrest the cell cycle in the S-phase than trans,trans-farnesol and α-humulene, being the last almost inactive. The relative hepatoprotection effect of fifteen sesquiterpenic compounds, presenting different chemical structures and commonly found in plants and plant-derived foods and beverages, was assessed. Endogenous lipid peroxidation and induced lipid peroxidation with tert-butyl hydroperoxide were evaluated in liver homogenates from Wistar rats. With the exception of α-humulene, all the sesquiterpenic compounds under study (1 mM) were effective in reducing the malonaldehyde levels in both endogenous and induced lipid peroxidation up to 35% and 70%, respectively. The developed 3D-QSAR models, relating the hepatoprotection activity with molecular properties, showed good fit (R2LOO > 0.819) with good prediction power (Q2 > 0.950 and SDEP < 2%) for both models. A network of effects associated with structural and chemical features of sesquiterpenic compounds such as shape, branching, symmetry, and presence of electronegative fragments, can modulate the hepatoprotective activity observed for these compounds. In conclusion, this study allowed the development of rapid and in-depth methods for the assessment of varietal volatile compounds that might have a positive impact on sensorial and health attributes related to Vitis vinifera L. These approaches can be extended to the analysis of other related food matrices, including grapes and musts, among others. In addition, the results of in vitro assays open a perspective for the promising use of the sesquiterpenic compounds, with similar chemical structures such as those studied in the present work, as antioxidants, hepatoprotective and antiproliferative agents, which meets the current challenges related to diseases of modern civilization.

Derivation of soil screening values for metals in Portuguese natural soil

The increasing human activity has been responsible by profound changes and a constinuos degradation of the soil compartment in all the European territory. Some European policies are appearing focusing soil’s protection and the management of contaminated sites, in order to recover land for other uses. To regulate the risk assessment and the management of contaminated soils, many European member states adopted soil guideline values, as for example soil screnning values (SSV).These values are particularly useful for the the first tier of the Ecological Risk Assessment (ERA) processes of contaminated sites,especially for a first screening of sites requiring a more site-specific evaluation. Hence, the approriate definition of regional SSVs will have relevant economic impacts in the management of contaminated sites. Portugal is one of European Member States that still lack these soil guideline values. In this context, this study gaves a remarkable contribution in the generation of ecotoxicological data for soil microbiological parameters, terrestrial plants and invertebrates for the derivation of SSVs for uranium (U), cadmium (Cd) and copper (Cu), using a Portuguese natural soil, representative of a dominant type of soil in the Portuguese territory. SSVs were derived based on two methods proposed by the the Technical Guidance Document for Risk Assessment of the European Commission; namely the assessment factor method (AF) and the species sensitivity distribution (SSD) method (with some adaptations). The outputs of both methods were compared and discussed. Further, this study laid the foundation for a deeper reflection about the cut-off (hazard concentration for a given percentage of species - HCps) to be estimated from the SSDs, and to be selected for the derivation of SSVs, with the adequate level of protection. It was proven that this selection may vary for different contaminants, however a clear justification should be given, in each case. The SSvs proposed in this study were for: U (151.4 mg U kg-1dw), Cd (5.6 mg Cd kg-1dw), and Cu (58.5 mg Cu kg-1dw) These values should now be tested for their descriminating power of soils with different levels of contamination. However, this studies clarifies the approach that should be followed for the derivation of SSVs for other metals and organic contaminants, and for other dominant types of Portuguese natural soils.

Evaluation of the matabolic response of osteosarcoma cells to conventional and new anticancer drugs by NMR metabolomics

The main scope of this work was to evaluate the metabolic effects of anticancer agents (three conventional and one new) in osteosarcoma (OS) cells and osteoblasts, by measuring alterations in the metabolic profile of cells by nuclear magnetic resonance (NMR) spectroscopy metabolomics. Chapter 1 gives a theoretical framework of this work, beginning with the main metabolic characteristics that globally describe cancer as well as the families and mechanisms of action of drugs used in chemotherapy. The drugs used nowadays to treat OS are also presented, together with the Palladium(II) complex with spermine, Pd2Spm, potentially active against cancer. Then, the global strategy for cell metabolomics is explained and the state of the art of metabolomic studies that analyze the effect of anticancer agents in cells is presented. In Chapter 2, the fundamentals of the analytical techniques used in this work, namely for biological assays, NMR spectroscopy and multivariate and statistical analysis of the results are described. A detailed description of the experimental procedures adopted throughout this work is given in Chapter 3. The biological and analytical reproducibility of the metabolic profile of MG-63 cells by high resolution magic angle spinning (HRMAS) NMR is evaluated in Chapter 4. The metabolic impact of several factors (cellular integrity, spinning rate, temperature, time and acquisition parameters) on the 1H HRMAS NMR spectral profile and quality is analysed, enabling the definition of the best acquisition parameters for further experiments. The metabolic consequences of increasing number of passages in MG-63 cells as well as the duration of storage are also investigated. Chapter 5 describes the metabolic impact of drugs conventionally used in OS chemotherapy, through NMR metabolomics studies of lysed cells and aqueous extracts analysis. The results show that MG-63 cells treated with cisplatin (cDDP) undergo a strong up-regulation of lipid contents, alterations in phospholipid constituents (choline compounds) and biomarkers of DNA degradation, all associated with cell death by apoptosis. Cells exposed to doxorubicin (DOX) or methotrexate (MTX) showed much slighter metabolic changes, without any relevant alteration in lipid contents. However, metabolic changes associated with altered Krebs cycle, oxidative stress and nucleotides metabolism were detected and were tentatively interpreted at the light of the known mechanisms of action of these drugs. The metabolic impact of the exposure of MG-63 cells and osteoblasts to cDDP and the Pd2Spm complex is described in Chapter 6. Results show that, despite the ability of the two agents to bind DNA, the metabolic consequences that arise from exposure to them are distinct, namely in what concerns to variation in lipid contents (absent for Pd2Spm). Apoptosis detection assays showed that, differently from what was seen for MG-63 cells treated with cDDP, the decreased number of living cells upon exposure to Pd2Spm was not due to cell death by apoptosis or necrosis. Moreover, the latter agent induces more marked alterations in osteoblasts than in cancer cells, while the opposite seemed to occur upon cDDP exposure. Nevertheless, the results from MG-63 cells exposure to combination regimens with cDDP- or Pd2Spm-based cocktails, described in Chapter 7, revealed that, in combination, the two agents induce similar metabolic responses, arising from synergy mechanisms between the tested drugs. Finally, the main conclusions of this thesis are summarized in Chapter 8, and future perspectives in the light of this work are presented.