|
||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
ARTICLE |
Interaction between Extracellular TEA and an Aromatic Residue in Potassium Channels
Correspondence to Richard Horn: Richard.Horn{at}jefferson.edu
|
|
|---|
interaction between TEA and such an aromatic residue contributes to TEA block using the in vivo nonsense suppression method to incorporate a series of increasingly fluorinated Phe side chains at position 449. Fluorination, which is known to decrease the cation
binding ability of an aromatic ring, progressively increased the inhibitory constant Ki for the TEA block of Shaker. A larger increase in Ki was observed when the benzene ring of Phe449 was substituted by nonaromatic cyclohexane. These results support a strong cation
component to the TEA block. The data provide an empirical basis for choosing between Shaker models that are based on two classes of reported crystal structures for the bacterial channel KcsA, showing residue Tyr82 in orientations either compatible or incompatible with a cation
mechanism. We propose that the aromatic residue at this position in Shaker is favorably oriented for a cation
interaction with the permeation pathway. This choice is supported by high level ab initio calculations of the predicted effects of Phe modifications on TEA binding energy.
| INTRODUCTION |
|---|
|
|
|---|
interaction was critical to the high affinity binding of TEA. An optimal cation
interaction is supported only if the cation (in this case TEA) interacts with the face, not the edge, of the aromatic ring, an arrangement we will refer to as en face.
Reported crystal structures of KcsA actually show two orientations for residue Tyr82. In support of the cation
hypothesis, a recent crystal structure of a noninactivating mutant of KcsA shows an en face orientation of Tyr82 (Fig. 1 B, orange side chain; Cordero-Morales et al., 2006
), suggesting that this orientation might occur in eukaryotic potassium channels.
However, other reported KcsA structures pose a possible problem. Specifically, the side chain of Tyr82 in many published structures of KcsA (e.g., 1BL8 [Doyle et al., 1998
], 1K4C [Zhou et al., 2001
], 1R3J [Zhou and MacKinnon, 2003
], 2A9H [Yu et al., 2005
], and 1ZWI [Cordero-Morales et al., 2006
]; in the Protein Data Bank [www.rcsb.org/pdb/]) is aligned such that a TEA positioned to block the channel would not interact with the face of the aromatic ring (Fig. 1, blue side chains). Furthermore, several molecular dynamics studies of KcsA complexed with TEA support this nonoptimal geometry at Tyr82 such that the edges of the tyrosine rings point toward the blocker (Crouzy et al., 2001
; Luzhkov and Aqvist, 2001
; Guidoni and Carloni, 2002
; Luzhkov et al., 2003
). Indeed, a subsequent crystal structure of KcsA complexed with the TEA analogue tetraethylarsonium (TEAs) also has this edge-on conformation (Lenaeus et al., 2005
). Thus, it was concluded that a cation
interaction was not involved in binding TEA to KcsA but rather that the crucial Tyr influenced the local hydration structure of the binding site.
|
interaction at the aromatic residue at the entrance to potassium channels. Simulation studies of the edge-on KcsA channel argue against a cation
interaction, but the binding studies of TEA to eukaryotic potassium channels and the recent en face KcsA structure allow for such an interaction. In this study, we aim to test directly for a cation
mechanism in Shaker potassium channels using a definitive probe developed for functional ion channels and receptors. Modifying the side chain of an aromatic amino acid with fluorine substantially diminishes the cation
binding ability of the aromatic ring, and multiple substitutions produce additive effects (Dougherty, 1996
interactions. In comparison with conventional mutagenesis, the steric perturbation introduced by fluorine substitution is universally considered to be minimal. Also, fluorine substitution does not substantially alter the hydrophobicity of the ring. For example, benzene and hexafluorobenzene have nearly identical logP (water/octanol partition) values (Leo et al., 1971
In this study, we provide an empirical rationale for the choice between the edge-on and en face orientations by applying the side chain fluorination approach to position 449 of the Shaker channel. We find a compelling correlation between the degree of fluorination of Phe449 and the loss of TEA affinity. The correlation is supported by ab initio quantum mechanical calculations of the two disparate models of the TEA-binding site. We conclude that there is a significant cation
interaction involved in the binding of TEA to the Shaker channel, suggesting that an aromatic residue at this position adopts an en face orientation.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Electrophysiology
Voltage-clamped potassium currents were recorded with two microelectrodes using either the OpusXpress system (Axon Instruments, Inc.), which includes an automated solution delivery system, or an OC-725C voltage clamp (Warner) in a standard Ringers solution (116 mM NaCl, 2 mM KCl, 1 mM MgCl2, 0.5 mM CaCl2, and 5 mM HEPES, pH 7.5). The chloride salt of TEAs was a gift from A. Gross and M. Lenaeus (Northwestern University Medical School, Chicago, IL).
Computations
For ab initio calculations, we used Gaussian 03 software (Gaussian). To ensure an appropriate level of theory for the examination of cation
interactions, we optimized the structure of benzene plus a Na+ ion (see Fig. 5 A) with Hartree-Fock (HF) self-consistent field theory using the 6-31G** polarized basis set (Mecozzi et al., 1996a
). The same level of theory was used to determine the single point energy of this structure and those containing derivatives of benzene. For binding energy determinations, we used counterpoise correction for basis set superposition errors (Simon et al., 1996
), which were always <1 kcal/mol.
The 117-atom reduced molecular model of Fig. 5 B is based on the coordinates of TEAs cocrystallized with KcsA (Lenaeus et al., 2005
). First, the hydroxyls of Tyr82 residues were removed to convert them into Phe. The arsenic of TEAs was then replaced by a nitrogen, and the structure of the resultant TEA molecule was optimized at the HF/6-31G* level while freezing the position of the central nitrogen and the four Phe82 residues. This optimization preserved both the D2d conformation of the blocker and the approximate direction of its four alkyl chains. The en face model of Fig. 5 C was generated from the reduced molecular model by a 60° rotation around the CßC
bond of Phe82. To replace Phe82 with cyclohexylalanine (Cha), we substituted the aromatic ring with cyclohexane and optimized the structure of the ring only (freezing the position of the C
carbon) in the context of the reduced model at the HF/STO-3G level. TEA binding energy was then calculated at the HF/6-31G** level as described above for Na+ binding to benzene. Structures are displayed with DS/ViewerPro software (Accelrys).
| RESULTS |
|---|
|
|
|---|
role in the interaction between extracellular TEA and an aromatic residue at position 449 of Shaker, we used Phe derivatives chosen to reduce the negative electrostatic potential on the face of the aromatic ring. Fig. 2 shows examples of inactivation-removed Shaker potassium currents obtained with this method.
In all cases, the cRNA encoding for Shaker was identical, containing a UAG codon at position 449. This cRNA was coinjected into oocytes with one of six different constructs of suppressor tRNAs, each containing the anticodon CUA and all but one synthetically acylated with an amino acid. Fig. 2 (AD) shows families of potassium currents in which the tRNAs were acylated with Phe and three of its fluorinated derivatives. Insets (Fig. 2) show color-coded electrostatic potential surfaces, which were determined by ab initio calculations, of benzene and these fluorinated derivatives. On this color scale, red and blue indicate negative and positive electrostatic potential, respectively. As shown, successive substitution of fluorine atoms for hydrogen atoms on a benzene ring produces a successive decrease in the negative electrostatic potential on the face of the aromatic ring, leading to a monotonic decrease in intrinsic cation binding ability (Mecozzi et al., 1996a
|
25 in Shaker T449F (unpublished data).
|
|
|
interaction involving the aromatic side chain at residue 449, one would expect increasing the substitution of fluorines to decrease the binding affinity monotonically (Mecozzi et al., 1996a
interaction involving the site undergoing substitution (Zhong et al., 1998
interaction between acetylcholine and a Trp residue in the nicotinic acetylcholine receptor (Zhong et al., 1998
The largest effect on TEA affinity, 640-fold, occurs with the nonaromatic residue Cha. As discussed in the next section, this observation provides further support for a cation
interaction in TEA binding and contributes valuable information on proposed alternative models underlying TEA affinity.
| DISCUSSION |
|---|
|
|
|---|
interaction between TEA and an aromatic amino acid at position 449 in Shaker. As in other studies of this type (Zhong et al., 1998
contribution to TEA block. Moreover, a previous study shows that fluorination has little effect on nonelectrostatic energetic components of cation binding to a benzene ring (e.g., donor acceptor, charge transfer, and induced dipoles in the aromatic; Ma and Dougherty, 1997
|
interaction, a comparison of Phe and 3,4,5F3-Phe is useful. Examination of the electrostatic potential surfaces of Fig. 2 shows that three fluorines effectively erase the electrostatic attractiveness of the aromatic ring, and earlier studies showed that trifluorobenzene is a good model for an aromatic that has little or no electrostatic binding ability (Williams, 1993
G values of Table I suggest that the cation
interaction is the dominant factor in distinguishing strong from weak TEA binders in that the energetic consequence of the Phe
F3-Phe mutation (3.5 kcal/mol) is comparable with that of the Phe
Thr mutation (2.6 kcal/mol).
The biggest decrease in TEA affinity is seen for the Phe
Cha mutation (3.8 kcal/mol). This is consistent with both theoretical and experimental studies showing that cyclohexane binds cations more poorly than any of the aromatics considered here (Shepodd et al., 1988
; Mecozzi et al., 1996a
). Note that the side chains of Phe and Cha, modeled here as benzene and cyclohexane, respectively, are similar in size, shape, and hydrophobicity (Fig. 2). The most noticeable difference, which is in the context of noncovalent binding interactions, is that the negative electrostatic potential that leads to cation binding in benzene disappears in cyclohexane. Thus, the Phe
Cha mutation further confirms that the
character and its associated negative electrostatic potential dominate strong TEA binding when Phe is present at position 449. Also note that cyclohexane is considerably more polarizable than benzene (Craven et al., 1989
), ruling out any special role for an induced dipole in the binding of TEA by Phe and derivatives.
Previous computational studies of the KcsA channel are based on structures with Tyr82 in the edge-on orientation and have led others to conclude that the cation
interaction does not contribute substantially to external TEA blockade (Crouzy et al., 2001
; Luzhkov and Aqvist, 2001
; Guidoni and Carloni, 2002
; Luzhkov et al., 2003
). These studies instead emphasize local hydration structure as a crucial component of the energetics of TEA binding, although electrostatic stabilization by the backbone carbonyls of Tyr78 and Gly79 (in KcsA) may also play a role (Guidoni and Carloni, 2002
). The distribution of K+ ions within the selectivity filter also appears to affect TEA binding energy (Crouzy et al., 2001
; Luzhkov and Aqvist, 2001
; Guidoni and Carloni, 2002
), another possible source of electrostatic energy. Energetic contributions of local hydration include both electrostatic and hydrophobic forces. However, an essential role of attractive hydrophobic forces between TEA and its binding site is contradicted by earlier data from Shaker, in which relatively hydrophobic residues at site 449 are not generally better at binding TEA (Heginbotham and MacKinnon, 1992
). For example, comparison of the essentially isosteric residues Val and Thr shows that polar Thr is associated with a higher TEA affinity (Heginbotham and MacKinnon, 1992
). Also, Leu, Ile, and Val are all substantially more hydrophobic than Phe (Radzicka et al., 1988
) but produce weak TEA-binding sites. The Phe/Cha pair described here presents another isosteric comparison, and, again, hydrophobics do not explain the results. In addition, the size but not hydrophobicity of the alkylammonium blocker appears to be a critical factor in binding affinity (Hille, 1967
; Heginbotham and MacKinnon, 1992
). Furthermore, examination of the temperature dependence of TEA block shows that enthalpy rather than entropy largely accounts for the free energy of binding (Heginbotham and MacKinnon, 1992
), again arguing against a fundamental role of hydrophobic forces in TEA binding.
Computational Examination of the TEAChannel Interaction
To further investigate the evidence for a cation
mechanism, we modeled the TEAside chain interaction computationally. We were constrained to use KcsA structures for this purpose because it is the only crystallized potassium channel with an aromatic residue (Tyr82) at this position. We specifically started with the structure 2BOC, a cocrystal of KcsA and TEAs (Fig. 1, blue side chain; Lenaeus et al., 2005
). KcsA is highly homologous to voltage-gated potassium channels in this region of the protein. In the vicinity of the selectivity filter, KcsA has the sequence TTVGYGDLY, whereas in KV1.2, it is TTVGYGDMV. In addition, the crystallographic distances across the tetramer between the
carbons of these residues agree within 1 Å when comparing KcsA structures with the mammalian KV1.2 (2A79). With respect to a cation
mechanism, the majority of published structures of KcsA display an unfavorable orientation of the Tyr82 aromatic rings. Moreover, the closest carboncarbon distance between the blocker and the aromatic side chains (4.1 Å) is larger than typically observed for simple cation
interactions (e.g., 2.4 Å between Na+ and the center of a benzene ring; Dougherty, 1996
; Ma and Dougherty, 1997
). To examine the energetic consequences of withdrawing
electrons by fluorinating the aromatic ring of this residue, we constructed a reduced molecular model of TEA and its four coordinating Phe residues based on the 2BOC structure (Lenaeus et al., 2005
). In this reduced system (117 atoms; Fig. 5 B), we performed high level ab initio quantum mechanical calculations of the binding energy of TEA. Instead of decreasing the binding, as we observed experimentally (Fig. 5 A, squares), successively fluorinating these Phe residues in silico had a nonmonotonic effect on TEA binding energy (Fig. 5 A, diamonds), the main effect being a substantial increase of TEA affinity in the di- and trifluorinated derivatives.
Our ab initio calculations can be rationalized as follows. The cationic blocker is directed more toward the edge than the face of the aromatic ring (Fig. 5 B). The electrostatic potential along the edge of a simple aromatic ring is positive (Fig. 2 A, inset) and would therefore electrostatically repel TEA. However, fluorination introduces negative electrostatic potential to the edge of the ring via the fluorines (Fig. 2), and this would tend to attract TEA, especially in the three or five position of the aromatic ring (Figs. 1 and 5 B). In fact, the 3,5-difluoro and 3,4,5-trifluoro derivatives of Phe show the most dramatic increases in TEA affinity in the ab initio calculations.
Because these self-consistent theoretical results disagree with our experimental data, we considered the possibility that the aromatic rings were rotated to face the central axis of the pore (Fig. 5 C), as proposed before the crystal structures of KcsA were available (Heginbotham and MacKinnon, 1992
; Kumpf and Dougherty, 1993
). In this en face orientation, TEA would be attracted to the negative electrostatic potential on the face of the aromatic rings. Fluorination would then decrease the negative electrostatic potential presented to TEA and destabilize its binding, which is a prediction consistent with our experimental data. To quantify this idea, we calculated the energetic consequences of fluorinating the aromatic rings of Phe82 on TEA binding in this conformation. The only modification of the structure of Fig. 5 B was a 60° rotation of the aromatic rings, maintaining all other coordinates of the four Phe residues. With this single change, the results of the ab initio calculations now resemble our experimental results, showing a monotonic destabilization of TEA binding as the aromatic rings are increasingly fluorinated (Fig. 5 A, triangles). Moreover, the rotation itself enhanced TEA affinity by
3 kcal/mol. Although we did not examine the optimal dihedral angle in detail, rotations of 40 or 80° were less effective at replicating our experimental data. Given the crudeness of this model and the fact that it is based on a crystal structure of KcsA rather than a eukaryotic voltage-gated potassium channel, the agreement is quite acceptable. Even the larger perturbation associated with the nonaromatic Cha residue is reproduced by the calculations.
Although the 144°
2 dihedral angle in our en face model is unusual for both Tyr and Phe residues in known protein structures (Dunbrack and Cohen, 1997
), an en face orientation of Tyr82 was reported for one of the crystal structures of the functional E71A mutant of KcsA (Cordero-Morales et al., 2006
). Fig. 1 B shows this structure (2ATK; orange Tyr82 side chains) aligned with that of KcsA cocrystallized with TEAs (blue Tyr82 side chains). The similarity to our en face model (Fig. 5 C) is apparent. There are other subtle structural changes in this region of the 2ATK structure that enable an en face geometry. We also note that the selectivity filter has conformational variants in wild-type KcsA depending on the concentration of permeant ions (Zhou et al., 2001
), supporting the possibility that the functional conformation of Shaker is compatible with an en face orientation of Phe449 residues. It also raises the possibility that the open states of potassium channels in general have conformational flexibility, including the orientation of an aromatic residue at this position.
The slope of the experimental fluorination plot as shown in Fig. 5 A is considered to be a crude indication of the energetic magnitude of the cation
interaction. The slope of this plot (0.19) is comparable with those seen in studies of agonists such as GABA and serotonin binding to neuroreceptors (Beene et al., 2002
; Lummis et al., 2005
). However, the receptor experiments measure the binding of an
ion to a single aromatic residue. Here, we are simultaneously modifying four aromatic residues. Thus, the slope in the present case suggests that any single TEAaromatic interaction is not very strong. This is consistent with the low affinity of TEA and the fact that TEA is a very diffuse cation and so is expected to experience a weaker cation
interaction than more focused cations such as
. This may also reflect the fact that TEA cannot make van der Waals contact with all four rings at once in either of the en face conformations we have considered (Fig. 1 B and 5 C), and this should weaken the cation
interaction. Nevertheless, the magnitude of the cation
interaction does not fall off steeply with distance (Dougherty, 1996
), and, therefore, four nonoptimal interactions can combine to create a moderate effect. It is worth noting that a closer approach of the aromatic side chains to the blocker may be impossible if the open-channel structure is maintained because this residue is only two amino acids downstream from the TVGYGD signature sequence of the selectivity filter. Functional data argue against a large-scale collapse of the aromatics onto the blocker because TEA has no effect on open dwell times of single channels (Spruce et al., 1987
).
Although the calculations we describe do not consider all energetic factors, this is intentional. Our calculations primarily evaluate one component of the binding interaction: the electrostatic attraction between a cation and the side chain (the cation
interaction). The fact that binding energies correlate so nicely with this term establishes an important role for electrostatic attractions in the binding. We neither rule out nor rule in additional effects such as hydrophobics; we simply note that electrostatics must be important. The functional data on eukaryotic potassium channels show that although there may well be a hydrophobic contribution to binding, it is not the factor that discriminates among residues. Another factor that might contribute to TEA binding in our experiments is the possibility that the electronegative fluorine atoms could influence the local structure of water molecules in the vicinity of the blocker. However, it is not clear how this effect would produce the systematic destabilization of TEA binding that we observe (Figs. 3 C and 5 A). To explore this possibility, we estimated the Gibbs free energy of TEA binding in our reduced edge-on model, accounting both for the gas phase ab initio energies as well as the energies of hydration using a self-consistent reaction field method with the Poisson-Boltzmann solver in the program Jaguar (Schrodinger). As in the gas phase calculations of Fig. 5 A, difluorinating the Phe residue in this conformation enhanced TEA binding (by 2.3 kcal/mol; not depicted), which is again inconsistent with our experimental results and, therefore, tends to rule out the edge-on orientation of Phe449 in Shaker.
How do we resolve our conclusions with earlier works disputing a role for a cation
mechanism? As noted, previous simulations were based on the analysis of two edge-on crystallographic structures (1BL8 and 2BOC) of the bacterial channel KcsA rather than on a eukaryotic voltage-gated potassium channel or on the en face structure of a KcsA mutant (2ATK). Note that although KcsA contains an aromatic residue at the aligned position (Tyr82), it does not bind TEA strongly. Its Ki for TEA block is 3.2 mM (Heginbotham et al., 1999
), which is almost 10-fold higher than for Shaker T449F. Also, KcsA Y82T shows a TEA-blocking constant of 143 mM, which is significantly higher than that obtained in Shaker with its wild-type Thr449. Finally, the TEA analogue TEAs used for the 2BOC crystal structure shown in Fig. 1 has a 25-fold larger inhibitory constant than TEA in Shaker T449F, allowing for the possibility that the structure may not be representative of TEA in its binding site. We conclude that inherent but unknown aspects of the KcsA channel, perhaps including a nonoptimal arrangement of Tyr82 residues, underlie a relatively weak TEA-binding site.
In conclusion, we show that in voltage-gated potassium channels, a cation
interaction makes a substantial contribution to high affinity blockade by extracellular TEA. Evidently, the orientation of an aromatic residue at this position is favorable for a cation
mechanism. This hypothesis is supported by our experimental and theoretical results, and it provides a simple explanation of the longstanding observation that aromatics at position 449 are necessary and sufficient for high affinity external TEA blockade in voltage-gated potassium channels. We do not yet understand the (presumably subtle) conditions that govern the crystallization of KcsA in one state or the other, but the functional form relevant to Shaker in this region is apparently the en face orientation.
| ACKNOWLEDGMENTS |
|---|
This research was sponsored by grants from the National Institutes of Health (GM079427, NS11756, and NS34407).
Lawrence G. Palmer served as editor.
Submitted: 21 August 2006
Accepted: 27 October 2006
| REFERENCES |
|---|
|
|
|---|
interactions. Science. 261:17081710.
binding interaction with a tyrosine in the binding site of the GABAc receptor. Chem. Biol. 12:993997.[CrossRef][Medline]
interaction. Chem. Rev. 97:13031324.[CrossRef][Medline]
interactions in simple aromatics: electrostatics provide a predictive tool. J. Am. Chem. Soc. 118:23072308.[CrossRef]
binding site in the nicotinic receptor. Proc. Natl. Acad. Sci. USA. 95:1208812093.Related Article
This article has been cited by other articles:
![]() |
X. Xu, M. Recanatini, M. Roberti, and G.-N. Tseng Probing the Binding Sites and Mechanisms of Action of Two Human Ether-a-go-go-Related Gene Channel Activators, 1,3-bis-(2-Hydroxy-5-trifluoromethyl-phenyl)-urea (NS1643) and 2-[2-(3,4-Dichloro-phenyl)-2,3-dihydro-1H-isoindol-5-ylamino]-nicotinic acid (PD307243) Mol. Pharmacol., June 1, 2008; 73(6): 1709 - 1721. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. A. Ahern, A. L. Eastwood, D. A. Dougherty, and R. Horn Electrostatic Contributions of Aromatic Residues in the Local Anesthetic Receptor of Voltage-Gated Sodium Channels Circ. Res., January 4, 2008; 102(1): 86 - 94. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Schroeder and U.-P. Hansen Saturation and Microsecond Gating of Current Indicate Depletion-induced Instability of the MaxiK Selectivity Filter J. Gen. Physiol., July 1, 2007; 130(1): 83 - 97. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. P. Santarelli, A. L. Eastwood, D. A. Dougherty, R. Horn, and C. A. Ahern A Cation-{pi} Interaction Discriminates among Sodium Channels That Are Either Sensitive or Resistant to Tetrodotoxin Block J. Biol. Chem., March 16, 2007; 282(11): 8044 - 8051. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Roux Extracellular Blockade of Potassium Channels by TEA+: The Tip of the Iceberg? J. Gen. Physiol., December 1, 2006; 128(6): 635 - 636. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|