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.

Longterm peripheral baroreflex and chemoreflex function after bilateral eversion carotid endarterectomy

Introduction The “eversion” technique for carotid endarterectomy (e-CEA), that involves the transection of the internal carotid artery at the carotid bulb and its eversion over the atherosclerotic plaque, has been associated with an increased risk of postoperative hypertension possibly due to a direct iatrogenic damage to the carotid sinus fibers. The aim of this study is to assess the long-term effect of the e-CEA on arterial baroreflex and peripheral chemoreflex function in humans. Methods A retrospective review was conducted on a prospectively compiled computerized database of 3128 CEAs performed on 2617 patients at our Center between January 2001 and March 2006. During this period, a total of 292 patients who had bilateral carotid stenosis ≥70% at the time of the first admission underwent staged bilateral CEAs. Of these, 93 patients had staged bilateral e-CEAs, 126 staged bilateral s- CEAs and 73 had different procedures on each carotid. CEAs were performed with either the eversion or the standard technique with routine Dacron patching in all cases. The study inclusion criteria were bilateral CEA with the same technique on both sides and an uneventful postoperative course after both procedures. We decided to enroll patients submitted to bilateral e-CEA to eliminate the background noise from contralateral carotid sinus fibers. Exclusion criteria were: age >70 years, diabetes mellitus, chronic pulmonary disease, symptomatic ischemic cardiac disease or medical therapy with b-blockers, cardiac arrhythmia, permanent neurologic deficits or an abnormal preoperative cerebral CT scan, carotid restenosis and previous neck or chest surgery or irradiation. Young and aged-matched healthy subjects were also recruited as controls. Patients were assessed by the 4 standard cardiovascular reflex tests, including Lying-to-standing, Orthostatic hypotension, Deep breathing, and Valsalva Maneuver. Indirect autonomic parameters were assessed with a non-invasive approach based on spectral analysis of EKG RR interval, systolic arterial pressure, and respiration variability, performed with an ad hoc software. From the analysis of these parameters the software provides the estimates of spontaneous baroreflex sensitivity (BRS). The ventilatory response to hypoxia was assessed in patients and controls by means of classic rebreathing tests. Results A total of 29 patients (16 males, age 62.4±8.0 years) were enrolled. Overall, 13 patients had undergone bilateral e-CEA (44.8%) and 16 bilateral s-CEA (55.2%) with a mean interval between the procedures of 62±56 days. No patient showed signs or symptoms of autonomic dysfunction, including labile hypertension, tachycardia, palpitations, headache, inappropriate diaphoresis, pallor or flushing. The results of standard cardiovascular autonomic tests showed no evidence of autonomic dysfunction in any of the enrolled patients. At spectral analysis, a residual baroreflex performance was shown in both patient groups, though reduced, as expected, compared to young controls. Notably, baroreflex function was better maintained in e-CEA, compared to standard CEA. (BRS at rest: young controls 19.93 ± 2.45 msec/mmHg; age-matched controls 7.75 ± 1.24; e-CEA 13.85 ± 5.14; s-CEA 4.93 ± 1.15; ANOVA P=0.001; BRS at stand: young controls 7.83 ± 0.66; age-matched controls 3.71 ± 0.35; e-CEA 7.04 ± 1.99; s-CEA 3.57 ± 1.20; ANOVA P=0.001). In all subjects ventilation (VÝ E) and oximetry data fitted a linear regression model with r values > 0.8. Oneway analysis of variance showed a significantly higher slope both for ΔVE/ΔSaO2 in controls compared with both patient groups which were not different from each other (-1.37 ± 0.33 compared with -0.33±0.08 and -0.29 ±0.13 l/min/%SaO2, p<0.05, Fig.). Similar results were observed for and ΔVE/ΔPetO2 (-0.20 ± 0.1 versus -0.01 ± 0.0 and -0.07 ± 0.02 l/min/mmHg, p<0.05). A regression model using treatment, age, baseline FiCO2 and minimum SaO2 achieved showed only treatment as a significant factor in explaining the variance in minute ventilation (R2= 25%). Conclusions Overall, we demonstrated that bilateral e-CEA does not imply a carotid sinus denervation. As a result of some expected degree of iatrogenic damage, such performance was lower than that of controls. Interestingly though, baroreflex performance appeared better maintained in e-CEA than in s-CEA. This may be related to the changes in the elastic properties of the carotid sinus vascular wall, as the patch is more rigid than the endarterectomized carotid wall that remains in the e-CEA. These data confirmed the safety of CEA irrespective of the surgical technique and have relevant clinical implication in the assessment of the frequent hemodynamic disturbances associated with carotid angioplasty stenting.

Computational methods for the analysis of protein structure and function

The vast majority of known proteins have not yet been experimentally characterized and little is known about their function. The design and implementation of computational tools can provide insight into the function of proteins based on their sequence, their structure, their evolutionary history and their association with other proteins. Knowledge of the three-dimensional (3D) structure of a protein can lead to a deep understanding of its mode of action and interaction, but currently the structures of <1% of sequences have been experimentally solved. For this reason, it became urgent to develop new methods that are able to computationally extract relevant information from protein sequence and structure. The starting point of my work has been the study of the properties of contacts between protein residues, since they constrain protein folding and characterize different protein structures. Prediction of residue contacts in proteins is an interesting problem whose solution may be useful in protein folding recognition and de novo design. The prediction of these contacts requires the study of the protein inter-residue distances related to the specific type of amino acid pair that are encoded in the so-called contact map. An interesting new way of analyzing those structures came out when network studies were introduced, with pivotal papers demonstrating that protein contact networks also exhibit small-world behavior. In order to highlight constraints for the prediction of protein contact maps and for applications in the field of protein structure prediction and/or reconstruction from experimentally determined contact maps, I studied to which extent the characteristic path length and clustering coefficient of the protein contacts network are values that reveal characteristic features of protein contact maps. Provided that residue contacts are known for a protein sequence, the major features of its 3D structure could be deduced by combining this knowledge with correctly predicted motifs of secondary structure. In the second part of my work I focused on a particular protein structural motif, the coiled-coil, known to mediate a variety of fundamental biological interactions. Coiled-coils are found in a variety of structural forms and in a wide range of proteins including, for example, small units such as leucine zippers that drive the dimerization of many transcription factors or more complex structures such as the family of viral proteins responsible for virus-host membrane fusion. The coiled-coil structural motif is estimated to account for 5-10% of the protein sequences in the various genomes. Given their biological importance, in my work I introduced a Hidden Markov Model (HMM) that exploits the evolutionary information derived from multiple sequence alignments, to predict coiled-coil regions and to discriminate coiled-coil sequences. The results indicate that the new HMM outperforms all the existing programs and can be adopted for the coiled-coil prediction and for large-scale genome annotation. Genome annotation is a key issue in modern computational biology, being the starting point towards the understanding of the complex processes involved in biological networks. The rapid growth in the number of protein sequences and structures available poses new fundamental problems that still deserve an interpretation. Nevertheless, these data are at the basis of the design of new strategies for tackling problems such as the prediction of protein structure and function. Experimental determination of the functions of all these proteins would be a hugely time-consuming and costly task and, in most instances, has not been carried out. As an example, currently, approximately only 20% of annotated proteins in the Homo sapiens genome have been experimentally characterized. A commonly adopted procedure for annotating protein sequences relies on the "inheritance through homology" based on the notion that similar sequences share similar functions and structures. This procedure consists in the assignment of sequences to a specific group of functionally related sequences which had been grouped through clustering techniques. The clustering procedure is based on suitable similarity rules, since predicting protein structure and function from sequence largely depends on the value of sequence identity. However, additional levels of complexity are due to multi-domain proteins, to proteins that share common domains but that do not necessarily share the same function, to the finding that different combinations of shared domains can lead to different biological roles. In the last part of this study I developed and validate a system that contributes to sequence annotation by taking advantage of a validated transfer through inheritance procedure of the molecular functions and of the structural templates. After a cross-genome comparison with the BLAST program, clusters were built on the basis of two stringent constraints on sequence identity and coverage of the alignment. The adopted measure explicity answers to the problem of multi-domain proteins annotation and allows a fine grain division of the whole set of proteomes used, that ensures cluster homogeneity in terms of sequence length. A high level of coverage of structure templates on the length of protein sequences within clusters ensures that multi-domain proteins when present can be templates for sequences of similar length. This annotation procedure includes the possibility of reliably transferring statistically validated functions and structures to sequences considering information available in the present data bases of molecular functions and structures.

Morphogenesis of follicular epithelium in Drosophila melanogaster: function of von Hippel-Lindau tumour suppressor gene

During my PhD I have been involved in several projects regarding the morphogenesis of the follicular epithelium, such as the analysis of the pathways that correlate follicular epithelium patterning and eggshell genes expression. Moreover, I used the follicular epithelium as a model system to analyze the function of the Drosophila homolog of the human von Hippel-Lindau (d-VHL) during oogenesis, in order to gain insight into the role of h-VHL for the pathogenesis of VHL disease. h-VHL is implicated in a variety of processes and there is now a greater appreciation of HIF-independent h-VHL functions that are relevant to tumour development, including maintenance and organization of the primary cilium, maintenance of the differentiated phenotype in renal cells and regulation of epithelial-mesenchymal transition. However, the function of h-VHL gene during development has not been fully understood. It was previously shown that d-VHL down-regulates the motility of tubular epithelial cells (tracheal cells) during embryogenesis. Epithelial morphogenesis is important for organogenesis and pivotal for carcinogenesis, but mechanisms that control it are poorly understood. The Drosophila follicular epithelium is a genetically tractable model to understand these mechanisms in vivo. Therefore, to examine whether d-VHL has a role in epithelial morphogenesis and maintenance, I performed genetic and molecular analyses by using in vivo and in vitro approaches. From my analysis, I determined that d-VHL binds to and stabilizes microtubules. Loss of d-VHL depolymerizes the microtubule network during oogenesis, leading to a possible deregulation in the subcellular trafficking transport of polarity markers from Golgi apparatus to the different domains in which follicle cells are divided. The analysis carried out has allowed to establish a significant role of d-VHL in the maintenance of the follicular epithelium integrity.

Complexity of N-Myc transcriptional function in childhood neuroblastoma

Myc is a transcription factor that can activate transcription of several hundreds genes by direct binding to their promoters at specific DNA sequences (E-box). However, recent studies have also shown that it can exert its biological role by repressing transcription. Such studies collectively support a model in which c-Myc-mediated repression occurs through interactions with transcription factors bound to promoter DNA regions but not through direct recognition of typical E-box sequences. Here, we investigated whether N-Myc can also repress gene transcription, and how this is mechanistically achieved. We used human neuroblastoma cells as a model system in that N-MYC amplification/over-expression represents a key prognostic marker of this tumour. By means of transcription profile analyses we could identify at least 5 genes (TRKA, p75NTR, ABCC3, TG2, p21) that are specifically repressed by N-Myc. Through a dual-step-ChIP assay and genetic dissection of gene promoters, we found that N-Myc is physically associated with gene promoters in vivo, in proximity of the transcription start site. N-Myc association with promoters requires interaction with other proteins, such as Sp1 and Miz1 transcription factors. Furthermore, we found that N-Myc may repress gene expression by interfering directly with Sp1 and/or with Miz1 activity (i.e. TRKA, p75NTR, ABCC3, p21) or by recruiting Histone Deacetylase 1 (Hdac1) (i.e. TG2). In vitro analyses show that distinct N-Myc domains can interact with Sp1, Miz1 and Hdac1, supporting the idea that Myc may participate in distinct repression complexes by interacting specifically with diverse proteins. Finally, results show that N-Myc, through repressed genes, affects important cellular functions, such as apoptosis, growth, differentiation and motility. Overall, our results support a model in which N-Myc, like c-Myc, can repress gene transcription by direct interaction with Sp1 and/or Miz1, and provide further lines of evidence on the importance of transcriptional repression by Myc factors in tumour biology.

Velocity structure and seismic anisotropy in the crust and upper mantle from Receiver Function analysis: three case studies in Italy

The research for this PhD project consisted in the application of the RFs analysis technique to different data-sets of teleseismic events recorded at temporary and permanent stations located in three distinct study regions: Colli Albani area, Northern Apennines and Southern Apennines. We found some velocity models to interpret the structures in these regions, which possess very different geologic and tectonics characteristics and therefore offer interesting case study to face. In the Colli Albani some of the features evidenced in the RFs are shared by all the analyzed stations: the Moho is almost flat and is located at about 23 km depth, and the presence of a relatively shallow limestone layer is a stable feature; contrariwise there are features which vary from station to station, indicating local complexities. Three seismic stations, close to the central part of the former volcanic edifice, display relevant anisotropic signatures­­­ with symmetry axes consistent with the emplacement of the magmatic chamber. Two further anisotropic layers are present at greater depth, in the lower crust and the upper mantle, respectively, with symmetry axes directions related to the evolution of the volcano complex. In Northern Apennines we defined the isotropic structure of the area, finding the depth of the Tyrrhenian (almost 25 km and flat) and Adriatic (40 km and dipping underneath the Apennines crests) Mohos. We determined a zone in which the two Mohos overlap, and identified an anisotropic body in between, involved in the subduction and going down with the Adiratic Moho. We interpreted the downgoing anisotropic layer as generated by post-subduction delamination of the top-slab layer, probably made of metamorphosed crustal rocks caught in the subduction channel and buoyantly rising toward the surface. In the Southern Apennines, we found the Moho depth for 16 seismic stations, and highlighted the presence of an anisotropic layer underneath each station, at about 15-20 km below the whole study area. The moho displays a dome-like geometry, as it is shallow (29 km) in the central part of the study area, whereas it deepens peripherally (down to 45 km); the symmetry axes of anisotropic layer, interpreted as a layer separating the upper and the lower crust, show a moho-related pattern, indicated by the foliation of the layer which is parallel to the Moho trend. Moreover, due to the exceptional seismic event occurred on April 6th next to L’Aquila town, we determined the Vs model for two station located next to the epicenter. An extremely high velocity body is found underneath AQU station at 4-10 km depth, reaching Vs of about 4 km/s, while this body is lacking underneath FAGN station. We compared the presence of this body with other recent works and found an anti-correlation between the high Vs body, the max slip patches and earthquakes distribution. The nature of this body is speculative since such high velocities are consistent with deep crust or upper mantle, but can be interpreted as a as high strength barrier of which the high Vs is a typical connotation.

Studies of OPA1 pathogenic mechanisms in Dominant Optic Atrophy and novel protein function in mitochondrial DNA stability

The mitochondrion is an essential cytoplasmic organelle that provides most of the energy necessary for eukaryotic cell physiology. Mitochondrial structure and functions are maintained by proteins of both mitochondrial and nuclear origin. These organelles are organized in an extended network that dynamically fuses and divides. Mitochondrial morphology results from the equilibrium between fusion and fission processes, controlled by a family of “mitochondria-shaping” proteins. It is becoming clear that defects in mitochondrial dynamics can impair mitochondrial respiration, morphology and motility, leading to apoptotic cell death in vitro and more or less severe neurodegenerative disorders in vivo in humans. Mutations in OPA1, a nuclear encoded mitochondrial protein, cause autosomal Dominant Optic Atrophy (DOA), a heterogeneous blinding disease characterized by retinal ganglion cell degeneration leading to optic neuropathy (Delettre et al., 2000; Alexander et al., 2000). OPA1 is a mitochondrial dynamin-related guanosine triphosphatase (GTPase) protein involved in mitochondrial network dynamics, cytochrome c storage and apoptosis. This protein is anchored or associated on the inner mitochondrial membrane facing the intermembrane space. Eight OPA1 isoforms resulting from alternative splicing combinations of exon 4, 4b and 5b have been described (Delettre et al., 2001). These variants greatly vary among diverse organs and the presence of specific isoforms has been associated with various mitochondrial functions. The different spliced exons encode domains included in the amino-terminal region and contribute to determine OPA1 functions (Olichon et al., 2006). It has been shown that exon 4, that is conserved throughout evolution, confers functions to OPA1 involved in maintenance of the mitochondrial membrane potential and in the fusion of the network. Conversely, exon 4b and exon 5b, which are vertebrate specific, are involved in regulation of cytochrome c release from mitochondria, and activation of apoptosis, a process restricted to vertebrates (Olichon et al., 2007). While Mgm1p has been identified thanks to its role in mtDNA maintenance, it is only recently that OPA1 has been linked to mtDNA stability. Missense mutations in OPA1 cause accumulation of multiple deletions in skeletal muscle. The syndrome associated to these mutations (DOA-1 plus) is complex, consisting of a combination of dominant optic atrophy, progressive external ophtalmoplegia, peripheral neuropathy, ataxia and deafness (Amati- Bonneau et al., 2008; Hudson et al., 2008). OPA1 is the fifth gene associated with mtDNA “breakage syndrome” together with ANT1, PolG1-2 and TYMP (Spinazzola et al., 2009). In this thesis we show for the first time that specific OPA1 isoforms associated to exon 4b are important for mtDNA stability, by anchoring the nucleoids to the inner mitochondrial membrane. Our results clearly demonstrate that OPA1 isoforms including exon 4b are intimately associated to the maintenance of the mitochondrial genome, as their silencing leads to mtDNA depletion. The mechanism leading to mtDNA loss is associated with replication inhibition in cells where exon 4b containing isoforms were down-regulated. Furthermore silencing of exon 4b associated isoforms is responsible for alteration in mtDNA-nucleoids distribution in the mitochondrial network. In this study it was evidenced that OPA1 exon 4b isoform is cleaved to provide a 10kd peptide embedded in the inner membrane by a second transmembrane domain, that seems to be crucial for mitochondrial genome maintenance and does correspond to the second transmembrane domain of the yeasts orthologue encoded by MGM1 or Msp1, which is also mandatory for this process (Diot et al., 2009; Herlan et al., 2003). Furthermore in this thesis we show that the NT-OPA1-exon 4b peptide co-immuno-precipitates with mtDNA and specifically interacts with two major components of the mitochondrial nucleoids: the polymerase gamma and Tfam. Thus, from these experiments the conclusion is that NT-OPA1- exon 4b peptide contributes to the nucleoid anchoring in the inner mitochondrial membrane, a process that is required for the initiation of mtDNA replication and for the distribution of nucleoids along the network. These data provide new crucial insights in understanding the mechanism involved in maintenance of mtDNA integrity, because they clearly demonstrate that, besides genes implicated in mtDNA replications (i.e. polymerase gamma, Tfam, twinkle and genes involved in the nucleotide pool metabolism), OPA1 and mitochondrial membrane dynamics play also an important role. Noticeably, the effect on mtDNA is different depending on the specific OPA1 isoforms down-regulated, suggesting the involvement of two different combined mechanisms. Over two hundred OPA1 mutations, spread throughout the coding region of the gene, have been described to date, including substitutions, deletions or insertions. Some mutations are predicted to generate a truncated protein inducing haploinsufficiency, whereas the missense nucleotide substitutions result in aminoacidic changes which affect conserved positions of the OPA1 protein. So far, the functional consequences of OPA1 mutations in cells from DOA patients are poorly understood. Phosphorus MR spectroscopy in patients with the c.2708delTTAG deletion revealed a defect in oxidative phosphorylation in muscles (Lodi et al., 2004). An energetic impairment has been also show in fibroblasts with the severe OPA1 R445H mutation (Amati-Bonneau et al., 2005). It has been previously reported by our group that OPA1 mutations leading to haploinsufficiency are associated in fibroblasts to an oxidative phosphorylation dysfunction, mainly involving the respiratory complex I (Zanna et al., 2008). In this study we have evaluated the energetic efficiency of a panel of skin fibroblasts derived from DOA patients, five fibroblast cell lines with OPA1 mutations causing haploinsufficiency (DOA-H) and two cell lines bearing mis-sense aminoacidic substitutions (DOA-AA), and compared with control fibroblasts. Although both types of DOA fibroblasts maintained a similar ATP content when incubated in a glucose-free medium, i.e. when forced to utilize the oxidative phosphorylation only to produce ATP, the mitochondrial ATP synthesis through complex I, measured in digitonin-permeabilized cells, was significantly reduced in cells with OPA1 haploinsufficiency only, whereas it was similar to controls in cells with the missense substitutions. Furthermore, evaluation of the mitochondrial membrane potential (DYm) in the two fibroblast lines DOA-AA and in two DOA-H fibroblasts, namely those bearing the c.2819-2A>C mutation and the c.2708delTTAG microdeletion, revealed an anomalous depolarizing response to oligomycin in DOA-H cell lines only. This finding clearly supports the hypothesis that these mutations cause a significant alteration in the respiratory chain function, which can be unmasked only when the operation of the ATP synthase is prevented. Noticeably, oligomycin-induced depolarization in these cells was almost completely prevented by preincubation with cyclosporin A, a well known inhibitor of the permeability transition pore (PTP). This results is very important because it suggests for the first time that the voltage threshold for PTP opening is altered in DOA-H fibroblasts. Although this issue has not yet been addressed in the present study, several are the mechanisms that have been proposed to lead to PTP deregulation, including in particular increased reactive oxygen species production and alteration of Ca2+ homeostasis, whose role in DOA fibroblasts PTP opening is currently under investigation. Identification of the mechanisms leading to altered threshold for PTP regulation will help our understanding of the pathophysiology of DOA, but also provide a strategy for therapeutic intervention.

Exercise associated factors affecting cardiovascular health and function: from myocyte signaling to genetic variations

The exact mechanisms of the exercise induced adaptations is not lucid, but recent studies have delineated two means of signaling by which the adaptations occur (1) substrate availability signaling (metabolic stress) (2) hormone-receptor signaling. We have decided to specifically investigate two metabolic signaling enzymes [AMP-activated kinase (AMPK) and Sirtuin 1(SIRT1)] and two hormones [Adiponectin and Adrenergic stimulation].Tis based on four papers with the following conclusions: (1)Increase in SIRT1 activity and expression in H9c2 cells treated with phenylephrine is an adaptive response to the hypertrophic stress, mediated by AMPK. (2)The lack of optimal nutritional conditions (energetic substrates) due to a prolonged activation of AMPK can contrast the establishment of hypertrophy, possibly also by means of the negative modulation of ODC activity. (3) Our findings offer a possibile hypothesis as to the fact the the G allele on site 45 could lead to the increasd risk of Type II diabetes through a decrease in lean body mass. (4) Our results suggest that there is an ADIPOQ gene effect in relation to bone parameters. Statistical analysis show that the presence of the T allele in position 45 favors an increase in lumbar spine bone mineral content (BMC) when compared to subjects with a G allele substitution, which can be do the the increase in lean body mass in this genotype group.

Effect of estrogen suppression on lung function in pulmonary Lymphangioleiomyomatosis

Background: Lymphangioleiomyomatosis (LAM), a rare progressive disease, is characterized by the proliferation of abnormal smooth muscle cells (LAM cells) in the lung, which leads to cystic parenchymal destruction and progressive respiratory failure. Estrogen receptors are present in LAM cells. LAM affects almost exclusively women of childbearing age. These findings, along with reports of disease progression during pregnancy or treatment with exogenous estrogens, have led to the assumption that hormonal factors play an important role in the pathogenesis of LAM. So, various therapies aim at preventing estrogen receptors (ER) by lowering circulating estrogen levels, by trying to block ER activity, or by attempting to lower ER expression in LAM. Prior experience have yielded conflicting results. Objective: The goal of this study was to evaluate, retrospectively, the effect of estrogen suppression in 21 patients with LAM. Design: We evaluated hormonal assays, pulmonary function tests and gas-exchange at baseline and after 12, 24 and 36 months after initiating hormonal manipulation. Results: The mean yearly rates of decline in FEV1 and DLCO are lower than those observed in prior studies and just DLCO decline was statistically significant. We also found an improvement of mean value of FVC and PaO2. Conclusions: Estrogen suppression appears to prevent decline in lung function in LAM.

Pathogenic role of IL-33-mediated eosinophil infiltration and function in experimental inflammatory bowel disease

IL-33/ST2 axis is known to promote Th2 immune responses and has been linked to several autoimmune and inflammatory disorders, including inflammatory bowel disease (IBD), and recent evidences show that it can regulate eosinophils (EOS) infiltration and function. Based also on the well documented relationship between EOS and IBD, we assessed the role of IL-33-mediated eosinophilia and ileal inflammation in SAMP1/YitFc (SAMP) murine model of Th1/Th2 chronic enteritis, and we found that IL-33 is related to inflammation progression and EOS infiltration as well as IL-5 and eotaxins increase. Administering IL-33 to SAMP and AKR mice augmented eosinophilia, eotaxins mRNA expression and Th2 molecules production, whereas blockade of ST2 and/or typical EOS molecules, such as IL-5 and CCR3, resulted in a marked decrease of inflammation, EOS infiltration, IL-5 and eotaxins mRNA expression and Th2 cytokines production. Human data supported mice’s showing an increased colocalization of IL-33 and EOS in the colon mucosa of UC patients, as well as an augmented IL-5 and eotaxins mRNA expression, when compared to non-UC. Lastly we analyzed SAMP raised in germ free (GF) condition to see the microbiota effect on IL-33 expression and Th2 responses leading to chronic intestinal inflammation. We found a remarkable decrease in ileal IL-33 and Th2 cytokines mRNA expression as well as EOS infiltration in GF versus normal SAMP with comparable inflammatory scores. Moreover, EOS depletion in normal SAMP didn’t affect IL-33 mRNA expression. These data demonstrate a pathogenic role of IL-33-mediated eosinophilia in chronic intestinal inflammation, and that blockade of IL-33 and/or downstream EOS activation may represent a novel therapeutic modality to treat patients with IBD. Also they highlight the gut microbiota role in IL-33 production, and the following EOS infiltration in the intestinal mucosa, confirming that the microbiota is essential in mounting potent Th2 response leading to chronic ileitis in SAMP.

Regulation of wake-sleep states and state-dependent cardiovascular function in diet-induced obesity rats

Obesity often predisposes to coronary heart disease, heart failure, and sudden death. Also, several studies suggest a reciprocal enhancing interaction between obesity and sleep curtailment. Aim of the present study was to go deeper in the understanding of sleep and cardiovascular regulation in an animal model of diet-induced obesity (DIO). According to this, Wake-Sleep (W-S) regulation, and W-S dependent regulation of cardiovascular and metabolic/thermoregulatory function was studied in DIO rats, under normal laboratory conditions and during sleep deprivation and the following recovery period, enhancing either wake or sleep, respectively. After 8 weeks of the delivery of a hypercaloric (HC) diet, treated animals were heavier than those fed a normocaloric (NC) diet (NC: 441 ±17g; HC: 557±17g). HC rats slept more than NC ones during the activity period (Dark) of the normal 12h:12h light-dark (LD) cycle (Wake: 67.3±1.2% and 57.2 ±1.6%; NREM sleep (NREMS): 26.8±1.0% and 34.0±1.4%; REM sleep (REMS): 5.7±0. 6% and 8.6±0.7%; for NC and HC, respectively; p<0.05 for all). HC rats were hypertensive throughout the W-S states, as shown by the mean arterial blood pressure values across the 24-h period (Wake: 90.0±5.3 and 97.3±1.3; NREMS: 85.1±5.5 and 92.2±1.2; REMS: 87.2±4.5 and 96.5±1.1, mmHg for NC and HC, respectively; p<0.05 for all). Also, HC rats appeared to be slightly bradycardic compared to NC ones (Wake: 359.8±9.3 and 352.4±7.7; NREMS: 332.5±10.1 and 328.9±5.4; REMS: 338.5±9.3 and 334.4±5.8; bpm for NC and HC, respectively; p<0.05 for Wake). In HC animals, sleep regulation was not apparently altered during the sleep rebound observed in the recovery period following sleep deprivation, although REMS rebound appeared to be quicker in NC animals. In conclusion, these results indicate that in the rat obesity interfere with W-S and cardiovascular regulation and that DIO rats are suitable for further studies aimed at a better understanding of obesity comorbidities.