Wall mediated transport in confined spaces: Exact theory for low density

We present a theory for the transport of molecules adsorbed in slit and cylindrical nanopores at low density, considering the axial momentum gain of molecules oscillating between diffuse wall reflections. Good agreement with molecular dynamics simulations is obtained over a wide range of pore sizes, including the regime of single-file diffusion where fluid-fluid interactions are shown to have a negligible effect on the collective transport coefficient. We show that dispersive fluid-wall interactions considerably attenuate transport compared to classical hard sphere theory.

Hydrodynamic origin of diffusion in nanopores

We study the transport of a subcritical Lennard-Jones fluid in a cylindrical nanopore, using a combination of equilibrium and nonequilibrium as well as dual control volume grand canonical molecular dynamics methods. We show that all three techniques yield the same value of the transport coefficient for diffusely reflecting pore walls, even in the presence of viscous transport. We also demonstrate that the classical Knudsen mechanism is not manifested, and that a combination of viscous flow and momentum exchange at the pore wall governs the transport over a wide range of densities.

Modelling the Mechanical Behavior of Polymer-Based Nanocomposites

Molecular dynamics simulations were employed to analyze the mechanical properties of polymer-based nanocomposites with varying nanofiber network parameters. The study was focused on nanofiber aspect ratio, concentration and initial orientation. The reinforcing phase affects the behavior of the polymeric nanocomposite. Simulations have shown that the fiber concentration has a significant effect on the properties, with higher loadings resulting in higher stress levels and higher stiffness, matching the general behavior from experimental knowledge in this field. The results also indicate that, within the studied range, the observed effect of the aspect ratio and initial orientation is smaller than that of the concentration, and that these two parameters are interrelated.

Biochemical and structural characterization of β-lactamases tem-180 and tem-201

With the constant development of new antibiotics, selective pressure is a force to reckon when investigating antibiotic resistance. Although advantageous for medical treatments, it leads to increasing resistance. It is essential to use more potent and toxic antibiotics. Enzymes capable of hydrolyzing antibiotics are among the most common ways of resistance and TEM variants have been detected in several resistant isolates. Due to the rapid evolution of these variants, complex phenotypes have emerged and the need to understand their biological activity becomes crucial. To investigate the biochemical properties of TEM-180 and TEM-201 several computational methodologies have been used, allowing the comprehension of their structure and catalytic activity, which translates into their biological phenotype. In this work we intent to characterize the interface between these proteins and the several antibiotics used as ligands. We performed explicit solvent molecular dynamics (MD) simulations of these complexes and studied a variety of structural and energetic features. The interfacial residues show a distinct behavior when in complex with different antibiotics. Nevertheless, it was possible to identify some common Hot Spots among several complexes – Lys73, Tyr105 and Glu166. The structural changes that occur during the Molecular Dynamic (MD) simulation lead to the conclusion that these variants have an inherent capacity of adapting to the various antibiotics. This capability might be the reason why they can hydrolyze antibiotics that have not been described until now to be degraded by TEM variants. The results obtained with computational and experimental methodologies for the complex with Imipenem have shown that in order to this type of enzymes be able to acylate the antibiotics, they need to be capable to protect the ligand from water molecules.

Mapping stability : real rational maps of degree zero

In memory of our beloved Professor José Rodrigues Santos de Sousa Ramos (1948-2007), who João Cabral, one of the authors of this paper, had the honor of being his student between 2000 and 2006, we wrote this paper following the research by experimentation, using the new technologies to capture a new insight about a problem, as him so much love to do it. His passion was to create new relations between different fields of mathematics. He was a builder of bridges of knowledge, encouraging the birth of new ways to understand this science. One of the areas that Sousa Ramos researched was the iteration of maps and the description of its behavior, using the symbolic dynamics. So, in this issue of this journal, honoring his memory, we use experimental results to find some stable regions of a specific family of real rational maps, the ones that he worked with João Cabral. In this paper we describe a parameter space (a,b) to the real rational maps fa,b(x) = (x2 −a)/(x2 −b), using some tools of dynamical systems, as the study of the critical point orbit and Lyapunov exponents. We give some results regarding the stability of these family of maps when we iterate it, specially the ones connected to the order 3 of iteration. We hope that our results would help to understand better the behavior of these maps, preparing the ground to a more efficient use of the Kneading Theory on these family of maps, using symbolic dynamics.

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

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

Reversible current power supply for fast-field cycling nuclear magnetic resonance equipment

The Fast Field-Cycling Nuclear Magnetic Resonance (FFC-NMR) is a technique used to study the molecular dynamics of different types of materials. The main elements of this equipment are a magnet and its power supply. The magnet used as reference in this work is basically a ferromagnetic core with two sets of coils and an air-gap where the materials' sample is placed. The power supply should supply the magnet being the magnet current controlled in order to perform cycles. One of the technical issues of this type of solution is the compensation of the non-linearities associated to the magnetic characteristic of the magnet and to parasitic magnetic fields. To overcome this problem, this paper describes and discusses a solution for the FFC-NMR power supply based on a four quadrant DC/DC converter.

A low-cost method to test cytotoxic effects of Crotalus vegrandis (Serpentes: Viperidae) venom on kidney cell cultures

The pathogenesis of the renal lesion upon envenomation by snakebite has been related to myolysis, hemolysis, hypotension and/or direct venom nephrotoxicity caused by the venom. Both primary and continuous cell culture systems provide an in vitro alternative for quantitative evaluation of the toxicity of snake venoms. Crude Crotalus vegrandis venom was fractionated by molecular exclusion chromatography. The toxicity of C. vegrandis crude venom, hemorrhagic, and neurotoxic fractions were evaluated on mouse primary renal cells and a continuous cell line of Vero cells maintained in vitro. Cells were isolated from murine renal cortex and were grown in 96 well plates with Dulbecco's Modified Essential Medium (DMEM) and challenged with crude and venom fractions. The murine renal cortex cells exhibited epithelial morphology and the majority showed smooth muscle actin determined by immune-staining. The cytotoxicity was evaluated by the tetrazolium colorimetric method. Cell viability was less for crude venom, followed by the hemorrhagic and neurotoxic fractions with a CT50 of 4.93, 18.41 and 50.22 µg/mL, respectively. The Vero cell cultures seemed to be more sensitive with a CT50 of 2.9 and 1.4 µg/mL for crude venom and the hemorrhagic peak, respectively. The results of this study show the potential of using cell culture system to evaluate venom toxicity.

Dielectric study of triton X100: a glass-forming liquid

Dissertação para obtenção do grau de mestre em Engenharia de Materiais

Development and modulation of magnetic responsive supported ionic liquid membranes

The work presented in this thesis aims at developing a new separation process based on the application of supported magnetic ionic liquid membranes, SMILMs, using magnetic ionic liquids, MILs. MILs have attracted growing interest due to their ability to change their physicochemical characteristics when exposed to variable magnetic field conditions. The magnetic responsive behavior of MILs is thus expected to contribute for the development of more efficient separation processes, such as supported liquid membranes, where MILs may be used as a selective carrier. Driven by the MILs behavior, these membranes are expected to switch reversibly their permeability and selectivity by in situ and non-invasive adjustment of the conditions (e.g. intensity, direction vector and uniformity) of an external applied magnetic field. The development of these magnetic responsive membrane processes were anticipated by studies, performed along the first stage of this PhD work, aiming at getting a deep knowledge on the influence of magnetic field on MILs properties. The influence of the magnetic field on the molecular dynamics and structural rearrangement of MILs ionic network was assessed through a 1H-NMR technique. Through the 1H-NMR relaxometry analysis it was possible to estimate the self-diffusion profiles of two different model MILs, [Aliquat][FeCl4] and [P66614][FeCl4]. A comparative analysis was established between the behavior of magnetic and non-magnetic ionic liquids, MILs and ILs, to facilitate the perception of the magnetic field impact on MILs properties. In contrast to ILs, MILs show a specific relaxation mechanism, characterized by the magnetic dependence of their self-diffusion coefficients. MILs self-diffusion coefficients increased in the presence of magnetic field whereas ILs self-diffusion was not affected. In order to understand the reasons underlying the magnetic dependence of MILs self-diffusion, studies were performed to investigate the influence of the magnetic field on MILs’ viscosity. It was observed that the MIL´s viscosity decreases with the increase of the magnetic field, explaining the increase of MILs self-diffusion according to the modified Stokes- Einstein equation. Different gas and liquid transport studies were therefore performed aiming to determine the influence of the magnetic behavior of MILs on solute transport through SMILMs. Gas permeation studies were performed using pure CO2 andN2 gas streams and air, using a series of phosphonium cation based MILs, containing different paramagnetic anions. Transport studies were conducted in the presence and absence of magnetic field at a maximum intensity of 1.5T. The results revealed that gas permeability increased in the presence of the magnetic field, however, without affecting the membrane selectivity. The increase of gas permeability through SMILMs was related to the decrease of the MILs viscosity under magnetic field conditions.(...)

A minimal version of a Protein Tyrosine Phosphatase

Phosphatase and tensin homologue (PTEN) protein belongs to the family of protein tyrosine phos-phatase. Mutations on the phosphatase and tensin homologue (PTEN) protein are highly observed in diverse types of human tumors, being mostly identified on the phosphatase domain of the protein. Although PTEN is a modular protein composed by a phosphatase domain and a C2 domain for mem-brane anchoring, this work aimed at developing a minimal version of PTEN´s phosphatase domain. The minimal version (Small Domain) comprises a 28 residue peptide, with the PTEN 8-mer catalytic peptide accommodated between a α-helix and β-turn as observed in PTEN native structure. Firstly, a de novo prediction of the Small Domain´s secondary structure was carried out by molecular modeling tools. The stability of the predicted structures were then evaluated by Molecular Dynamics. Automated molecular docking of PTEN natural substrate PIP3, its analogue (Inositol) and a PTEN inhibitor (L-tar-tare) were performed with the modeled structure, and PTEN used as a positive control. The gene en-coding for Small Domain was designed and cloned into an expression vector at N-terminal of Green Fluorescence Protein (GFP) encoding gene. The fusion protein was then expressed in Escherichia coli cells. Different expression conditions have been explored for the production of the fusion protein to minimize the formation of inclusion bodies.

A breakthrough on Amanita phalloides poisoning: an effective antidotal effect by polymyxin B

Amanita phalloides is responsible for more than 90 % of mushroom-related fatalities, and no effective antidote is available. a-Amanitin, the main toxin of A. phalloides, inhibits RNA polymerase II (RNAP II), causing hepatic and kidney failure. In silico studies included docking and molecular dynamics simulation coupled to molecular mechanics with generalized Born and surface area method energy decomposition on RNAP II. They were performed with a clinical drug that shares chemical similarities to a-amanitin, polymyxin B. The results show that polymyxin B potentially binds to RNAP II in the same interface of a-amanitin, preventing the toxin from binding to RNAP II. In vivo, the inhibition of the mRNA transcripts elicited by a-amanitin was efficiently reverted by polymyxin B in the kidneys. Moreover, polymyxin B significantly decreased the hepatic and renal a-amanitin-induced injury as seen by the histology and hepatic aminotransferases plasma data. In the survival assay, all animals exposed to a-amanitin died within 5 days, whereas 50 % survived up to 30 days when polymyxin B was administered 4, 8, and 12 h post-a-amanitin. Moreover, a single dose of polymyxin B administered concomitantly with a-amanitin was able to guarantee 100 % survival. Polymyxin B protects RNAP II from inactivation leading to an effective prevention of organ damage and increasing survival in a-amanitin-treated animals. The present use of clinically relevant concentrations of an already human-use-approved drug prompts the use of polymyxin B as an antidote for A. phalloides poisoning in humans.

New self-assembled supramolecular hydrogels based on dehydropeptides

Supramolecular hydrogels rely on small molecules that self-assemble in water as a result of the cooperative effect of several relatively weak intermolecular interactions. Peptide-based low molecular weight hydrogelators have attracted enormous interest owing to the simplicity of small molecules combined with the versatility and biocompatibility of peptides. In this work, naproxen, a well known non-steroidal anti-inflammatory drug, was N-conjugated with various dehydrodipeptides to give aromatic peptide amphiphiles that resist proteolysis. Molecular dynamics simulations were used to obtain insight into the underlying molecular mechanism of self-assembly and to rationalize the design of this type of hydrogelators. The results obtained were in excellent agreement with the experimental observations. Only dehydrodipeptides having at least one aromatic amino acid gave hydrogels. The characterization of the hydrogels was carried out using transmission electron microscopy (TEM), circular dichroism (CD), fluorescence spectroscopy and also rheological assays.

Synthesis of iminosugars through diastereo/enantioselective diels-alder cycloaddition

Tese de Doutoramento em Ciências (Especialidade em Química)

Nanocomposites: the contribution of multiscale modelling to improve the dispersion levels

Tese de Doutoramento em Ciência e Engenharia de Polímeros e Compósitos