The implementation of slab geometry for membrane-channel molecular dynamics simulations

Exterior site occupancy infers chloride-induced proton gating in a prokaryotic homolog of the ClC chloride channel

The ClC family of anion channels mediates the efficient, selective permeation of Cl(-) across the biological membranes of living cells under the driving force of an electrochemical gradient. In some eukaryotes, these channels are known to exhibit a unique gating mechanism, which appears to be triggered by the permeant Cl(-) anion. We infer details of this gating mechanism by studying the free energetics of Cl(-) occupancy in the pore of a prokaryotic ClC homolog. These free energetics were gleaned from 30 ns of molecular dynamics simulation on an approximately 133,000-atom system consisting of a hydrated membrane embedded StClC transporter. The binding sites for Cl(-) in the transporter were determined for the cases where the putative gating residue, Glu(148), was protonated and unprotonated. When the glutamate gate is protonated, Cl(-) favorably occupies an exterior site, S(ext), to form a queue of anions in the pore. However, when the glutamate gate is unprotonated, Cl(-) cannot occupy this site nor, consequently, pass through the pore. An additional, previously undetected, site was found in the pore near the outer membrane that exists regardless of the protonation state of Glu(148). Although this suggests that, for the prokaryotic homolog, protonation of Glu(148) may be the first step in transporting Cl(-) at the expense of H(+) transport in the opposite direction, an evolutionary argument might suggest that Cl(-) opens the ClC gate in eukaryotic channels by inducing the conserved glutamate's protonation. During an additional 20 ns free dynamics simulation, the newly discovered outermost site, S(out), and the innermost site, S(int), were seen to allow spontaneous exchange of Cl(-) ions with the bulk electrolyte while under depolarization conditions.

Filter flexibility and distortion in a bacterial inward rectifier K+ channel: simulation studies of KirBac1.1

The bacterial channel KirBac1.1 provides a structural homolog of mammalian inward rectifier potassium (Kir) channels. The conformational dynamics of the selectivity filter of Kir channels are of some interest in the context of possible permeation and gating mechanisms for this channel. Molecular dynamics simulations of KirBac have been performed on a 10-ns timescale, i.e., comparable to that of ion permeation. The results of five simulations (total simulation time 50 ns) based on three different initial ion configurations and two different model membranes are reported. These simulation data provide evidence for limited (<0.1 nm) filter flexibility during the concerted motion of ions and water molecules within the filter, such local changes in conformation occurring on an approximately 1-ns timescale. In the absence of K(+) ions, the KirBac selectivity filter undergoes more substantial distortions. These resemble those seen in comparable simulations of other channels (e.g., KcsA and KcsA-based homology models) and are likely to lead to functional closure of the channel. This suggests filter distortions may provide a mechanism of K-channel gating in addition to changes in the hydrophobic gate formed at the intracellular crossing point of the M2 helices. The simulation data also provide evidence for interactions of the "slide" (pre-M1) helix of KirBac with phospholipid headgroups.

Ionic permeation free energy in gramicidin: a semimicroscopic perspective

Ion permeation through the gramicidin channel is studied using a model that circumvents two major difficulties inherent to standard simulational methods. It exploits the timescale separation between electronic and structural contributions to dielectric stabilization, accounting for the influence of electronic polarization by embedding the channel in a dielectric milieu that describes this polarization in a mean sense. The explicit mobile moieties are the ion, multipolar waters, and the carbonyls and amides of the peptide backbone. The model treats the influence of aromatic residues and the membrane dipole potential. A new electrical geometry is introduced that treats long-range electrostatics exactly and avoids problems related to periodic boundary conditions. It permits the translocating ion to make a seamless transition from nearby electrolyte to the channel interior. Other degrees of freedom (more distant bulk electrolyte and nonpolar lipid) are treated as dielectric continua. Reasonable permeation free energy profiles are obtained for potassium, rubidium, and cesium; binding wells are shallow and the central barrier is small. Estimated cationic single-channel conductances are smaller than experiment, but only by factors between 2 (rubidium) and 50 (potassium). When applied to chloride the internal barrier is large, with a corresponding miniscule single-channel conductance. The estimated relative single-channel conductances of gramicidin A, B, and C agree well with experiment.

Computer simulations of membrane proteins

Computer simulations are rapidly becoming a standard tool to study the structure and dynamics of lipids and membrane proteins. Increasing computer capacity allows unbiased simulations of lipid and membrane-active peptides. With the increasing number of high-resolution structures of membrane proteins, which also enables homology modelling of more structures, a wide range of membrane proteins can now be simulated over time spans that capture essential biological processes. Longer time scales are accessible by special computational methods. We review recent progress in simulations of membrane proteins.

Lipid-protein interactions of integral membrane proteins: a comparative simulation study

The interactions between membrane proteins and their lipid bilayer environment play important roles in the stability and function of such proteins. Extended (15-20 ns) molecular dynamics simulations have been used to explore the interactions of two membrane proteins with phosphatidylcholine bilayers. One protein (KcsA) is an alpha-helix bundle and embedded in a palmitoyl oleoyl phosphatidylcholine bilayer; the other (OmpA) is a beta-barrel outer-membrane protein and is in a dimyristoyl phosphatidylcholine bilayer. The simulations enable analysis in detail of a number of aspects of lipid-protein interactions. In particular, the interactions of aromatic amphipathic side chains (i.e., Trp, Tyr) with lipid headgroups, and "snorkeling" interactions of basic side chains (i.e., Lys, Arg) with phosphate groups are explored. Analysis of the number of contacts and of H-bonds reveal fluctuations on an approximately 1- to 5-ns timescale. There are two clear bands of interacting residues on the surface of KcsA, whereas there are three such bands on OmpA. A large number of Arg-phosphate interactions are seen for KcsA; for OmpA, the number of basic-phosphate interactions is smaller and shows more marked fluctuations with respect to time. Both classes of interaction occur in clearly defined interfacial regions of width approximately 1 nm. Analysis of lateral diffusion of lipid molecules reveals that "boundary" lipid molecules diffuse at about half the rate of bulk lipid. Overall, these simulations present a dynamic picture of lipid-protein interactions: there are a number of more specific interactions but even these fluctuate on an approximately 1- to 5-ns timescale.

OmpT: molecular dynamics simulations of an outer membrane enzyme

Five molecular dynamics simulations (total duration >25 ns) have been performed on the Escherichia coli outer membrane protease OmpT embedded in a dimyristoylphosphatidylcholine lipid bilayer. Globally the protein is conformationally stable. Some degree of tilt of the beta-barrel is observed relative to the bilayer plane. The greatest degree of conformational flexibility is seen in the extracellular loops. A complex network of fluctuating H-bonds is formed between the active site residues, such that the Asp210-His212 interaction is maintained throughout, whereas His212 and Asp83 are often bridged by a water molecule. This supports a catalytic mechanism whereby Asp83 and His212 bind a water molecule that attacks the peptide carbonyl. A configuration yielded by docking calculations of OmpT simulation snapshots and a model substrate peptide Ala-Arg-Arg-Ala was used as the starting point for an extended Huckel calculation on the docked peptide. These placed the lowest unoccupied molecular orbital mainly on the carbon atom of the central C=O in the scissile peptide bond, thus favoring attack on the central peptide by the water held by residues Asp83 and His212. The trajectories of water molecules reveal exchange of waters between the intracellular face of the membrane and the interior of the barrel but no exchange at the extracellular mouth. This suggests that the pore-like region in the center of OmpT may enable access of water to the active site from below. The simulations appear to reveal the presence of specific lipid interaction sites on the surface of the OmpT barrel. This reveals the ability of extended MD simulations to provide meaningful information on protein-lipid interactions.

Atomistic simulations of biologically realistic transmembrane potential gradients

We present all-atom molecular dynamics simulations of biologically realistic transmembrane potential gradients across a DMPC bilayer. These simulations are the first to model this gradient in all-atom detail, with the field generated solely by explicit ion dynamics. Unlike traditional bilayer simulations that have one bilayer per unit cell, we simulate a 170 mV potential gradient by using a unit cell consisting of three salt-water baths separated by two bilayers, with full three-dimensional periodicity. The study shows that current computational resources are powerful enough to generate a truly electrified interface, as we show the predicted effect of the field on the overall charge distribution. Additionally, starting from Poisson's equation, we show a new derivation of the double integral equation for calculating the potential profile in systems with this type of periodicity.

Kv channel S6 helix as a molecular switch: simulation studies

Ion channels form pores of nanoscopic dimensions in biological membranes and play a key role in the physiology of cells. The majority of ion channels are gated, i.e. they contain a molecular switch that allows a transition between a closed (functionally 'off') and open (functionally 'on') state. Comparison of crystal structures of potassium channels suggest that the gating mechanism of voltage-gated potassium (Kv) channels involves a key role for the pore-lining S6 helix. There is a conserved PVP sequence motif in the S6 helix. Molecular dynamics simulations are used here to explore the conformational dynamics of the S6 helix hinge in models of fragments of a Kv channel, namely an S5-P-S6 monomer and an (S5-P-S6)4 tetramer. The latter is a model of the complete pore-forming domain of a Kv channel. All models were simulated embedded in an octane slab (a simple membrane mimetic). The results of these simulations indicate that the PVP motif may form a molecular hinge, even when the S6 helix forms part of a more complex model. The conformational dynamics of S6 are modulated by the remainder of protein, but it remains flexible. These simulation results are compatible with a channel gating model in which S6 bends in the vicinity of the PVP motif in addition to the region around the conserved glycine (G466) that is N-terminal to the PVP motif. This model is supported by comparison of the Kv S6 models with the S6 helix of the bacterial KvAP channel crystal structure. Thus, K channel gating may depend on a complex nanoswitch with three rigid helical sections linked by two molecular hinges.

KcsA closed and open: modelling and simulation studies

Bacterial homologues of mammalian potassium channels provide structures of two states of a gated K channel. Thus, the crystal structure of KcsA represents a closed state whilst that of MthK represents an open state. Using homology modelling and molecular dynamics simulations we have built a model of the transmembrane domain of KcsA in an open state and have compared its conformational stability with that of the same domain of KcsA in a closed state. Approximate Born energy calculations of monovalent cations within the two KcsA channel states suggest that the intracellular hydrophobic gate in the closed state provides a barrier of height approximately 5 kT to ion permeation, whilst in the open state the barrier is absent. Simulations (10 ns duration) in an octane slab (a simple membrane mimetic) suggest that closed- and open-state models are of comparable conformational stability, both exhibiting conformational drifts of approximately 3.3 A Calpha RMSD relative to the respective starting models. Substantial conformational fluctuations are observed in the intracellular gate region during both simulations (closed state and open state). In the simulation of open-state KcsA, rapid ( < 5 ns) exit of all three K+ ions occurs through the intracellular mouth of the channel. Helix kink and swivel motion is observed at the molecular hinge formed by residue G99 of the M2 helix. This motion is more substantial for the open- than for the closed-state model of the channel.