La Sarcolipine (SLN) est un hélice transmembranaire de 31 résidus dont la function est de reguler l’ATPase-Ca2+ SERCA1a. Ce régulateur peut subir une modification post-traductionnelle chez certaines espèces. Par exemple, chez le lapin, il est palmitoyle ou oleoylé sur le résidu Cys9. Pour comprendre au niveau moléculaire l’effet de cette modification post-traductionnelle sur la SLN, nous avons réalisé des simulations de dynamique moléculaire de la SLN de lapin insérée dans une bicouche de 1-palmitoyl-2-oléoyl-sn-glycéro-3-phosphocholine (POPC), non acylée et palmitoylée. L’analyse de ces simulations démontre que la palmitoylation n’affecte pas la structure secondaire, l’orientation (tilt et azimut), ainsi que l’enfouissemnt de la SLN dans la membrane. De plus, l’analyse de simulations tout atome de la SLN humaine insérée dans une bicouche de POPC montre que la SLN humaine a la même structure secondaire et orientation que la SLN de lapin mais est plus enfouie dans la membraneque celle de lapin, du fait de sa sequence en acides amines N-terminale plus hydrophobe.L’ATPase-Ca2+ SERCA1a, une ATPase de type P, est localisée dans la membrane du reticulum sarcoplasmique des cellules du muscle squelettique. Elle est impliquée dans les processus de contraction/relaxation musculaire en transportant rapidement le Ca2+ cytosolique dans le lumen du reticulum sarcoplasmique grâce à l’énergie fournie par l’hydrolyse de l’ATP. D’importants changements conformationnels de SERCA1a ont lieu durant son cycle catalytique comme le montrent les nombreuses structures cristallines de SERCA1a. En particulier, à l’état E1, la cavité contenant les sites de fixation du Ca2+ est ouverte vers le cytoplasme, tandis qu’à l’état E2, cette cavité est ouverte vers le lumen. La transition de l’état E1 à E2 nécessite la phosphorylation du résidu Asp351. Des structures 3D du complexe SERCA1a-SLN ont été déterminées par diffraction aux rayons X avec SERCA1a dans un état E1-Mg2+. Pour comprendre le mécanisme détaillé de la regulation de SERCA1a par la SLN, des simulations de dynamique moléculaire et des analyses des modes normaux (NMA) ont été réalisées en utilisant la structure 3D du complexe SERCA1a-SLN inséré dans une bicouche de POPC. Les résultats principaux de ces analyses sont les suivants : 1) la SLN régule les transitions E1.Mg2+ → E1.2Ca2+ et E1.Mg2+ → E2 ; 2) l’interaction de la SLN influe sur la structure et la dynamique de SERCA1a et modifie la position de l’hélice transmembranaire TM1 de sorte à ce que la cavité contenant les sites de fixation du Ca2+soit plus ouverte et que les sites soient plus accessibles ; 3) l’interaction de la SLN avec TM6 affecte deux regions de SERCA1a indispensables à sa fonction : en modifiant la structure et la dynamique de TM6, la SLN perturbe la position et la fluctuation des résidus des sites de fixation du Ca2+, leur conférant une conformation inapte à fixer le Ca2+. De même, l’interaction avec TM6 induit la courbure de TM5, ce qui affecte de façon indirecte le site de phosphorylation (éloigné de plus de 35 Å de la SLN) et conduit à l’inhibition de la phosphorylation du résidu Asp351.Nos résultats de cette étude in silico fournissent de nouveaux éléments concernant le mécanisme par lequel la SLN régule SERCA1a et qui pourrait être complétés par des travaux expérimentaux. / Sarcolipin (SLN), a transmembrane helix of 31 residues, binds to and regulates the Ca2+-ATPase SERCA1a. This regulator is post-translationnally modified in some species. For example, in rabbit, it is palmitoylated or oleoylated on its Cys9 residue. To understand at a molecular level, the effect of this post-translationnal modification on SLN, all-atom molecular dynamics simulations of unacylated and palmitoylated rabbit SLN embedded in a POPC bilayer were performed. Analysis of the simulations demonstrates that palmitoylation does not affect the secondary structure, the orientation (tilt and azimuth) as well as the burying of SLN within the membrane. In addition, the analyses of all-atom simulations of human SLN embedded in a POPC bilayer show that human SLN has the same secondary structure and orientation as rabbit SLN but is more buried within the membrane than rabbit SLN as a result of its more hydrophobic N-terminal amino acids sequence.The Ca2+ pump SERCA1a, a P-type ATPase, is localized in the sarcoplasmic reticulum membrane of striated muscle cells. It is involved in the contraction/relaxation process by fast pumping the cytoplasmic Ca2+ from the cytosol to the lumen of the sarcoplasmic reticulum using the energy of ATP hydrolysis. Large conformational changes of SERCA1a occur during its catalytic cycle as evidenced by the various crystal structures of SERCA1a. In particular, in the E1 state, the cavity that contains the Ca2+ binding sites is open toward the cytoplasm while in the E2 state, this cavity is open toward the lumen. The transition from the E1 to the E2 state involves the phosphorylation of Asp351 residue. 3D structures of SERCA1a-SLN complex have been determined by X-Ray diffraction, with SERCA1a in a E1-Mg2+ state. To understand the detailed mechanisms of SERCA1a regulation by SLN, molecular dynamics (MD) simulations and normal mode analysis (NMA) were performed using the 3D structures of SERCA1a-SLN complex embedded in a POPC bilayer. Main results from these analyses are the followings: 1) SLN regulates the E1-Mg2+ → E1-2Ca2+ and E1-Mg2+ → E2 state transitions; 2) interaction of SLN with SERCA1a impact the structure and dynamic of SERCA1a and modifies the position of the transmembrane helix TM1 such that the cavity that contains the Ca2+ binding sites is more widely opened and the Ca2+ binding sites more accessible; 3) SLN interaction with affects two regions essential to its function. By changing the structure and dynamic of TM6, SLN alters the position and fluctuations of residues involved in the Ca2+ binding sites, such that those sites are unable to bind Ca2+. This interaction with TM6 also induces TM5 bending and thus, indirectly modifies the phosphorylation site conformation, leading to the inhibition of Asp351 phosphorylation.Our results from these in silico studies provide new insights into the mechanism by which SLN regulates SERCA1a activity and could be completed by experimental work.

Gelöster Wasserstoff in Metallen führt in vielen Fällen zu einer Reduzierung der Güte von mechanischen Eigenschaften. Dies äußert sich auf vielfältige Weise und wird unter dem Begriff Wasserstoffversprödung zusammengefasst. Für ein grundlegendes Verständnis dieses Phänomens müssen die Vorgänge im Metall auf mikroskopischer Skala ergründet werden. Im Rahmen dieser Arbeit wurde daher ein Aspekt der Wasserstoffversprödung, die Interaktion von Wasserstoff mit Versetzungen, näher untersucht. Zur Untersuchung des Einflusses von Wasserstoff auf die Versetzungsbildung wurden verschiedene Verformungsexperimente an Palladium und Vanadium durchgeführt. Prinzipielle Vorgänge der Defektbildung wurden durch Versuche an einzelnen Versetzungen unter Verwendung von Nanoindentation und Zugexperimenten im ETEM durchgeführt, um einen breiten Überblick zu erlangen. Zusätzlich wurden zum besseren Verständnis der Vorgänge Molekulardynamiksimulationen von derartigen Versuchen ausgeführt. Zur Untersuchung der Interaktion von Versetzungen miteinander wurden Säulen im Mikrometerbereich verformt und Blech durch Kaltwalzen verformt. Des Weiteren wurde durch Hochdruck-Torsion maximale Verformungen realisiert. Die verwendeten Modellmaterialien erlauben es verschiedene prinzipielle Vorgänge der Defektbildung zu untersuchen und so einen breiten Überblick über prinzipielle Vorgänge im kfz Gitter (Palladium) bzw. krz Gitter (Vanadium) zu erhalten.

Le Nafion dans la couche catalytique des PEMFCs peut former un revêtement de film ultra-mince enrobant la surface du catalyseur et ses supports. La morphologie du Nafion se révèle être très sensibles à la nature du matériau sur lequel le film est déposé, et en particulier le caractère hydrophobe/phile de ces matériaux. Notre travail consiste à effectuer une enquête complète sur les effets hydrophiles du substrat sur les propriétés physiques du film ultra-mince de Nafion à des niveaux différents d'hydratation. Par conséquent, nous étudions selon un cadre unique une variété d'environnements spécifiques de la couche du catalyseur de la PEMFC, pouvant aller d'un substrat hydrophobe (carbone) à hydrophile (platine). La méthode numérique choisie pour ce travail est une simulation par Dynamique Moléculaire classique. Les configurations de films ultra-minces correctement thermalisés ont été décrites en détail en fonction de leurs propriétés structurales et dynamiques. / Nafion inside Polymer Electrolyte Membrane Fuel Cells (PEMFC) catalyst layers can be found as an ultra-thin film coating the catalyst and the catalyst support surfaces. Nafion morphology shows to be strongly sensitive to the type of material where the film is deposited, especially the hydrophobic/philic character of these materials. Our work consists in performing a complete investigation of the substrate hydrophilicity effects on the physical properties of Nafion ultra-thin film at different hydration levels. We investigate in a unique framework a variety of environments peculiar of the PEMFC catalyst layer, ranging from hydrophobic (carbon) to hydrophilic (platinum) substrates. The numerical method chosen for this work is classical Molecular Dynamics simulations. The well-thermalized thin-film configurations were described in details in terms of their structural and dynamical properties.

Microbial cell disruption using pressurized gas is a promising approach to improve the lipid extraction yield directly from the wet biomass by eliminating the energy-intensive drying process, which is an integral part of traditional methods. As the process starts with the solubilization of the gas in lipid-rich microbial cells, it is important to understand the solubility of different potential gases in both lipid (triglyceride) and lipid-rich microbial cell culture to design efficient cell disruption processes. In this study, we determined the solubility of different gases (e.g., CO2, CH4, N2, and Ar) in canola oil (triglyceride) using a pressure drop gas apparatus developed in our laboratory. The solubility of different gases in triglyceride followed the trend CO2 > CH4 > Ar > N2. Since the solubility of CO2 was found to be higher compared to other gases, the solubility of CO2 in lipid rich cell culture, cell culture media, and spent media was also determined. It was found that CO2 is more soluble in triglycerides, but less soluble in lipid-rich cell culture compared to CO2 in water. From both thermodynamic models and Monte Carlo simulations, the correlated solubility was found to be in good agreement with the experimental results. CO2 was found to be the most suitable gas for microbial cell disruption because almost 100% cell death occurred when using CO2 whereas more than 85% cells were found to be active after treatment with CH4, N2, and Ar. The optimization of microbial cell disruption was conducted using the combination of Box-Behnken design of experiment (DOE) technique and response surface methodology. The optimized cell disruption conditions were found to be 3900 kPa, 296.5 K, 360 min, and 325 rpm where almost 100% cell death was predicted from the statistical modeling. Finally, it was found that 86% of the total lipid content can be recovered from the wet biomass after treatment with pressurized CO2 under optimized conditions compared to control where up to 74% of the total lipid content can be recovered resulting in 12% increase in the lipid extraction yield using pressurized CO2.

The knowledge of solid-state second phase heterogeneous nucleation process is limited due to the experimental difficulty, such as tiny length scale, short time period, and high temperature condition. In recent years, some significant breakthroughs in nucleation studies have been achieved by aid of computational techniques. In this study, we apply molecular dynamics (MD) simulations to perform with heterogeneous nucleation occurring at grain boundaries (GB) during the austenite (FCC) phase to ferrite (BCC) phase transformation in a pure Fe polycrystalline system. A neighbor vector analysis (NVA) method has been introduced and it is shown how the NVA can be used to determine the misorientation of grain or interphase boundaries, which allow a further investigation of the boundary structure correlated to interfacial energy and mobility during the nucleation and early grain growth stage. Meanwhile, benefited from the MD technique, the bulk energy, grain boundary energy, and interfacial energy can be individually captured during the simulations, which allow a detail analyze of the shape, critical size and nucleation energy of specific nuclei, through the classical nucleation theory (CNT) and according to a faceted-spherical cap geometric model (FSC). In addition, we also compared the results from the classical approach with a new algorithm that combination of the multi-phase field model (MPFM) and the nudged elastic band (NEB) method to demonstrate the CNT in the solid-state conduction. Finally, we extend our simulation method to a more complex triple GB junction nucleation event, and investigate the non-classical barrier-free nucleation behaviors. The results support the critical informations to clarify the initial state of austenite to ferrite transition, and improve our knowledge of the heterogeneous nucleation process, which help to bridge the gap between the experimental measurements and the theoretical calculations. The simulation method also provided a new approach for studying the complicate heterogeneous nucleation phenomenon in solid-state for a wide variety of polycrystalline material systems. / Thesis / Doctor of Philosophy (PhD)

Molekulardynamik (MD) Simulationen wurden an zwei verschiedenen Systemen durchgeführt: 1. Grenzfläche zwischen Polyethylen und isotaktischem Polypropylen (PE-iPP) und 2. Zylindrische Mizellen, bestehend aus Tetradecyltrimethylammoniumbromid (C14TAB), in wässriger Lösung und an Fest-Flüssig-Grenzflächen. Die allgemeinen Schwierigkeiten bei der Simulation von Grenzflächen kristalliner Polymere wurden diskutiert und eine Methode für solche Simulationen vorgeschlagen. Diese Methode wurde zur epitaxialen Kristallisation von PE auf iPP benutzt. Experimentelle Ergebnisse der epitaxialen Kristallisation konnten durch die Simulation bestätigt werden. Ferner konnte vorhergesagt werden, dass PE bevorzugt auf einer iPP-Oberfläche mit hoher Methylgruppenkonzentration kristallisiert. Ebenso wurde durch die MD Simulation vorhergesagt, dass PE in der Grenzflächenregion von einer orthorhombischen zur monoklinischen Kristallstruktur wechselt. Die Simulationsdauer für die Mizellen betrug einige Nanosekunden. Die Ergebnisse für die Mizellen in wässriger Lösung stehen hierbei in guter Übereinstimmung mit experimentellen Werten. Im Widerspruch zur allgemein üblichen Vorstellung führte die Simulation der Mizellen zur Ausbildung eines Hohlraums in ihrer Mitte sowie zu einer inhomogenen Dichte des hydrophoben Mizellkerns. Dies wurde zum Teil der inhomogenen Verteilung der terminalen Methylgruppen im Mizellkern zugeschrieben. Zylindrische und halbzylindrische Mizellen wurden an den Paraffin/Wasser- und Gold/Wasser-Grenzflächen simuliert. / Molecular Dynamic (MD) simulation of two different systems was performed: 1) Polyethylene- isotactic Polypropylene (PE-iPP) interfaces and 2) cylindrical micelles formed by tetradecyl trimethylammonium bromide (C14TAB) molecules in aqueous solution and at solid liquid interfaces. The general difficulties of simulation of polymer crystalline interfaces were discussed and one method was proposed for such simulations. Thise method was used to simulate epitaxial crystallisation of PE on iPP. The experimental results on epitaxial crystallisation were confirmed by MD simulation and in addition epitaxial crystallisation of PE on iPP surface with high dencity of methyl groups was predicted. MD simulation also predicted that PE should change at the interfacial region from the orthorhombic to monoclinic crystalline structure. Several nanoseconds of life of cylindrical micelles were simulated. The simulation results for the micelle in aqueous solution were favourably compared with experimental results. In contradiction to the standard picture of an ionic micelle the simulated micelle formed hole in its centre and the density of the hydrophobic micelle core was inhomogeneous. This effect partially was explained by the inhomogeneous distribution of the terminal methyl groups in the micelle core. Cylindrical and half cylindrical micelles of C14TAB molecules were simulated at the paraffin- and gold-aqueous interfaces.

No description available.

A multiscale model is proposed to study the influence of interfacial interactions at the nanoscale in carbon nanotube(CNT)-polymer nanocomposites on the macroscale bulk elastic material properties. The efficiency of CNT reinforcement in terms of interfacial load transferring is assessed for the non-functionalized and functionalized interfaces between the CNTs and polymer matrix using force field based molecular dynamic simulations at the nanoscale. Polyethylene (PE) as a thermoplastic material is adopted and studied first because of its simplicity. Characterization of the nanoscale load transfer has been done through the identification of representative nanoscale interface elements for unfunctionalized CNT-PE interface models which are studied parametrically in terms of the length of the PE chains, the number of the PE chains and the "grip" position. Referring to the non-functionalized interface, CNTs interact with surrounding polymer only through weakly nonbonded van der Waals (vdW) forces in our study. Once appropriate values of these parameters are deemed to yield sufficiently converged results, the representative interface elements are subjected to normal and sliding mode simulations in order to obtain the force-separation responses at 100K and 300K for unfunctionalized CNT-PE interfaces. To study the functionalization effects, atomistic interface representative elements for functionalized CNT-PE interface are built based on non-functionalized interface models by grafting functional groups between the PE matrix and the graphene sheet. This introduces covalent bonding forces in addition to the non-bonded vdW forces. A modified consistent covalent force field (CVFF) and adaptive intermolecular reactive empirical bond order (AIREBO) potentials, both of which account for bond breaking, are applied to investigate the interfacial characteristic of functionalized CNT-PE interface in terms of the force-separation responses at 100K in both normal opening and sliding mode separations. In these studies, the focus has been on the influence of the functionalization density on the load transfer at the nanoscale interface. As an important engineering material, Epon 862/DETDA epoxy polymer,a thermoset plastic, has also been used as the polymer matrix material in order to see the difference in interfacial load transfer between a network structured polymer and the amorphous entangled structure of the PE matrix. As for thermoset epoxy polymer, emphasis has been put on investigating the effects of the crosslink density of the epoxy network on the interfacial load transfer ability for both non-functionalized and functionalized CNT-Epoxy interface at different temperatures(100K and 300K) and on the functionalization effect influenceing the interfacial interactions at the functionalized CNT-Epoxy interface. Cohesive zone traction-displacement laws are developed based on the force-separation responses obtained from the MD simulations for both non-functionalied and functionalized CNT-PE/epoxy interfaces. Using the cohesive zone laws, the influence of the interface on the effective elastic material properties of the nanocomposites are observed and determined in continuum level models using analytic and computational micromechanics approaches, allowing for the assessment of the improvement in reinforcement efficiency of CNTs due to the functionalization. It is found that the inclusion of the nanoscale interface in place of the perfectly bonded interface results in effective elastic properties which are dependent on the applied strain and temperature in accordance with the interface sensitivity to those effects, and which are significantly diminished from those obtained under the perfect interface assumption for non-functionalized nanocomposites. Better reinforcement efficiency of CNTs are also observed for the nanocomposites with the functionalized interface between CNTs and polymer matrix, which results in large increasing for the effective elastic material properties relative to the non-functionalized nanocomposites with pristine CNTs. Such observations indicates that trough controlling the degree of functionalization, i.e. the number and distribution of covalent bonds between the embedded CNTs and the enveloping polymer, one can tailor to some degree the interfacial load transfer and hence, the effective mechanical properties. The multiscale model developed in this study bridges the atomistic modeling and micromechanics approaches with cohesive zone models, which demonstrates to deepen the understanding of the nanoscale load transfer mechanism at the interface and its effects on the effective mechanical properties of the nanocomposites. It is anticipated that the results can offer insights about how to engineer the interface and improve the design of nanocomposites. / Ph. D.

L’ADN, support de l’information génétique, est constamment altéré par des agents physiques ou chimiques d’origines endogènes (métabolisme) et exogènes (UV, radiations ionisantes, produits chimiques) dont les effets sont génotoxiques. Ces modifications structurales délétères de l’ADN sont éliminées par de nombreux mécanismes de réparation. Parmi eux, le système de réparation par excision de bases (BER) est initié par les ADN glycosylases qui reconnaissent et éliminent les bases endommagées. Dans certaines stratégies anti-cancéreuses, l’utilisation de la chimiothérapie et la radiothérapie ont pour but la destruction des cellules cancéreuses en altérant leur ADN. Dans ce contexte, les ADN glycosylases réparent l’ADN des cellules traitées et induisent une résistance non désirée au traitement, faisant de ces enzymes des cibles thérapeutiques intéressantes. Le but de ces travaux est d’approfondir la compréhension des mécanismes de réparation des ADN glycosylases de la superfamille structurale Fpg/Nei grâce à la modélisation moléculaire et de pouvoir identifier et concevoir des inhibiteurs de ces enzymes. Les simulations de dynamique moléculaire (DM) nous ont permis d’étudier la « Lesion Capping Loop » (LCL) et de l’associer à la stabilisation de la base endommagée positionnée dans le site actif. Nous avons également étudié les chemins de sortie possibles de la base après coupure par l’enzyme et l’implication de la boucle LCL dans ce phénomène grâce à des simulations de DM ciblée (TMD-1). De plus, les simulations de DM couplées à un protocole d’amarrage moléculaire « aveugle » nous ont permis d’identifier 2 sites de fixations possibles majoritaires pour des petites molécules potentiellement inhibitrices. Un de ces sites correspondant au site actif de hNEIL1 a fait l’objet d’un criblage virtuel d’une partie de la base de molécules Ambinter. Ceci nous a permis d’identifier des molécules potentiellement inhibitrices dont les effets seront prochainement testés in vitro dans l’équipe sur la protéine humaine hNeil1. / The DNA, genetic information support, is frequently damaged by physical or chemical agents from endogenous (cell metabolism) and exogenous (UV, ionizing radiations, chemicals) factors whose effects are genotoxic. These deleterious DNA structural alterations are removed by many DNA repair mechanisms. Among them, the base excision repair (BER) is initiated by DNA glycosylases which recognize and remove damaged bases. In some anti-cancer strategies, the use of chemo- and radiotherapy is aimed to cancerous cells destruction by altering their DNA. In that specific context, DNA glycosylases repair the DNA of treated cells and induce unwanted resistance to treatments, making these enzymes interesting therapeutic targets. The purpose of this work is to deepen the repair mechanism knowledge of Fpg/Nei structural superfamily of DNA glycosylases using molecular modeling and designing inhibitors of these enzymes. Molecular dynamic simulations allowed us to study the « Lesion Capping Loop » (LCL) and to associate its role to substrate stabilization in the enzyme active site. We also studied some possible excision’s product release pathways and LCL implication in this phenomena by targeted molecular dynamic simulations (TMD-1). Furthermore, molecular dynamic simulations coupled to a blind molecular docking protocol allowed us to identify 2 possible main binding sites of potential inhibitiors. One of these binding sites corresponding to the hNEIL1 active site has been the object of a virtual screening of the Greenpharma database. This allowed us to identify potential inhibitors whom effects will be soon tested in vitro on the humain protein hNEIL1.

Orientador: Sergio Persio Ravagnani / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-08-07T03:49:14Z (GMT). No. of bitstreams: 1 Trochmann_JoseLuizLino_D.pdf: 1070584 bytes, checksum: 3407aee7ad6d88d9de0a1326aaf3d29d (MD5) Previous issue date: 2006 / Resumo: Neste trabalho de tese foi realizado um estudo do potencial de predição de propriedades de transporte em matrizes poliméricas de poli - imidas, utilizando a simulação dinâmica molecular de gases simples como Oxigênio, Nitrogênio e Dióxido de Carbono. A propriedade de transporte de interesse prático, a permeabilidade de uma membrana polimérica a um dado penetrante, envolve a determinação de propriedades de ordem cinética e termodinâmica, respectivamente a determinação do coeficiente de difusão e da solubilidade deste penetrante na matriz polimérica. Atenção especial foi conferida à propriedade cinética, pela predição do coeficiente de autodifusão dos penetrantes. Num procedimento experimental clássico é de vital importância para significância das conclusões derivadas dos experimentos, o uso de amostras de membranas poliméricas adequadamente preparadas quanto à composição química, estrutura física e morfologia. Analogamente, quando se utiliza a simulação molecular para a predição de propriedades, tais como o coeficiente de autodifusão, também é de fundamental relevância para os resultados obtidos, a qualidade dos modelos moleculares das matrizes poliméricas, que serão usados como base. Assim para a preparação de modelos moleculares com o adequado empacotamento, um procedimento para a obtenção de modelos bem equilibrados foi desenvolvido neste trabalho. Os modelos moleculares desenvolvidos foram usados para a obtenção dos valores de massa específica em função da temperatura, e comparados aos valores experimentais disponíveis e quando necessário a, valores preditos por meio da expressão de massa específica em função da temperatura, acima e abaixo da temperatura de transição. A capacidade do modelo molecular desenvolvido em predizer a massa especifica e temperatura de transição vítrea foi usada como critério para a validação da adequação do empacotamento proposto para o referido modelo molecular da matriz polimérica. Os modelos validados de empacotamento, células amorfas, foram utilizados para o cálculo do coeficiente de autodifusão dos gases acima mencionados, através do da simulação dinâmica molecular. A comparação dos coeficientes de autodifusão obtidos das poli-imidas aromáticas e éster imidas, BAAF, 6FDA-ODA, PMDA-ODA e BA-20DA, para os gases O2, N2 e CO2, com os dados experimentais, permitiu concluir a adequação das células amorfas e do esquema de simulação dinâmica molecular para a predição do coeficiente de autodifusão.. A versão preditiva de Vrentas e Duda, baseada na teoria do volume livre, foi utilizada para a predição dos coeficientes de autodifusão da água e do etanol para as poli-imidas acima. , Estes valores, quando comparados com os valores obtidos através da simulação dinâmica molecular mostram a validade de ambas as teorias para a predição da cinética de difusão de penetrantes em matrizes poliméricas complexas / Abstract: In this thesis a study of the predictive potential of the molecular dynamic simulation was performed for transport properties of light gases in polyimide matrix. From de practical point of view permeability is the property of most interest, and involves kinetics as well as thermodynamics properties, diffusion coefficient and solubility of the penetrants molecule in the bulk polymeric matrix, this work will be focus in the former. As important as is in as experimental work, a well prepared polymeric membrane is essential for the significance of the draw conclusions. Therefore a special attention was take in the preparation of the bulk molecular polymeric model, the so called amorphous cell, in order to obtain well-equilibrated molecular packing models for the polyimide matrixes. The amorphous cells were prepared throughout thermodynamic transforms, using one or more of the statistical ensembles and cell specific volume obtained as a function of temperature, this data was compared against the experimental data available, and when necessary to data obtained via predictive methods. The molecular packing model ability to predict the glass transition temperature was used as criteria to validate de amorphous cell, to be used in the molecular dynamic' simulations allow the matrix to be locally flexible and coupled to the classic molecular dynamics simulation. The resulting self diffusion coefficients for the polyimide, BAAF, 6FDA-ODA, PMDA­ODA and BA-20DA for the gases O2, N2 e CO2 were compared to the experimental data. The lack of quality experimental diffusion data available for polyimide membranes for larger penetrants as water and ethanol, showed up as a good opportunity to assess the predictive capability of the molecular dynamic simulation for self diffusion coefficients, considering the relevant technological relevance of polyimide membranes for pervaporation process. The data of self diffusion coefficient produced by the predictive version of free-volume theory after Vrentas and Duda, was compared with the data produced via coupled molecular dynamic simulation for the water and ethanol penetrants, showing the relevance of both theories for the prediction of penetrants kinetic in complex polymeric matrixes / Doutorado / Ciencia e Tecnologia de Materiais / Doutor em Engenharia Química