Uma perspectiva computacional sobre catálise enzimática

Enzymes are extremely efficient catalysts. Here, part of the mechanisms proposed to explain this catalytic power will be compared to quantitative experimental results and computer simulations. Influence of the enzymatic environment over species along the reaction coordinate will be analysed. Concepts of transition state stabilisation and reactant destabilisation will be confronted. Divided site model and near-attack conformation hypotheses will also be discussed. Molecular interactions such as covalent catalysis, general acid-base catalysis, electrostatics, entropic effects, steric hindrance, quantum and dynamical effects will also be analysed as sources of catalysis. Reaction mechanisms, in particular that catalysed by protein tyrosine phosphatases, illustrate the concepts.

Avaliação do destino e bioacumulação de benzo(a)pireno através de simulação computacional

The objective of this work was to evaluate the environmental distribution of benzo(a)pirene, a polycyclic aromatic hydrocarbon, by the EQC model. The modeling of the contaminant distribution was accomplished by means of the fugacity model applied to a hypothetical scenario constituted by air, water, soil and sediment. The modeling and simulations revealed that the soil is the preferential compartment. We also discuss the implications of the results about fate and ecological risks associated with benzo(a)pirene. We concluded that the emissions of HPAs can not be ignored and bioaccumulation among others risks can be induced.

Implementação computacional do modelo carga-fluxo de carga-fluxo de dipolo para cálculo e interpretação das intensidades do espectro infravermelho

The first computational implementation that automates the procedures involved in the calculation of infrared intensities using the charge-charge flux-dipole flux model is presented. The atomic charges and dipoles from the Quantum Theory of Atoms in Molecules (QTAIM) model was programmed for Morphy98, Gaussian98 and Gaussian03 programs outputs, but for the ChelpG parameters only the Gaussian programs are supported. Results of illustrative but new calculations for the water, ammonia and methane molecules at the MP2/6-311++G(3d,3p) theoretical level, using the ChelpG and QTAIM/Morphy charges and dipoles are presented. These results showed excellent agreement with analytical results obtained directly at the MP2/6-311++G(3d,3p) level of theory.

Modelagem da combustão da dinitramida de amônio por simulação computacional

In this work, the combustion process of ammonium dinitramide, ADN, has been modeled in two different situations: decomposition in open environment, with abundant air and decomposition in a rocket motor internal environmental conditions. The profiles of the two processes were achieved, based on molar fractions of the species that compose the products of ADN combustion. The velocity of formation and quantity of species in the open environment was bigger than the ones in the rocket motor environment, showing the effect of the different atmosphere in the reactions kinetics.

Propriedades dinâmicas de fluidos por simulação computacional: métodos híbridos atomístico-contínuo

Computational methods for the calculation of dynamical properties of fluids might consider the system as a continuum or as an assembly of molecules. Molecular dynamics (MD) simulation includes molecular resolution, whereas computational fluid dynamics (CFD) considers the fluid as a continuum. This work provides a review of hybrid methods MD/CFD recently proposed in the literature. Theoretical foundations, basic approaches of computational methods, and dynamical properties typically calculated by MD and CFD are first presented in order to appreciate the similarities and differences between these two methods. Then, methods for coupling MD and CFD, and applications of hybrid simulations MD/CFD, are presented.

Comparação da labilidade de metais empregando voltametria, difusão em filmes finos por gradiente de concentração (DGT) e modelo computacional

Speciation of metals in a synthetic freshwater was comparatively evaluated using Anodic Stripping Voltammetry, Diffusive Gradients in Thin Films and a Chemical Equilibrium Model. The labile fractions of Cu and Zn quantified by DGT were similar to the ones measured by ASV. The labile species of Cd and Pb could not be determined by both experimental methods due to the formation of inert complexes with organic ligands in the sample. Despite the differences among the methods, the speciation results obtained by the use of DGT and ASV agreed well with predictions made by the chemical equilibrium model.

Validação computacional de métodos compostos no estudo de propriedades moleculares

Composite methods using ONIOM and different basis sets have been used to calculate proton and electron affinities for a set of alcohols at QCISD(T)/6-311++G(2df,p) level of theory. The study was carried out considering HF, MP2 and DFT (25 exchange correlation functional) methods. The calculation performed at ONIOM2(QCISD(T)/6-311++G(2df,p):HF/6-31G(d))//ONIOM2(O3LYP/6-31G(d):HF/6-31G(d)) resulted in the smallest average absolute deviation for AP and AE, 4,75 kJ/mol e 0,43 eV, respectively.

Caracterização composicional e transesterificação de óleo de microalga: uma abordagem computacional

A procedure for compositional characterization of a microalgae oil is presented and applied to investigate a microalgae based biodiesel production process through process simulation. The methodology consists of: proposing a set of triacylglycerides (TAG) present in the oil; assuming an initial TAG composition and simulating the transesterification reaction (UNISIM Design, Honeywell) to obtain FAME characterization values (methyl ester composition); evaluating deviations of experimental from calculated values; minimizing the sum of squared deviations by a non-linear optimization algorithm, with TAG molar fractions as decision variables. Biodiesel from the characterized oil is compared to a rapeseed based biodiesel.

QSAR modeling: um novo pacote computacional open source para gerar e validar modelos QSAR

QSAR modeling is a novel computer program developed to generate and validate QSAR or QSPR (quantitative structure- activity or property relationships) models. With QSAR modeling, users can build partial least squares (PLS) regression models, perform variable selection with the ordered predictors selection (OPS) algorithm, and validate models by using y-randomization and leave-N-out cross validation. An additional new feature is outlier detection carried out by simultaneous comparison of sample leverage with the respective Studentized residuals. The program was developed using Java version 6, and runs on any operating system that supports Java Runtime Environment version 6. The use of the program is illustrated. This program is available for download at lqta.iqm.unicamp.br.

Eletroquímica computacional - princípios e diretrizes para utilização

In this paper, a detailed guide for the application of computational electrochemistry is presented. The basic framework of the electrochemical models and their computational solutions are described. We highlighted that the availability of commercial software allows application of the technique by experimentalists with minimal mathematical and computational expertise. The most used packages are indicated. Simulations of typical examples are presented and some references cited to illustrate the wide applicability of computational electrochemistry.

Índexos de localització i deslocalització derivats de la densitat bielectrònica: anàlisi i aplicacions en estructura molecular, reactivitat química i aromaticitat

La present tesi doctoral s'emmarca dins el camp de la química teòrica i computacional. El primer objectiu, que va ser el de partida, consistia a veure si els índexos bielectrònics derivats de la teoria AIM podien ser útils per estudiar la reorganització electrònica al llarg d'una reacció. De forma paral·lela es va plantejar un objectiu metodològic, el càlcul dels índexos mencionats a altres nivells de teoria, com podia ser DFT o CI. CI. Aquests dos objectius varen portar a molts d'altres, tot i que els principals són: Objectiu 1. Estudi de la reorganització de parells d'electrons al llarg d'una reacció. Objectiu 2. Efecte de la solvatació sobre l'estructura de parells d'electrons d'una molècula. Objectiu 3. Càlcul dels índexos de localització i deslocalització als nivells de teoria del funcional de la densitat i post-Hartree-Fock. Objectiu 4. Caracterització de l'enllaç per pont d'hidrogen per mitjà de l'estructura dels parells d'electrons. Objectiu 5. L'índex de deslocalització: nou criteri electrònic d'aromaticitat.

Modelagem da composição química do leite através de indicadores metabólicos em vacas leiteiras de alta produção

A lactação é um processo fisiológico complexo que ainda não foi compreendido na sua totalidade. Inúmeros fatores intervêm na síntese e secreção do leite, sendo os mais importantes a nutrição e o metabolismo endógeno dos nutrientes. A qualidade do leite é valorizada tanto pela sua composição química, como pelo conteúdo de células somáticas. No entanto, visando a comercialização do leite, as maiores mudanças e melhoras na qualidade podem ser atingidas através da manipulação da dieta dos animais, em especial em vacas leiteiras de alta produção. Avaliar os processos de absorção de alimentos, bem como o metabolismo catabólico e anabólico direcionado para a síntese do leite, têm sido uma grande preocupação na pesquisa de nutrição e bioquímica da produção animal. O principal objetivo da presente pesquisa foi gerar modelos matemáticos que pudessem explicar a participação de diferentes metabólitos sobre a composição química do leite. Neste intuito foram coletadas amostras de fluído ruminal, sangue, urina e leite de 140 vacas da raça Holandesa nas primeiras semanas de lactação e mantidas sob sistema semi-intensivo de produção e dieta controlada. Os animais foram selecionados de sistemas de produção no ecossistema do Planalto Médio de Rio Grande do Sul e foram amostrados em dois períodos climáticos críticos. No fluido ruminal foram avaliados o pH e o tempo de redução do azul de metileno. No sangue foram determinados os metabólitos: glicose, colesterol, β-hidroxibutirato (BHB), triglicerídeos, fructosamina, ácidos graxos não esterificados (NEFA), proteínas totais, albumina, globulina, uréia, creatinina, cálcio, fósforo e magnésio. As enzimas: aspartato amino transferase (AST), gama glutamil transferase (GGT) e creatina kinase (CK). Os hormônios: cortisol, insulina, triiodotironina (T3), tiroxina (T4), e leptina. Foi efetuado hemograma, para conhecer: hematócrito, hemoglobina, e contagem total e diferencial de células brancas. Na urina foram dosados: corpos cetônicos, pH e densidade. No leite foi determinada: proteína, gordura, lactose, sólidos totais, sólidos não gordurosos, contagem de células somáticas e uréia. Para a determinação de cada um dos metabólitos ou compostos foram usadas técnicas específicas validadas internacionalmente. Os diferentes valores obtidos constituíram os parâmetros básicos de entrada para a construção dos diversos modelos matemáticos executados para predizer a composição do leite. Mediante procedimentos de regressão linear múltipla algoritmo Stepwise, procedimentos de correlação linear simples de Pearson e procedimentos de análise computacional através de redes neurais, foram gerados diferentes modelos para identificar os parâmetros endógenos de maior relevância na predição dos diferentes componentes do leite. A parametrização das principais rotas bioquímicas, do controle endócrino, do estado de funcionamento hepático, da dinâmica ruminal e da excreção de corpos cetônicos aportou informação suficiente para predizer com diferente grau de precisão o conteúdo dos diferentes sólidos no leite. O presente trabalho é apresentado na forma de quatro artigos correspondentes aos parâmetros energéticos, de controle endócrino, modelagem matemática linear múltipla e predição através de Artificial Neural Networks (ANN).

Implementação de um modelo computacional para estudo do processo Fischer-Tropsch em reator de leito de lama

This work aims at the implementation and adaptation of a computational model for the study of the Fischer-Tropsch reaction in a slurry bed reactor from synthesis gas (CO+H2) for the selective production of hydrocarbons (CnHm), with emphasis on evaluation of the influence of operating conditions on the distribution of products formed during the reaction.The present model takes into account effects of rigorous phase equilibrium in a reactive flash drum, a detailed kinetic model able of predicting the formation of each chemical species of the reaction system, as well as control loops of the process variables for pressure and level of slurry phase. As a result, a system of Differential Algebraic Equations was solved using the computational code DASSL (Petzold, 1982). The consistent initialization for the problem was based on phase equilibrium formed by the existing components in the reactor. In addition, the index of the system was reduced to 1 by the introduction of control laws that govern the output of the reactor products. The results were compared qualitatively with experimental data collected in the Fischer-Tropsch Synthesis plant installed at Laboratório de Processamento de Gás - CTGÁS-ER-Natal/RN

Modelagem Matemática e Computacional de fenômenos eletrocinéticos em meios porosos carregados eletricamente

In this work we present a mathematical and computational modeling of electrokinetic phenomena in electrically charged porous medium. We consider the porous medium composed of three different scales (nanoscopic, microscopic and macroscopic). On the microscopic scale the domain is composed by a porous matrix and a solid phase. The pores are filled with an aqueous phase consisting of ionic solutes fully diluted, and the solid matrix consists of electrically charged particles. Initially we present the mathematical model that governs the electrical double layer in order to quantify the electric potential, electric charge density, ion adsorption and chemical adsorption in nanoscopic scale. Then, we derive the microscopic model, where the adsorption of ions due to the electric double layer and the reactions of protonation/ deprotanaç~ao and zeta potential obtained in modeling nanoscopic arise in microscopic scale through interface conditions in the problem of Stokes and Nerst-Planck equations respectively governing the movement of the aqueous solution and transport of ions. We developed the process of upscaling the problem nano/microscopic using the homogenization technique of periodic structures by deducing the macroscopic model with their respectives cell problems for effective parameters of the macroscopic equations. Considering a clayey porous medium consisting of kaolinite clay plates distributed parallel, we rewrite the macroscopic model in a one-dimensional version. Finally, using a sequential algorithm, we discretize the macroscopic model via the finite element method, along with the interactive method of Picard for the nonlinear terms. Numerical simulations on transient regime with variable pH in one-dimensional case are obtained, aiming computational modeling of the electroremediation process of clay soils contaminated

Estudo por simulação computacional da interação entre polieletrólitos fracos e macroíons esféricos

Pós-graduação em Biofísica Molecular - IBILCE