Modelação molecular de associações entre aniões e receptores sintéticos

A aplicação de simulações de mecânica e dinâmica molecular ao estudo de sistemas supramoleculares tem adquirido, ao longo dos últimos anos, enorme relevância. A sua utilização não só tem levado a uma melhor compreensão dos mecanismos de formação desses mesmos sistemas, como também tem fornecido um meio para o desenvolvimento de novas arquitecturas supramoleculares. Nesta tese são descritos os trabalhos de mecânica e dinâmica molecular desenvolvidos no âmbito do estudo de associações supramoleculares entre aniões e receptores sintéticos do tipo [2]catenano, [2]rotaxano e pseudorotaxano. São ainda estudados complexos supramoleculares envolvendo receptores heteroditópicos do tipo calix[4]diquinona e pares iónicos formados por aniões halogeneto e catiões alcalinos e amónio. Os estudos aqui apresentados assentam essencialmente em duas vertentes: no estudo das propriedades dinâmicas em solução dos vários complexos supramoleculares considerados e no cálculo das energias livres de Gibbs de associação relativas dos vários iões aos receptores sintéticos. As metodologias utilizadas passaram por dinâmica molecular convencional e REMD (Replica Exchange Molecular Dynamics), para o estudo das propriedades em solução, e por cálculos de integração termodinâmica e MMPBSA (Molecular Mechanics – Poisson Boltzmann Surface Area), para a computação das energias livres de associação relativas. Os resultados obtidos, além de terem permitido uma visão mais detalhada dos mecanismos envolvidos no reconhecimento e associação dos vários receptores aos aniões e pares iónicos abordados, encontram-se, globalmente, de acordo com os análogos determinados experimentalmente, validando assim as metodologias empregadas. Em jeito de conclusão, investigou-se ainda a capacidade de um dos receptores heteroditópicos estudados para assistir favoravelmente na migração do par iónico KCl através da interface água-clorofórmio. Para tal, foram utilizadas simulações SMD (Steered Molecular Dynamics) para a computação do perfil de energia livre de Gibbs associada à migração do par iónico através da interface.

Estudo de receptores sintéticos no reconhecimento molecular de substratos aniónicos e iões metálicos

O trabalho descrito insere-se no âmbito da Química Supramolecular e consistiu no desenvolvimento de receptores artificiais, na forma protonada ou complexada, para o reconhecimento molecular de substratos aniónicos derivados de ácidos carboxílicos incluindo os herbicidas PMG2-, ATCP- e 2,4-D-. Foram investigadas duas séries de aniões, uma alifática (ox2-, mal2-, suc2-, glu2-, adip2-, cit3- e cta3-) e outra aromática (bzc-, naphc-, anthc-, pyrc-, ph2-, iph2-, tph2-, btc3-, dihyac2-, 4,4-dibzc2-, 3-nitrobzc- e 4-nitrobzc-). Foram sintetizados sete novos ligandos macrocíclicos simétricos constituídos por anéis aromáticos piridina ou fenantrolina ligados por cadeias de poliaminas saturadas. O comportamento ácido-base destes macrociclos foi investigado em solução aquosa e as constantes de protonação correspondentes determinadas por métodos potenciométricos e de RMN de 1H. As propriedades de complexação destes ligandos com os iões metálicos Ni2+, Cu2+, Zn2+, Cd2+ e Pb2+ foram também estudadas por métodos potenciométricos nas mesmas condições experimentais, tendo revelado que os macrociclos de dimensão média são capazes de acomodar um ou dois iões metálicos. O complexo dinuclear de Cu(II) derivado do macrociclo com dois grupos piridina foi utilizado como receptor de aniões carboxilato originando complexos ternários. Todos os complexos foram caracterizados em solução por espectroscopias de UV/vis/IVpróx e de RMN. As espécies paramagnéticas foram também caracterizadas por espectroscopia de RPE. A formação de espécies ternárias foi ainda investigada por espectrometria de massa ESI-MS e ESI-MS/MS. As estruturas cristalinas de alguns dos complexos foram determinadas por difracção de raios-X. As formas protonadas dos macrociclos foram utilizadas como receptores de uma grande variedade de aniões carboxilato. O reconhecimento molecular entre os receptores e os substratos aniónicos foi investigado em solução por métodos potenciométricos e de espectroscopia de RMN com determinação das constantes de associação. Os agregados supramoleculares foram caracterizados no estado sólido por difracção de raios-X. Finalmente as associações supramoleculares foram estudadas em solução por métodos de dinâmica molecular com determinação dos termos entrópicos e entálpicos das energias livres de ligação. Em resumo, nesta tese apresentam-se os resultados de estudos realizados com duas famílias de macrociclos: desde a síntese dos compostos, passando por estudos em solução e finalizando com simulação molecular. Este estudo sistemático através da conjugação de metodologias complementares permitiu caracterizar ao nível macroscópico e microscópico as associações moleculares.

Síntese de azacalixarenos para o reconhecimento molecular de aniões

No âmbito desta dissertação são apresentados novos recetores sintéticos baseados na plataforma macrocíclica tetraazacalix[2]areno[2]triazina incorporando bases de purina ou unidades de aminoácidos naturais. Estes recetores foram preparados tendo como objectivo o reconhecimento seletivo de fármacos e aniões biologicamente relevantes. No primeiro capítulo é apresentada uma revisão bibliográfica sobre recetores derivados de tetraazacalix[2]areno[2]triazina demonstrando-se que a funcionalização desta plataforma, nos anéis de triazina ou nos anéis benzénicos, encontra-se ainda na infância bem como a investigação das suas afinidades para aniões. Subsequentemente, considerando os aniões principais estudados nesta tese, aniões carboxilatos alifáticos, aromáticos e aminoácidos é apresentada uma revisão bibliográfica sobre os estudos de reconhecimento molecular reportados na literatura, entre derivados de calix[4]areno e este tipo de aniões. No segundo capítulo descrevem-se as sínteses de recetores com uma ou duas unidades de aminoácido, L-alanina (AC1A, AC2A) e L-triptofano (AC2T), ancoradas nos anéis benzénicos da plataforma tetraazacalix[2]areno[2]triazina. Como demostrado no capítulo 3, os grupos N-H dos azotos em ponte são locais de reconhecimento de aniões. Neste contexto, estes macrociclos foram também metilados nos azotos em ponte (Me4AC2A e Me4AC2T) de modo a direcionar o reconhecimento cooperativo de aniões exclusivamente através dos grupos N-H das unidades de aminoácidos. A lipofilicidade destas moléculas foi alterada por substituição dos átomos de cloro das triazinas por dihexilamina. Os compostos sintetizados foram caracterizados através de técnicas espetroscópicas complementadas em alguns compostos por difração de raios X de cristal único. No capítulo 3 apresentam-se os estudos de reconhecimento molecular entre os recetores sintetizados e aniões derivados de ácidos carboxílicos (mono, di e tricarboxílicos) alifáticos, aromáticos, isoméricos e aniões inorgânicos. Estes estudos foram efetuados em DMSO-d6 ou em CDCl3 por titulações de RMN de 1H com determinação das respetivas constantes de associação. Para os recetores AC1A e AC2A o reconhecimento da maioria dos aniões estudados ocorre simultaneamente através dos grupos N-H dos azotos em ponte e do N-H da amida do braço da L-alanina. Contudo, no caso do isoftalato e tricarboxilato, com dois e três grupos carboxilato em posição meta, ocorre preferencialmente através dos dois braços de L-alanina como sugerido por estudos de dinâmica molecular em DMSO. Os estudos de associação realizados para os macrociclos que contêm unidades de L-triptofano, AC2T e Me4AC2T, mostraram que o reconhecimento de aniões é efetuado preferencialmente através de ligações de hidrogénio estabelecidas com os grupos amida e amina do grupo indole dos braços de aminoácido em detrimento das aminas em ponte como em AC1A e AC2A. Para a série dos carboxilatos alifáticos verifica-se que o recetor Me4AC2A tem maior afinidade com o anião glutarato (K= 389 M-1, DMSO-d6) enquanto o recetor Me4AC2T associa-se mais fortemente com o anião oxalato (K= 776 M-1, CDCl3). Todas as Job plots efetuadas confirmaram que as entidades estudadas obedeciam a uma estequiometria recetor/anião de 1:1.

Diffusion coefficients in dense and supercritical fluids

Os coeficientes de difusão (D 12) são propriedades fundamentais na investigação e na indústria, mas a falta de dados experimentais e a inexistência de equações que os estimem com precisão e confiança em fases comprimidas ou condensadas constituem limitações importantes. Os objetivos principais deste trabalho compreendem: i) a compilação de uma grande base de dados para valores de D 12 de sistemas gasosos, líquidos e supercríticos; ii) o desenvolvimento e validação de novos modelos de coeficientes de difusão a diluição infinita, aplicáveis em amplas gamas de temperatura e densidade, para sistemas contendo componentes muito distintos em termos de polaridade, tamanho e simetria; iii) a montagem e teste de uma instalação experimental para medir coeficientes de difusão em líquidos e fluidos supercríticos. Relativamente à modelação, uma nova expressão para coeficientes de difusão a diluição infinita de esferas rígidas foi desenvolvida e validada usando dados de dinâmica molecular (desvio relativo absoluto médio, AARD = 4.44%) Foram também estudados os coeficientes de difusão binários de sistemas reais. Para tal, foi compilada uma extensa base de dados de difusividades de sistemas reais em gases e solventes densos (622 sistemas binários num total de 9407 pontos experimentais e 358 moléculas) e a mesma foi usada na validação dos novos modelos desenvolvidos nesta tese. Um conjunto de novos modelos foi proposto para o cálculo de coeficientes de difusão a diluição infinita usando diferentes abordagens: i) dois modelos de base molecular com um parâmetro específico para cada sistema, aplicáveis em sistemas gasosos, líquidos e supercríticos, em que natureza do solvente se encontra limitada a apolar ou fracamente polar (AARDs globais na gama 4.26-4.40%); ii) dois modelos de base molecular biparamétricos, aplicáveis em todos os estados físicos, para qualquer tipo de soluto diluído em qualquer solvente (apolar, fracamente polar e polar). Ambos os modelos dão origem a erros globais entre 2.74% e 3.65%; iii) uma correlação com um parâmetro, específica para coeficientes de difusão em dióxido de carbono supercrítico (SC-CO2) e água líquida (AARD = 3.56%); iv) nove correlações empíricas e semi-empíricas que envolvem dois parâmetros, dependentes apenas da temperatura e/ou densidade do solvente e/ou viscosidade do solvente. Estes últimos modelos são muito simples e exibem excelentes resultados (AARDs entre 2.78% e 4.44%) em sistemas líquidos e supercríticos; e v) duas equações preditivas para difusividades de solutos em SC-CO2, em que os erros globais de ambas são inferiores a 6.80%. No global, deve realçar-se o facto de os novos modelos abrangerem a grande variedade de sistemas e moléculas geralmente encontrados. Os resultados obtidos são consistentemente melhores do que os obtidos com os modelos e abordagens encontrados na literatura. No caso das correlações com um ou dois parâmetros, mostrou-se que estes mesmos parâmetros podem ser ajustados usando um conjunto muito pequeno de dados, e posteriormente serem utilizados na previsão de valores de D 12 longe do conjunto original de pontos. Uma nova instalação experimental para medir coeficientes de difusão binários por técnicas cromatográficas foi montada e testada. O equipamento, o procedimento experimental e os cálculos analíticos necessários à obtenção dos valores de D 12 pelo método de abertura do pico cromatográfico, foram avaliados através da medição de difusividades de tolueno e acetona em SC-CO2. Seguidamente, foram medidos coeficientes de difusão de eucaliptol em SC-CO2 nas gamas de 202 – 252 bar e 313.15 – 333.15 K. Os resultados experimentais foram analisados através de correlações e modelos preditivos para D12.

A computational study of ionic liquids used in the production of fuels and biofuels

For the past decades it has been a worldwide concern to reduce the emission of harmful gases released during the combustion of fossil fuels. This goal has been addressed through the reduction of sulfur-containing compounds, and the replacement of fossil fuels by biofuels, such as bioethanol, produced in large scale from biomass. For this purpose, a new class of solvents, the Ionic Liquids (ILs), has been applied, aiming at developing new processes and replacing common organic solvents in the current processes. ILs can be composed by a large number of different combinations of cations and anions, which confer unique but desired properties to ILs. The ability of fine-tuning the properties of ILs to meet the requirements of a specific application range by mixing different cations and anions arises as the most relevant aspect for rendering ILs so attractive to researchers. Nonetheless, due to the huge number of possible combinations between the ions it is required the use of cheap predictive approaches for anticipating how they will act in a given situation. Molecular dynamics (MD) simulation is a statistical mechanics computational approach, based on Newton’s equations of motion, which can be used to study macroscopic systems at the atomic level, through the prediction of their properties, and other structural information. In the case of ILs, MD simulations have been extensively applied. The slow dynamics associated to ILs constitutes a challenge for their correct description that requires improvements and developments of existent force fields, as well as larger computational efforts (longer times of simulation). The present document reports studies based on MD simulations devoted to disclose the mechanisms of interaction established by ILs in systems representative of fuel and biofuels streams, and at biomass pre-treatment process. Hence, MD simulations were used to evaluate different systems composed of ILs and thiophene, benzene, water, ethanol and also glucose molecules. For the latter molecules, it was carried out a study aiming to ascertain the performance of a recently proposed force field (GROMOS 56ACARBO) to reproduce the dynamic behavior of such molecules in aqueous solution. The results here reported reveal that the interactions established by ILs are dependent on the individual characteristics of each IL. Generally, the polar character of ILs is deterministic in their propensity to interact with the other molecules. Although it is unquestionable the advantage of using MD simulations, it is necessary to recognize the need for improvements and developments of force fields, not only for a successful description of ILs, but also for other relevant compounds such as the carbohydrates.