Membrane Transport and Ion Channels

Introduction

The lipid bilayer is a superb barrier. Its hydrophobic interior effectively blocks the passage of ions and most polar molecules, creating the enclosed chemical environments that cells need to function. But a perfect barrier would be useless — cells must constantly import nutrients, export waste, maintain ion gradients, and transmit electrical signals. The solution lies in an enormous family of membrane transport proteins that provide controlled pathways across the bilayer.

This lesson examines membrane transport at textbook depth: the fundamental physics of permeability, the two major classes of transport proteins (carriers and channels), the ATP-driven pumps that build ion gradients, and the ion channels that exploit those gradients to generate electrical signals. We follow the logic from single pumps and channels through to the integrative physiology of neurons, action potentials, and synaptic transmission. Along the way, we introduce the bioinformatics tools for classifying transporters, predicting drug-transporter interactions, and interpreting channelopathy variants.

11.1 — Principles of Membrane Transport

Protein-Free Lipid Bilayers Are Impermeable to Ions

The permeability of a pure lipid bilayer depends on the size and polarity of the molecule attempting to cross. Small nonpolar molecules (O₂, CO₂, N₂) and very small uncharged polar molecules (water, ethanol) diffuse across readily because they can dissolve in the hydrophobic interior. Larger polar molecules (glucose, amino acids) cross far more slowly, and ions are essentially impermeant — the energetic cost of moving a charged species out of water and into a hydrocarbon environment is prohibitive.

The energy barrier for an ion crossing a pure lipid bilayer is on the order of 40 kcal/mol, giving permeability coefficients of approximately 10−12; cm/s for Na⁺ and K⁺. In biological terms, this means that without transport proteins, it would take hours to days for a biologically meaningful number of ions to cross the membrane. Cells solve this by embedding transport proteins that provide hydrophilic pathways.

Two Classes of Transport Proteins: Carriers and Channels

All membrane transport proteins fall into one of two classes:

Channels are essentially gated pores. When open, selected ions flow through at near-diffusion-limited rates. They can only mediate passive transport — movement down an electrochemical gradient.

Carriers (also called transporters) bind their substrates on one side of the membrane and undergo conformational changes to release them on the other side. This alternating-access mechanism is inherently slower than channel-mediated transport but allows carriers to couple transport to an energy source and drive molecules against their concentration gradient.

Active Transport Must Be Coupled to an Energy Source

Passive transport (facilitated diffusion) moves solutes down their electrochemical gradient and requires no energy input. Active transport moves solutes against their gradient and requires energy. There are two forms:

• Primary active transport — energy comes directly from ATP hydrolysis (e.g., the Na⁺/K⁺-ATPase)
• Secondary active transport — energy comes from the electrochemical gradient of another ion, typically Na⁺ or H⁺, which was itself established by a primary pump

Secondary active transport can be further classified as symport (two solutes move in the same direction) or antiport (two solutes move in opposite directions). For example, the Na⁺/glucose symporter in intestinal epithelial cells uses the Na⁺ gradient to drive glucose uptake against its concentration gradient.

Bioinformatics of Transporter Classification

The Transporter Classification Database (TCDB)

The Transporter Classification Database (TCDB; tcdb.org) provides a comprehensive, evolution-based classification system for membrane transport proteins. Analogous to the Enzyme Commission (EC) system for enzymes, the TC system assigns every transporter a hierarchical identifier with five levels:

1. Class — mechanism of transport (channel, carrier, primary pump, group translocator, or accessory factor)
2. Subclass — energy source or driving force
3. Family — phylogenetic grouping
4. Subfamily — function within the family
5. Individual transporter — specific protein

For example, the human glucose transporter GLUT1 is classified as TC 2.A.1.1.28 — a carrier (class 2), active or facilitated (subclass A), of the Major Facilitator Superfamily (family 1).

Transporter Substrate Prediction and SLC Family Annotation

Predicting the substrate specificity of an uncharacterized transporter is a key bioinformatics challenge. Approaches include:

• Sequence homology — BLAST or HMM searches against TCDB or UniProt to identify close relatives with known substrates
• Machine learning classifiers — trained on structural and sequence features to predict substrate class (amino acid, sugar, ion, drug)
• Structural modeling — homology models or AlphaFold predictions analyzed for binding pocket shape and electrostatics

The Solute Carrier (SLC) superfamily is the largest group of membrane transporters in the human genome, comprising over 450 members in 66 families (SLC1 through SLC66). SLC transporters move amino acids, sugars, nucleotides, inorganic ions, and drugs across membranes. Many are secondary active transporters driven by Na⁺ or H⁺ gradients. The SLC tables maintained by the HUGO Gene Nomenclature Committee (HGNC) provide systematic annotation of substrate specificity, tissue expression, and disease associations for every member.

ABC Transporter Identification and Pharmacogenomics

ATP-Binding Cassette (ABC) transporters constitute a second major superfamily with 48 members in 7 subfamilies (ABCA through ABCG) in humans. Unlike SLC transporters, ABC transporters are primary active pumps driven by ATP hydrolysis. They are identified computationally by the presence of highly conserved nucleotide-binding domains (NBDs) containing the Walker A, Walker B, and ABC signature motifs. Profile HMMs from Pfam (PF00005) provide sensitive detection of ABC transporter genes in any genome.

11.2 — Carrier Proteins and Active Transport

Active Transport Driven by Ion-Concentration Gradients

The Na⁺ gradient established by the Na⁺/K⁺ pump is the primary energy currency for secondary active transport in animal cells. In bacteria and plants, H⁺ gradients serve the same purpose. The free energy stored in an ion gradient can be calculated from the concentration ratio and the membrane potential:

ΔG = RT ln([ion]_in/[ion]_out) + zFΔV

where R is the gas constant, T is temperature, z is the ion charge, F is Faraday’s constant, and ΔV is the membrane potential. For Na⁺ entering a typical animal cell, both terms are negative (concentration gradient and electrical gradient both favor entry), making the total driving force approximately −12 to −14 kJ/mol — sufficient to power the uphill transport of sugars, amino acids, and other solutes.

Carrier Proteins Regulate Cytosolic pH

Cytosolic pH is maintained near 7.2 despite metabolic acid production. Two families of carrier proteins are primarily responsible:

• Na⁺/H⁺ exchangers (NHE) — antiporters that export H⁺ from the cell using the Na⁺ gradient, raising cytosolic pH
• Cl⁻/HCO₃⁻ exchangers — antiporters that import Cl⁻ and export HCO₃⁻, lowering cytosolic pH (also called band 3 or AE1 in red blood cells)

Together, these transporters form a pH buffering system that responds dynamically to acid or base loads. In addition, Na⁺-driven HCO₃⁻ transporters (NBC family) contribute to pH regulation in many cell types. The importance of pH regulation is underscored by the fact that even small deviations from normal pH — as little as 0.1 pH units — can significantly alter enzyme activity and protein folding.

Asymmetric Distribution in Epithelial Cells — Transepithelial Transport

Epithelial cells line the intestine, kidney tubules, and other body cavities. They are polarized, with an apical surface facing the lumen and a basolateral surface facing the blood. Different transport proteins are targeted to each surface, enabling transepithelial transport — the directional movement of solutes across an entire cell layer.

A classic example is glucose absorption in the small intestine:

1. The Na⁺/glucose symporter (SGLT1) on the apical surface uses the Na⁺ gradient to import glucose from the intestinal lumen against its concentration gradient
2. GLUT2 on the basolateral surface allows glucose to exit the cell down its concentration gradient into the blood
3. The Na⁺/K⁺-ATPase on the basolateral surface maintains the Na⁺ gradient that drives step 1

Tight junctions between cells prevent transport proteins from diffusing between apical and basolateral domains, maintaining the asymmetry essential for vectorial transport.

Three Classes of ATP-Driven Pumps

Primary active transporters that use ATP directly fall into three families:

The P-type Ca²⁺ Pump in the Sarcoplasmic Reticulum

The SERCA pump (Sarco/Endoplasmic Reticulum Ca²⁺-ATPase) is a P-type ATPase that maintains the low cytosolic Ca²⁺ concentration (~100 nM) by pumping Ca²⁺ into the ER or sarcoplasmic reticulum (SR) lumen, where the concentration reaches ~0.5 mM — a 5,000-fold gradient. In muscle cells, the SR serves as a Ca²⁺ reservoir; release of stored Ca²⁺ through ryanodine receptor channels triggers contraction, and SERCA pumps it back to allow relaxation.

The SERCA cycle involves phosphorylation of a conserved aspartate residue, which drives large conformational changes between an inward-open state (E1, high Ca²⁺ affinity) and an outward-open state (E2, low Ca²⁺ affinity). Two Ca²⁺ ions are transported per ATP hydrolyzed. Thapsigargin, a plant-derived toxin, is a specific SERCA inhibitor used widely in research.

The Na⁺/K⁺ Pump Establishes the Na⁺ Gradient

The Na⁺/K⁺-ATPase is arguably the most important pump in animal cells. Each cycle hydrolyzes one ATP and transports 3 Na⁺ out and 2 K⁺ in, maintaining the steep gradients that are essential for:

• Electrical excitability (the basis of the resting membrane potential)
• Secondary active transport (Na⁺-driven symporters and antiporters)
• Cell volume regulation (preventing osmotic swelling)

The pump consumes roughly 30% of a cell’s total ATP — rising to 70% in neurons, where intense electrical activity constantly dissipates the Na⁺ and K⁺ gradients. The pump has binding sites for both ouabain (a cardiac glycoside that inhibits the pump) and digitalis (used clinically to treat heart failure by indirectly increasing cardiac contractility through elevated intracellular Na⁺ and, consequently, Ca²⁺).

The typical ion concentrations maintained by the pump are:

ABC Transporters — The Largest Transporter Family

ABC transporters form the largest family of membrane transport proteins, found in all domains of life. Each functional transporter contains two transmembrane domains (TMDs) that form the translocation pathway and two nucleotide-binding domains (NBDs) that hydrolyze ATP. The NBDs dimerize upon ATP binding, driving a conformational change that opens the pathway to the opposite side of the membrane.

In humans, ABC transporters have diverse physiological roles:

• ABCB1 (MDR1/P-glycoprotein) — exports hydrophobic drugs and toxins from cells; a major cause of multidrug resistance in cancer
• ABCC7 (CFTR) — a Cl⁻ channel (unique among ABC transporters); mutations cause cystic fibrosis
• ABCB2/B3 (TAP1/TAP2) — transport peptides into the ER for MHC class I antigen presentation
• ABCA1 — exports cholesterol and phospholipids; mutations cause Tangier disease
• ABCG2 (BCRP) — breast cancer resistance protein; drug efflux pump

ABC Transporters Pump Drugs and Toxins

The role of ABC transporters in multidrug resistance (MDR) is a major clinical problem. Cancer cells that overexpress P-glycoprotein (ABCB1) can pump chemotherapy drugs out of the cell before they reach their intracellular targets, rendering the cancer resistant to multiple structurally unrelated drugs. P-glycoprotein has remarkably broad substrate specificity — its binding pocket can accommodate hundreds of hydrophobic compounds, including vinca alkaloids, anthracyclines, taxanes, and many targeted therapies.

In normal tissues, ABC transporters at the blood-brain barrier, intestinal epithelium, and renal tubules protect the body by limiting drug absorption and promoting drug excretion. This is why the pharmacokinetics of many drugs are critically dependent on transporter activity.

Pharmacogenomics of Transporters

ADME Gene Polymorphism Analysis

Pharmacogenomics studies how genetic variation affects drug response. Transporters are a key component of the ADME pathway (Absorption, Distribution, Metabolism, Excretion). Common polymorphisms in transporter genes can dramatically alter drug efficacy and toxicity:

• ABCB1 (P-glycoprotein) — the synonymous SNP C3435T alters mRNA stability and protein folding, affecting the pharmacokinetics of digoxin, cyclosporine, and HIV protease inhibitors
• SLC22A1 (OCT1) — reduced-function variants decrease the hepatic uptake of metformin, the most widely prescribed diabetes drug
• SLCO1B1 (OATP1B1) — the *5 variant (Val174Ala) impairs statin uptake into hepatocytes, increasing plasma statin levels and the risk of statin-induced myopathy
• ABCG2 (BCRP) — the Q141K variant reduces uric acid excretion and increases gout risk

The Clinical Pharmacogenetics Implementation Consortium (CPIC) and the PharmGKB database provide curated guidelines linking transporter genotypes to drug dosing recommendations.

Drug-Transporter Interaction Prediction

Computational prediction of drug-transporter interactions is essential for drug development. Approaches include:

• Molecular docking — predicting binding poses of drugs in transporter structures (crystal structures or AlphaFold models)
• Quantitative structure-activity relationship (QSAR) models — machine learning on molecular descriptors to predict whether a compound is a substrate, inhibitor, or non-interactor of a specific transporter
• Physiologically based pharmacokinetic (PBPK) modeling — integrating transporter kinetics, tissue expression, and blood flow to predict whole-body drug distribution

The International Transporter Consortium (ITC) has published recommendations identifying key transporters that should be evaluated during drug development, including P-glycoprotein, BCRP, OATP1B1/1B3, OAT1/3, OCT2, MATE1/2-K, and BSEP.

Transporter Expression Across Tissues

The Human Protein Atlas and GTEx (Genotype-Tissue Expression) projects provide transcriptomic and proteomic data for transporter expression across human tissues. These resources reveal that many transporters have highly tissue-specific expression patterns — for example, SGLT2 is almost exclusively expressed in the kidney proximal tubule (explaining why SGLT2 inhibitors, a major class of diabetes drugs, act specifically there), while P-glycoprotein is abundant at barrier tissues (intestine, blood-brain barrier, placenta).

11.3 — Ion Channels

Ion Channels Are Ion-Selective and Fluctuate Between Open and Closed States

Ion channels are transmembrane proteins that form aqueous pores selective for specific ions. Unlike carriers, they do not bind and release their substrates — ions flow through the open pore at near-diffusion rates of up to 10⁸ ions per second. Channels fluctuate stochastically between open and closed conformations. The fraction of time a channel spends open (the open probability, P_o) is regulated by gating signals.

Three major gating mechanisms exist:

• Voltage-gated — respond to changes in membrane potential (Na⁺, K⁺, Ca²⁺ channels in excitable cells)
• Ligand-gated — open when a specific molecule binds (neurotransmitter receptors)
• Mechanically gated — respond to membrane tension or physical force (touch receptors, hair cells)

The Membrane Potential Depends on K⁺ Leak Channels

The resting membrane potential of a typical animal cell is approximately −70 mV (inside negative). It arises primarily from K⁺ leak channels (also called two-pore domain K⁺ channels, or K2P channels) that are constitutively open.

The logic is straightforward: the Na⁺/K⁺ pump maintains a high intracellular K⁺ concentration (~140 mM inside vs. ~5 mM outside). K⁺ leak channels allow K⁺ to flow out of the cell down its concentration gradient. As positive charges leave, the cell interior becomes negative, creating an electrical force that opposes further K⁺ efflux. At equilibrium, the electrical driving force exactly balances the concentration gradient — this is the K⁺ equilibrium potential (E_K), approximately −90 mV. The actual resting potential (−70 mV) is slightly less negative than E_K because the membrane has a small permeability to Na⁺, which tends to depolarize the cell.

The equilibrium potential for any ion can be calculated using the Nernst equation:

E_ion = (RT/zF) × ln([ion]_outside/[ion]_inside)

For K⁺ at 37°C: E_K = (61.5 mV / +1) × log(5/140) ≈ −90 mV.

The Resting Potential Decays Slowly Without the Pump

If the Na⁺/K⁺ pump is poisoned (e.g., with ouabain), the resting potential does not collapse immediately. Because the membrane is far more permeable to K⁺ than to Na⁺, the K⁺ gradient dissipates slowly, and the resting potential is maintained for many minutes. Over time, however, the gradients run down as K⁺ leaks out and Na⁺ leaks in, eventually depolarizing the cell. This demonstrates that the pump is not directly responsible for the resting potential on short timescales — its role is to maintain the ion gradients that the K⁺ leak channels exploit.

K⁺ Channel Structure

The crystal structure of the bacterial KcsA K⁺ channel (solved by Roderick MacKinnon, Nobel Prize 2003) revealed the molecular basis of ion selectivity. The channel is a tetramer, with each subunit contributing one α-helix to a central pore. The selectivity filter — a narrow region near the extracellular surface — is lined with backbone carbonyl oxygen atoms spaced at precisely the distances needed to coordinate a dehydrated K⁺ ion. Na⁺, which is smaller, cannot make energetically favorable contacts with these carbonyls, explaining the channel’s remarkable ~10,000-fold selectivity for K⁺ over Na⁺.

The selectivity filter contains four K⁺ binding sites in single file. Electrostatic repulsion between adjacent ions drives rapid throughput: as a new K⁺ enters the filter, it pushes the ion ahead of it out the other side, achieving transport rates of ~10⁸ ions per second.

Mechanically and Thermally Gated Channels

Mechanosensitive channels open in response to membrane stretch or mechanical force. The Piezo channels (Piezo1, Piezo2) are the primary mechanotransducers in mammals — Piezo1 mediates vascular blood flow sensing, while Piezo2 is essential for light touch sensation, proprioception, and baroreception. Piezo channels are remarkably large, with each subunit containing 38 transmembrane helices arranged in a propeller-like structure.

Thermosensitive channels of the TRP (Transient Receptor Potential) family respond to temperature. TRPV1 is activated by heat (>43°C) and by capsaicin (the active ingredient in chili peppers); TRPM8 is activated by cold (<25°C) and by menthol. These channels allow the nervous system to sense temperature across a wide range.

Aquaporins

While water can diffuse slowly through the lipid bilayer, many cells require much higher water permeability. Aquaporins are tetrameric water channels that allow water to flow at rates of ~3 × 10⁹ molecules per second per channel while completely excluding ions and protons. The exclusion of protons (H₃O⁺) is critical and is achieved by a combination of electrostatic repulsion (positive charges in the channel interior) and reorientation of water molecules as they pass through the pore, breaking the hydrogen-bonded chain that would otherwise conduct protons.

There are 13 aquaporins in humans. AQP1 is abundant in the kidney proximal tubule and red blood cells; AQP2 in the kidney collecting duct is regulated by vasopressin (antidiuretic hormone) — vasopressin triggers insertion of AQP2-containing vesicles into the apical membrane, increasing water reabsorption. Mutations in AQP2 cause nephrogenic diabetes insipidus. Peter Agre received the Nobel Prize in Chemistry in 2003 for the discovery of aquaporins.

Nerve Cells Have an Elongated Structure

Neurons are specialized for rapid, long-distance electrical signaling. A typical motor neuron has a cell body (soma) containing the nucleus, branching dendrites that receive signals, and a long axon (up to 1 meter in humans) that transmits signals to distant targets. The axon terminates in synaptic boutons that form connections (synapses) with other neurons or muscle cells. The highly elongated structure of neurons means that simple diffusion is hopelessly slow for communication over distances — a chemical signal would take years to diffuse from the spinal cord to the foot. Instead, neurons use electrical signals (action potentials) that propagate at speeds up to 120 m/s.

Voltage-Gated Channels Generate Action Potentials

Action potentials are rapid, all-or-nothing electrical signals generated by the coordinated opening and closing of voltage-gated Na⁺ and K⁺ channels:

1. A stimulus depolarizes the membrane to threshold (~−55 mV)
2. Voltage-gated Na⁺ channels open rapidly (within 0.1 ms) → Na⁺ floods in → membrane depolarizes to ~+40 mV
3. Na⁺ channels undergo inactivation (a separate gate blocks the pore from the cytoplasmic side within ~1 ms)
4. Voltage-gated K⁺ channels open (with a slight delay) → K⁺ flows out → membrane repolarizes
5. K⁺ channels close slowly → brief hyperpolarization (undershoot below −70 mV)
6. Resting potential is restored as K⁺ channels close completely

The entire action potential lasts about 1–2 milliseconds. The refractory period (during which Na⁺ channels are inactivated and cannot reopen) ensures that action potentials propagate in one direction and limits their maximum frequency to about 500–1,000 per second.

Only a tiny fraction of the total Na⁺ and K⁺ in the cell crosses the membrane during a single action potential — approximately 1 in 100,000 ions — so the gradients are barely perturbed and thousands of action potentials can fire before the pump must fully restore them.

Myelination Increases Speed and Efficiency

In vertebrates, many axons are wrapped in myelin — an insulating sheath of tightly packed membrane formed by Schwann cells (peripheral nervous system) or oligodendrocytes (central nervous system). Myelin prevents ion flow across the axonal membrane except at the nodes of Ranvier — small gaps (~1 μm) between myelin segments where voltage-gated Na⁺ channels are concentrated.

The action potential “jumps” from node to node in a process called saltatory conduction, which increases speed from ~1 m/s (unmyelinated) to up to 120 m/s (myelinated) while also reducing the metabolic cost by limiting the membrane area over which ions must be pumped back. Demyelinating diseases such as multiple sclerosis disrupt saltatory conduction, causing neurological dysfunction.

Patch-Clamp Recording

The patch-clamp technique, developed by Erwin Neher and Bert Sakmann (Nobel Prize 1991), revolutionized electrophysiology by allowing the recording of current through single ion channels in real time. A glass micropipette with a polished tip (~1 μm) is pressed against the cell membrane to form a tight seal (gigaohm seal, resistance >10 GΩ). In the cell-attached configuration, one can record the opening and closing of individual channels in the patch of membrane under the pipette. Other configurations (whole-cell, inside-out, outside-out) allow different experimental questions to be addressed.

Patch-clamp recordings reveal that individual channels switch between discrete conductance states — they are either fully open or fully closed, with no intermediate states. The time spent in each state follows stochastic kinetics. The macroscopic current recorded from many channels is the sum of these individual stochastic events.

Voltage-Gated Na⁺, K⁺, and Ca²⁺ Channels Are Structurally Related

The major voltage-gated ion channels (Na_v, K_v, Ca_v) share a common evolutionary origin and structural plan:

All contain a voltage sensor (the S4 helix, with positively charged arginine residues at every third position) that moves in response to changes in membrane potential, mechanically opening or closing the channel gate.

Transmitter-Gated Ion Channels at Synapses

At chemical synapses, the presynaptic neuron releases neurotransmitters from synaptic vesicles into the synaptic cleft (~20 nm wide). Neurotransmitters bind ligand-gated ion channels (also called ionotropic receptors) on the postsynaptic cell, directly opening an ion channel and producing a rapid electrical response.

Major transmitter-gated channels include:

Excitatory and Inhibitory Synapses

Whether a synapse is excitatory or inhibitory depends on the ion channels it activates:

• Excitatory synapses open channels permeable to Na⁺ (and sometimes Ca²⁺), depolarizing the postsynaptic cell toward threshold. The resulting depolarization is called an excitatory postsynaptic potential (EPSP).
• Inhibitory synapses open channels permeable to Cl⁻ or K⁺, hyperpolarizing the cell or stabilizing it near the resting potential. The resulting change is an inhibitory postsynaptic potential (IPSP).

The postsynaptic cell integrates EPSPs and IPSPs from thousands of synapses — a typical neuron in the brain receives input from 1,000 to 10,000 other neurons. If the sum of all inputs at the axon hillock exceeds threshold, an action potential is initiated. This process of synaptic summation is the fundamental computation of the nervous system.

The Neuromuscular Junction

The neuromuscular junction (NMJ) is the synapse between a motor neuron and a skeletal muscle fiber. It is the best-characterized synapse in biology. Acetylcholine released from the motor nerve terminal binds nicotinic acetylcholine receptors on the muscle fiber, opening channels that admit Na⁺, depolarizing the muscle cell (the end-plate potential). The end-plate potential is normally large enough to exceed threshold reliably, ensuring faithful transmission.

Myasthenia gravis is an autoimmune disease in which antibodies attack the nicotinic acetylcholine receptor, reducing the number of functional receptors and causing muscle weakness. Curare and α-bungarotoxin (from snake venom) block the receptor, causing paralysis.

Neurotransmitter Receptors as Drug Targets

Neurotransmitter receptors are among the most important drug targets in medicine:

• Benzodiazepines (diazepam, lorazepam) enhance GABA_A receptor function → anti-anxiety, sedation
• Barbiturates also potentiate GABA_A receptors → anesthesia, anticonvulsant
• Ketamine blocks NMDA receptors → dissociative anesthesia, rapid-acting antidepressant
• Memantine blocks NMDA receptors → treatment of Alzheimer’s disease
• Nicotine activates nicotinic acetylcholine receptors → addiction, cognitive enhancement
• SSRIs (fluoxetine/Prozac) block serotonin reuptake, prolonging serotonergic signaling

Approximately 19% of all FDA-approved drugs target ion channels or ionotropic receptors, making this protein class second only to GPCRs as a drug target family.

Synaptic Summation and Neural Computation

Each neuron acts as an integrator, summing excitatory and inhibitory inputs both in space (spatial summation — inputs arriving simultaneously from different synapses) and time (temporal summation — inputs arriving in rapid succession at the same synapse). The decision to fire an action potential depends on whether the sum of all inputs depolarizes the axon hillock past threshold.

This summation is not simple addition. The geometry of the dendritic tree, the distance of each synapse from the axon hillock, the time constants of the postsynaptic potentials, and the presence of voltage-gated channels in the dendrites all shape how inputs are integrated. Understanding these computations is central to neuroscience.

Long-Term Potentiation Depends on Ca²⁺ Entry Through NMDA Receptors

Long-term potentiation (LTP) is a persistent strengthening of synaptic transmission that is widely regarded as a cellular mechanism of learning and memory. LTP at glutamatergic synapses in the hippocampus depends on the NMDA receptor, which has unique properties:

1. It requires both glutamate binding and postsynaptic depolarization to open (the channel is blocked by Mg²⁺ at resting potential)
2. When open, it is permeable to Ca²⁺ in addition to Na⁺ and K⁺

The NMDA receptor thus acts as a coincidence detector — it opens only when the presynaptic neuron is active (releasing glutamate) and the postsynaptic neuron is depolarized (removing the Mg²⁺ block). The resulting Ca²⁺ influx activates CaMKII (calcium/calmodulin-dependent protein kinase II) and other signaling cascades that increase the number and conductance of AMPA receptors at the synapse, strengthening transmission.

Ion Channel and Neuroscience Informatics

Ion Channel Databases

Several curated databases support ion channel research:

• IUPHAR/BPS Guide to Pharmacology (guidetopharmacology.org) — the authoritative resource for ion channel nomenclature, pharmacology, and ligand interactions, maintained by the International Union of Basic and Clinical Pharmacology
• ChannelPedia (channelpedia.epfl.ch) — detailed kinetic models and electrophysiological data for individual ion channels
• Ion Channel Genealogy (icg.neurotheory.ox.ac.uk) — traces the evolutionary relationships among computational models of ion channels

Ion Channel Homology Modeling and Structure Prediction

With the explosion of cryo-EM structures and AlphaFold predictions, high-quality structural models are now available for most human ion channel families. Homology modeling and molecular dynamics simulations are used to:

• Predict the binding sites of drugs and toxins
• Understand the structural basis of channelopathy mutations
• Model ion permeation and selectivity at atomic resolution
• Design novel channel modulators for drug development

Electrophysiology Data Analysis

Modern electrophysiology generates large datasets — automated patch-clamp platforms can record from hundreds of cells per day. Computational pipelines include:

• Single-channel analysis — detecting opening and closing events, fitting dwell-time distributions to Markov models
• Spike sorting — separating action potentials from different neurons in multi-electrode recordings
• Current-voltage (I-V) curve fitting — characterizing channel conductance, rectification, and voltage dependence

Neural Network Modeling

Computational neuroscience uses detailed biophysical simulations to understand neural circuit function:

• NEURON (neuron.yale.edu) — the most widely used simulator for biophysically detailed neuron models, incorporating Hodgkin-Huxley-type channel kinetics and realistic morphologies
• Brian (briansimulator.org) — a Python-based simulator for spiking neural networks, optimized for large-scale network simulations
• The Human Brain Project and Allen Brain Atlas provide connectomics data (wiring diagrams) and gene expression maps that constrain and validate computational models

Channelopathy Variant Interpretation

Channelopathies are diseases caused by mutations in ion channel genes. Examples include long QT syndrome (K⁺ or Na⁺ channel mutations causing cardiac arrhythmia), epilepsy (Na⁺ or GABA_A channel mutations), and cystic fibrosis (CFTR Cl⁻ channel). Variant interpretation combines:

• ClinVar and HGMD annotations for known pathogenic variants
• In silico prediction tools (SIFT, PolyPhen-2, CADD) for novel variants
• Functional domain mapping — variants in the selectivity filter, voltage sensor, or pore-lining helices are more likely to be pathogenic
• Electrophysiological characterization of mutant channels in heterologous expression systems (the gold standard)

The growing availability of population-scale sequencing (gnomAD) allows assessment of variant frequency — common variants are less likely to cause severe disease, while ultra-rare variants in functionally critical domains are strong candidates.

Interactive Exploration

Ion channels achieve extraordinary selectivity through short stretches of amino acids called selectivity filters. Let us compare the properties of the K+ and Na+ channel selectivity filters to understand how their chemistry determines which ions can pass.

let k_channel_filter = "TVGYG"
let na_channel_filter = "DEKA"
print("K+ channel selectivity filter:")
print(Struct.protein_props(k_channel_filter))
print("Na+ channel selectivity filter:")
print(Struct.protein_props(na_channel_filter))

The K+ filter (TVGYG) is largely uncharged and moderately hydrophobic, using backbone carbonyls to coordinate dehydrated K+ ions. The Na+ filter (DEKA) contains charged residues (aspartate, glutamate, lysine) that create an electrostatic environment favoring the smaller Na+ ion.

The Na⁺/K⁺ pump and other primary active transporters maintain steep ion concentration gradients across the membrane. These gradients store energy that drives secondary active transport and electrical signaling. Let us visualize the dramatic asymmetry:

let ion_data = '[{"label": "K+ (inside)", "value": 140}, {"label": "K+ (outside)", "value": 5}, {"label": "Na+ (inside)", "value": 12}, {"label": "Na+ (outside)", "value": 145}, {"label": "Ca2+ (inside)", "value": 0.0001}, {"label": "Ca2+ (outside)", "value": 1.8}]'
let chart = Viz.bar(ion_data, '{"title": "Ion Concentrations (mM)", "color": "#06B6D4"}')
print(chart)

Note the massive gradients: K+ is ~28-fold higher inside cells, Na+ is ~12-fold higher outside, and Ca2+ is maintained at a staggering 15,000-fold gradient. These gradients represent a substantial energy investment — the Na+/K+ pump alone consumes ~30% of cellular ATP.

Membrane transport proteins must recognize and move specific molecules. The physical and chemical properties of transported substrates determine whether they can cross the lipid bilayer freely or require dedicated transporters. Let us compare glucose (which requires a transporter) with ATP (the energy currency that powers active transport):

let glucose = "OC[C@H]1OC(O)[C@H](O)[C@@H](O)[C@@H]1O"
let atp = "c1nc(N)c2ncn(C3OC(COP(=O)(O)OP(=O)(O)OP(=O)(O)O)C(O)C3O)c2n1"
print("Glucose properties:")
print(Chem.properties(glucose))
print("ATP properties:")
print(Chem.properties(atp))

Glucose is a polar, hydrophilic molecule that cannot cross the lipid bilayer unaided — it requires GLUT transporters (passive) or SGLT symporters (active, Na+-driven). ATP is even larger and more polar, with multiple phosphate groups carrying negative charges at physiological pH, making it essentially impermeant to membranes.

Exercise: Predict Channel Selectivity from Filter Sequence

Ion channels achieve extraordinary selectivity through their selectivity filter — a short stretch of amino acids that determines which ions can pass. Compare the properties of K+ and Na+ channel filters:

let k_filter = "TVGYG"
let na_filter = "DEKA"
let k_props = Struct.protein_props(k_filter)
let na_props = Struct.protein_props(na_filter)
print("K+ channel filter:")
print(k_props)
print("Na+ channel filter:")
print(na_props)
// Which channel's filter is more hydrophobic?
let answer = "K+"
print(answer)

Exercise: Visualize Transport Rates

Different transport mechanisms move solutes at vastly different rates. Create a visualization comparing the transport rates of channels, carriers, and pumps:

let rates = '[{"label": "Ion channels", "value": 10000000}, {"label": "Carriers", "value": 1000}, {"label": "ATP pumps", "value": 100}]'
let chart = Viz.bar(rates, '{"title": "Transport Rate (ions/sec)", "color": "#10B981"}')
print(chart)
// Which mechanism has the highest transport rate?
let answer = "channels"
print(answer)

Exercise: Molecular Properties and Membrane Permeability

A molecule’s ability to cross the lipid bilayer depends on its size, charge, and hydrophobicity. Compare glucose (needs a transporter) with ethanol (crosses freely):

let glucose = "OC[C@H]1OC(O)[C@H](O)[C@@H](O)[C@@H]1O"
let ethanol = "CCO"
print("Glucose:")
print(Chem.properties(glucose))
print("Ethanol:")
print(Chem.properties(ethanol))
// Which molecule can freely cross the lipid bilayer?
let answer = "ethanol"
print(answer)

Knowledge Check

Summary

In this lesson you covered membrane transport and ion channels in depth:

• Protein-free lipid bilayers are essentially impermeable to ions — transport proteins provide the only physiologically relevant pathways
• Channels form gated aqueous pores for rapid passive transport (~10⁸ ions/sec); carriers undergo conformational changes and can couple to energy sources for active transport
• Primary active transport (P-type ATPases, V-type ATPases, ABC transporters) is driven directly by ATP hydrolysis; secondary active transport uses ion gradients established by primary pumps
• The Na⁺/K⁺-ATPase pumps 3 Na⁺ out and 2 K⁺ in per ATP, consuming ~30% of cellular energy and establishing the gradients essential for electrical signaling, secondary transport, and volume regulation
• ABC transporters are the largest transporter family; P-glycoprotein (MDR1) causes multidrug resistance in cancer; CFTR mutations cause cystic fibrosis
• Epithelial cells use asymmetric transporter distribution for directional (transepithelial) transport of nutrients
• The resting membrane potential (~−70 mV) is set primarily by K⁺ leak channels and the K⁺ concentration gradient
• Voltage-gated Na⁺ and K⁺ channels generate action potentials — rapid, all-or-nothing electrical signals that propagate along axons
• Myelination enables saltatory conduction, increasing speed up to 120 m/s and reducing metabolic cost
• Ligand-gated ion channels at synapses mediate fast excitatory (glutamate, acetylcholine) and inhibitory (GABA, glycine) neurotransmission
• Synaptic summation of thousands of EPSPs and IPSPs is the fundamental computation of the nervous system
• NMDA receptors act as coincidence detectors; their Ca²⁺ permeability drives long-term potentiation (LTP), a cellular basis of learning and memory
• Bioinformatics resources (TCDB, PharmGKB, IUPHAR, ChannelPedia) support transporter classification, pharmacogenomics, and channelopathy variant interpretation
• Pharmacogenomic variation in transporter genes (ABCB1, SLCO1B1, SLC22A1) significantly affects drug efficacy and toxicity