Investigating the size and dynamics of voltage-gated sodium channel fenestrations

Eukaryotic voltage-gated sodium channels (VGSCs) are essential for the initiation and propagation of action potentials in electrically excitable cells, and are important pharmaceutical targets for the treatment of neurological disorders such as epilepsy, cardiac arrhythmias, and chronic pain. Evidence suggests that small, hydrophobic, VGSC-blocking drugs can gain access to binding residues within the central cavity of these channels by passing through lateral, lipid-filled “fenestrations” which run between the exterior of the protein and its central pore. Here, we use molecular dynamics simulations to investigate how the size and shape of fenestrations change over time in several bacterial VGSC models and a homology model of Nav1.4. We show that over the course of the simulations, the size of the fenestrations is primarily influenced by rapid protein motions, such as amino acid side-chain rotation, and highlight that differences between fenestration bottleneck-contributing residues are the primary cause of variations in fenestration size between the 6 bacterial models. In the eukaryotic channel model, 2 fenestrations are wide, but 2 are narrow due to differences in the amino acid sequence in the 4 domains. Lipid molecules are found to influence the size of the fenestrations by protruding acyl chains into the fenestrations and displacing amino acid side-chains. Together, the results suggest that fenestrations provide viable pathways for small, flexible, hydrophobic drugs.     

Characterizing fenestration size in sodium channel subtypes and their accessibility to inhibitors

Voltage-gated sodium channels (Nav) underlie the electrical activity of nerve and muscle cells. Humans have nine different subtypes of these channels, which are the target of small-molecule inhibitors commonly used to treat a range of conditions. Structural studies have identified four lateral fenestrations within the Nav pore module that have been shown to influence Nav pore blocker access during resting-state inhibition. However, the structural differences among the nine subtypes are still unclear. In particular, the dimensions of the four individual fenestrations across the Nav subtypes and their differential accessibility to pore blockers is yet to be characterized. To address this, we applied classical molecular dynamics simulations to study the recently published structures of Nav1.1, Nav1.2, Nav1.4, Nav1.5, and Nav1.7. Although there is significant variability in the bottleneck sizes of the Nav fenestrations, the subtypes follow a common pattern, with wider DI-II and DIII-IV fenestrations, a more restricted DII-III fenestration, and the most restricted DI-IV fenestration. We further identify the key bottleneck residues in each fenestration and show that the motions of aromatic residue sidechains govern the bottleneck radii. Well-tempered metadynamics simulations of Nav1.4 and Nav1.5 in the presence of the pore blocker lidocaine also support the DI-II fenestration being the most likely access route for drugs. Our computational results provide a foundation for future in vitro experiments examining the route of drug access to sodium channels. Understanding the fenestrations and their accessibility to drugs is critical for future analyses of diseases mutations across different sodium channel subtypes, with the potential to inform pharmacological development of resting-state inhibitors and subtype-selective drug design.     

Influence of lipids on the hydrophobic barrier within the pore of the TWIK-1 K2P channel

Several recent ion channel structures have revealed large side portals, or ‘fenestrations’ at the interface between their transmembrane helices that potentially expose the ion conduction pathway to the lipid core of the bilayer. In a recent study we demonstrated that functional activity of the TWIK-1 K2P channel is influenced by the presence of hydrophobic residues deep within the inner pore. These residues are located near the fenestrations in the TWIK-1 structure and promote dewetting of the pore by forming a hydrophobic barrier to ion conduction. During our previous MD simulations, lipid tails were observed to enter these fenestrations. In this addendum to that study, we investigate lipid contribution to the dewetting process. Our results demonstrate that lipid tails from both the upper and lower leaflets can occupy the fenestrations and partially penetrate into the pore. The lipid tails do not sterically occlude the pore, but there is an inverse correlation between the presence of water within the hydrophobic barrier and the number of lipids tails within the lining of the pore. However, dewetting still occurs in the absence of lipids tails, and pore hydration appears to be determined primarily by those side-chains lining the narrowest part of the pore cavity.     

Influence of lipids on the hydrophobic barrier within the pore of the TWIK-1 K2P channel

Several recent ion channel structures have revealed large side portals, or ‘fenestrations’ at the interface between their transmembrane helices that potentially expose the ion conduction pathway to the lipid core of the bilayer. In a recent study we demonstrated that functional activity of the TWIK-1 K2P channel is influenced by the presence of hydrophobic residues deep within the inner pore. These residues are located near the fenestrations in the TWIK-1 structure and promote dewetting of the pore by forming a hydrophobic barrier to ion conduction. During our previous MD simulations, lipid tails were observed to enter these fenestrations. In this addendum to that study, we investigate lipid contribution to the dewetting process. Our results demonstrate that lipid tails from both the upper and lower leaflets can occupy the fenestrations and partially penetrate into the pore. The lipid tails do not sterically occlude the pore, but there is an inverse correlation between the presence of water within the hydrophobic barrier and the number of lipids tails within the lining of the pore. However, dewetting still occurs in the absence of lipids tails, and pore hydration appears to be determined primarily by those side-chains lining the narrowest part of the pore cavity.     

Structure of Membrane-active Toxin from Crab Spider Heriaeus melloteei Suggests Parallel Evolution of Sodium Channel Gating Modifiers in Araneomorphae and Mygalomorphae

We present a structural and functional study of a sodium channel activation inhibitor from crab spider venom. Hm-3 is an insecticidal peptide toxin consisting of 35 amino acid residues from the spider Heriaeus melloteei (Thomisidae). We produced Hm-3 recombinantly in Escherichia coli and determined its structure by NMR spectroscopy. Typical for spider toxins, Hm-3 was found to adopt the so-called “inhibitor cystine knot” or “knottin” fold stabilized by three disulfide bonds. Its molecule is amphiphilic with a hydrophobic ridge on the surface enriched in aromatic residues and surrounded by positive charges. Correspondingly, Hm-3 binds to both neutral and negatively charged lipid vesicles. Electrophysiological studies showed that at a concentration of 1 μm Hm-3 effectively inhibited a number of mammalian and insect sodium channels. Importantly, Hm-3 shifted the dependence of channel activation to more positive voltages. Moreover, the inhibition was voltage-dependent, and strong depolarizing prepulses attenuated Hm-3 activity. The toxin is therefore concluded to represent the first sodium channel gating modifier from an araneomorph spider and features a “membrane access” mechanism of action. Its amino acid sequence and position of the hydrophobic cluster are notably different from other known gating modifiers from spider venom, all of which are described from mygalomorph species. We hypothesize parallel evolution of inhibitor cystine knot toxins from Araneomorphae and Mygalomorphae suborders.     

Molecular Dynamic Simulations Reveal the Structural Determinants of Fatty Acid Binding to Oxy-Myoglobin

The mechanism(s) by which fatty acids are sequestered and transported in muscle have not been fully elucidated. A potential key player in this process is the protein myoglobin (Mb). Indeed, there is a catalogue of empirical evidence supporting direct interaction of globins with fatty acid metabolites; however, the binding pocket and regulation of the interaction remains to be established. In this study, we employed a computational strategy to elucidate the structural determinants of fatty acids (palmitic & oleic acid) binding to Mb. Sequence analysis and docking simulations with a horse (Equus caballus) structural Mb reference reveals a fatty acid-binding site in the hydrophobic cleft near the heme region in Mb. Both palmitic acid and oleic acid attain a “U” shaped structure similar to their conformation in pockets of other fatty acid-binding proteins. Specifically, we found that the carboxyl head group of palmitic acid coordinates with the amino group of Lys45, whereas the carboxyl group of oleic acid coordinates with both the amino groups of Lys45 and Lys63. The alkyl tails of both fatty acids are supported by surrounding hydrophobic residues Leu29, Leu32, Phe33, Phe43, Phe46, Val67, Val68 and Ile107. In the saturated palmitic acid, the hydrophobic tail moves freely and occasionally penetrates deeper inside the hydrophobic cleft, making additional contacts with Val28, Leu69, Leu72 and Ile111. Our simulations reveal a dynamic and stable binding pocket in which the oxygen molecule and heme group in Mb are required for additional hydrophobic interactions. Taken together, these findings support a mechanism in which Mb acts as a muscle transporter for fatty acid when it is in the oxygenated state and releases fatty acid when Mb converts to deoxygenated state.     

Water transport in human aquaporin-4: molecular dynamics (MD) simulations

Aquaporin-4 (AQP4) is the predominant water channel in the central nervous system, where it has been reported to be involved in many pathophysiological roles including water transport. In this paper, the AQP4 tetramer was modeled from its PDB structure file, embedded in a palmitoyl-oleoyl-phosphatidyl-choline (POPC) lipid bilayer, solvated in water, then minimized and equilibrated by means of molecular dynamics simulations. Analysis of the equilibrated structure showed that the central pore along the fourfold axis of the tetramers is formed with hydrophobic amino acid residues. In particular, Phe-195, Leu-191 and Leu-75, form the narrowest part of the pore. Therefore water molecules are not expected to transport through the central pore, which was confirmed by MD simulations. Each monomer of the AQP4 tetramers forms a channel whose walls consist mostly of hydrophilic residues. There are eight water molecules in single file observed in each of the four channels, transporting through the selectivity filter containing Arg-216, His-201, Phe-77, Ala-210, and the two conserved Asn-Pro-Ala (NPA) motifs containing Asn-213 and Asn-97. By using Brownian dynamics fluctuation–dissipation-theorem (BD-FDT), the overall free-energy profile was obtained for water transporting through AQP4 for the first time, which gives a complete map of the entire channel of water permeation.     

Block of Brain Sodium Channels by Peptide Mimetics of the Isoleucine,Phenylalanine, and Methionine (IFM) Motif from the Inactivation Gate

Inactivation of sodium channels is thought to be mediated by an inactivation gate formed by the intracellular loop connecting domains III and IV. A hydrophobic motif containing the amino acid sequence isoleucine, phenylalanine, and methionine (IFM) is required for the inactivation process. Peptides containing the IFM motif, when applied to the cytoplasmic side of these channels, produce two types of block: fast block, which resembles the inactivation process, and slow, use-dependent block stimulated by strong depolarizing pulses. Fast block by the peptide ac-KIFMK-NH₂, measured on sodium channels whose inactivation was slowed by the α-scorpion toxin from Leiurus quinquestriatus (LqTx), was reversed with a time constant of 0.9 ms upon repolarization. In contrast, control and LqTx-modified sodium channels were slower to recover from use-dependent block. For fast block, linear peptides of three to six amino acid residues containing the IFM motif and two positive charges were more effective than peptides with one positive charge, whereas uncharged IFM peptides were ineffective. Substitution of the IFM residues in the peptide ac-KIFMK-NH₂ with smaller, less hydrophobic residues prevented fast block. The positively charged tripeptide IFM-NH₂ did not cause appreciable fast block, but the divalent cation IFM-NH(CH₂)₂NH₂ was as effective as the pentapeptide ac-KIFMK-NH₂. The constrained peptide cyclic KIFMK containing two positive charges did not cause fast block. These results indicate that the position of the positive charges is unimportant, but flexibility or conformation of the IFM-containing peptide is important to allow fast block. Slow, use-dependent block was observed with IFM-containing peptides of three to six residues having one or two positive charges, but not with dipeptides or phenylalanine-amide. In contrast to its lack of fast block, cyclic KIFMK was an effective use-dependent blocker. Substitutions of amino acid residues in the tripeptide IFM-NH₂ showed that large hydrophobic residues are preferred in all three positions for slow, use-dependent block. However, substitution of the large hydrophobic residue diphenylalanine or the constrained residues phenylglycine or tetrahydroisoquinoline for phe decreased potency, suggesting that this phe residue must be able to enter a restricted hydrophobic pocket during the binding of IFM peptides. Together, the results on fast block and slow, use-dependent block indicate that IFM peptides form two distinct complexes of different stability and structural specificity with receptor site(s) on the sodium channel. It is proposed that fast block represents binding of these peptides to the inactivation gate receptor, while slow, use-dependent block represents deeper binding of the IFM peptides in the pore.     

Homology Models of the Tetramerization Domain of Six Eukaryotic Voltage-gated Potassium Channels Kv1.1-Kv1.6

Abstract

The homology models of the tetramerization (T1) domain of six eukaryotic potassium channels, Kv1.1-Kv1.6, were constructed based on the crystal structure of the Shaker T1 domain. The results of amino acid sequence alignment indicate that the T1 domains of these K⁺ channels are highly conserved, with the similarities varying from 77% between Shaker and Kv1.6 to 93% between Kv1.2 and Kv1.3. The homology models reveal that the T1 domains of these Kv channels exhibit similar folds as those of Shaker K⁺ channel. These models also show that each T1 monomer consists of three distinct layers, with N-terminal layer 1 and C- terminal layer 3 facing the cytoplasm and the membrane, respectively. Layer 2 exhibits the highest structural conservation because it is located around the central hydrophobic core. For each Kv channel, four identical subunits assemble into the homotetramer architecture around a four-fold axis through the hydrogen bonds and salt bridges formed by 15 highly conserved polar residues. The narrowest opening of the pore is formed by the four conserved residues corresponding to R115 of the Shaker T1 domain. The homology models of these Kv T1 domains provide particularly attractive targets for further structure-based studies.     

The coiled-coil domain of Escherichia coli FtsLB is a structurally detuned element critical for modulating its activation in bacterial cell division

The FtsLB complex is a key regulator of bacterial cell division, existing in either an off state or an on state, which supports the activation of septal peptidoglycan synthesis. In Escherichia coli, residues known to be critical for this activation are located in a region near the C-terminal end of the periplasmic coiled-coil domain of FtsLB, raising questions about the precise role of this conserved domain in the activation mechanism. Here, we investigate an unusual cluster of polar amino acids found within the core of the FtsLB coiled coil. We hypothesized that these amino acids likely reduce the structural stability of the domain and thus may be important for governing conformational changes. We found that mutating these positions to hydrophobic residues increased the thermal stability of FtsLB but caused cell division defects, suggesting that the coiled-coil domain is a “detuned” structural element. In addition, we identified suppressor mutations within the polar cluster, indicating that the precise identity of the polar amino acids is important for fine-tuning the structural balance between the off and on states. We propose a revised structural model of the tetrameric FtsLB (named the “Y-model”) in which the periplasmic domain splits into a pair of coiled-coil branches. In this configuration, the hydrophilic terminal moieties of the polar amino acids remain more favorably exposed to water than in the original four-helix bundle model (“I-model”). We propose that a shift in this architecture, dependent on its marginal stability, is involved in activating the FtsLB complex and triggering septal cell wall reconstruction.     

A Conserved Residue Cluster That Governs Kinetics of ATP-dependent Gating of Kir6.2 Potassium Channels

ATP-sensitive potassium (K_ATP) channels are heteromultimeric complexes of an inwardly rectifying Kir channel (Kir6.x) and sulfonylurea receptors. Their regulation by intracellular ATP and ADP generates electrical signals in response to changes in cellular metabolism. We investigated channel elements that control the kinetics of ATP-dependent regulation of K_ATP (Kir6.2 + SUR1) channels using rapid concentration jumps. WT Kir6.2 channels re-open after rapid washout of ATP with a time constant of ∼60 ms. Extending similar kinetic measurements to numerous mutants revealed fairly modest effects on gating kinetics despite significant changes in ATP sensitivity and open probability. However, we identified a pair of highly conserved neighboring amino acids (Trp-68 and Lys-170) that control the rate of channel opening and inhibition in response to ATP. Paradoxically, mutations of Trp-68 or Lys-170 markedly slow the kinetics of channel opening (500 and 700 ms for W68L and K170N, respectively), while increasing channel open probability. Examining the functional effects of these residues using φ value analysis revealed a steep negative slope. This finding implies that these residues play a role in lowering the transition state energy barrier between open and closed channel states. Using unnatural amino acid incorporation, we demonstrate the requirement for a planar amino acid at Kir6.2 position 68 for normal channel gating, which is potentially necessary to localize the ϵ-amine of Lys-170 in the phosphatidylinositol 4,5-bisphosphate-binding site. Overall, our findings identify a discrete pair of highly conserved residues with an essential role for controlling gating kinetics of Kir channels.     

RNA editing modulates the binding of drugs and highly unsaturated fatty acids to the open pore of Kv potassium channels

The time course of inactivation of voltage‐activated potassium (Kv) channels is an important determinant of the firing rate of neurons. In many Kv channels highly unsaturated lipids as arachidonic acid, docosahexaenoic acid and anandamide can induce fast inactivation. We found that these lipids interact with hydrophobic residues lining the inner cavity of the pore. We analysed the effects of these lipids on Kv1.1 current kinetics and their competition with intracellular tetraethylammonium and Kvβ subunits. Our data suggest that inactivation most likely represents occlusion of the permeation pathway, similar to drugs that produce ‘open‐channel block’. Open‐channel block by drugs and lipids was strongly reduced in Kv1.1 channels whose amino acid sequence was altered by RNA editing in the pore cavity, and in Kv1.x heteromeric channels containing edited Kv1.1 subunits. We show that differential editing of Kv1.1 channels in different regions of the brain can profoundly alter the pharmacology of Kv1.x channels. Our findings provide a mechanistic understanding of lipid‐induced inactivation and establish RNA editing as a mechanism to induce drug and lipid resistance in Kv channels.     

Molecular Analysis of the Putative Inactivation Particle in the Inactivation Gate of Brain Type IIA Na+ Channels

Fast Na⁺ channel inactivation is thought to involve binding of phenylalanine 1489 in the hydrophobic cluster IFM in L_III-IV of the rat brain type IIA Na⁺ channel. We have analyzed macroscopic and single channel currents from Na⁺ channels with mutations within and adjacent to hydrophobic clusters in L_III-IV. Substitution of F1489 by a series of amino acids disrupted inactivation to different extents. The degree of disruption was closely correlated with the hydrophilicity of the amino acid at position 1489. These mutations dramatically destabilized the inactivated state and also significantly slowed the entry into the inactivated state, consistent with the idea that F1489 forms a hydrophobic interaction with a putative receptor during the fast inactivation process. Substitution of a phe residue at position 1488 or 1490 in mutants lacking F1489 did not restore normal inactivation, indicating that precise location of F1489 is critical for its function. Mutations of T1491 disrupted inactivation substantially, with large effects on the stability of the inactivated state and smaller effects on the rate of entry into the inactivated state. Mutations of several other hydrophobic residues did not destabilize the inactivated state at depolarized potentials, indicating that the effects of mutations at F1489 and T1491 are specific. The double mutant YY1497/8QQ slowed macroscopic inactivation at all potentials and accelerated recovery from inactivation at negative membrane potentials. Some of these mutations in L_III-IV also affected the latency to first opening, indicating coupling between L_III-IV and channel activation. Our results show that the amino acid residues of the IFM hydrophobic cluster and the adjacent T1491 are unique in contributing to the stability of the inactivated state, consistent with the designation of these residues as components of the inactivation particle responsible for fast inactivation of Na⁺ channels.     

A prokaryotic glutamate receptor: homology modelling and molecular dynamics simulations of GluR0

GluR0 is a prokaryotic homologue of mammalian glutamate receptors that forms glutamate-activated, potassium-selective ion channels. The topology of its transmembrane (TM) domain is similar to that of simple potassium channels such as KcsA. Two plausible alignments of the sequence of the TM domain of GluR0 with KcsA are possible, differing in the region of the P helix. We have constructed homology models based on both alignments and evaluated them using 6 ns duration molecular dynamics simulations in a membrane-mimetic environment. One model, in which an insertion in GluR0 relative to KcsA is located in the loop between the M1 and P helices, is preferred on the basis of lower structural drift and maintenance of the P helix conformation during simulation. This model also exhibits inter-subunit salt bridges that help to stabilise the TM domain tetramer. During the simulation, concerted K(+) ion-water movement along the selectivity filter is observed, as is the case in simulations of KcsA. K(+) ion exit from the central cavity is associated with opening of the hydrophobic gate formed by the C-termini of the M2 helices. In the intact receptor the opening of this gate will be controlled by interactions with the extramembranous ligand-binding domains.     

Exploration of the chlorpyrifos escape pathway from acylpeptide hydrolases using steered molecular dynamics simulations

Acylpeptide hydrolases (APH) catalyze the removal of an N-acylated amino acid from blocked peptides. APH is significantly more sensitive than acetylcholinesterase, a target of Alzheimer’s disease, to inhibition by organophosphorus (OP) compounds. Thus, OP compounds can be used as a tool to probe the physiological functions of APH. Here, we report the results of a computational study of molecular dynamics simulations of APH bound to the OP compounds and an exploration of the chlorpyrifos escape pathway using steered molecular dynamics (SMD) simulations. In addition, we apply SMD simulations to identify potential escape routes of chlorpyrifos from hydrolase hydrophobic cavities in the APH-inhibitor complex. Two previously proposed APH pathways were reliably identified by CAVER 3.0, with the estimated relative importance of P1 > P2 for its size. We identify the major pathway, P2, using SMD simulations, and Arg526, Glu88, Gly86, and Asn65 are identified as important residues for the ligand leaving via P2. These results may help in the design of APH-targeting drugs with improved efficacy, as well as in understanding APH selectivity of the inhibitor binding in the prolyl oligopeptidase family.     

Positive and negative design in stability and thermal adaptation of natural proteins

The aim of this work is to elucidate how physical principles of protein design are reflected in natural sequences that evolved in response to the thermal conditions of the environment. Using an exactly solvable lattice model, we design sequences with selected thermal properties. Compositional analysis of designed model sequences and natural proteomes reveals a specific trend in amino acid compositions in response to the requirement of stability at elevated environmental temperature: the increase of fractions of hydrophobic and charged amino acid residues at the expense of polar ones. We show that this “from both ends of the hydrophobicity scale” trend is due to positive (to stabilize the native state) and negative (to destabilize misfolded states) components of protein design. Negative design strengthens specific repulsive non-native interactions that appear in misfolded structures. A pressure to preserve specific repulsive interactions in non-native conformations may result in correlated mutations between amino acids that are far apart in the native state but may be in contact in misfolded conformations. Such correlated mutations are indeed found in TIM barrel and other proteins.     

Exploring Volatile General Anesthetic Binding to a Closed Membrane-Bound Bacterial Voltage-Gated Sodium Channel via Computation

Despite the clinical ubiquity of anesthesia, the molecular basis of anesthetic action is poorly understood. Amongst the many molecular targets proposed to contribute to anesthetic effects, the voltage gated sodium channels (VGSCs) should also be considered relevant, as they have been shown to be sensitive to all general anesthetics tested thus far. However, binding sites for VGSCs have not been identified. Moreover, the mechanism of inhibition is still largely unknown. The recently reported atomic structures of several members of the bacterial VGSC family offer the opportunity to shed light on the mechanism of action of anesthetics on these important ion channels. To this end, we have performed a molecular dynamics “flooding” simulation on a membrane-bound structural model of the archetypal bacterial VGSC, NaChBac in a closed pore conformation. This computation allowed us to identify binding sites and access pathways for the commonly used volatile general anesthetic, isoflurane. Three sites have been characterized with binding affinities in a physiologically relevant range. Interestingly, one of the most favorable sites is in the pore of the channel, suggesting that the binding sites of local and general anesthetics may overlap. Surprisingly, even though the activation gate of the channel is closed, and therefore the pore and the aqueous compartment at the intracellular side are disconnected, we observe binding of isoflurane in the central cavity. Several sampled association and dissociation events in the central cavity provide consistent support to the hypothesis that the “fenestrations” present in the membrane-embedded region of the channel act as the long-hypothesized hydrophobic drug access pathway.     

An energy-efficient gating mechanism in the acetylcholine receptor channel suggested by molecular and Brownian dynamics

Acetylcholine receptors mediate electrical signaling between nerve and muscle by opening and closing a transmembrane ion conductive pore. Molecular and Brownian dynamics simulations are used to shed light on the location and mechanism of the channel gate. Four separate 5 ns molecular dynamics simulations are carried out on the imaged structure of the channel, a hypothetical open structure with a slightly wider pore and a mutant structure in which a central ring of hydrophobic residues is replaced by polar groups. Water is found to partially evacuate the pore during molecular simulations of the imaged structure, whereas ions face a large energy barrier and do not conduct through the channel in Brownian dynamics simulations. The pore appears to be in a closed configuration despite containing an unobstructed pathway across the membrane as a series of hydrophobic residues in the center of the channel provide an unfavorable home to water and ions. When the channel is widened slightly, water floods into the channel and ions conduct at a rate comparable to the currents measured experimentally in open channels. The pore remains permeable to ions provided the extracellular end of the pore-lining helix is restrained near the putative open configuration to mimic the presence of the ligand binding domain. Replacing some of the hydrophobic residues with polar ones decreases the barrier for ion permeation but does not result in significant currents. The channel is posited to utilize an energy efficient gating mechanism in which only minor conformational changes of the hydrophobic region of the pore are required to create macroscopic changes in conductance.     

BATMAS30: amino acid substitution matrix for alignment of bacterial transporters

Aligned amino acid sequences of three functionally independent samples of transmembrane (TM) transport proteins have been analyzed. The concept of TM-kernel is proposed as the most probable transmembrane region of a sequence. The average amino acid composition of TM-kernels differs from the published amino acid composition of transmembrane segments. TM-kernels contain more alanines, glycines, and less polar, charged, and aromatic residues in contrast to non-TM-proteins. There are also differences between TM-kernels of bacterial and eukaryotic proteins. We have constructed amino acid substitution matrices for bacterial TM-kernels, named the BATMAS (BActerial Transmembrane MAtrix of Substitutions) series. In TM-kernels, polar and charged residues, as well as proline and tyrosine, are highly conserved, whereas there are more substitutions within the group of hydrophobic residues, in contrast to non-TM-proteins that have fewer, relatively more conserved, hydrophobic residues. These results demonstrate that alignment of transmembrane proteins should be based on at least two amino acid substitution matrices, one for loops (e.g., the BLOSUM series) and one for TM-segments (the BATMAS series), and the choice of the TM-matrix should be different for eukaryotic and bacterial proteins.     

Brownian dynamics simulations of the recognition of the scorpion toxin maurotoxin with the voltage-gated potassium ion channels

The recognition of the scorpion toxin maurotoxin (MTX) by the voltage-gated potassium (Kv1) channels, Kv1.1, Kv1.2, and Kv1.3, has been studied by means of Brownian dynamics (BD) simulations. All of the 35 available structures of MTX in the Protein Data Bank (http://www.rcsb.org/pdb) determined by nuclear magnetic resonance were considered during the simulations, which indicated that the conformation of MTX significantly affected both the recognition and the binding between MTX and the Kv1 channels. Comparing the top five highest-frequency structures of MTX binding to the Kv1 channels, we found that the Kv1.2 channel, with the highest docking frequencies and the lowest electrostatic interaction energies, was the most favorable for MTX binding, whereas Kv1.1 was intermediate, and Kv1.3 was the least favorable one. Among the 35 structures of MTX, the 10th structure docked into the binding site of the Kv1.2 channel with the highest probability and the most favorable electrostatic interactions. From the MTX-Kv1.2 binding model, we identified the critical residues for the recognition of these two proteins through triplet contact analyses. MTX locates around the extracellular mouth of the Kv1 channels, making contacts with its beta-sheets. Lys23, a conserved amino acid in the scorpion toxins, protrudes into the pore of the Kv1.2 channel and forms two hydrogen bonds with the conserved residues Gly401(D) and Tyr400(C) and one hydrophobic contact with Gly401(C) of the Kv1.2 channel. The critical triplet contacts for recognition between MTX and the Kv1.2 channel are Lys23(MTX)-Asp402(C)(Kv1), Lys27(MTX)-Asp378(D)(Kv1), and Lys30(MTX)-Asp402(A)(Kv1). In addition, six hydrogen-bonding interactions are formed between residues Lys23, Lys27, Lys30, and Tyr32 of MTX and residues Gly401, Tyr400, Asp402, Asp378, and Thr406 of Kv1.2. Many of them are formed by side chains of residues of MTX and backbone atoms of the Kv1.2 channel. Five hydrophobic contacts exist between residues Pro20, Lys23, Lys30 and Tyr32 of MTX and residues Asp402, Val404, Gly401, and Arg377 of the Kv1.2 channel. The simulation results are in agreement with the previous molecular biology experiments and explain the binding phenomena between MTX and Kv1 channels at the molecular level. The consistency between the results of the BD simulations and the experimental data indicated that our three-dimensional model of the MTX-Kv1.2 channel complex is reasonable and can be used in additional biological studies, such as rational design of novel therapeutic agents blocking the voltage-gated channels and in mutagenesis studies in both the toxins and the Kv1 channels. In particular, both the BD simulations and the molecular mechanics refinements indicate that residue Asp378 of the Kv1.2 channel is critical for its recognition and binding functionality toward MTX. This phenomenon has not been appreciated in the previous mutagenesis experiments, indicating this might be a new clue for additional functional study of Kv1 channels.