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.

Immunohistochemical and Molecular Prognostic/Predictive Markers in Neoplastic Diseases

Traditional morphological examinations are not anymore sufficient for a complete evaluation of tumoral tissue and the use of neoplastic markers is of utmost importance. Neoplastic markers can be classified in: diagnostic, prognostic and predictive markers. Three markers were analyzed. 1) Insulin-like growth factor binding protein 2 (IGFBP2) was immunohistochemically examined in prostatic tissues: 40 radical prostatectomies from hormonally untreated patients with their preoperative biopsies, 10 radical prostatectomies from patients under complete androgen ablation before surgery and 10 simple prostatectomies from patients with bladder outlet obstruction. Results were compared with α-methylacyl-CoA racemase (AMACR). IGFBP2 was expressed in the cytoplasm of untreated adenocarcinomas and, to a lesser extent, in HG-PIN; the expression was markedly lower in patients after complete androgen ablation. AMACR was similarly expressed in both adenocarcinoma and HG-PIN, the level being similar in both lesions; the expression was slightly lower in patients after complete androgen ablation. IGFBP2 may be used a diagnostic marker of prostatic adenocarcinomas. 2) Heparan surface proteoglycan immunohistochemical expression was examined in 150 oral squamous cell carcinomas. Follow up information was available in 93 patients (range: 6-34 months, mean: 19±7). After surgery, chemotherapy was performed in 8 patients and radiotherapy in 61 patients. Multivariate and univariate overall survival analyses showed that high expression of syndecan-1 (SYN-1) was associated with a poor prognosis. In patients treated with radiotherapy, such association was higher. SYN-1 is a prognostic marker in oral squamous cell carcinomas; it may also represent a predictive factor for responsiveness to radiotherapy. 3) EGFR was studied in 33 pulmonary adenocarcinomas with traditional DNA sequencing methods and with two mutation-specific antibodies. Overall, the two antibodies had 61.1% sensitivity and 100% specificity in detecting EGFR mutations. EGFR mutation-specific antibodies may represent a predictive marker to identify patients candidate to tyrosine kinase inhibitors therapy.

Mechanisms of cell senescence: p21 and DNA damage response in aging and longevity

Theory of aging postulates that aging is a remodeling process where the body of survivors progressively adapts to internal and external damaging agents they are exposed to during several decades. Thus , stress response and adaptation mechanisms play a fundamental role in the aging process where the capability of adaptating effects, certainly, also is related the lifespan of each individual. A key gene linking aging to stress response is indeed p21, an induction of cyclin-dependent kinase inhibitor which triggers cell growth arrest associated with senescence and damage response and notably is involved in the up-regulation of multiple genes that have been associated with senescence or implicated in age-related . This PhD thesis project that has been performed in collaboration with the Roninson Lab at Ordway Research Institute in Albany, NY had two main aims: -the testing the hypothesis that p21 polymorphisms are involved in longevity -Evaluating age-associated differences in gene expression and transcriptional response to p21 and DNA damage In the first project, trough PCR-sequencing and Sequenom strategies, we we found out that there are about 30 polymorphic variants in the p21 gene. In addition, we found an haplotpype located in -5kb region of the p21 promoter whose frequency is ~ 2 fold higher in centenarians than in the general population (Large-scale analysis of haplotype frequencies is currently in progress). Functional studies I carried out on the promoter highilighted that the ―centenarian‖ haplotype doesn’t affect the basal p21 promoter activity or its response to p53. However, there are many other possible physiological conditions in which the centenarian allele of the p21 promoter may potentially show a different response (IL6, IFN,progesterone, vitamin E, Vitamin D etc). In the second part, project #2, trough Microarrays we seeked to evaluate the differences in gene expression between centenarians, elderly, young in dermal fibroblast cultures and their response to p21 and DNA damage. Microarray analysis of gene expression in dermal fibroblast cultures of individuals of different ages yielded a tentative "centenarian signature". A subset of genes that were up- or downregulated in centenarians showed the same response to ectopic expression of p21, yielding a putative "p21-centenarian" signature. Trough RQ-PCR (as well Microarrays studies whose analysis is in progress) we tested the DNA damage response of the p21-centenarian signature genes showing a correlation stress/aging in additional sets of young and old samples treated with p21-inducing drug doxorubicin thus finding for a subset of of them , a response to stress age-related.

Down Syndrome: Neuropsychological phenotype and mitochondrial DNA

Introduction. Down Syndrome (DS) is the most known autosomal trisomy, due to the presence in three copies of chromosome 21. Many studies were designed to identify phenotypic and clinical consequences related to the triple gene dosage. However, the general conclusion is a senescent phenotype; in particular, the most features of physiological aging, such as skin and hair changes, vision and hearing impairments, thyroid dysfunction, Alzheimer-like dementia, congenital heart defects, gastrointestinal malformations, immune system changes, appear in DS earlier than in normal age-matched subjects. The only established risk factor for the DS is advanced maternal age, responsible for changes in the meiosis of oocytes, in particular the meiotic nondisjunction of chromosome 21. In this process mitochondria play an important role since mitochondrial dysfunction, due to a variety of extrinsic and intrinsic influences, can profoundly influence the level of ATP generation in oocytes, required for a correct chromosomal segregation. Aim. The aim of this study is to investigate an integrated set of molecular genetic parameters (sequencing of complete mtDNA, heteroplasmy of the mtDNA control region, genotypes of APOE gene) in order to identify a possible association with the early neurocognitive decline observed in DS. Results. MtDNA point mutations do not accumulate with age in our study sample and do not correlate with early neurocognitive decline of DS subjects. It seems that D-loop heteroplasmy is largely not inherited and tends to accumulate somatically. Furthermore, in our study sample no association of cognitive impairment and ApoE genotype is found. Conclusions. Overall, our data cast some doubts on the involvement of these mutations in the decline of cognitive functions observed in DS.