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.

Pattern recognition analysis on heavy ion reaction data

One of the problems in the analysis of nucleus-nucleus collisions is to get information on the value of the impact parameter b. This work consists in the application of pattern recognition techniques aimed at associating values of b to groups of events. To this end, a support vec- tor machine (SVM) classifier is adopted to analyze multifragmentation reactions. This method allows to backtracing the values of b through a particular multidimensional analysis. The SVM classification con- sists of two main phase. In the first one, known as training phase, the classifier learns to discriminate the events that are generated by two different model:Classical Molecular Dynamics (CMD) and Heavy- Ion Phase-Space Exploration (HIPSE) for the reaction: 58Ni +48 Ca at 25 AMeV. To check the classification of events in the second one, known as test phase, what has been learned is tested on new events generated by the same models. These new results have been com- pared to the ones obtained through others techniques of backtracing the impact parameter. Our tests show that, following this approach, the central collisions and peripheral collisions, for the CMD events, are always better classified with respect to the classification by the others techniques of backtracing. We have finally performed the SVM classification on the experimental data measured by NUCL-EX col- laboration with CHIMERA apparatus for the previous reaction.

Salt stress responses in pear and quince: physiological and molecular aspects

The productivity of agricultural crops is seriously limited by salinity. This problem is rapidly increasing, particularly in irrigated lands. Like almost all the fruit tree species, Pyrus communis is generally considered a salt sensitive species, but only little information is available on its behavior under saline conditions. Previous studies, carried out in the Department of Fruit Tree and Woody Plant Science (University of Bologna), focused their attention on pear and quince salt stress responses to understand which rootstock would be the most suitable for pear in order to tolerate a salt stress condition. It has been reported that pear and quince have different ability in the uptake, translocation and accumulation of chloride (Cl-) and sodium (Na+) ions, when plants were irrigated for one season with saline water (5 dS/m). The aim of the present work was to deepen these aspects and investigate salt stress responses in pear and quince. Two different experiments have been performed: a “short-term” trial in a growth chamber and a “long-term” experiment in the open field. In the short-term experiment, three different genotypes usually adopted as pear rootstocks (MC, BA29 and Farold®40) and the pear variety Abbé Fétel own rooted have been compared under salt stress conditions. The trial was performed in a hydroponic culture system, applying a 90 mM NaCl stress to half of the plants, after five weeks of normal growth in Hoagland’s solution. During the three-weeks of salt stress treatment, physiological, mineral and molecular analyses were performed in order to monitor, for each genotype, the development of the salt stress responses in comparison with the corresponding “unstressed” plants. Farold®40 and Abbé Fétel own rooted showed the onset of leaf necrosis, due to salt toxicity, one week before quinces. Moreover, quinces displayed a significant delay in premature senescence of old leaves, while pears emerged for their ability to regenerate new leaves from apparently dead foliage with the salt stress still running. Physiological measurements, such as shoots length, chlorophyll (Chl) content, and photosynthesis, have been carried out and revealed that pears exhibited a significant reduction in water content and a wilting aspect, while for quinces a decrease in Chl content and a growth slowdown were observed. At the end of the trial, all plants were collected and organs separated for dry weight estimation and mineral analyses (Cu, Fe, Mn, Zn Mg, Ca, K, Na and Cl). Mineral contents have been affected by salinity; same macro/micro nutrients were altered in some organs or relocated within the plant. This plant response could have partially contributed to face the salt stress. Leaves and roots have been harvested for molecular analyses at four different times during stress conditions. Molecular analyses consisted of the gene expression study of three main ion transporters, well known in Arabidopsis thaliana as salt-tolerance determinants in the “SOS” pathway: NHX1 (tonoplast Na+/H+ antiporter), SOS1 (plasmalemma Na+/H+ antiporter) and HKT1 (K+ high-affinity and Na+ low-affinity transporter). These studies showed that two quince rootstocks adopted different responsive mechanisms to NaCl stress. BA29 increased its Na+ sequestration activity into leaf vacuoles, while MC enhanced temporarily the same ability, but in roots. Farold®40, instead, exhibited increases in SOS1 and HKT1 expression mainly at leaf level in the attempt to retrieve Na+ from xylem, while Abbé Fétel differently altered the expression of these genes in roots. Finally, each genotype showed a peculiar response to salt stress that was the sum of its ability in Na+ exclusion, osmotic tolerance and tissue tolerance. In the long-term experiment, potted trees of the pear variety Abbé Fétel grafted on different rootstocks (MC, BA29 and Farold®40), or own rooted and also rootstocks only were subjected to a salt stress through saline water irrigation with an electrical conductivity of 5 dS/m for two years. The purposes of this study were to evaluate salinity effects on physiological (shoot length, number of buds, photosynthesis, etc.) and yield parameters of cultivar Abbé Fétel in the different combinations and to determine the salt amount that pear is able to tolerate over the years. With this work, we confirmed the previous hypothesis that pear, despite being classified as a salt-sensitive fruit tree, can be cultivated for two years under saline water irrigation, without showing any salt toxicity symptoms or severe drawbacks on plant development and production. Among different combinations, Abbé Fétel grafted on MC resulted interesting for its peculiar behaviors under salt stress conditions. In the near future, further investigations on physiological and molecular aspects will be necessary to enrich and broaden the knowledge of salt stress responses in pear.

A hybrid technology for parallel recording of single ion channels

Hybrid technologies, thanks to the convergence of integrated microelectronic devices and new class of microfluidic structures could open new perspectives to the way how nanoscale events are discovered, monitored and controlled. The key point of this thesis is to evaluate the impact of such an approach into applications of ion-channel High Throughput Screening (HTS)platforms. This approach offers promising opportunities for the development of new classes of sensitive, reliable and cheap sensors. There are numerous advantages of embedding microelectronic readout structures strictly coupled to sensing elements. On the one hand the signal-to-noise-ratio is increased as a result of scaling. On the other, the readout miniaturization allows organization of sensors into arrays, increasing the capability of the platform in terms of number of acquired data, as required in the HTS approach, to improve sensing accuracy and reliabiity. However, accurate interface design is required to establish efficient communication between ionic-based and electronic-based signals. The work made in this thesis will show a first example of a complete parallel readout system with single ion channel resolution, using a compact and scalable hybrid architecture suitable to be interfaced to large array of sensors, ensuring simultaneous signal recording and smart control of the signal-to-noise ratio and bandwidth trade off. More specifically, an array of microfluidic polymer structures, hosting artificial lipid bilayers blocks where single ion channel pores are embededed, is coupled with an array of ultra-low noise current amplifiers for signal amplification and data processing. As demonstrating working example, the platform was used to acquire ultra small currents derived by single non-covalent molecular binding between alpha-hemolysin pores and beta-cyclodextrin molecules in artificial lipid membranes.