Computational analyses on the structure-function relation in ion channels

Ion channels are protein molecules, embedded in the lipid bilayer of the cell membranes. They act as powerful sensing elements switching chemicalphysical stimuli into ion-fluxes. At a glance, ion channels are water-filled pores, which can open and close in response to different stimuli (gating), and one once open select the permeating ion species (selectivity). They play a crucial role in several physiological functions, like nerve transmission, muscular contraction, and secretion. Besides, ion channels can be used in technological applications for different purpose (sensing of organic molecules, DNA sequencing). As a result, there is remarkable interest in understanding the molecular determinants of the channel functioning. Nowadays, both the functional and the structural characteristics of ion channels can be experimentally solved. The purpose of this thesis was to investigate the structure-function relation in ion channels, by computational techniques. Most of the analyses focused on the mechanisms of ion conduction, and the numerical methodologies to compute the channel conductance. The standard techniques for atomistic simulation of complex molecular systems (Molecular Dynamics) cannot be routinely used to calculate ion fluxes in membrane channels, because of the high computational resources needed. The main step forward of the PhD research activity was the development of a computational algorithm for the calculation of ion fluxes in protein channels. The algorithm - based on the electrodiffusion theory - is computational inexpensive, and was used for an extensive analysis on the molecular determinants of the channel conductance. The first record of ion-fluxes through a single protein channel dates back to 1976, and since then measuring the single channel conductance has become a standard experimental procedure. Chapter 1 introduces ion channels, and the experimental techniques used to measure the channel currents. The abundance of functional data (channel currents) does not match with an equal abundance of structural data. The bacterial potassium channel KcsA was the first selective ion channels to be experimentally solved (1998), and after KcsA the structures of four different potassium channels were revealed. These experimental data inspired a new era in ion channel modeling. Once the atomic structures of channels are known, it is possible to define mathematical models based on physical descriptions of the molecular systems. These physically based models can provide an atomic description of ion channel functioning, and predict the effect of structural changes. Chapter 2 introduces the computation methods used throughout the thesis to model ion channels functioning at the atomic level. In Chapter 3 and Chapter 4 the ion conduction through potassium channels is analyzed, by an approach based on the Poisson-Nernst-Planck electrodiffusion theory. In the electrodiffusion theory ion conduction is modeled by the drift-diffusion equations, thus describing the ion distributions by continuum functions. The numerical solver of the Poisson- Nernst-Planck equations was tested in the KcsA potassium channel (Chapter 3), and then used to analyze how the atomic structure of the intracellular vestibule of potassium channels affects the conductance (Chapter 4). As a major result, a correlation between the channel conductance and the potassium concentration in the intracellular vestibule emerged. The atomic structure of the channel modulates the potassium concentration in the vestibule, thus its conductance. This mechanism explains the phenotype of the BK potassium channels, a sub-family of potassium channels with high single channel conductance. The functional role of the intracellular vestibule is also the subject of Chapter 5, where the affinity of the potassium channels hEag1 (involved in tumour-cell proliferation) and hErg (important in the cardiac cycle) for several pharmaceutical drugs was compared. Both experimental measurements and molecular modeling were used in order to identify differences in the blocking mechanism of the two channels, which could be exploited in the synthesis of selective blockers. The experimental data pointed out the different role of residue mutations in the blockage of hEag1 and hErg, and the molecular modeling provided a possible explanation based on different binding sites in the intracellular vestibule. Modeling ion channels at the molecular levels relates the functioning of a channel to its atomic structure (Chapters 3-5), and can also be useful to predict the structure of ion channels (Chapter 6-7). In Chapter 6 the structure of the KcsA potassium channel depleted from potassium ions is analyzed by molecular dynamics simulations. Recently, a surprisingly high osmotic permeability of the KcsA channel was experimentally measured. All the available crystallographic structure of KcsA refers to a channel occupied by potassium ions. To conduct water molecules potassium ions must be expelled from KcsA. The structure of the potassium-depleted KcsA channel and the mechanism of water permeation are still unknown, and have been investigated by numerical simulations. Molecular dynamics of KcsA identified a possible atomic structure of the potassium-depleted KcsA channel, and a mechanism for water permeation. The depletion from potassium ions is an extreme situation for potassium channels, unlikely in physiological conditions. However, the simulation of such an extreme condition could help to identify the structural conformations, so the functional states, accessible to potassium ion channels. The last chapter of the thesis deals with the atomic structure of the !- Hemolysin channel. !-Hemolysin is the major determinant of the Staphylococcus Aureus toxicity, and is also the prototype channel for a possible usage in technological applications. The atomic structure of !- Hemolysin was revealed by X-Ray crystallography, but several experimental evidences suggest the presence of an alternative atomic structure. This alternative structure was predicted, combining experimental measurements of single channel currents and numerical simulations. This thesis is organized in two parts, in the first part an overview on ion channels and on the numerical methods adopted throughout the thesis is provided, while the second part describes the research projects tackled in the course of the PhD programme. The aim of the research activity was to relate the functional characteristics of ion channels to their atomic structure. In presenting the different research projects, the role of numerical simulations to analyze the structure-function relation in ion channels is highlighted.

Caratterizzazione di membrane da nanofiltrazione

Nanofiltration (NF) is a pressure-driven membrane process, intermediate between reverse osmosis and ultrafiltration. Commercially available polymeric membranes have been used in a wide range of applications, such as drinking, process industry and waste water treatment. For all the applications requiring high stability and harsh washing procedures inorganic membranes are preferred due to their high chemical inertia. Typically, γ – Al2O3 as well as TiO2 and ZrO2 selective layers are used; the latter show higher chemical stability in a wide range of pH and temperatures. In this work the experimental characterization of two different type of membrane has been performed in order to investigate permeation properties, separation performance and efficiency with aqueous solutions containing strong inorganic electrolytes. The influence of salt concentration and feed pH as well as the role of concentration polarization and electrolyte type on the membrane behavior are investigated. Experimentation was performed testing a multi–layer structured NF membrane in α-Al2O3, TiO2 and ZrO2, and a polymeric membrane, in polyamide supported on polysulfone, with binary aqueous solutions containing NaCl, Na2SO4 or CaCl2; the effect of salt composition and pH in the feed side was studied both on flux and salt rejection. All the NF experimental data available for the two membranes were used to evaluate the volumetric membrane charge (X) corresponding to each operative conditions investigated, through the Donnan Steric Pore Model and Dielectric Exclusion (DSPM&DE). The results obtained allow to understand which are the main phenomena at the basis of the different behaviors observed.

Caratterizzazione di membrane inorganiche per la separazione di idrogeno da gas di reforming

The work of this thesis has been focused on the characterisation of inorganic membranes for the hydrogen purification from steam reforming gas. Composite membranes based on porous inorganic supports coated with palladium silver alloys and ceramic membranes have been analysed. A brief resume of theoretical laws governing transport of gases through dense and porous inorganic membranes and an overview on different methods to prepare inorganic membranes has been also reported. A description of the experimental apparatus used for the characterisation of gas permeability properties has been reported. The device used permits to evaluate transport properties in a wide range of temperatures (till 500°C) and pressures (till 15 bar). Data obtained from experimental campaigns reveal a good agreement with Sievert law for hydrogen transport through dense palladium based membranes while different transport mechanisms, such as Knudsen diffusion and Hagen-Poiseuille flow, have been observed for porous membranes and for palladium silver alloy ones with pinholes in the metal layer. Mixtures permeation experiments reveal also concentration polarisation phenomena and hydrogen permeability reduction due to carbon monoxide adsorption on metal surface.

Characterization and modeling of the barrier properties in nanostructured systems

The object of the present study is the process of gas transport in nano-sized materials, i.e. systems having structural elements of the order of nanometers. The aim of this work is to advance the understanding of the gas transport mechanism in such materials, for which traditional models are not often suitable, by providing a correct interpretation of the relationship between diffusive phenomena and structural features. This result would allow the development new materials with permeation properties tailored on the specific application, especially in packaging systems. The methods used to achieve this goal were a detailed experimental characterization and different simulation methods. The experimental campaign regarded the determination of oxygen permeability and diffusivity in different sets of organic-inorganic hybrid coatings prepared via sol-gel technique. The polymeric samples coated with these hybrid layers experienced a remarkable enhancement of the barrier properties, which was explained by the strong interconnection at the nano-scale between the organic moiety and silica domains. An analogous characterization was performed on microfibrillated cellulose films, which presented remarkable barrier effect toward oxygen when it is dry, while in the presence of water the performance significantly drops. The very low value of water diffusivity at low activities is also an interesting characteristic which deals with its structural properties. Two different approaches of simulation were then considered: the diffusion of oxygen through polymer-layered silicates was modeled on a continuum scale with a CFD software, while the properties of n-alkanthiolate self assembled monolayers on gold were analyzed from a molecular point of view by means of a molecular dynamics algorithm. Modeling transport properties in layered nanocomposites, resulting from the ordered dispersion of impermeable flakes in a 2-D matrix, allowed the calculation of the enhancement of barrier effect in relation with platelets structural parameters leading to derive a new expression. On this basis, randomly distributed systems were simulated and the results were analyzed to evaluate the different contributions to the overall effect. The study of more realistic three-dimensional geometries revealed a prefect correspondence with the 2-D approximation. A completely different approach was applied to simulate the effect of temperature on the oxygen transport through self assembled monolayers; the structural information obtained from equilibrium MD simulations showed that raising the temperature, makes the monolayer less ordered and consequently less crystalline. This disorder produces a decrease in the barrier free energy and it lowers the overall resistance to oxygen diffusion, making the monolayer more permeable to small molecules.

Caratterizzazione di membrane metalliche e ionomeriche per applicazioni energetiche

The work of this thesis has been focused on the characterization of metallic membranes for the hydrogen purification from steam reforming process and also of perfluorosulphonic acid ionomeric (PFSI) membranes suitable as electrolytes in fuel cell applications. The experimental study of metallic membranes was divided in three sections: synthesis of palladium and silver palladium coatings on porous ceramic support via electroless deposition (ELD), solubility and diffusivity analysis of hydrogen in palladium based alloys (temperature range between 200 and 400 °C up to 12 bar of pressure) and permeation experiments of pure hydrogen and mixtures containing, besides hydrogen, also nitrogen and methane at high temperatures (up to 600 °C) and pressures (up to 10 bar). Sequential deposition of palladium and silver on to porous alumina tubes by ELD technique was carried out using two different procedures: a stirred batch and a continuous flux method. Pure palladium as well as Pd-Ag membranes were produced: the Pd-Ag membranes’ composition is calculated to be close to 77% Pd and 23% Ag by weight which was the target value that correspond to the best performance of the palladium-based alloys. One of the membranes produced showed an infinite selectivity through hydrogen and relatively high permeability value and is suitable for the potential use as a hydrogen separator. The hydrogen sorption in silver palladium alloys was carried out in a gravimetric system on films produced by ELD technique. In the temperature range inspected, up to 400°C, there is still a lack in literature. The experimental data were analyzed with rigorous equations allowing to calculate the enthalpy and entropy values of the Sieverts’ constant; the results were in very good agreement with the extrapolation made with literature data obtained a lower temperature (up to 150 °C). The information obtained in this study would be directly usable in the modeling of hydrogen permeation in Pd-based systems. Pure and mixed gas permeation tests were performed on Pd-based hydrogen selective membranes at operative conditions close to steam-reforming ones. Two membranes (one produced in this work and another produced by NGK Insulators Japan) showed a virtually infinite selectivity and good permeability. Mixture data revealed the existence of non negligible resistances to hydrogen transport in the gas phase. Even if the decrease of the driving force due to polarization concentration phenomena occurs, in principle, in all membrane-based separation systems endowed with high perm-selectivity, an extensive experimental analysis lack, at the moment, in the palladium-based membrane process in literature. Moreover a new procedure has been introduced for the proper comparison of the mass transport resistance in the gas phase and in the membrane. Another object of study was the water vapor sorption and permeation in PFSI membranes with short and long side chains was also studied; moreover the permeation of gases (i.e. He, N2 and O2) in dry and humid conditions was considered. The water vapor sorption showed strong interactions between the hydrophilic groups and the water as revealed from the hysteresis in the sorption-desorption isotherms and thermo gravimetric analysis. The data obtained were used in the modeling of water vapor permeation, that was described as diffusion-reaction of water molecules, and in the humid gases permeation experiments. In the dry gas experiments the permeability and diffusivity was found to increase with temperature and with the equivalent weight (EW) of the membrane. A linear correlation was drawn between the dry gas permeability and the opposite of the equivalent weight of PFSI membranes, based on which the permeability of pure PTFE is retrieved in the limit of high EW. In the other hand O2 ,N2 and He permeability values was found to increase significantly, and in a similar fashion, with water activity. A model that considers the PFSI membrane as a composite matrix with a hydrophilic and a hydrophobic phase was considered allowing to estimate the variation of gas permeability with relative humidity on the basis of the permeability in the dry PFSI membrane and in pure liquid water.

Fluid dynamic analysis and modelling of industrial chemical equipment

The research is aimed at contributing to the identification of reliable fully predictive Computational Fluid Dynamics (CFD) methods for the numerical simulation of equipment typically adopted in the chemical and process industries. The apparatuses selected for the investigation, specifically membrane modules, stirred vessels and fluidized beds, were characterized by a different and often complex fluid dynamic behaviour and in some cases the momentum transfer phenomena were coupled with mass transfer or multiphase interactions. Firs of all, a novel modelling approach based on CFD for the prediction of the gas separation process in membrane modules for hydrogen purification is developed. The reliability of the gas velocity field calculated numerically is assessed by comparison of the predictions with experimental velocity data collected by Particle Image Velocimetry, while the applicability of the model to properly predict the separation process under a wide range of operating conditions is assessed through a strict comparison with permeation experimental data. Then, the effect of numerical issues on the RANS-based predictions of single phase stirred tanks is analysed. The homogenisation process of a scalar tracer is also investigated and simulation results are compared to original passive tracer homogenisation curves determined with Planar Laser Induced Fluorescence. The capability of a CFD approach based on the solution of RANS equations is also investigated for describing the fluid dynamic characteristics of the dispersion of organics in water. Finally, an Eulerian-Eulerian fluid-dynamic model is used to simulate mono-disperse suspensions of Geldart A Group particles fluidized by a Newtonian incompressible fluid as well as binary segregating fluidized beds of particles differing in size and density. The results obtained under a number of different operating conditions are compared with literature experimental data and the effect of numerical uncertainties on axial segregation is also discussed.

Numerical simulation of ion transport through ion channels and solid-state nanopores

Ion channels are pore-forming proteins that regulate the flow of ions across biological cell membranes. Ion channels are fundamental in generating and regulating the electrical activity of cells in the nervous system and the contraction of muscolar cells. Solid-state nanopores are nanometer-scale pores located in electrically insulating membranes. They can be adopted as detectors of specific molecules in electrolytic solutions. Permeation of ions from one electrolytic solution to another, through a protein channel or a synthetic pore is a process of considerable importance and realistic analysis of the main dependencies of ion current on the geometrical and compositional characteristics of these structures are highly required. The project described by this thesis is an effort to improve the understanding of ion channels by devising methods for computer simulation that can predict channel conductance from channel structure. This project describes theory, algorithms and implementation techniques used to develop a novel 3-D numerical simulator of ion channels and synthetic nanopores based on the Brownian Dynamics technique. This numerical simulator could represent a valid tool for the study of protein ion channel and synthetic nanopores, allowing to investigate at the atomic-level the complex electrostatic interactions that determine channel conductance and ion selectivity. Moreover it will provide insights on how parameters like temperature, applied voltage, and pore shape could influence ion translocation dynamics. Furthermore it will help making predictions of conductance of given channel structures and it will add information like electrostatic potential or ionic concentrations throughout the simulation domain helping the understanding of ion flow through membrane pores.

Chitosan based hydrogels for transmucosal drug delivery

The aim of this thesis was the formulation of new chitosan based delivery systems for transmucosal drug administration. Transmucosal routes, such as buccal, vaginal and nasal routes, allow the circumvention of the hepatic first pass metabolism and avoid the gastrointestinal chemical and enzymatic degradations. Moreover, transmucosal drug administration can allow to avoid pain or discomfort caused by injections, when drugs are administered through parenteral routes, thus increasing patient compliance. On the other side, the major disadvantage of transmucosal drug administration is represented by the presence of biological fluids and mucus that can remove drug systems from the application site, thus reducing the contact time between drug and mucosa and consequently, decreasing drug bioavailability. For this reason, in this study, the investigation of chitosan delivery systems as mucoadhesive formulations able to increase drugs residence time and to improve their bioavailability, was taken into account. In the paper 1, buccal films based on chitosan-gelatin complexes were prepared and loaded with propranolol hydrochloride. The complexes were characterized and studied in order to evaluate their physical- chemical properties and their ability to release the drug and to allow its permeation through buccal mucosa. In the paper 2, vaginal inserts based on chitosan/alginate complexes were formulated for local delivery of chlorhexidine digluconate. Tests to evaluate the interaction between the polymers and to study drug release properties were performed, as well as the determination of antimicrobial activity against the patogens responsible of vaginitis and candidosis. In the project 3, chitosan based nanoparticles containing cyclodextrin and other excipients, with the capacity to modify insulin bioavailabity were formulated for insulin nasal delivery. Nanoparticles were characterized in terms of size, stability and drug release. Moreover, in vivo tests were performed in order to study the hypoglycemic reduction in rats blood samples.

Polymeric membranes for CO2 separation: effect of aging, humidity and facilitated transport

Polymeric membranes represent a promising technology for gas separation processes, thanks to low costs, reduced energy consumption and limited waste production. The present thesis aims at studying the transport properties of two membrane materials, suitable for CO2 purification applications. In the first part, a polyimide, Matrimid 5218, has been throughout investigated, with particular reference to the effect of thermal treatment, aging and the presence of water vapor in the gas transport process. Permeability measurements showed that thermal history affects relevantly the diffusion of gas molecules across the membrane, influencing also the stability of the separation performances. Subsequently, the effect of water on Matrimid transport properties has been characterized for a wide set of incondensable penetrants. A monotonous reduction of permeability took place at increasing the water concentration within the polymer matrix, affecting the investigated gaseous species to the same extent, despite the different thermodynamic and kinetic features. In this view, a novel empirical model, based on the Free Volume Theory, has been proposed to qualitatively describe the phenomenon. Moreover, according to the accurate representation of the experimental data, the suggested approach has been combined with a more rigorous thermodynamic tool (NELF Model), allowing an exhaustive description of water influence on the single parameters contributing to the gas permeation across the membrane. In the second part, the study has focused on the synthesis and characterization of facilitated transport membranes, able to achieving outstanding separation performances thanks to the chemical enhancement of CO2 permeability. In particular, the transport properties have been investigated for high pressure CO2 separation applications and specific solutions have been proposed to solve stability issues, frequently arising under such severe conditions. Finally, the effect of different process parameters have been investigated, aiming at the identification of the optimal conditions capable to maximize the separation performance.