Ligand-gated ion channels couple intra- and extracellular chemical signals to cellular excitability. In response to a specific ligand, these channels change their permeability to certain ions by opening or closing their ion conductive pathway, a controlling mechanism known as channel gating. Although recent studies with X-ray crystallography and site-directed mutagenesis have revealed several structures potentially important for channel gating, the gating mechanism is still elusive. Ligand-dependent channel gating involves a series of transient events and asymmetric movements of individual subunits. Understanding of these events appears to be a challenge to current approaches in gating studies by using the homomeric wild-type or mutant channels. I therefore took an alternative approach by constructing heteromeric channels. Subunit stoichiometric studies of the Kir1.1 channel showed that a minimum of one functional subunit was required for the pH-dependent gating of the channel. Four subunits in this channel were coordinated as dynamic functional dimers. In Kir6.2 channel, stoichiometry for proton-binding was almost identical to that for channel gating in the M2 helix, suggesting a one-to-one direct coupling of proton binding in C-terminus to channel gating in M2 helix. Positive cooperativity was suggested among subunits in both the proton binding and channel gating. Ligand binding can be differentiated from channel gating by studying the ATP-dependent gating of Kir6.2 channel. Disruptions in ATP binding were found to change both the potency and efficacy of the concentration-dependent curves, while the baseline activity instead of maximum inhibition was affected by disruptions of channel gating. Four subunits in the Kir6.2 channel undergo negative cooperativity in ATP binding and positive cooperativity in channel gating. The ligand binding was coupled to the gating mechanism in the same subunit and neighboring subunits, although the intrasubunit coupling was more effective. These results are well described with the operational model which we have applied to ion channel studies for the first time. By manipulating the relative distance and the interaction of two transmembrane helices, the inner helix bundle of crossing was found to not only serve as a gate but also determine the consequence of ligand binding.

No description available.

Thesis (Ph. D.)--University of Missouri--Columbia, 1998. / Typescript. Vita. Includes bibliographical references (l. 140-148). Also available on the Internet.

KCNQ1 is a homotetrameric voltage-gated potassium channel expressed in cardiomyocytes and epithelial tissues. However, currents arising from KCNQ1 have never been physiologically observed. KCNQ1 is able to provide the diverse potassium conductances required by these distinct cell types through coassembly with and modulation by type I transmembrane β-subunits of the KCNE gene family. KCNQ1-KCNE K+ channels play important physiological roles. In cardiac tissues the association of KCNQ1 with KCNE1 gives rise to IKs, the slow delayed outwardly rectifying potassium current. IKs is in part responsible for repolarizing heart muscle, and is therefore crucial in maintaining normal heart rhymicity. IKschannels help terminate each action potential and provide cardiac repolarization reserve. As such, mutations in either subunit can lead to Romano-Ward Syndrome or Jervell and Lange-Nielsen Syndrome, two forms of Q-T prolongation. In epithelial cells, KCNQ1-KCNE1, KCNQ1-KCNE2 and KCNQ1-KCNE3 give rise to potassium currents required for potassium recycling and secretion. These functions arise because the biophysical properties of KCNQ1 are always dramatically altered by KCNE co-expression. We wanted to understand how KCNE peptides are able to modulate KCNQ1. In Chapter II, we produce partial truncations of KCNE3 and demonstrate the transmembrane domain is necessary and sufficient for both assembly with and modulation of KCNQ1. Comparing these results with published results obtained from chimeric KCNE peptides and partial deletion mutants of KCNE1, we propose a bipartite modulation residing in KCNE peptides. Transmembrane modulation is either active (KCNE3) or permissive (KCNE1). Active transmembrane KCNE modulation masks juxtamembranous C-terminal modulation of KCNQ1, while permissive modulation allows C-terminal modulation of KCNQ1 to express. We test our hypothesis, and demonstrate C-terminal Long QT point mutants in KCNE1 can be masked by active trasnsmembrane modulation. Having confirmed the importance the C-terminus of KCNE1, we continue with two projects designed to elucidate KCNE1 C-terminal structure. In Chapter III we conduct an alanine-perturbation scan within the C-terminus. C-terminal KCNE1 alanine point mutations result in changes in the free energy for the KCNQ1-KCNE1 channel complex. High-impact point mutants cluster in an arrangement consistent with an alphahelical secondary structure, "kinked" by a single proline residue. In Chapter IV, we use oxidant-mediated disulfide bond formation between non-native cysteine residues to demonstrate amino acid side chains residing within the C-terminal domain of KCNE1 are close and juxtaposed to amino acid side chains on the cytoplasmic face of the KCNQ1 pore domain. Many of the amino acids identified as high impact through alanine perturbation correspond with residues identified as able to form disulfide bonds with KCNQ1. Taken together, we demonstrate that the interaction between the C-terminus of KCNE1 and the pore domain of KCNQ1 is required for the proper modulation of KCNQ1 by KCNE1, and by extension, normal IKs function and heart rhymicity.

GLIC är en pH-känslig pentamerisk ligandstyrd jonkanal (pLGIC) som finns i cellmembranet hos prokaryoten Gloeobacter violaceus. GLIC är en bakteriell homolog till flera receptorer som är viktiga i nervsystemet hos de flesta eukaryotiska organismer. Dessa receptorer fungerar som mallar för utvecklingen av målstyrda bedövnings- och stimulerande läkemedel som påverkar nervsystemet. Förståelsen av ett proteins mekanismer har därför hög prioritet inför läkemedelsutvecklingen. Eukaryota pLGICs är dock mycket komplexa eftersom några av de är heteromera, har flera domäner, och de pågår eftertranslationella ändringar. GLIC, å andra sidan, har en enklare struktur och det räcker att analysera strukturen av en subenhet - eftersom alla subenheter är helt lika. Flertalet möjliga grindmekanismer föreslogs av vetenskapen men riktiga öppningsmekanismen av GLIC är fortfarande oklar. Projektets mål är att genomföra maskininlärning (ML) för att upptäcka nya grindmekanismer med hjälp av datormetoder. Urspungsdatan togs från tidigare forskning där andra ML-redskap såsom molekyldynamik (MD), elastisk nätverksstyrd Brownsk dynamik (eBDIMS) och Markovstillståndsmodeller (MSM) användes. Utifrån dessa redskap simulerades proteinet som vildtyp samt med funktionsförstärkt mutation vid två olika pH värden. Fem makrotillstånd byggdes: två öppna, två stängda och ett mellanliggande. I projektet användes ett annat ML redskap: KL-divergens. Detta redskap användes för att hitta skillnader i avståndfördelning mellan öppet och stängt makrotillstånd. Utifrån ursprungsdatan byggdes en tensor som lagrade alla parvisa aminosyrornas avstånd. Varje aminosyrapar hade sin egen metadata som i sin tur användes för att frambringa alla fem avståndsfördelningar fråm MSMs som byggdes i förväg. Sedan bräknades medel-KL-divergens mellan två avståndfördelningar av intresse för att filtrera bort aminosyropar med överlappande avståndsfördelningar. För att se till att aminosyror inom aminosyrapar som låg kvar kan påverka varandra, filtrerades bort alla par vars minsta och medelavstånd var stora. De kvarvarande aminosyroparen utvärderades i förhållande till alla fem makrotillstånd Viktiga nya grindmekanismer som hittades genom både KL-divergens och makrotillståndsfördelningar innefattade loopen mellan M2-M3 helixarna av en subenhet och både loopen mellan sträckor β8 och β9 (Loop F)/N-terminal β9-sträckan och pre-M1/N-terminal M1 av närliggande subenheten. Loopen mellan sträckor β8 och β9 (Loop F) visade höga KL-värden också med loopen mellan sträckor β1 och β2 loop samt med loopen mellan sträckor β6 och β7 (Pro-loop) och avståndet mellan aminosyror minskade vid kanalens grind. Övriga intressanta grindmekanismer innefattade parning av aminosyror från loopen β4-β5 (Loop A) med aminosyror från sträckor β1 och β6 samt böjning av kanalen porangränsande helix. KL-divergens påvisades vara ett viktigt redskap för att filtrera tillgänglig data och de nya grindmekanismer kan bli användbara både för akademin, som vill reda ut GLIC:s fullständiga grindmekanismer, och läkemedelsföretag, som letar efter bindningsställen inom molekylen för att utveckla nya läkemedel. / GLIC is a transmembrane proton-gated pentameric ligand-gated ion channel (pLGIC) that is found in the prokaryote Gloeobacter violaceus. GLIC is the prokaryotic homolog to several receptors that are found in the nervous system of many eukaryotic organisms. These receptors are targets for the development of pharmaceutical drugs that interfere with the gating of these channels - such drugs involve anesthetics and stimulants. Understanding the mechanism of a drug’s target is a high priority for the development of a novel medicine. However, eukaryotic pLGICs are complex to analyse, because some of them are heteromeric, have more domains, and because of their post-translational modifications (PTMs). GLIC, on the other hand, has a simpler structure and it is enough to study the structure of only one subunit - since all subunits are identical. Several possible gating mechanisms have been proposed by the scientific community, but the complete gating of GLIC remains unclear. The goal of this project is to implement machine learning (ML) to discover novel gating mechanisms by computational approaches. The starting data was extracted from a previous research where computational tools like unbiased molecular dynamics (MD), elastic network-driven Brownian Dynamics (eBDIMS), and Markov state models (MSMs) were used. From those tools, the protein was simulated in wild-type and in a gain-of-function mutation at two different pH values. Five macrostates were constructed: two open, two closed, and an intermediate. In this project another ML tool was used: KL divergence. This tool was used to score the difference between the distance distributions of one open and one closed macrostate. The starting data was used to create a tensor that stored all residue-residue distances. Each residue pair had its own metadata, which in turn was used to yield the distance distributions of all five pre-build MSMs. Then the average KL scores between two states of interest were calculated and were used to filter out the residue pairs with overlapping distance distributions. To make sure that the residues within a pair can interact with each other, all residue pairs with very high minimum and average distance were filtered out as well. The residue pairs that remained were later evaluated across all five macrostates for further studies. Important novel mechanisms discovered in this project through both the KL divergence and the macrostate distributions involved the M2-M3 loop of one subunit and both the β8-β9 loop/N-terminal β9 strand and the preM1/N-terminal M1 region of the neighboring subunit. The β8-β9 loop (Loop F) showed high KL scores with the β1-β2 and β6-β7 (Pro-loop) loops as well with decreasing distances upon the channel’s opening. Other notable gating mechanisms involved are the pairing of residues from the β1-β2 loop (Loop A) with residues from the strands β1 and β6, as well as the kink of the pore-lining helix. KL divergence proved a valuable tool to filter available data and the novel mechanisms can prove useful both to the academic community that seeks to unravel the complete gating mechanism of GLIC and to the pharmaceutical companies that search for new binding sites within the molecule for new drugs.

Computational methods, including homology modelling, in-silico dockings, and molecular dynamics simulations have been used to study the functional dynamics and interactions of K⁺ channels. Molecular models were built of the inwardly rectifying K⁺ channel Kir2.2, the bacterial homolog K⁺ channel KirBac3.1, and the twin pore (K2P) K⁺ channels TREK-1 and TRESK. To investigate the electrostatic energy profile of K⁺ permeating through these homology models, continuum electrostatic calculations were performed. The primary mechanism of KirBac3.1 gating is believed to involve an opening at the helix bundle crossing (HBC). However, simulations of Kir channels have not yet revealed opening at the HBC. Here, in simulations of the new KirBac3.1-S129R X-ray crystal structure, in which the HBC was trapped open by the S129R mutation in the inner pore-lining helix (TM2), the HBC was found to exhibit considerable mobility. In a simulation of the new KirBac3.1-S129R-S205L double mutant structure, if the S129R and the S205L mutations were converted back to the wild-type serine, the HBC would close faster than in the simulations of the KirBac3.1-S129R single mutant structure. The double mutant structure KirBac3.1-S129R-S205L therefore likely represents a higher-energy state than the single mutant KirBac3.1-S129R structure, and these simulations indicate a staged pathway of gating in KirBac channels. Molecular modelling and MD simulations of the Kir2.2 channel structure demonstrated that the HBC would tend to open if the C-linker between the transmembrane and cytoplasmic domain was modelled helical. The electrostatic energy barrier for K⁺ permeation at the helix bundle crossing was found to be sensitive to subtle structural changes in the C-linker. Charge neutralization or charge reversal of the PIP2-binding residue R186 on the C-linker decreased the electrostatic barrier for K⁺ permeation through the HBC, suggesting an electrostatic contribution to the PIP2-dependent gating mechanism. Multi-scale simulations determined the PIP2 binding site in Kir2.2, in good agreement with crystallographic predictions. A TREK-1 homology model was built, based on the TRAAK structure. Two PIP2 binding sites were found in this TREK-1 model, at the C-terminal end, in line with existing functional data, and between transmembrane helices TM2 and TM3. The TM2-TM3 site is in reasonably good agreement with electron density attributed to an acyl tail in a recently deposited TREK-2 structure.

Ion channels are fundamental molecules in the nervous system that catalyze the flux of ions across the cell membrane. There are mounting evidences suggesting that the kinetic properties of ion channels undergo activity-dependent changes in various pathophysiological conditions. Here such activity-dependent changes were studied in case of two different ion channels; the rat brain derived voltage-gated Na+ channel, rNav1.2 and the human background leak K+ channel, hTREK1 using the single channel patch-clamp technique. Our results on the voltage-gated Na+ channel (Chapter III) illustrated that sustained membrane depolarization, as seen in pathophysiological conditions like epilepsy, induced a defined non-linear variation in the unitary conductance, activation, inactivation and recovery kinetic properties of the channel. Signal processing tools attributed a pseudo-oscillatory nature to the non-linearity observed in the channel properties. Prolonged membrane depolarization also induced a “molecular memory” phenomenon, characterized by clustering of dwell time events and strong autocorrelation in the dwell time series. The persistence of such molecular memory was found to be dependent on the duration of depolarization. Similar plastic changes were observed in case of the hTREK1 channel in presence of saturating concentrations of agonist, trichloroethanol (TCE) (Chapter IV). TREK1 channel behaves similar to single enzyme molecules with a single binding site for the substrate K+ ion whereas TCE acts as an allosteric activator of the channel. We observed that with increasing concentration of TCE (10 M to 10 mM) the catalytic turnover rate exhibited progressive departure from monoexponential to multi-exponential distribution suggesting the presence of ‘dynamic disorder’ analogous to single enzyme molecules. In addition, we observed the induction of strong correlation in successive waiting times and flux intensities, exemplified by distinct mode switching between high and low flux activity, which implied the induction of memory in single ion channel. Our observation of such molecular memory in two different ion channels in different experimental conditions highlights the importance and generality of the phenomenon which is normally hidden under the ensemble behaviour of ion channels. In the final part of the work (chapter V) we observed strong negative cooperativity and half-of-sites saturation kinetics in the interaction of local anesthetic, lidocaine with hTREK1 channel. We also mapped the specific anesthetic binding site in the c-terminal domain of the channel. Further, single channel analysis and the heterodimer studies enabled us to propose a model for this interaction and provide a plausible paradigm for the inhibitory action of lidocaine on hTREK1.

Thesis (Ph. D.)--University of Virginia, 2008. / Title from title page. Includes bibliographical references. Also available online through Digital Dissertations.