Peptide fingerprinting of the neurotoxic fractions isolated from the secretions of sea anemones Stichodactyla helianthus and Bunodosoma granulifera. New members of the APETx-like family identified by a 454 pyrosequencing approach

Sea anemones are known to contain a wide diversity of biologically active peptides, mostly unexplored according to recent peptidomic and transcriptomic studies. In the present work, the neurotoxic fractions from the exudates of Stichodactyla helianthus and Bunodosoma granulifera were analyzed by reversed-phase chromatography and mass spectrometry. The first peptide fingerprints of these sea anemones were assessed, revealing the largest number of peptide components (156) so far found in sea anemone species, as well as the richer peptide diversity of B. granulifera in relation to S. helianthus. The transcriptomic analysis of B. granulifera, performed by massive cDNA sequencing with 454 pyrosequencing approach allowed the discovery of five new APETx-like peptides (U-AITX-Bg1a-e - including the full sequences of their precursors for four of them), which together with type 1 sea anemone sodium channel toxins constitute a very distinguishable feature of studied sea anemone species belonging to genus Bunodosoma. The molecular modeling of these new APETx-like peptides showed a distribution of positively charged and aromatic residues in putative contact surfaces as observed in other animal toxins. On the other hand, they also showed variable electrostatic potentials, thus suggesting a docking onto their targeted channels in different spatial orientations. Moreover several crab paralyzing toxins (other than U-AITX-Bg1a-e), which induce a variety of symptoms in crabs, were isolated. Some of them presumably belong to new classes of crab-paralyzing peptide toxins, especially those with molecular masses below 2 kDa, which represent the smallest peptide toxins found in sea anemones. (C) 2011 Elsevier Inc. All rights reserved.

Hyperbolic-like estimates for higher order equations

The main goal of this paper is to derive long time estimates of the energy for the higher order hyperbolic equations with time-dependent coefficients. in particular, we estimate the energy in the hyperbolic zone of the extended phase space by means of a function f (t) which depends on the principal part and on the coefficients of the terms of order m - 1. Then we look for sufficient conditions that guarantee the same energy estimate from above in all the extended phase space. We call this class of estimates hyperbolic-like since the energy behavior is deeply depending on the hyperbolic structure of the equation. In some cases, these estimates produce a dissipative effect on the energy. (C) 2012 Elsevier Inc. All rights reserved.

An estimate on the fractal dimension of attractors of gradient-like dynamical systems

This paper is dedicated to estimate the fractal dimension of exponential global attractors of some generalized gradient-like semigroups in a general Banach space in terms of the maximum of the dimension of the local unstable manifolds of the isolated invariant sets, Lipschitz properties of the semigroup and the rate of exponential attraction. We also generalize this result for some special evolution processes, introducing a concept of Morse decomposition with pullback attractivity. Under suitable assumptions, if (A, A*) is an attractor-repeller pair for the attractor A of a semigroup {T(t) : t >= 0}, then the fractal dimension of A can be estimated in terms of the fractal dimension of the local unstable manifold of A*, the fractal dimension of A, the Lipschitz properties of the semigroup and the rate of the exponential attraction. The ingredients of the proof are the notion of generalized gradient-like semigroups and their regular attractors, Morse decomposition and a fine analysis of the structure of the attractors. As we said previously, we generalize this result for some evolution processes using the same basic ideas. (C) 2012 Elsevier Ltd. All rights reserved.

Form III-like conformation and Form I-like packing in a chloroform channel solvate of the diuretic drug chlortalidone

Chlortalidone (CTD) is an antihypertensive drug for which only two solid state phases have been structurally elucidated thus far. Here, we have prepared a chloroform solvate thereof, namely, CTD Form IV, and its structure was compared to those of Form I and Form III. Its two conformers exhibit a dual structural feature in relation to the antecedent polymorphs. Both CTD molecules of Form IV adopt a Form III-like conformation, which is featured, if the conformation of CTD Form I is used as a reference, by a rotation of about 90 degrees on the axis of the C-C bond bridging the substituted benzene and isoindolinyl rings. However, CTD Form IV assembles as in the Form I crystal packing despite the different stacking fashion of their centrosymmetric dimers. In contrast to Form I, there is no offset stacking in Form IV, which forces a bend of ca. 24 degrees between the planes passing through the isoindolinyl moieties of two [100]-stacked dimers. Chloroform molecules at a maximum stoichiometry of 0.25 mol per mol of the drug play a stabilizing role in the assembly of Form IV by filling the channels formed on the crystals.

Lanthanum carbonate, like sevelamer-HCl, retards the progression of vascular calcification and atherosclerosis in uremic apolipoprotein E-deficient mice

Atherosclerosis and vascular calcification (VC) progression in chronic kidney disease is favored by disturbances of mineral metabolism. We compared the effect of phosphate binder lanthanum (La) carbonate with sevelamer-HCl on atherosclerosis, VC and bone structure and function in mice with chronic renal failure (CRF). Apolipoprotein E-deficient (apoE(-/-)) mice were randomized to one non-CRF and three CRF groups, fed with standard diet (one non-CRF and one CRF) or diet supplemented with either 3% lanthanum carbonate (La3%) or 3% sevelamer-HCl (Sev3%). Both La3% and Sev3% supplemented CRF mice displayed a decrease of serum phosphorus, calcification at both intimal and medial aortic sites and atherosclerosis. This was associated with a reduction of plaque Type I collagen expression by both binders and of positive nitrotyrosine staining in response to sevelamer-HCl only. Increased mineral apposition and bone formation rates in unsupplemented CRF mice were reduced by Sev3% but not by La3%. The beneficial effects of La carbonate and sevelamer-HCl on the progression of VC and atherosclerosis in CRF mice could be mainly due to a decrease in phosphate retention and likewise a reduction of arterial Type I collagen expression. The effect of La carbonate differed from that of sevelamer-HCl in that it did not appear to exert its vascular effects via changes in oxidative stress or bone remodeling in the present model.

Carlosbarbosaite, ideally (UO2)(2)Nb2O6(OH)(2)center dot 2H(2)O, a new hydrated uranyl niobate mineral with tunnels from Jaguaracu, Minas Gerais, Brazil: description and crystal structure

Carlosbarbosaite, ideally (UO2)(2)Nb2O6(OH)(2)center dot 2H(2)O, is a new mineral which occurs as a late cavity filling in albite in the Jaguaracu pegmatite, Jaguaracu municipality, Minas Gerais, Brazil. The name honours Carlos do Prado Barbosa (1917-2003). Carlosbarbosaite forms long flattened lath-like crystals with a very simple orthorhombic morphology. The crystals are elongated along [001] and flattened on (100); they are up to 120 mu m long and 2-5 mu m thick. The colour is cream to pale yellow, the streak yellowish white and the lustre vitreous. The mineral is transparent (as individual crystals) to translucent (massive). It is not fluorescent under either long-wave or short-wave ultraviolet radiation. Carlosbarbosaite is biaxial(+) with alpha = 1.760(5), beta = 1.775(5), gamma = 1.795(5), 2V(meas) = 70(1)degrees, 2V(calc) = 83 degrees. The orientation is X parallel to a, Y parallel to b, Z parallel to c. Pleochroism is weak, in yellowish green shades, which are most intense in the Z direction. Two samples were analysed. For sample I, the composition is: UO3 54.52, CaO 2.07, Ce2O3 0.33, Nd2O3 0.49, Nb2O5 14.11, Ta2O5 15.25, TiO2 2.20, SiO2 2.14, Fe2O3 1.08, Al2O3 0.73, H2O (calc.) 11.49, total 104.41 wt.%; the empirical formula is (square 0.68Ca0.28Nd0.02Ce0.02)(Sigma=1.00)[U-1.44 square O-0.56(2.88)(H2O)(1.12)](Nb0.80Ta0.52Si0.27Ti0.21Al0.11Fe0.10)(Sigma=2.01) O-4.72(OH)(3.20)(H2O)(2.08). For sample 2, the composition is: UO3 41.83, CaO 2.10, Ce2O3 0.31, Nd2O3 1.12, Nb2O5 14.64, Ta2O5 16.34, TiO2 0.95, SiO2 3.55, Fe2O3 0.89, Al2O3 0.71, H2O (calc.) 14.99, total 97.43 wt.%; the empirical formula is (square 0.67Ca0.27Nd0.05Ce0.01)(Sigma=1.00)[U-1.04 square O-0.96(2.08)(H2O)(1.92)] (Nb0.79Ta0.53Si0.42Ti0.08Al0.10Fe0.08)(Sigma=2.00)O-4.00(OH)(3.96)(H2O)(2.04). The ideal endmember formula is (UO2)(2)Nb2O6(OH)(2)center dot 2H(2)O. Calculated densities are 4.713 g cm(-3) (sample 1) and 4.172 g cm(-3) (sample 2). Infrared spectra show that both (OH) and H2O are present. The strongest eight X-ray powder-diffraction lines [listed as d in angstrom(I)(hkl)] are: 8.405(8)(110), 7.081(10)(200), 4.201(9)(220), 3.333(6)(202), 3.053(8)(022), 2.931(7)(420), 2.803(6)(222) and 2.589(5)(040,402). The crystal structure was solved using single-crystal X-ray diffraction (R = 0.037) which gave the following data: orthorhombic, Cmem, a = 14.150(6), b = 10.395(4), c = 7.529(3) angstrom, V = 1107(1) angstrom(3), Z = 4. The crystal structure contains a single U site with an appreciable deficiency in electron scattering, which is populated by U atoms and vacancies. The U site is surrounded by seven 0 atoms in a pentagonal bipyramidal arrangemet. The Nb site is coordinated by four 0 atoms and two OH groups in an octahedral arrangement. The half-occupied tunnel Ca site is coordinated by four 0 atoms and four H2O groups. Octahedrally coordinated Nb polyhedra share edges and comers to form Nb2O6(OH)(2) double chains, and edge-sharing pentagonal bipyramidal U polyhedra form UO5 chains. The Nb2O6(OH)(2) and UO5 chains share edges to form an open U-Nb-phi framework with tunnels along [001] that contain Ca(H2O)(4) clusters. Carlosbarbosaite is closely related to a family of synthetic U-Nb-O framework tunnel structures, it differs in that is has an (OH)-bearing framework and Ca(H2O)(4) tunnel occupant. The structure of carlosbarbosaite resembles that of holfertite.

Crystal structure of Leishmania tarentolae hypoxanthine-guanine phosphoribosyltransferase

Abstract Background Hypoxanthine-guanine phosphoribosyltransferase (HGPRT) (EC 2.4.2.8) is a central enzyme in the purine recycling pathway. Parasitic protozoa of the order Kinetoplastida cannot synthesize purines de novo and use the salvage pathway to synthesize purine bases, making this an attractive target for antiparasitic drug design. Results The glycosomal HGPRT from Leishmania tarentolae in a catalytically active form purified and co-crystallized with a guanosine monophosphate (GMP) in the active site. The dimeric structure of HGPRT has been solved by molecular replacement and refined against data extending to 2.1 Å resolution. The structure reveals the contacts of the active site residues with GMP. Conclusion Comparative analysis of the active sites of Leishmania and human HGPRT revealed subtle differences in the position of the ligand and its interaction with the active site residues, which could be responsible for the different reactivities of the enzymes to allopurinol reported in the literature. The solution and analysis of the structure of Leishmania HGPRT may contribute to further investigations leading to a full understanding of this important enzyme family in protozoan parasites.

Toll-like receptor 1 N248S single-nucleotide polymorphism is associated with leprosy risk and regulates immune activation during mycobacterial infection

Conflicting findings about the association between leprosy and TLR1 variants N248S and I602S have been reported. Here, we performed case-control and family based studies, followed by replication in 2 case-control populations from Brazil, involving 3162 individuals. Results indicated an association between TLR1 248S and leprosy in the case-control study (SS genotype odds ratio [OR], 1.81; P = .004) and the family based study (z = 2.02; P = .05). This association was consistently replicated in other populations (combined OR, 1.51; P < .001), corroborating the finding that 248S is a susceptibility factor for leprosy. Additionally, we demonstrated that peripheral blood mononuclear cells (PBMCs) carrying 248S produce a lower tumor necrosis factor/interleukin-10 ratio when stimulated with Mycobacterium leprae but not with lipopolysaccharide or PAM3cysK4. The same effect was observed after infection of PBMCs with the Moreau strain of bacillus Calmette-Guerin but not after infection with other strains. Finally, molecular dynamics simulations indicated that the Toll-like receptor 1 structure containing 248S amino acid is different from the structure containing 248N. Our results suggest that TLR1 248S is associated with an increased risk for leprosy, consistent with its hypoimmune regulatory function.

Purification and characterization of a surfactin-like molecule produced by Bacillus sp. H2O-1 and its antagonistic effect against sulfate reducing bacteria

Abstract Background Bacillus sp. H2O-1, isolated from the connate water of a Brazilian reservoir, produces an antimicrobial substance (denoted as AMS H2O-1) that is active against sulfate reducing bacteria, which are the major bacterial group responsible for biogenic souring and biocorrosion in petroleum reservoirs. Thus, the use of AMS H2O-1 for sulfate reducing bacteria control in the petroleum industry is a promising alternative to chemical biocides. However, prior to the large-scale production of AMS H2O-1 for industrial applications, its chemical structure must be elucidated. This study also analyzed the changes in the wetting properties of different surfaces conditioned with AMS H2O-1 and demonstrated the effect of AMS H2O-1 on sulfate reducing bacteria cells. Results A lipopeptide mixture from AMS H2O-1 was partially purified on a silica gel column and identified via mass spectrometry (ESI-MS). It comprises four major components that range in size from 1007 to 1049 Da. The lipid moiety contains linear and branched β-hydroxy fatty acids that range in length from C13 to C16. The peptide moiety contains seven amino acids identified as Glu-Leu-Leu-Val-Asp-Leu-Leu. Transmission electron microscopy revealed cell membrane alteration of sulfate reducing bacteria after AMS H2O-1 treatment at the minimum inhibitory concentration (5 μg/ml). Cytoplasmic electron dense inclusions were observed in treated cells but not in untreated cells. AMS H2O-1 enhanced the osmosis of sulfate reducing bacteria cells and caused the leakage of the intracellular contents. In addition, contact angle measurements indicated that different surfaces conditioned by AMS H2O-1 were less hydrophobic and more electron-donor than untreated surfaces. Conclusion AMS H2O-1 is a mixture of four surfactin-like homologues, and its biocidal activity and surfactant properties suggest that this compound may be a good candidate for sulfate reducing bacteria control. Thus, it is a potential alternative to the chemical biocides or surface coating agents currently used to prevent SRB growth in petroleum industries.

Panic-like defensive behavior but not fear-induced antinociception is differently organized by dorsomedial and posterior hypothalamic nuclei of Rattus norvegicus (Rodentia, Muridae)

The hypothalamus is a forebrain structure critically involved in the organization of defensive responses to aversive stimuli. Gamma-aminobutyric acid (GABA)ergic dysfunction in dorsomedial and posterior hypothalamic nuclei is implicated in the origin of panic-like defensive behavior, as well as in pain modulation. The present study was conducted to test the difference between these two hypothalamic nuclei regarding defensive and antinociceptive mechanisms. Thus, the GABA A antagonist bicuculline (40 ng/0.2 µL) or saline (0.9% NaCl) was microinjected into the dorsomedial or posterior hypothalamus in independent groups. Innate fear-induced responses characterized by defensive attention, defensive immobility and elaborate escape behavior were evoked by hypothalamic blockade of GABA A receptors. Fear-induced defensive behavior organized by the posterior hypothalamus was more intense than that organized by dorsomedial hypothalamic nuclei. Escape behavior elicited by GABA A receptor blockade in both the dorsomedial and posterior hypothalamus was followed by an increase in nociceptive threshold. Interestingly, there was no difference in the intensity or in the duration of fear-induced antinociception shown by each hypothalamic division presently investigated. The present study showed that GABAergic dysfunction in nuclei of both the dorsomedial and posterior hypothalamus elicit panic attack-like defensive responses followed by fear-induced antinociception, although the innate fear-induced behavior originates differently in the posterior hypothalamus in comparison to the activity of medial hypothalamic subdivisions.

Structure and peroxidase activity of ferric Streptomyces clavuligerus orf10-encoded protein P450CLA: UV-visible, CD, MCD and EPR spectroscopic characterization

The present study reports the spectroscopic characterization by UV-visible absorption spectroscopy, magnetic circular dichroism (MCD) and electron paramagnetic resonance (EPR) of the recombinant orf10-encoded P450-camphor like protein (P450CLA)of Streptomyces clavuligerus expressed in Escherichia coli Rosetta in the native form and associated to external ligands containing the β-lactam, oxazole and alkylamine-derived (alcohol) moieties of the clavulamic acid. Considering the diversity of potential applications for the enzyme, the reactivity with tert-butylhydroperoxide (tert-BuOOH) was also characterized. P450CLA presents a covalently bound heme group and exhibited the UV-visible, CD and MCD spectral features of P450CAM including the fingerprint Soret band at 450 nm generated by the ferrous CO-complex. P450CLA was converted to high valence species by tert-BuOOH and promoted homolytic scission of the O-O bond. The radical profile of the reaction was tert-butyloxyl as primary and methyl and butylperoxyl as secondary radicals. The secondary methyl and butylperoxyl radicals resulted respectively from the β-scission of the alkoxyl radical and from the reaction of methyl radical with molecular oxygen.

On the three-finger protein domain fold and CD59-like proteins in Schistosoma mansoni

Background: It is believed that schistosomes evade complement-mediated killing by expressing regulatory proteins on their surface. Recently, six homologues of human CD59, an important inhibitor of the complement system membrane attack complex, were identified in the schistosome genome. Therefore, it is important to investigate whether these molecules could act as CD59-like complement inhibitors in schistosomes as part of an immune evasion strategy. Methodology/Principal Findings: Herein, we describe the molecular characterization of seven putative SmCD59-like genes and attempt to address the putative biological function of two isoforms. Superimposition analysis of the 3D structure of hCD59 and schistosome sequences revealed that they contain the three-fingered protein domain (TFPD). However, the conserved amino acid residues involved in complement recognition in mammals could not be identified. Real-time RT-PCR and Western blot analysis determined that most of these genes are up-regulated in the transition from free-living cercaria to adult worm stage. Immunolocalization experiments and tegument preparations confirm that at least some of the SmCD59-like proteins are surface-localized; however, significant expression was also detected in internal tissues of adult worms. Finally, the involvement of two SmCD59 proteins in complement inhibition was evaluated by three different approaches: (i) a hemolytic assay using recombinant soluble forms expressed in Pichia pastoris and E. coli; (ii) complement-resistance of CHO cells expressing the respective membrane-anchored proteins; and (iii) the complement killing of schistosomula after gene suppression by RNAi. Our data indicated that these proteins are not involved in the regulation of complement activation. Conclusions: Our results suggest that this group of proteins belongs to the TFPD superfamily. Their expression is associated to intra-host stages, present in the tegument surface, and also in intra-parasite tissues. Three distinct approaches using SmCD59 proteins to inhibit complement strongly suggested that these proteins are not complement inhibitors and their function in schistosomes remains to be determined.

ScPdZn and ScPtZn with YAlGe type structure: group-subgroup relation and 45Sc solid state NMR spectroscopy

The intermetallic compounds ScPdZn and ScPtZn were prepared from the elements by high-frequency melting in sealed tantalum ampoules. Both structures were refined from single crystal X-ray diffractometer data: YAlGe type, Cmcm, a = 429.53(8), b = 907.7(1), c = 527.86(1) pm, wR2 = 0.0375, 231 F2 values, for ScPdZn and a = 425.3(1), b = 918.4(2), c = 523.3(1) pm, wR2 = 0.0399, 213 F2 values for ScPtZn with 14 variables per refinement. The structures are orthorhombically distorted variants of the AlB2 type. The scandium and palladium (platinum atoms) build up ordered networks Sc3Pd3 and Sc3Pt3 (boron networks) which are slightly shifted with respect to each other. These networks are penetrated by chains of zinc atoms (262 pm in ScPtZn) which correspond to the aluminum positions, i.e. Zn(ScPd) and Zn(ScPt). The corresponding group-subgroup scheme and the differences in chemical bonding with respect to other AlB2-derived REPdZn and REPtZn compounds are discussed. 45Sc solid state NMR spectra confirm the single crystallographic scandium sites. From electronic band structure calculations the two compounds are found metallic with free electron like behavior at the Fermi level. A larger cohesive energy for ScPtZn suggests a more strongly bonded intermetallic than ScPdZn. Electron localization and overlap population analyses identify the largest bonding for scandium with the transition metal (Pd, Pt).

Analyzing the dependence structure of microarray data: a copula–based approach

The main aim of this Ph.D. dissertation is the study of clustering dependent data by means of copula functions with particular emphasis on microarray data. Copula functions are a popular multivariate modeling tool in each field where the multivariate dependence is of great interest and their use in clustering has not been still investigated. The first part of this work contains the review of the literature of clustering methods, copula functions and microarray experiments. The attention focuses on the K–means (Hartigan, 1975; Hartigan and Wong, 1979), the hierarchical (Everitt, 1974) and the model–based (Fraley and Raftery, 1998, 1999, 2000, 2007) clustering techniques because their performance is compared. Then, the probabilistic interpretation of the Sklar’s theorem (Sklar’s, 1959), the estimation methods for copulas like the Inference for Margins (Joe and Xu, 1996) and the Archimedean and Elliptical copula families are presented. In the end, applications of clustering methods and copulas to the genetic and microarray experiments are highlighted. The second part contains the original contribution proposed. A simulation study is performed in order to evaluate the performance of the K–means and the hierarchical bottom–up clustering methods in identifying clusters according to the dependence structure of the data generating process. Different simulations are performed by varying different conditions (e.g., the kind of margins (distinct, overlapping and nested) and the value of the dependence parameter ) and the results are evaluated by means of different measures of performance. In light of the simulation results and of the limits of the two investigated clustering methods, a new clustering algorithm based on copula functions (‘CoClust’ in brief) is proposed. The basic idea, the iterative procedure of the CoClust and the description of the written R functions with their output are given. The CoClust algorithm is tested on simulated data (by varying the number of clusters, the copula models, the dependence parameter value and the degree of overlap of margins) and is compared with the performance of model–based clustering by using different measures of performance, like the percentage of well–identified number of clusters and the not rejection percentage of H0 on . It is shown that the CoClust algorithm allows to overcome all observed limits of the other investigated clustering techniques and is able to identify clusters according to the dependence structure of the data independently of the degree of overlap of margins and the strength of the dependence. The CoClust uses a criterion based on the maximized log–likelihood function of the copula and can virtually account for any possible dependence relationship between observations. Many peculiar characteristics are shown for the CoClust, e.g. its capability of identifying the true number of clusters and the fact that it does not require a starting classification. Finally, the CoClust algorithm is applied to the real microarray data of Hedenfalk et al. (2001) both to the gene expressions observed in three different cancer samples and to the columns (tumor samples) of the whole data matrix.

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.