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.

Development of a comprehensive on-line multidimensional high performance column liquid chromatography system for protein and peptide mapping with integrated size selective sample fractionation

Aufbau einer kontinuierlichen, mehrdimensionalen Hochleistungs-flüssigchromatographie-Anlage für die Trennung von Proteinen und Peptiden mit integrierter größenselektiver ProbenfraktionierungEs wurde eine mehrdimensionale HPLC-Trennmethode für Proteine und Peptide mit einem Molekulargewicht von <15 kDa entwickelt.Im ersten Schritt werden die Zielanalyte von höhermolekularen sowie nicht ionischen Bestandteilen mit Hilfe von 'Restricted Access Materialien' (RAM) mit Ionenaustauscher-Funktionalität getrennt. Anschließend werden die Proteine auf einer analytischen Ionenaustauscher-Säule sowie auf Reversed-Phase-Säulen getrennt. Zur Vermeidung von Probenverlusten wurde ein kontinuierlich arbeitendes, voll automatisiertes System auf Basis unterschiedlicher Trenngeschwindigkeiten und vier parallelen RP-Säulen aufgebaut.Es werden jeweils zwei RP-Säulen gleichzeitig, jedoch mit zeitlich versetztem Beginn eluiert, um durch flache Gradienten ausreichende Trennleistungen zu erhalten. Während die dritte Säule regeneriert wird, erfolgt das Beladen der vierte Säule durch Anreicherung der Proteine und Peptide am Säulenkopf. Während der Gesamtanalysenzeit von 96 Minuten werden in Intervallen von 4 Minuten Fraktionen aus der 1. Dimension auf die RP-Säulen überführt und innerhalb von 8 Minuten getrennt, wobei 24 RP-Chromatogramme resultieren.Als Testsubstanzen wurden u.a. Standardproteine, Proteine und Peptide aus humanem Hämofiltrat sowie aus Lungenfibroblast-Zellkulturüberständen eingesetzt. Weiterhin wurden Fraktionen gesammelt und mittels MALDI-TOF Massenspektrometrie untersucht. Bei einer Injektion wurden in den 24 RP-Chromatogrammen mehr als 1000 Peaks aufgelöst. Der theoretische Wert der Peakkapazität liegt bei ungefähr 3000.

Cellulose based Lithium ion polymer electrolytes for Lithium batteries

The separator membrane in batteries and fuel cells is of crucial importance for the function of these devices. In lithium ion batteries the separator membrane as well as the polymer matrix of the electrodes consists of polymer electrolytes which are lithium ion conductors. To overcome the disadvantage of currently used polymer electrolytes which are highly swollen with liquids and thus mechanically and electrochemically unstable, the goal of this work is a new generation of solid polymer electrolytes with a rigid backbone and a soft side chain structure. Moreover the novel material should be based on cheap substrates and its synthesis should not be complicated aiming at low overall costs. The new materials are based on hydroxypropylcellulose and oligoethyleneoxide derivatives as starting materials. The grafting of the oligoethyleneoxide side chains onto the cellulose was carried out following two synthetic methods. One is based on a bromide derivative and another based on p-toluolsulfonyl as a leaving group. The side chain reagents were prepared form tri(ethylene glycol) monoethyl ether. In order to improve the mechanical properties the materials were crosslinked. Two different conceptions have been engaged based on either urethane chemistry or photosensitive dimethyl-maleinimide derivatives. PEO - graft - cellulose derivatives with a high degree of substitution between 2,9 and 3,0 were blended with lithium trifluoromethane-sulfonate, lithium bis(trifluorosulfone)imide and lithium tetrafluoroborate. The molar ratios were in the range from 0,02 to 0,2 [Li]/[O]. The products have been characterized with nuclear magnetic resonance (NMR), gel permeation chromatography (GPC) and laserlight scattering (LS) with respect to their degree of substitution and molecular weight. The effect of salt concentration on ionic conductivity, thermal behaviour and morphology has been investiga-ted with impedance spectroscopy, differential scanning calorimetry (DSC) and thermal gravimetric analysis (TGA). The crosslinking reactions were controlled with dynamic mechanical analysis (DMS). The degree of substitution of our products is varying between 2,8 and 3,0 as determined by NMR. PEO - graft - cellulose derivatives are highly viscous liquids at room temperature with glass transition temperatures around 215 K. The glass transition temperature for the Lithium salt complexes of PEO - graft - cellulose deri-vatives increase with increasing salt content. The maximum conductivity at room temperature is about 10-4 and at 100°C around 10-3 Scm-1. The presence of lithium salt decreases the thermal stability of the complexes in comparison to pure PEO - graft - cellulose derivatives. Complexes heated over 140 – 150°C completely lose their ionic conductivity. The temperature dependence of the conductivity presented as Arrhenius-type plots for all samples is similar in shape and follows a VTF behaviour. This proofs that the ionic transport is closely related to the segmental motions of the polymer chains. Novel cellulose derivatives with grafted oligoethylen-oxide side chains with well-defined chemical structure and high side chain grafting density have been synthesized. Cellulose was chosen as stiff, rod like macromolecule for the backbone while oligoethylen-oxides are chosen as flexible side chains. A maximum grafting density of 3.0 have been obtained. The best conductivity reaches 10-3 Scm-1 at 100°C for a Li-triflate salt complex with a [Li]/[O] ratio of 0.8. The cross-linked complexes containing the lithium salts form elastomeric films with convenient mechanical stability. Our method of cellulose modification is based on relatively cheap and commercially available substrates and as such appears to be a promising alternative for industrial applications.

Development and validation of a liquid chromatography-tandem mass spectrometry method for the simultaneous determination of nine corticosteroid residues in bovine liver samples

A liquid chromatography tandem mass spectrometry (LC-MS/MS) confirmatory method for the simultaneous determination of nine corticosteroids in liver, including the four MRL compounds listed in Council Regulation 37/2010, was developed. After an enzymatic deconjugation and a solvent extraction of the liver tissue, the resulting solution was cleaned up through an SPE Oasis HLB cartridge. The analytes were then detected by liquid chromatography-negative-ion electrospray tandem mass spectrometry, using deuterium-labelled internal standards. The procedure was validated as a quantitative confirmatory method according to the Commission Decision 2002/657/EC criteria. The results showed that the method was suitable for statutory residue testing regarding the following performance characteristics: instrumental linearity, specificity, precision (repeatability and intra-laboratory reproducibility), recovery, decision limit (CCα), detection capability (CCβ) and ruggedness. All the corticosteroids can be detected at a concentration around 1 μg kg(-1); the recoveries were above 62% for all the analytes. Repeatability and reproducibility (within-laboratory reproducibility) for all the analytes were below 7.65% and 15.5%, respectively.

A closer look at the high affinity benzodiazepine binding site on GABAA receptors

Ligands of the benzodiazepine binding site of the GABA(A) receptor come in three flavors: positive allosteric modulators, negative allosteric modulators and antagonists all of which can bind with high affinity. The GABA(A) receptor is a pentameric protein which forms a chloride selective ion channel and ligands of the benzodiazepine binding site stabilize three different conformations of this protein. Classical benzodiazepines exert a positive allosteric effect by increasing the apparent affinity of channel opening by the agonist γ-aminobutyric acid (GABA). We concentrate here on the major adult isoform, the α(1)β(2)γ(2) GABA(A) receptor. The classical binding pocket for benzodiazepines is located in a subunit cleft between α(1) and γ(2) subunits in a position homologous to the agonist binding site for GABA that is located between β(2) and α(1) subunits. We review here approaches to this picture. In particular, point mutations were performed in combination with subsequent analysis of the expressed mutant proteins using either electrophysiological techniques or radioactive ligand binding assays. The predictive power of these methods is assessed by comparing the results with the predictions that can be made on the basis of the recently published crystal structure of the acetylcholine binding protein that shows homology to the N-terminal, extracellular domain of the GABA(A) receptor. In addition, we review an approach to the question of how the benzodiazepine ligands are positioned in their binding pocket. We also discuss a newly postulated modulatory site for benzodiazepines at the α(1)/β(2) subunit interface, homologous to the classical benzodiazepine binding pocket.

Quantification of plasma carnitine and acylcarnitines by high-performance liquid chromatography-tandem mass spectrometry using online solid-phase extraction

Carnitine is an amino acid derivative that plays a key role in energy metabolism. Endogenous carnitine is found in its free form or esterified with acyl groups of several chain lengths. Quantification of carnitine and acylcarnitines is of particular interest for screening for research and metabolic disorders. We developed a method with online solid-phase extraction coupled to high-performance liquid chromatography and tandem mass spectrometry to quantify carnitine and three acylcarnitines with different polarity (acetylcarnitine, octanoylcarnitine, and palmitoylcarnitine). Plasma samples were deproteinized with methanol, loaded on a cation exchange trapping column and separated on a reversed-phase C8 column using heptafluorobutyric acid as an ion-pairing reagent. Considering the endogenous nature of the analytes, we quantified with the standard addition method and with external deuterated standards. Solid-phase extraction and separation were achieved within 8 min. Recoveries of carnitine and acylcarnitines were between 98 and 105 %. Both quantification methods were equally accurate (all values within 84 to 116 % of target concentrations) and precise (day-to-day variation of less than 18 %) for all carnitine species and concentrations analyzed. The method was used successfully for determination of carnitine and acylcarnitines in different human samples. In conclusion, we present a method for simultaneous quantification of carnitine and acylcarnitines with a rapid sample work-up. This approach requires small sample volumes and a short analysis time, and it can be applied for the determination of other acylcarnitines than the acylcarnitines tested. The method is useful for applications in research and clinical routine.

Divalent metal-ion transporter DMT1 mediates both H+ -coupled Fe2+ transport and uncoupled fluxes

The H(+) -coupled divalent metal-ion transporter DMT1 serves as both the primary entry point for iron into the body (intestinal brush-border uptake) and the route by which transferrin-associated iron is mobilized from endosomes to cytosol in erythroid precursors and other cells. Elucidating the molecular mechanisms of DMT1 will therefore increase our understanding of iron metabolism and the etiology of iron overload disorders. We expressed wild type and mutant DMT1 in Xenopus oocytes and monitored metal-ion uptake, currents and intracellular pH. DMT1 was activated in the presence of an inwardly directed H(+) electrochemical gradient. At low extracellular pH (pH(o)), H(+) binding preceded binding of Fe(2+) and its simultaneous translocation. However, DMT1 did not behave like a typical ion-coupled transporter at higher pH(o), and at pH(o) 7.4 we observed Fe(2+) transport that was not associated with H(+) influx. His(272) --> Ala substitution uncoupled the Fe(2+) and H(+) fluxes. At low pH(o), H272A mediated H(+) uniport that was inhibited by Fe(2+). Meanwhile H272A-mediated Fe(2+) transport was independent of pH(o). Our data indicate (i) that H(+) coupling in DMT1 serves to increase affinity for Fe(2+) and provide a thermodynamic driving force for Fe(2+) transport and (ii) that His-272 is critical in transducing the effects of H(+) coupling. Notably, our data also indicate that DMT1 can mediate facilitative Fe(2+) transport in the absence of a H(+) gradient. Since plasma membrane expression of DMT1 is upregulated in liver of hemochromatosis patients, this H(+) -uncoupled facilitative Fe(2+) transport via DMT1 can account for the uptake of nontransferrin-bound plasma iron characteristic of iron overload disorders.

Functional properties of multiple isoforms of human divalent metal-ion transporter 1 (DMT1)

DMT1 (divalent metal-ion transporter 1) is a widely expressed metal-ion transporter that is vital for intestinal iron absorption and iron utilization by most cell types throughout the body, including erythroid precursors. Mutations in DMT1 cause severe microcytic anaemia in animal models. Four DMT1 isoforms that differ in their N- and C-termini arise from mRNA transcripts that vary both at their 5'-ends (starting in exon 1A or exon 1B) and at their 3'-ends giving rise to mRNAs containing (+) or lacking (-) the 3'-IRE (iron-responsive element) and resulting in altered C-terminal coding sequences. To determine whether these variations result in functional differences between isoforms, we explored the functional properties of each isoform using the voltage clamp and radiotracer assays in cRNA-injected Xenopus oocytes. 1A/IRE+-DMT1 mediated Fe2+-evoked currents that were saturable (K(0.5)(Fe) approximately 1-2 microM), temperature-dependent (Q10 approximately 2), H+-dependent (K(0.5)(H) approximately 1 muM) and voltage-dependent. 1A/IRE+-DMT1 exhibited the provisional substrate profile (ranked on currents) Cd2+, Co2+, Fe2+, Mn2+>Ni2+, V3+>>Pb2+. Zn2+ also evoked large currents; however, the zinc-evoked current was accounted for by H+ and Cl- conductances and was not associated with significant Zn2+ transport. 1B/IRE+-DMT1 exhibited the same substrate profile, Fe2+ affinity and dependence on the H+ electrochemical gradient. Each isoform mediated 55Fe2+ uptake and Fe2+-evoked currents at low extracellular pH. Whereas iron transport activity varied markedly between the four isoforms, the activity for each correlated with the density of anti-DMT1 immunostaining in the plasma membrane, and the turnover rate of the Fe2+ transport cycle did not differ between isoforms. Therefore all four isoforms of human DMT1 function as metal-ion transporters of equivalent efficiency. Our results reveal that the N- and C-terminal sequence variations among the DMT1 isoforms do not alter DMT1 functional properties. We therefore propose that these variations serve as tissue-specific signals or cues to direct DMT1 to the appropriate subcellular compartments (e.g. in erythroid cells) or the plasma membrane (e.g. in intestine).

Metal-ion-dependent biological properties of a chelator-derived somatostatin analogue for tumour targeting

Somatostatin-based radioligands have been shown to have sensitive imaging properties for neuroendocrine tumours and their metastases. The potential of [(55)Co(dotatoc)] (dotatoc =4,7,10-tricarboxymethyl-1,4,7,10-tetraazacyclododecane-1-ylacetyl-D-Phe-(Cys-Tyr-D-Trp-Lys-Thr-Cys)-threoninol (disulfide bond)) as a new radiopharmaceutical agent for PET has been evaluated. (57)Co was used as a surrogate of the positron emitter (55)Co and the pharmacokinetics of [(57)Co(dotatoc)] were investigated by using two nude mouse models. The somatostatin receptor subtype (sst1-sst5) affinity profile of [(nat)Co(dotatoc)] on membranes transfected with human somatostatin receptor subtypes was assessed by using autoradiographic methods. These studies revealed that [(57)Co(dotatoc)] is an sst2-specific radiopeptide which presents the highest affinity ever found for the sst2 receptor subtype. The rate of internalisation into the AR4-2J cell line also was the highest found for any somatostatin-based radiopeptide. Biodistribution studies, performed in nude mice bearing an AR4-2J tumour or a transfected HEK-sst2 cell-based tumour, showed high and specific uptake in the tumour and in other sst-receptor-expressing tissues, which reflects the high receptor binding affinity and the high rate of internalisation. The pharmacologic differences between [(57)Co(dotatoc)] and [(67)Ga(dotatoc)] are discussed in terms of the structural parameters found for the chelate models [Co(II)(dota)](2-) and [Ga(III)(dota)](-) whose X-ray structures have been determined. Both chelates show six-fold coordination in pseudo-octahedral arrangements.

Fenofibrate metabolism in the cynomolgus monkey using ultraperformance liquid chromatography-quadrupole time-of-flight mass spectrometry-based metabolomics

Fenofibrate, widely used for the treatment of dyslipidemia, activates the nuclear receptor, peroxisome proliferator-activated receptor alpha. However, liver toxicity, including liver cancer, occurs in rodents treated with fibrate drugs. Marked species differences occur in response to fibrate drugs, especially between rodents and humans, the latter of which are resistant to fibrate-induced cancer. Fenofibrate metabolism, which also shows species differences, has not been fully determined in humans and surrogate primates. In the present study, the metabolism of fenofibrate was investigated in cynomolgus monkeys by ultraperformance liquid chromatography-quadrupole time-of-flight mass spectrometry (UPLC-QTOFMS)-based metabolomics. Urine samples were collected before and after oral doses of fenofibrate. The samples were analyzed in both positive-ion and negative-ion modes by UPLC-QTOFMS, and after data deconvolution, the resulting data matrices were subjected to multivariate data analysis. Pattern recognition was performed on the retention time, mass/charge ratio, and other metabolite-related variables. Synthesized or purchased authentic compounds were used for metabolite identification and structure elucidation by liquid chromatographytandem mass spectrometry. Several metabolites were identified, including fenofibric acid, reduced fenofibric acid, fenofibric acid ester glucuronide, reduced fenofibric acid ester glucuronide, and compound X. Another two metabolites (compound B and compound AR), not previously reported in other species, were characterized in cynomolgus monkeys. More importantly, previously unknown metabolites, fenofibric acid taurine conjugate and reduced fenofibric acid taurine conjugate were identified, revealing a previously unrecognized conjugation pathway for fenofibrate.

Determination of menthol in plasma and urine of rats and humans by headspace solid phase microextraction and gas chromatography--mass spectrometry

A method for the determination of menthol and menthol glucuronide (M-G) after enzymatic hydrolysis in plasma and urine of rats and humans was developed using headspace solid phase microextraction and gas chromatography-mass spectrometry in the selected ion monitoring mode (HS-SPME/GC-MS). The assay linearity for plasma ranged from 5 to 1000 ng/ml. The limit of quantification (LOQ) in plasma was 5 ng/ml. The intra- and inter-day precision for menthol and M-G were < or = 18.1% R.S.D. at the LOQ and < or = 4.0% at higher concentrations. Menthol and M-G were determined in rat and human plasma and urine after administration of menthol.

The SLC1 high-affinity glutamate and neutral amino acid transporter family

Glutamate transporters play important roles in the termination of excitatory neurotransmission and in providing cells throughout the body with glutamate for metabolic purposes. The high-affinity glutamate transporters EAAC1 (SLC1A1), GLT1 (SLC1A2), GLAST (SLC1A3), EAAT4 (SLC1A6), and EAAT5 (SLC1A7) mediate the cellular uptake of glutamate by the co-transport of three sodium ions (Na(+)) and one proton (H(+)), with the counter-transport of one potassium ion (K(+)). Thereby, they protect the CNS from glutamate-induced neurotoxicity. Loss of function of glutamate transporters has been implicated in the pathogenesis of several diseases, including amyotrophic lateral sclerosis and Alzheimer's disease. In addition, glutamate transporters play a role in glutamate excitotoxicity following an ischemic stroke, due to reversed glutamate transport. Besides glutamate transporters, the SLC1 family encompasses two transporters of neutral amino acids, ASCT1 (SLC1A4) and ASCT2 (SLC1A5). Both transporters facilitate electroneutral exchange of amino acids in neurons and/or cells of the peripheral tissues. Some years ago, a high resolution structure of an archaeal homologue of the SLC1 family was determined, followed by the elucidation of its structure in the presence of the substrate aspartate and the inhibitor d,l-threo-benzyloxy aspartate (d,l-TBOA). Historically, the first few known inhibitors of SLC1 transporters were based on constrained glutamate analogs which were active in the high micromolar range but often also showed off-target activity at glutamate receptors. Further development led to the discovery of l-threo-β-hydroxyaspartate derivatives, some of which effectively inhibited SLC1 transporters at nanomolar concentrations. More recently, small molecule inhibitors have been identified whose structures are not based on amino acids. Activators of SLC1 family members have also been discovered but there are only a few examples known.

Ion concentration and deduced chemical compounds of the GRIP and Dome Fuji ice cores, Greenland and Antarctica

We have proposed a method of deducing the chemical compounds found in deep polar ice cores by analyzing the balance between six major ions (Cl-, NO3-, SO4**2-, Na+, Mg2+, and Ca2+). The method is demonstrated for the Holocene and last glacial maximum regions of the Dome Fuji and GRIP ice cores. The dominant compounds depend only on the ion balance and the sequence of chemical reactions. In priority order, the principle salts are calcium sulfate, other sulfates, nitrate, chloride, and carbonate. The chemical abundances deduced by this method agree well with the results of Raman spectroscopy on individual salt inclusions. The abundances in the ice cores are shown to reflect differences in climatic periods (the acidic environment of the Holocene versus the reductive environment of the last glacial maximum) and regional conditions (the marine environment of Antarctica versus the continental environment of Greenland).