Efficacy of external tetraethylammonium block of the KcsA potassium channel: molecular and Brownian dynamics studies.

Blockade of the KcsA potassium channel by externally applied tetraethylammonium is investigated using molecular dynamics calculations and Brownian dynamics simulations. In KcsA, the aromatic rings of four tyrosine residues located just external to the selectivity filter create an attractive energy well or a binding cage for a tetraethylammonium molecule. We first investigate the effects of re-orienting the four tyrosine residues such that the centers of the aromatic rings face the tetraethylammonium molecule directly. Then, we systematically move the residues inward in both orientations so that the radius of the binding cage formed by them becomes smaller. For each configuration, we construct a one-dimensional free energy profile by bringing in a tetraethylammonium molecule from the external reservoir toward the selectivity filter. The free energy profile is then converted to a one-dimensional potential energy profile, taking the available space between the tyrosine residues and the tetraethylammonium molecule into account. Incorporating this potential energy profile into the Brownian dynamics algorithm, we determine the conductance properties of the channel under various conditions, construct the current-tetraethylammonium-concentration curve and compare it with the experimentally determined inhibitory constant k(i) for externally applied tetraethylammonium. We show that the experimentally determined binding affinity for externally applied tetraethylammonium can be replicated when each of the four tyrosine residues is moved inward by about 0.7 angstroms, irrespective of orientation of their aromatic rings.

Conduction of Na+ and K+ through the NaK channel: molecular and Brownian dynamics studies.

Conduction of ions through the NaK channel, with M0 helix removed, was studied using both Brownian dynamics and molecular dynamics. Brownian dynamics simulations predict that the truncated NaK has approximately a third of the conductance of the related KcsA K+ channel, is outwardly rectifying, and has a Michaelis-Menten current-concentration relationship. Current magnitude increases when the glutamine residue located near the intracellular gate is replaced with a glutamate residue. The channel is blocked by extracellular Ca2+. Molecular dynamics simulations show that, under the influence of a strong applied potential, both Na+ and K+ move across the selectivity filter, although conduction rates for Na+ ions are somewhat lower. The mechanism of conduction of Na+ differs significantly from that of K+ in that Na+ is preferentially coordinated by single planes of pore-lining carbonyl oxygens, instead of two planes as in the usual K+ binding sites. The water-containing filter pocket resulting from a single change in the selectivity filter sequence (compared to potassium channels) disrupts several of the planes of carbonyl oxygens, and thus reduces the filter's ability to discriminate against sodium.

The fast gating mechanism in ClC-0 channels.

We investigate and then modify the hypothesis that a glutamate side chain acts as the fast gate in ClC-0 channels. We first create a putative open-state configuration of the prokaryotic ClC Cl- channel using its crystallographic structure as a basis. Then, retaining the same pore shape, the prokaryotic ClC channel is converted to ClC-0 by replacing all the nonconserved polar and charged residues. Using this open-state channel model, we carry out molecular dynamics simulations to study how the glutamate side chain can move between open and closed configurations. When the side chain extends toward the extracellular end of the channel, it presents an electrostatic barrier to Cl- conduction. However, external Cl- ions can push the side chain into a more central position where, pressed against the channel wall, it does not impede the motion of Cl- ions. Additionally, a proton from a low-pH external solution can neutralize the extended glutamate side chain, which also removes the barrier to conduction. Finally, we use Brownian dynamics simulations to demonstrate the influence of membrane potential and external Cl- concentration on channel open probability.