Exploring conformational states and helical packings in the P2X receptor transmembrane domain by molecular dynamics simulation

The P2X receptor is a trimeric transmembrane protein that acts as an ATP-gated ion channel. Its transmembrane domain (TMD) contains only six helices and three of them, the M2 helices, line the ion conduction pathway. Here, using molecular dynamics simulation, I identify four conformational states of the TMD that are associated with four types of packing between M2 helices. Packing in the extracellular half of the M2 helix produces closed conformations, while packing in the intracellular half produces both open and closed conformations. State transition is observed and supports a mechanism where iris-like twisting of the M2 helices switches the location of helical packing between the extracellular and the intracellular halves of the helices. In addition, this twisting motion alters the position and orientation of residue side-chains relative to the pore and therefore influences the pore geometry and possibly ion permeation. Helical packing, on the other hand, may restrict the twisting motion and generate discrete conformational states.     

Open-state conformation of the KcsA K+ channel: Monte Carlo normal mode following simulations

Potassium channels fluctuate between closed and open states. The detailed mechanism of the conformational changes opening the intracellular pore in the K+ channel from Streptomyces lividans (KcsA) is unknown. Applying Monte Carlo normal mode following, we find that gating involves rotation and unwinding of the TM2 bundle, lateral movement of the TM2 helices away from the channel axis, and disappearance of the TM2 bundle. The open-state conformation of KcsA exhibits a very wide inner vestibule, with a radius approximately 5-7 A and inner helices bent at the A98-G99 hinge. Computed conformational changes demonstrate that spin labeling and X-ray experiments illuminate different stages in gating: transition begins with clockwise rotation of the TM2 helices ending at a final state with the TM2 bend hinged near residues A98-G99. The concordance between the computational and experimental results provides atomic-level insights into the structural rearrangements of the channel's inner pore.     

Coupling between the voltage-sensing and pore domains in a voltage-gated potassium channel

Voltage-dependent potassium (Kv), sodium (Nav), and calcium channels open and close in response to changes in transmembrane (TM) potential, thus regulating cell excitability by controlling ion flow across the membrane. An outstanding question concerning voltage gating is how voltage-induced conformational changes of the channel voltage-sensing domains (VSDs) are coupled through the S4-S5 interfacial linking helices to the opening and closing of the pore domain (PD). To investigate the coupling between the VSDs and the PD, we generated a closed Kv channel configuration from Aeropyrum pernix (KvAP) using atomistic simulations with experiment-based restraints on the VSDs. Full closure of the channel required, in addition to the experimentally determined TM displacement, that the VSDs be displaced both inwardly and laterally around the PD. This twisting motion generates a tight hydrophobic interface between the S4-S5 linkers and the C-terminal ends of the pore domain S6 helices in agreement with available experimental evidence.     

The closed structure of the MscS mechanosensitive channel. Cross-linking of single cysteine mutants

Mechanosensitive channels must make a large conformational change during the transition from the closed to the open state. The crystal structure of the open form of the Escherichia coli MscS channel was recently solved and depicts a homoheptamer (1). In this study, cross-linking of site-specific cysteine substitutions demonstrates that residues up to 10-33 A apart in the crystal structure readily form disulfide bridges in the closed form and can also be cross-linked by a 10-A linker. Cross-linking between adjacent subunits stabilizes the heptameric form of the channel providing biochemical evidence to support the crystal structure. The data are consistent with the published model (1) in that the membrane domain is highly flexible and that the closed to open transition may involve a significant displacement of transmembrane helices 1 and 2, possibly by as much as 30 A. The data are also consistent with significant flexibility of the cytoplasmic domain.     

Structural models of the MscL gating mechanism.

Three-dimensional structural models of the mechanosensitive channel of large conductance, MscL, from the bacteria Mycobacterium tuberculosis and Escherichia coli were developed for closed, intermediate, and open conformations. The modeling began with the crystal structure of M. tuberculosis MscL, a homopentamer with two transmembrane alpha-helices, M1 and M2, per subunit. The first 12 N-terminal residues, not resolved in the crystal structure, were modeled as an amphipathic alpha-helix, called S1. A bundle of five parallel S1 helices are postulated to form a cytoplasmic gate. As membrane tension induces expansion, the tilts of M1 and M2 are postulated to increase as they move away from the axis of the pore. Substantial expansion is postulated to occur before the increased stress in the S1 to M1 linkers pulls the S1 bundle apart. During the opening transition, the S1 helices and C-terminus amphipathic alpha-helices, S3, are postulated to dock parallel to the membrane surface on the perimeter of the complex. The proposed gating mechanism reveals critical spatial relationships between the expandable transmembrane barrel formed by M1 and M2, the gate formed by S1 helices, and "strings" that link S1s to M1s. These models are consistent with numerous experimental results and modeling criteria.     

Ion channel gating: insights via molecular simulations

Ion channels are gated, i.e. they can switch conformation between a closed and an open state. Molecular dynamics simulations may be used to study the conformational dynamics of ion channels and of simple channel models. Simulations on model nanopores reveal that a narrow (<4 A) hydrophobic region can form a functionally closed gate in the channel and can be opened by either a small (approximately 1 A) increase in pore radius or an increase in polarity. Modelling and simulation studies confirm the importance of hydrophobic gating in K channels, and support a model in which hinge-bending of the pore-lining M2 (or S6 in Kv channels) helices underlies channel gating. Simulations of a simple outer membrane protein, OmpA, indicate that a gate may also be formed by interactions of charged side chains within a pore, as is also the case in ClC channels.     

Asymmetric configurations and N-terminal rearrangements in connexin26 gap junction channels

Gap junction channels are unique in that they possess multiple mechanisms for channel closure, several of which involve the N terminus as a key component in gating, and possibly assembly. Here, we present electron crystallographic structures of a mutant human connexin26 (Cx26M34A) and an N-terminal deletion of this mutant (Cx26M34Adel2-7) at 6-Å and 10-Å resolutions, respectively. The three-dimensional map of Cx26M34A was improved by data from 60° tilt images and revealed a breakdown of the hexagonal symmetry in a connexin hemichannel, particularly in the cytoplasmic domain regions at the ends of the transmembrane helices. The Cx26M34A structure contained an asymmetric density in the channel vestibule ("plug") that was decreased in the Cx26M34Adel2-7 structure, indicating that the N terminus significantly contributes to form this plug feature. Functional analysis of the Cx26M34A channels revealed that these channels are predominantly closed, with the residual electrical conductance showing normal voltage gating. N-terminal deletion mutants with and without the M34A mutation showed no electrical activity in paired Xenopus oocytes and significantly decreased dye permeability in HeLa cells. Comparing this closed structure with the recently published X-ray structure of wild-type Cx26, which is proposed to be in an open state, revealed a radial outward shift in the transmembrane helices in the closed state, presumably to accommodate the N-terminal plug occluding the pore. Because both Cx26del2-7 and Cx26M34Adel2-7 channels are closed, the N terminus appears to have a prominent role in stabilizing the open configuration.     

Analysis of Hyperekplexia Mutations Identifies Transmembrane Domain Rearrangements That Mediate Glycine Receptor Activation

Pentameric ligand-gated ion channels (pLGICs) mediate numerous physiological processes and are therapeutic targets for a wide range of clinical indications. Elucidating the structural differences between their closed and open states may help in designing improved drugs that bias receptors toward the desired conformational state. We recently showed that two new hyperekplexia mutations, Q226E and V280M, induced spontaneous activity in α1 glycine receptors. Gln-226, located near the top of transmembrane (TM) 1, is closely apposed to Arg-271 at the top of TM2 in the neighboring subunit. Using mutant cycle analysis, we inferred that Q226E induces activation via an enhanced electrostatic attraction to Arg-271. This would tilt the top of TM2 toward TM1 and hence away from the pore axis to open the channel. We also concluded that the increased side chain volume of V280M, in the TM2-TM3 loop, exerts a steric repulsion against Ile-225 at the top of TM1 in the neighboring subunit. We infer that this steric repulsion would tilt the top of TM3 radially outwards against the stationary TM1 and thus provide space for TM2 to relax away from the pore axis to create an open channel. Because the transmembrane domain movements inferred from this functional analysis are consistent with the structural differences evident in the x-ray atomic structures of closed and open state bacterial pLGICs, we propyose that the model of pLGIC activation as outlined here may be broadly applicable across the eukaryotic pLGIC receptor family.     

Two different conformational states of the KirBac3.1 potassium channel revealed by electron crystallography

Potassium channels allow the selective flow of K(+) ions across membranes. In response to external gating signals, the potassium channel can move reversibly through a series of structural conformations from a closed to an open state. 2D crystals of the inwardly rectifying K(+) channel KirBac3.1 from Magnetospirillum magnetotacticum have been captured in two distinct conformations, providing "snap shots" of the gating process. Analysis by electron cryomicroscopy of these KirBac3.1 crystals has resulted in reconstructed images in projection at 9 A resolution. Kir channels are tetramers of four subunits arranged as dimers of dimers. Each subunit has two transmembrane helices (inner and outer). In one crystal form, the pore is blocked; in the other crystal form, the pore appears open. Modeling based on the KirBac1.1 (closed) crystal structure shows that opening of the ion conduction pathway could be achieved by bending of the inner helices and significant movements of the outer helices.     

Free-energy landscapes of ion-channel gating are malleable: changes in the number of bound ligands are accompanied by changes in the location of the transition state in acetylcholine-receptor channels

Acetylcholine-receptor channels (AChRs) are allosteric membrane proteins that mediate synaptic transmission by alternatively opening and closing ("gating") a cation-selective transmembrane pore. Although ligand binding is not required for the channel to open, the binding of agonists (for example, acetylcholine) increases the closed right harpoon over left harpoon open equilibrium constant because the ion-impermeable --> ion-permeable transition of the ion pathway is accompanied by a low-affinity --> high-affinity change at the agonist-binding sites. The fact that the gating conformational change of muscle AChRs can be kinetically modeled as a two-state reaction has paved the way to the experimental characterization of the corresponding transition state, which represents a snapshot of the continuous sequence of molecular events separating the closed and open states. Previous studies of fully (di) liganded AChRs, combining single-channel kinetic measurements, site-directed mutagenesis, and data analysis in the framework of the linear free-energy relationships of physical organic chemistry, have suggested a transition-state structure that is consistent with channel opening being an asynchronous conformational change that starts at the extracellular agonist-binding sites and propagates toward the intracellular end of the pore. In this paper, I characterize the gating transition state of unliganded AChRs, and report a remarkable difference: unlike that of diliganded gating, the unliganded transition state is not a hybrid of the closed- and open-state structures but, rather, is almost indistinguishable from the open state itself. This displacement of the transition state along the reaction coordinate obscures the mechanism underlying the unliganded closed right harpoon over left harpoon open reaction but brings to light the malleable nature of free-energy landscapes of ion-channel gating.     

Aromatic–aromatic interactions between residues in KCa3.1 pore helix and S5 transmembrane segment control the channel gating process

The Ca²⁺-activated potassium channel KCa3.1 is emerging as a therapeutic target for a large variety of health disorders. One distinguishing feature of KCa3.1 is that the channel open probability at saturating Ca²⁺ concentrations (Pomax) is low, typically 0.1–0.2 for KCa3.1 wild type. This observation argues for the binding of Ca²⁺ to the calmodulin (CaM)–KCa3.1 complex, promoting the formation of a preopen closed-state configuration leading to channel opening. We have previously shown that the KCa3.1 active gate is most likely located at the level of the selectivity filter. As Ca²⁺-dependent gating of KCa3.1 originates from the binding of Ca²⁺ to CaM in the C terminus, the hypothesis of a gate located at the level of the selectivity filter requires that the conformational change initiated in the C terminus be transmitted to the S5 and S6 transmembrane helices, with a resulting effect on the channel pore helix directly connected to the selectivity filter. A study was thus undertaken to determine to what extent the interactions between the channel pore helix with the S5 and S6 transmembrane segments contribute to KCa3.1 gating. Molecular dynamics simulations first revealed that the largest contact area between the pore helix and the S5 plus S6 transmembrane helices involves residue F248 at the C-terminal end of the pore helix. Unitary current recordings next confirmed that modulating aromatic–aromatic interactions between F248 and W216 of the S5 transmembrane helical segment and/or perturbing the interactions between F248 and residues in S6 surrounding the glycine hinge G274 cause important changes in Pomax. This work thus provides the first evidence for a key contribution of the pore helix in setting Pomax by stabilizing the channel closed configuration through aromatic–aromatic interactions involving F248 of the pore helix. We propose that the interface pore helix/S5 constitutes a promising site for designing KCa3.1 potentiators.     

Toward elucidating the heat activation mechanism of the TRPV1 channel gating by molecular dynamics simulation

As a key cellular sensor, the TRPV1 cation channel undergoes a gating transition from a closed state to an open state in response to various physical and chemical stimuli including noxious heat. Despite years of study, the heat activation mechanism of TRPV1 gating remains enigmatic at the molecular level. Toward elucidating the structural and energetic basis of TRPV1 gating, we have performed molecular dynamics (MD) simulations (with cumulative simulation time of 3 μs), starting from the high‐resolution closed and open structures of TRPV1 solved by cryo‐electron microscopy. In the closed‐state simulations at 30°C, we observed a stably closed channel constricted at the lower gate (near residue I679), while the upper gate (near residues G643 and M644) is dynamic and undergoes flickery opening/closing. In the open‐state simulations at 60°C, we found higher conformational variation consistent with a large entropy increase required for the heat activation, and both the lower and upper gates are dynamic with transient opening/closing. Through ensemble‐based structural analyses of the closed state versus the open state, we revealed pronounced closed‐to‐open conformational changes involving the membrane proximal domain (MPD) linker, the outer pore, and the TRP helix, which are accompanied by breaking/forming of a network of closed/open‐state specific hydrogen bonds. By comparing the closed‐state simulations at 30°C and 60°C, we observed heat‐activated conformational changes in the MPD linker, the outer pore, and the TRP helix that resemble the closed‐to‐open conformational changes, along with partial formation of the open‐state specific hydrogen bonds. Some of the residues involved in the above key hydrogen bonds were validated by previous mutational studies. Taken together, our MD simulations have offered rich structural and dynamic details beyond the static structures of TRPV1, and promising targets for future mutagenesis and functional studies of the TRPV1 channel. Proteins 2016; 84:1938–1949. © 2016 Wiley Periodicals, Inc.     

Conformational changes of the histidine ATP-binding cassette transporter studied by double electron–electron resonance spectroscopy

The conformational dynamics of the histidine ABC transporter HisQMP₂ from Salmonella enterica serovar Typhimurium, reconstituted into liposomes, is studied by site-directed spin labeling and double electron–electron resonance spectroscopy in the absence of nucleotides, in the ATP-bound, and in the post-hydrolysis state. The results show that the inter-dimer distances as measured between the Q-loops of HisP₂ in the intact transporter resemble those determined for the maltose transporter in all three states of the hydrolysis cycle. Only in the presence of liganded HisJ the closed conformation of the nucleotide binding sites is achieved revealing the transmembrane communication of the presence of substrate. Two conformational states can be distinguished for the periplasmic moiety of HisQMP₂ as detected by differences in distributions of interspin distances between positions 86 and 96 or 104 and 197. The observed conformational changes are correlated to proposed open, semi-open and closed conformations of the nucleotide binding domains HisP₂. Our results are in line with a rearrangement of transmembrane helices 4 and 4′ of HisQM during the closed to the semi-open transition of HisP₂ driven by the reorientation of the coupled helices 3a and 3b to occur upon hydrolysis.     

Conformational dynamics of M2 helices in KirBac channels: helix flexibility in relation to gating via molecular dynamics simulations

KirBac1.1 and 3.1 are bacterial homologues of mammalian inward rectifier K channels. We have performed extended molecular dynamics simulations (five simulations, each of >20 ns duration) of the transmembrane domain of KirBac in two membrane environments, a palmitoyl oleoyl phosphatidylcholine bilayer and an octane slab. Analysis of these simulations has focused on the conformational dynamics of the pore-lining M2 helices, which form the cytoplasmic hydrophobic gate of the channel. Principal components analysis reveals bending of M2, with a molecular hinge at the conserved glycine (Gly134 in KirBac1.1, Gly120 in KirBac3.1). More detailed analysis reveals a dimer-of-dimers type motion. The first two eigenvectors describing the motions of M2 correspond to helix kink and swivel motions. The conformational flexibility of M2 seen in these simulations correlates with differences in M2 conformation between that seen in the X-ray structures of closed channels (KcsA and KirBac) in which the helix is undistorted, and in open channels (e.g. MthK) in which the M2 helix is kinked. Thus, the simulations, albeit on a time scale substantially shorter than that required for channel gating, suggest a gating model in which the intrinsic flexibility of M2 about a molecular hinge is coupled to conformational transitions of an intracellular 'gatekeeper' domain, the latter changing conformation in response to ligand binding.     

Gating movement of acetylcholine receptor caught by plunge-freezing

The nicotinic acetylcholine (ACh) receptor converts transiently to an open-channel form when activated by ACh released into the synaptic cleft. We describe here the conformational change underlying this event, determined by electron microscopy of ACh-sprayed and freeze-trapped postsynaptic membranes. ACh binding to the α subunits triggers a concerted rearrangement in the ligand-binding domain, involving an ~ 1‐Å outward displacement of the extracellular portion of the β subunit where it interacts with the juxtaposed ends of α-helices shaping the narrow membrane-spanning pore. The β-subunit helices tilt outward to accommodate this displacement, destabilising the arrangement of pore-lining helices, which in the closed channel bend inward symmetrically to form a central hydrophobic gate. Straightening and tangential motion of the pore-lining helices effect channel opening by widening the pore asymmetrically and increasing its polarity in the region of the gate. The pore-lining helices of the α_γ and δ subunits, by flexing between alternative bent and straight conformations, undergo the greatest movements. This coupled allosteric transition shifts the structure from a tense (closed) state toward a more relaxed (open) state.     

A structural model of the acetylcholine receptor channel based on partition energy and helix packing calculations.

A structural model of the transmembrane portion of the acetylcholine receptor was developed from sequences of all its subunits by using transfer energy calculations to locate transmembrane alpha-helices and to calculate which helical side chains should be in contact with water inside the channel, with portions of other transmembrane helices, or with lipid hydrocarbon chains. "Knobs-into-holes" side chain packing calculations were used with other factors to stack the transmembrane alpha-helices together. In the model each subunit has the following structures in order along the sequence from the NH2 terminus: a large extracellular domain of undetermined structure, a short apolar alpha-helix that lies on the extracellular lipid surface of the membrane; three apolar transmembrane alpha-helices (I, II, and III), a cytoplasmic domain of undetermined structure, an amphipathic transmembrane alpha-helix (L) that forms the channel lining, a short extracellular alpha-helix, another apolar transmembrane alpha-helix (IV), and a small cytoplasmic domain formed by the COOH-terminal end of the chain. Three concentric layers form the pore. A bundle of five amphipathic L helices forms the channel lining. This bundle is surrounded by a bundle of 10 alternating II and III helices. Helices I and IV cover portions of the outer surface of the bundle formed by helices II and III. Positions of disulfide bridges are predicted and a mechanism for opening and closing conformational changes is proposed that requires tilting transmembrane helices and possibly a thiol-disulfide interchange reaction.     

Functional architecture of the CFTR chloride channel

Cystic fibrosis is caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR), a member of the ATP-binding cassette (ABC) family of membrane transport proteins. CFTR is unique among ABC proteins in that it functions not as an active transporter but as an ATP-gated Cl^? channel. As an ion channel, the function of the CFTR transmembrane channel pore that mediates Cl^? movement has been studied in great detail. On the other hand, only low resolution structural data is available on the transmembrane parts of the protein. The structure of the channel pore has, however, been modeled on the known structure of active transporter ABC proteins. Currently, significant barriers exist to building a unified view of CFTR pore structure and function. Reconciling functional data on the channel with indirect structural data based on other proteins with very different transport functions and substrates has proven problematic. This review summarizes current structural and functional models of the CFTR Cl^? channel pore, including a comprehensive review of previous electrophysiological investigations of channel structure and function. In addition, functional data on the three-dimensional arrangement of pore-lining helices, as well as contemporary hypotheses concerning conformational changes in the pore that occur during channel opening and closing, are discussed. Important similarities and differences between different models of the pore highlight current gaps in our knowledge of CFTR structure and function. In order to fill these gaps, structural and functional models of the membrane-spanning pore need to become better integrated.     

Local constraints in either the GluN1 or GluN2 subunit equally impair NMDA receptor pore opening

The defining functional feature of N-methyl-d-aspartate (NMDA) receptors is activation gating, the energetic coupling of ligand binding into opening of the associated ion channel pore. NMDA receptors are obligate heterotetramers typically composed of glycine-binding GluN1 and glutamate-binding GluN2 subunits that gate in a concerted fashion, requiring all four ligands to bind for subsequent opening of the channel pore. In an individual subunit, the extracellular ligand-binding domain, composed of discontinuous polypeptide segments S1 and S2, and the transmembrane channel-forming domain, composed of M1-M4 segments, are connected by three linkers: S1-M1, M3-S2, and S2-M4. To study subunit-specific events during pore opening in NMDA receptors, we impaired activation gating via intrasubunit disulfide bonds connecting the M3-S2 and S2-M4 in either the GluN1 or GluN2A subunit, thereby interfering with the movement of the M3 segment, the major pore-lining and channel-gating element. NMDA receptors with gating impairments in either the GluN1 or GluN2A subunit were dramatically resistant to channel opening, but when they did open, they showed only a single-conductance level indistinguishable from wild type. Importantly, the late gating steps comprising pore opening to its main long-duration open state were equivalently affected regardless of which subunit was constrained. Thus, the NMDA receptor ion channel undergoes a pore-opening mechanism in which the intrasubunit conformational dynamics at the level of the ligand-binding/transmembrane domain (TMD) linkers are tightly coupled across the four subunits. Our results further indicate that conformational freedom of the linkers between the ligand-binding and TMDs is critical to the activation gating process.     

New insights of membrane environment effects on MscL channel mechanics from theoretical approaches

The prokaryotic mechanosensitive channel of large conductance (MscL) is a remarkable integral membrane protein. During hypo-osmotic shock, it responses to membrane tension through large conformational changes, that lead to an open state of the pore. The structure of the channel from Mycobacterium tuberculosis has been resolved in the closed state. Numerous experiments have attempted to trap the channel in its open state but they did not succeed in obtaining a structure. A gating mechanism has been proposed based on different experimental data but there is no experimental technique available to follow this process in atomic details. In addition, it has been shown that a decrease of the lipid bilayer thickness lowered MscL activation energy and stabilized a structurally distinct closed channel intermediate. Here, we use atomistic molecular dynamics simulations to investigate the effect of the lipid bilayer thinning on our model of the structure of the Escherichia coli. We thoroughly analyze simulations of the channel embedded in two pre-equilibrated membranes differing by their hydrophobic tail length (DMPE and POPE). The MscL structure remains stable in POPE, whereas a distinct structural state is obtained in DMPE in response to hydrophobic mismatch. This latter is obtained by tilts and kinks of the transmembrane helices, leading to a widening and a diminution of the channel height. Part of these motions is guided by a competition between solvent and lipids for the interaction with the periplasmic loops. We finally conduct a principal component analysis of the simulation and compare anharmonic motions with harmonic ones, previously obtained from a coarse-grained normal mode analysis performed on the same structural model. Significant similarities exist between low-frequency harmonic motions and those observed with essential dynamics in DMPE. In summary, change in membrane thickness permits to accelerate the conformational changes involved in the mechanics of the E. coli channel, providing a closed structural intermediate en route to the open state. These results give clues for better understanding why the channel activation energy is lowered in a thinner membrane.     

Computational modelling of the open-state Kv1.5 ion channel block by bupivacaine

Binding of R(+)-bupivacaine to open-state homology models of the mammalian K_v1.5 membrane ion channel is studied using automated docking and molecular dynamics (MD) methods. Homology models of K_v1.5 are built using the 3D structures of the KcsA and MthK channels as a template. The packing of transmembrane (TM) α-helices in the KcsA structure corresponds to a closed channel state. Opening of the channel may be reached by a conformational transition yielding a bent structure of the internal S6 helices. Our first model of the K_v open state involves a PVP-type of bending hinge in the internal helices, while the second model corresponds to a Gly-type of bending hinge as found in the MthK channel. Ligand binding to these models is probed using the common local anaesthetic bupivacaine, where blocker binding from the intracellular side of the channel is considered. Conformational properties and partial atomic charges of bupivacaine are determined from quantum mechanical HF/6-31G* calculations with inclusion of solvent effects. The automated docking and MD calculations for the PVP-bend model predict that bupivacaine could bind either in the central cavity or in the PVP region of the channel pore. Linear interaction energy (LIE) estimates of the binding free energies for bupivacaine predict strongest binding to the PVP region. Surprisingly, no binding is predicted for the Gly-bend model. These results are discussed in light of electrophysiological data which show that the K_v1.5 channel is unable to close when bupivacaine is bound.