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.

Large-scale coupled-cluster calculations

Coupled-cluster theory provides one of the most successful concepts in electronic-structure theory. This work covers the parallelization of coupled-cluster energies, gradients, and second derivatives and its application to selected large-scale chemical problems, beside the more practical aspects such as the publication and support of the quantum-chemistry package ACES II MAB and the design and development of a computational environment optimized for coupled-cluster calculations. The main objective of this thesis was to extend the range of applicability of coupled-cluster models to larger molecular systems and their properties and therefore to bring large-scale coupled-cluster calculations into day-to-day routine of computational chemistry. A straightforward strategy for the parallelization of CCSD and CCSD(T) energies, gradients, and second derivatives has been outlined and implemented for closed-shell and open-shell references. Starting from the highly efficient serial implementation of the ACES II MAB computer code an adaptation for affordable workstation clusters has been obtained by parallelizing the most time-consuming steps of the algorithms. Benchmark calculations for systems with up to 1300 basis functions and the presented applications show that the resulting algorithm for energies, gradients and second derivatives at the CCSD and CCSD(T) level of theory exhibits good scaling with the number of processors and substantially extends the range of applicability. Within the framework of the ’High accuracy Extrapolated Ab initio Thermochemistry’ (HEAT) protocols effects of increased basis-set size and higher excitations in the coupled- cluster expansion were investigated. The HEAT scheme was generalized for molecules containing second-row atoms in the case of vinyl chloride. This allowed the different experimental reported values to be discriminated. In the case of the benzene molecule it was shown that even for molecules of this size chemical accuracy can be achieved. Near-quantitative agreement with experiment (about 2 ppm deviation) for the prediction of fluorine-19 nuclear magnetic shielding constants can be achieved by employing the CCSD(T) model together with large basis sets at accurate equilibrium geometries if vibrational averaging and temperature corrections via second-order vibrational perturbation theory are considered. Applying a very similar level of theory for the calculation of the carbon-13 NMR chemical shifts of benzene resulted in quantitative agreement with experimental gas-phase data. The NMR chemical shift study for the bridgehead 1-adamantyl cation at the CCSD(T) level resolved earlier discrepancies of lower-level theoretical treatment. The equilibrium structure of diacetylene has been determined based on the combination of experimental rotational constants of thirteen isotopic species and zero-point vibrational corrections calculated at various quantum-chemical levels. These empirical equilibrium structures agree to within 0.1 pm irrespective of the theoretical level employed. High-level quantum-chemical calculations on the hyperfine structure parameters of the cyanopolyynes were found to be in excellent agreement with experiment. Finally, the theoretically most accurate determination of the molecular equilibrium structure of ferrocene to date is presented.

Photoactive nanodevices for potential biological applications

Biological systems are complex and highly organized architectures governed by noncovalent interactions, which are responsible for molecular recognition, self-assembly, self-organization, adaptation and evolution processes. These systems provided the inspiration for the development of supramolecular chemistry, that aimed at the design of artificial multicomponent molecular assemblies, namely supramolecular systems, properly designed to perform different operations: each constituting unit performs a single act, whereas the entire supramolecular system is able to execute a more complex function, resulting from the cooperation of the constituting components. Supramolecular chemistry deals with the development of molecular systems able to mimic naturally occurring events, for example complexation and self-assembly through the establishment of noncovalent interactions. Moreover, the application of external stimuli, such as light, allows to perform these operations in a time- and space-controlled manner. These systems can interact with biological systems and, thus, can be applied for bioimaging, therapeutic and drug delivery purposes. In this work the study of biocompatible supramolecular species able to interact with light is presented. The first part deals with the photophysical, photochemical and electrochemical characterization of water-soluble blue emitting triazoloquinolinium and triazolopyridinium salts. Moreover, their interaction with DNA has been explored, in the perspective of developing water-soluble systems for bioimaging applications. In the second part, the effect exerted by the presence of azobenzene-bearing supramolecular species in liposomes, inserted both in the phospholipid bilayer and in the in the aqueous core of vesicles has been studied, in order to develop systems able to deliver small molecules and ions in a photocontrolled manner. Moreover, the versatility of azobenzene and its broad range of applications have been highlighted, since conjugated oligoazobenzene derivatives proved not to be adequate to be inserted in the phospholipid bilayer of liposomes, but their electrochemical properties made them interesting candidates as electron acceptor materials for photovoltaic applications.

Internally contracted multireference coupled-cluster methods

This thesis details the development of quantum chemical methods for the accurate theoretical description of molecular systems with a complicated electronic structure. In simple cases, a single Slater determinant, in which the electrons occupy a number of energetically lowest molecular orbitals, offers a qualitatively correct model. The widely used coupled-cluster method CCSD(T) efficiently includes electron correlation effects starting from this determinant and provides reaction energies in error by only a few kJ/mol. However, the method often fails when several electronic configurations are important, as, for instance, in the course of many chemical reactions or in transition metal compounds. Internally contracted multireference coupled-cluster methods (ic-MRCC methods) cure this deficiency by using a linear combination of determinants as a reference function. Despite their theoretical elegance, the ic-MRCC equations involve thousands of terms and are therefore derived by the computer. Calculations of energy surfaces of BeH2, HF, LiF, H2O, N2 and Be3 unveil the theory's high accuracy compared to other approaches and the quality of various hierarchies of approximations. New theoretical advances include size-extensive techniques for removing linear dependencies in the ic-MRCC equations and a multireference analog of CCSD(T). Applications of the latter method to O3, Ni2O2, benzynes, C6H7NO and Cr2 underscore its potential to become a new standard method in quantum chemistry.

Broadband Fourier transform rotational spectroscopy for structure determination: The water heptamer

Over the recent years chirped-pulse, Fourier-transform microwave (CP-FTMW) spectrometers have chan- ged the scope of rotational spectroscopy. The broad frequency and large dynamic range make possible structural determinations in molecular systems of increasingly larger size from measurements of heavy atom (13C, 15N, 18O) isotopes recorded in natural abundance in the same spectrum as that of the parent isotopic species. The design of a broadband spectrometer operating in the 2–8 GHz frequency range with further improvements in sensitivity is presented. The current CP-FTMW spectrometer performance is benchmarked in the analyses of the rotational spectrum of the water heptamer, (H2O)7, in both 2– 8 GHz and 6–18 GHz frequency ranges. Two isomers of the water heptamer have been observed in a pulsed supersonic molecular expansion. High level ab initio structural searches were performed to pro- vide plausible low-energy candidates which were directly compared with accurate structures provided from broadband rotational spectra. The full substitution structure of the most stable species has been obtained through the analysis of all possible singly-substituted isotopologues (H218O and HDO), and a least-squares rm(1) geometry of the oxygen framework determined from 16 different isotopic species compares with the calculated O–O equilibrium distances at the 0.01 Å level.

Broadband Fourier transform rotational spectroscopy for structure determination: The water heptamer

Over the recent years chirped-pulse, Fourier-transform microwave (CP-FTMW) spectrometers have changed the scope of rotational spectroscopy. The broad frequency and large dynamic range make possible structural determinations in molecular systems of increasingly larger size from measurements of heavy atom (C-13, N-15, O-18) isotopes recorded in natural abundance in the same spectrum as that of the parent isotopic species. The design of a broadband spectrometer operating in the 2-8 GHz frequency range with further improvements in sensitivity is presented. The current CP-FTMW spectrometer performance is benchmarked in the analyses of the rotational spectrum of the water heptamer, (H2O)(7), in both 2-8 GHz and 6-18 GHz frequency ranges. Two isomers of the water heptamer have been observed in a pulsed supersonic molecular expansion. High level ab initio structural searches were performed to provide plausible low-energy candidates which were directly compared with accurate structures provided from broadband rotational spectra. The full substitution structure of the most stable species has been obtained through the analysis of all possible singly-substituted isotopologues ((H2O)-O-18 and HDO), and a least-squares r(m)((1)) geometry of the oxygen framework determined from 16 different isotopic species compares with the calculated O-O equilibrium distances at the 0.01 angstrom level. (C) 2013 Elsevier B.V. All rights reserved.

Promising anchoring groups for single-molecule conductance measurements

The understanding of the charge transport through single molecule junctions is a prerequisite for the design and building of electronic circuits based on single molecule junctions. However, reliable and robust formation of such junctions is a challenging task to achieve. In this topical review, we present a systematic investigation of the anchoring group effect on single molecule junction conductance by employing two complementary techniques, namely scanning tunneling microscopy break junction (STM-BJ) and mechanically controllable break junction (MCBJ) techniques, based on the studies published in the literature and important results from our own work. We compared conductance studies for conventional anchoring groups described earlier with the molecular junctions formed through π-interactions with the electrode surface (Au, Pt, Ag) and we also summarized recent developments in the formation of highly conducting covalent Au–C σ-bonds using oligophenyleneethynylene (OPE) and an alkane molecular backbone. Specifically, we focus on the electron transport properties of diaryloligoyne, oligophenyleneethynylene (OPE) and/or alkane molecular junctions composed of several traditional anchoring groups, (dihydrobenzo[b]thiophene (BT), 5-benzothienyl analogue (BTh), thiol (SH), pyridyl (PY), amine (NH2), cyano (CN), methyl sulphide (SMe), nitro (NO2)) and other anchoring groups at the solid/liquid interface. The qualitative and quantitative comparison of the results obtained with different anchoring groups reveals structural and mechanistic details of the different types of single molecular junctions. The results reported in this prospective may serve as a guideline for the design and synthesis of molecular systems to be used in molecule-based electronic devices.

Pulse-shaping assisted multidimensional coherent electronic spectroscopy

Understanding nuclear and electronic dynamics of molecular systems has advanced considerably by probing their nonlinear responses with a suitable sequence of pulses. Moreover, the ability to control crucial parameters of the excitation pulses, such as duration, sequence, frequency, polarization, slowly varying envelope, or carrier phase, has led to a variety of advanced time-resolved spectroscopic methodologies. Recently, two-dimensional electronic spectroscopy with ultrashort pulses has become a more and more popular tool since it allows to obtain information on energy and coherence transfer phenomena, line broadening mechanisms, or the presence of quantum coherences in molecular complexes. Here, we present a high fidelity two-dimensional electronic spectroscopy setup designed for molecular systems in solution. It incorporates the versatility of pulse-shaping methods to achieve full control on the amplitude and phase of the individual exciting and probing pulses. Selective and precise amplitude- and phase-modulation is shown and applied to investigate electronic dynamics in several reference molecular systems.

Frequency modulation of ERK activation dynamics rewires cell fate.

Transient versus sustained ERK MAP kinase (MAPK) activation dynamics induce proliferation versus differentiation in response to epidermal (EGF) or nerve (NGF) growth factors in PC-12 cells. Duration of ERK activation has therefore been proposed to specify cell fate decisions. Using a biosensor to measure ERK activation dynamics in single living cells reveals that sustained EGF/NGF application leads to a heterogeneous mix of transient and sustained ERK activation dynamics in distinct cells of the population, different than the population average. EGF biases toward transient, while NGF biases toward sustained ERK activation responses. In contrast, pulsed growth factor application can repeatedly and homogeneously trigger ERK activity transients across the cell population. These datasets enable mathematical modeling to reveal salient features inherent to the MAPK network. Ultimately, this predicts pulsed growth factor stimulation regimes that can bypass the typical feedback activation to rewire the system toward cell differentiation irrespective of growth factor identity.

Investigation of cell fate decision in Schizosaccharomyces pombe by ncRNA annotation and statistical analyses

Studies in the fission yeast Schizosaccharomyces pombe (S. pombe) have done much to inform the view of heterochromatin and its control by the RNA interference (RNAi) machinery. Using cDNA synthesised from poly(A)-enriched RNA samples, numerous novel ncRNA loci were discovered, and the 50 and 30 ends of many other genes were refined in previous studies. Although some of these transcripts may encode novel proteins the function of the majority is yet to be determined. The authors have used strand-specific deep sequencing of RNA, irrespective of poly(A) status, to reveal a highly structured antisense programme that modulates gene expression to dictate cell fate decisions during sexual differentiation. They show that an extensive and elaborate array of ncRNA production accompanies sexual differentiation in the fission yeast S. pombe. Experimental manipulation suggests that these transcripts specifically regulate the function of the target genes.

Percentage of cultured cell-populations that stained positively and/or negatively for apoptotic markers cleaved caspase-3 and cleaved PARP, following DNA damage treatments induced by doxorubicin and TNF-alpha co-treatments

The present dataset contains the source data for Figure 2B of Tentner et al. (2012). The data shows the percentage of cultured cell-populations that stained positively and/or negatively for apoptotic markers cleaved caspase-3 and cleaved PARP, following DNA damage treatments induced by various doses of doxorubicin (0, 2 and 10 µmole/L) in the presence (100 ng/mL) or absence (0 ng/mL) of TNF-alpha co-treatment. For the six treatment conditions investigated, cell counts were made by flow cytometry at times 6, 12, 24, and 48 h following treatment; CULTURE DETAILS: U2OS cells were obtained from ATCC were maintained at 21% oxygen and 5% CO2 in Dulbecco's modified Eagle medium supplemented with 10% fetal bovine serum, penicillin, streptomycin, 2mM L-glutamine, and used within 15-20 passages. The first thymidine block was released by washing the plates three times with PBS, and incubating them in fresh thymidine-free media for 12 h. A second thymidine block was then performed by re-addition of thymidine to 2.5 mM followed by incubation for an additional 18 h. Media was aspirated, plates were washed 3 with PBS, and replaced with fresh media in the presence or absence of 10 mM aphidicolin; ANALYSIS DETAILS: See supplementary journal publication; RESULT: The authors of the supplementary journal publication conclude that TNF enhances dose-dependent cell death following doxorubicin-induced DNA damage with minimal affect on dose-dependent cell-cycle arrest.

Phosphorylation levels of JAK2 after addition of increasing Epo-concentrations: model approach (Figure 5a)

The present dataset contain source data for Figure 5a from Schilling et al., 2009. Cell fate decisions are regulated by the coordinated activation of signalling pathways such as the extracellular signal-regulated kinase (ERK) cascade, but contributions of individual kinase isoforms are mostly unknown. The authors combined quantitative data from erythropoietin-induced pathway activation in primary erythroid progenitor (colony-forming unit erythroid stage, CFU-E) cells with mathematical modelling, in order to predict and experimentally confirmed a distributive ERK phosphorylation mechanism in CFU-E cells. The authors found evidences that double-phosphorylated ERK1 attenuates proliferation beyond a certain activation level, whereas activated ERK2 enhances proliferation with saturation kinetics. Phosphorylation levels of JAK2 at 7 min after stimulation for Epo concentrations ranging from 0.1 to 1000 U/ml were simulated.

Modeled integrated responses of ppERK1 and ppERK2 of normal and elevated ERK1 and ERK2 levels (Figure 5b)

The present dataset contain source data for Figure 5b from Schilling et al., 2009. Cell fate decisions are regulated by the coordinated activation of signalling pathways such as the extracellular signal-regulated kinase (ERK) cascade, but contributions of individual kinase isoforms are mostly unknown. The authors combined quantitative data from erythropoietin-induced pathway activation in primary erythroid progenitor (colony-forming unit erythroid stage, CFU-E) cells with mathematical modelling, in order to predict and experimentally confirmed a distributive ERK phosphorylation mechanism in CFU-E cells. The authors found evidences that double-phosphorylated ERK1 attenuates proliferation beyond a certain activation level, whereas activated ERK2 enhances proliferation with saturation kinetics. They show integrated responses of double-phosphorylated ERK1 and ERK2 that were calculated for different Epo concentrations for the original model as well as for models with elevated ERK1 or ERK2 levels.

Phosphorylation levels of JAK2 after addition of increasing Epo-concentrations: experimental data (Figure 5a)

The present dataset data contain source data for Figure 5a from Schilling et al., 2009. Cell fate decisions are regulated by the coordinated activation of signalling pathways such as the extracellular signal-regulated kinase (ERK) cascade, but contributions of individual kinase isoforms are mostly unknown. The authors combined quantitative data from erythropoietin-induced pathway activation in primary erythroid progenitor (colony-forming unit erythroid stage, CFU-E) cells with mathematical modelling, in order to predict and experimentally confirmed a distributive ERK phosphorylation mechanism in CFU-E cells. The authors found evidences that double-phosphorylated ERK1 attenuates proliferation beyond a certain activation level, whereas activated ERK2 enhances proliferation with saturation kinetics. CFU-E cells were stimulated with the indicated Epo concentrations for 7 min and phosphorylation levels were determined by quantitative immunoblotting.

Proliferation of retrovirally transduced CFU-E cells after incubation with increasing Epo-concentration (Figure 5c)

Data contain source data for Figure 5c from Schilling et al., 2009. Cell fate decisions are regulated by the coordinated activation of signalling pathways such as the extracellular signal-regulated kinase (ERK) cascade, but contributions of individual kinase isoforms are mostly unknown. The authors combined quantitative data from erythropoietin-induced pathway activation in primary erythroid progenitor (colony-forming unit erythroid stage, CFU-E) cells with mathematical modelling, in order to predict and experimentally confirmed a distributive ERK phosphorylation mechanism in CFU-E cells. The authors found evidences that double-phosphorylated ERK1 attenuates proliferation beyond a certain activation level, whereas activated ERK2 enhances proliferation with saturation kinetics. Retrovirally transduced CFU-E cells were incubated with increasing Epo concentrations for 14 h and proliferation was measured by [3H]-thymidine incorporation.