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.

Template-based fluoroethylenepropylene piezoelectrets with tubular channels for transducer applications

We describe the concept, the fabrication, and the most relevant properties of a piezoelectric-polymer system: Two fluoroethylenepropylene (FEP) films with good electret properties are laminated around a specifically designed and prepared polytetrafluoroethylene (PTFE) template at 300 degrees C. After removing the PTFE template, a two-layer FEP film with open tubular channels is obtained. For electric charging, the two-layer FEP system is subjected to a high electric field. The resulting dielectric barrier discharges inside the tubular channels yield a ferroelectret with high piezoelectricity. d(33) coefficients of up to 160 pC/N have already been achieved on the ferroelectret films. After charging at suitable elevated temperatures, the piezoelectricity is stable at temperatures of at least 130 degrees C. Advantages of the transducer films include ease of fabrication at laboratory or industrial scales, a wide range of possible geometrical and processing parameters, straightforward control of the uniformity of the polymer system, flexibility, and versatility of the soft ferroelectrets, and a large potential for device applications e.g., in the areas of biomedicine, communications, production engineering, sensor systems, environmental monitoring, etc.

Synchronization of Discrete-Time Chaotic Systems in Bandlimited Channels

Over the last couple of decades, many methods for synchronizing chaotic systems have been proposed with communications applications in view. Yet their performance has proved disappointing in face of the nonideal character of usual channels linking transmitter and receiver, that is, due to both noise and signal propagation distortion. Here we consider a discrete-time master-slave system that synchronizes despite channel bandwidth limitations and an allied communication system. Synchronization is achieved introducing a digital filter that limits the spectral content of the feedback loop responsible for producing the transmitted signal. Copyright (C) 2009 Marcio Eisencraft et al.

Two formation channels of ultra-compact dwarf galaxies in Hickson compact groups

Context. The formation of ultra-compact dwarf galaxies (UCDs) is believed to be driven by interaction, and UCDs are abundant in the cores of galaxy clusters, environments that mark the end-point of galaxy evolution. Nothing is known about the properties of UCDs in compact groups of galaxies, environments where most of galaxy evolution and interaction is believed to occur and where UCDs in an intermediate stage in their evolution may be expected. Aims. The main goal of this study is to detect and characterize, for the first time, the UCD population of compact groups of galaxies. For that, two nearby groups in different evolutionary stages, HCG22 and HCG90, were targeted. Methods. We selected about 40 UCD candidates from pre-existing photometry of both groups, and obtained spectra of these candidates using the VLT FORS2 instrument in MXU mode. Archival HST/ACS imaging was used to measure their structural parameters. Results. We detect 16 and 5 objects belonging to HCG22 and HCG90, respectively, covering the magnitude range -10.0 > M(R) > -11.5 mag. Their integrated colours are consistent with old ages covering a broad range in metallicities (metallicities confirmed by the spectroscopic measurements). Photometric mass estimates put 4 objects in HCG90 and 9 in HCG22 in the mass range of UCDs (> 2 x 10(6) M(circle dot)) for an assumed age of 12Gyr. These UCDs are on average 2-3 times larger than the typical size of Galactic GCs, covering a range of 2 less than or similar to r(h) less than or similar to 21 pc. The UCDs in HCG22 are more concentrated around the central galaxy than in HCG90, at the 99% confidence level. They cover a broad range in [alpha/Fe] abundances from sub-to super-solar. The spectra of 3 UCDs (2 in HCG22, 1 in HCG90) show tentative evidence of intermediate age stellar populations. The clearest example is the largest and most massive UCD (similar to 10(7) M(circle dot)) in our sample, which is detected in HCG22. Its properties are most consistent with a stripped dwarf galaxy nucleus. We calculate the specific frequency (S(N)) of UCDs for both groups, finding that HCG22 has about three times higher S(N) than HCG90. Conclusions. The ensemble properties of the detected UCDs supports two co-existing formation channels: a star cluster origin (low-luminosity, compact sizes, old ages, super-solar alpha/Fe), and an origin as tidally stripped dwarf nuclei (more extended and younger stellar populations). Our results imply that the UCDs detected in both groups do not, in their majority, originate from relatively recent galaxy interactions. Most of the detected UCDs have likely been brought into the group along with their host galaxies.

ABNORMAL EXPRESSION OF VOLTAGE-GATED SODIUM CHANNELS Nav1.7, Nav1.3 AND Nav1.8 IN TRIGEMINAL NEURALGIA

Voltage-gated sodium channels have been implicated in acute and chronic neuropathic pain. Among subtypes, Nav1.7 single mutations can cause congenital indifference to pain or chronic neuropathic pain syndromes, including paroxysmal ones. This channel is co-expressed with Nav1.8, which sustains the initial action potential; Nav1.3 is an embrionary channel which is expressed in neurons after injury, as in neuropathic conditions. Few studies are focused on the expression of these molecules in human tissues having chronic pain. Trigeminal neuralgia (TN) is an idiopathic paroxysmal pain treated with sodium channel blockers. The aim of this study was to investigate the expression of Nav1.3, Nav1.7 and Nav1.8 by RT-PCR in patients with TN, compared to controls. The gingival tissue was removed from the correspondent trigeminal area affected. We found that Nav1.7 was downregulated in TN (P=0.017) and Nav1.3 was upregulated in these patients (P=0.043). We propose a physiopathological mechanism for these findings. Besides vascular compression of TN, this disease might be also a channelopathy. (C) 2009 IBRO. Published by Elsevier Ltd. All rights reserved.

A Class of Channels Resulting in Ill-Convergence for CMA in Decision Feedback Equalizers

This paper analyzes the convergence of the constant modulus algorithm (CMA) in a decision feedback equalizer using only a feedback filter. Several works had already observed that the CMA presented a better performance than decision directed algorithm in the adaptation of the decision feedback equalizer, but theoretical analysis always showed to be difficult specially due to the analytical difficulties presented by the constant modulus criterion. In this paper, we surmount such obstacle by using a recent result concerning the CM analysis, first obtained in a linear finite impulse response context with the objective of comparing its solutions to the ones obtained through the Wiener criterion. The theoretical analysis presented here confirms the robustness of the CMA when applied to the adaptation of the decision feedback equalizer and also defines a class of channels for which the algorithm will suffer from ill-convergence when initialized at the origin.

Chaotic Synchronization in Discrete-Time Systems Connected by Bandlimited Channels

Due to the broadband characteristic of chaotic signals, many of the methods that have been proposed for synchronizing chaotic systems do not usually present a satisfactory performance when applied to bandlimited communication channels. Here, the effects of bandwidth limitations imposed by the channel on the synchronous solution of a discrete-time chaotic master-slave network are investigated. The discrete-time system considered in this study is the Henon map. It is analytically shown that synchronism can be achieved in such a network by introducing a digital filter in the feedback loop responsible for generating the chaotic signal that will be sent to the slave node. Numerical simulations relating the filter parameters, such as its order and cut-off frequency, to the maximum Lyapunov exponent of the master node, which determines if the transmitted signal is chaotic or not, are also presented. These results can be useful for practical communication schemes based on chaos.

Acidosis induces relaxation mediated by nitric oxide and potassium channels in rat thoracic aorta

We investigated the mechanism by which extracellular acidification promotes relaxation in rat thoracic aorta. The relaxation response to HCl-induced extracellular acidification (7.4 to 6.5) was measured in aortic rings pre-contracted with phenylephrine (Phe, 10(-6) M) or KCl (45 mM). The vascular reactivity experiments were performed in endothelium-intact and denuded rings, in the presence or absence of indomethacin (10(-5) M), L-NAME (10(-4) M), apamin (10(-6) M), and glibenclamide (10(-5) M). The effect of extracellular acidosis (pH 7.0 and 6.5) on nitric oxide (NO) production was evaluated in isolated endothelial cells loaded with diaminofluorescein-FM diacetate (DAF-FM DA, 5 mu M). The extracellular acidosis failed to induce any changes in the vascular tone of aortic rings pre-contracted with KCl, however, it caused endothelium-dependent and independent relaxation in rings pre-contracted with Phe. This acidosis induced-relaxation was inhibited by L-NAME, apamin, and glibenclamide, but not by indomethacin. The acidosis (pH 7.0 and 6.5) also promoted a time-dependent increase in the NO production by the isolated endothelial cells. These results suggest that extracellular acidosis promotes vasodilation mediated by NO, K(ATP) and SK(Ca), and maybe other K(+) channels in isolated rat thoracic aorta. (C) 2011 Elsevier B.V. All rights reserved.

Characterization of calcium-activated chloride channels in patches excised from the dendritic knob of mammalian olfactory receptor neurons

We investigated the properties of calcium-activated chloride channels in inside-out membrane patches from the dendritic knobs of acutely dissociated rat olfactory receptor neurons. Patches typically contained large calcium-activated currents, with total conductances in the range 30-75 nS. The dose response curve for calcium exhibited an EC50 of about 26 mu M. In symmetrical NaCl solutions, the current-voltage relationship reversed at 0 mV and was linear between -80 and +70 mV. When the intracellular NaCl concentration was progressively reduced from 150 to 25 mM, the reversal potential changed in a manner consistent with a chloride-selective conductance. Indeed, modeling these data with the Goldman-Hodgkin-Katz equation revealed a P-Na/P-Cl of 0.034. The halide permeability sequence was P-Cl > P-F > P-I > P-Br indicating that permeation through the channel was dominated by ion binding sites with a high field strength. The channels were also permeable to the large organic anions, SCN-, acetate(-), and gluconate(-), with the permeability sequence P-Cl > P-SCN > gluconaie. Significant permeation to gluconate ions suggested that the channel pore had a minimum diameter of at least 5.8 Angstrom.

NMR structure of C-terminally tagged gramicidin channels

A biotin group was covalently attached to the C terminus of gramicidin A (gA) through a linker arm comprising a glycine residue with either one (gAXB) or two caproyl groups (gAXXB). High-resolution two-dimensional NMR spectroscopy was used to determine the structure of these modified gA analogues and [Lys(16)]gramicidin A (gA-Lys) in sodium dodecyl-d(25) sulphate micelles. Gated gA ion channels based on linking a receptor group to these gA analogues have been used recently as a component in a sensing device. The conformations of the gA backbones and amino acid side chains of lysinated gA and biotinylated gA in detergent micelles were found to be almost identical to that of native gA, i.e. that of an N-terminal to N-terminal (head to head) dimer formed by two right-handed, single-stranded beta(6.3) helices. The biotin tail of the gAXB and gAXXB and the lysine extremity of gA-Lys appeared to lie outside the micelle. Thus it appears that the covalent attachment of functional groups to the C terminus of gA does not disrupt the peptide's helical configuration. Further, single channel measurements of all three gA analogues showed that functioning ion channels were preserved within a membrane environment. (C) 1999 Elsevier Science B.V. All rights reserved.

Photolytic manipulation of [Ca2+](i) reveals slow kinetics of potassium channels underlying the afterhyperpolarization in hipppocampal pyramidal neurons

The identity of the potassium channel underlying the slow, apamin-insensitive component of the afterhyperpolarization current (sl(AHP)) remains unknown. We studied sl(AHP) in CA1 pyramidal neurons using simultaneous whole-cell recording, calcium fluorescence imaging, and flash photolysis of caged compounds. Intracellular calcium concentration ([Ca2+](i)) peaked earlier and decayed more rapidly than sl(AHP). Loading cells with low concentrations of the calcium chelator EGTA slowed the activation and decay of sl(AHP). In the presence of EGTA, intracellular calcium decayed with two time constants. When [Ca2+](i) was increased rapidly after photolysis of DM-Nitrophen, both apamin-sensitive and apamin-insensitive outward currents were activated. The apamin-sensitive current activated rapidly (<20 msec), whereas the apamin-insensitive current activated more slowly (180 msec). The apamin-insensitive current was reduced by application of serotonin and carbachol, confirming that it was caused by sl(AHP) channels. When [Ca2+](i) was decreased rapidly via photolysis of diazo-2, the decay of sl(AHP) was similar to control (1.7 sec). All results could be reproduced by a model potassium channel gated by calcium, suggesting that the channels underlying sl(AHP) have intrinsically slow kinetics because of their high affinity for calcium.

Structure-activity relationships of omega-conotoxins MVIIA, MVIIC and 14 loop splice hybrids at N and P/Q-type calcium channels

The omega-conotoxins are a set of structurally related, four-loop, six cysteine containing peptides, that have a range of selectivities for different subtypes of the voltage-sensitive calcium channel (VSCC). To investigate the basis of the selectivity displayed by these peptides, we have studied the binding affinities of two naturally occurring omega-conotoxins, MVIIA and MVIIC and a series of 14 MVIIA/MVIIC loop hybrids using radioligand binding assays for N and P/Q-type Ca2+ channels in rat brain tissue. A selectivity profile was developed from the ratio of relative potencies at N-type VSCCs (using [I-125]GVIA radioligand binding assays) and P/Q-type VSCCs (using [I-125]MVIIC radioligand binding assays). in these peptides, loops 2 and 4 make the greatest contribution to VSCC subtype selectivity, while the effects of loops 1 and 3 are negligible. Peptides with homogenous combinations of loop 2 and 4 display clear selectivity preferences, while those with heterogeneous combinations of loops 2 and 4 are less discriminatory. H-1 NMR spectroscopy revealed that the global folds of MVIIA, MVIIC and the 14 loop hybrid peptides were similar; however, several differences in local structure were identified. Based on the binding data and the 3D structures of MVIIA, GVIA and MVIIC, we have developed a preliminary pharmacophore based on the omega-conotoxin residues most Likely to interact with the N-type VSCC. (C) 1999 Academic Press.