Efficiency measures the conversion of agonist binding energy into receptor conformational change.

Receptors alternate between resting↔active conformations that bind agonists with low↔high affinity. Here, we define a new agonist attribute, energy efficiency (η), as the fraction of ligand-binding energy converted into the mechanical work of the activation conformational change. η depends only on the resting/active agonist-binding energy ratio. In a plot of activation energy versus binding energy (an "efficiency" plot), the slope gives η and the y intercept gives the receptor's intrinsic activation energy (without agonists; ΔG0). We used single-channel electrophysiology to estimate η for eight different agonists and ΔG0 in human endplate acetylcholine receptors (AChRs). From published equilibrium constants, we also estimated η for agonists of KCa1.1 (BK channels) and muscarinic, γ-aminobutyric acid, glutamate, glycine, and aryl-hydrocarbon receptors, and ΔG0 for all of these except KCa1.1. Regarding AChRs, η is 48-56% for agonists related structurally to acetylcholine but is only ∼39% for agonists related to epibatidine; ΔG0 is 8.4 kcal/mol in adult and 9.6 kcal/mol in fetal receptors. Efficiency plots for all of the above receptors are approximately linear, with η values between 12% and 57% and ΔG0 values between 2 and 12 kcal/mol. Efficiency appears to be a general attribute of agonist action at receptor binding sites that is useful for understanding binding mechanisms, categorizing agonists, and estimating concentration-response relationships.

Molecular recognition at cholinergic synapses: acetylcholine versus choline.

KEY POINTS: Neuromuscular acetylcholine (ACh) receptors have a high affinity for the neurotransmitter ACh and a low affinity for its metabolic product choline. At each transmitter binding site three aromatic groups determine affinity, and together provide ∼50% more binding energy for ACh than for choline. Deprotonation of αY190 by a nearby lysine strengthens the interaction between this aromatic ring and both ACh and choline. H-bonds position ACh and choline differently in the aromatic cage to generate the different affinities. ABSTRACT: Acetylcholine (ACh) released at the vertebrate nerve-muscle synapse is hydrolysed rapidly to choline (Cho), so endplate receptors (AChRs) are exposed to high concentrations of both of these structurally related ligands. To understand how these receptors distinguish ACh and Cho, we used single-channel electrophysiology to measure resting affinities (binding free energies) of these and other agonists in adult-type mouse AChRs having a mutation(s) at the transmitter-binding sites. The aromatic rings of αY190, αW149 and αY198 each provide ∼50% less binding energy for Cho compared to ACh. At αY198 a phenylalanine substitution had no effect, but at αY190 this substitution caused a large, agonist-independent loss in binding energy that depended on the presence of αK145. The results suggest that (1) αY190 is deprotonated by αK145 to strengthen the interaction between this benzene ring and the agonist's quaternary ammonium (QA) and (2) AChRs respond strongly to ACh because an H-bond positions the QA to interact optimally with the rings, and weakly to Cho because a different H-bond tethers the ligand to misalign the QA and form weaker interactions with the aromatic groups. The results suggest that the difference in ACh versus Cho binding energies is determined by different ligand positions within a fixed protein structure.

State-dependent inter-repeat contacts of exceptionally conserved asparagines in the inner helices of sodium and calcium channels.

Voltage-gated sodium and calcium channels play key roles in the physiology of excitable cells. The alpha-1 subunit of these channels folds from a polypeptide chain of four homologous repeats. In each repeat, the cytoplasmic halves of the pore-lining helices contain exceptionally conserved asparagines. Such conservation implies important roles, which are unknown. Mutations of the asparagines affect activation and inactivation gating as well as the action of pore-targeting ligands, including local anesthetics and steroidal agonists batrachotoxin and veratridine. In the absence of the open-channel structures, underlying mechanisms are unclear. Here, we modeled the pore module of Cav1.2 and Nav1.4 channels and their mutants in the open and closed states using the X-ray structures of potassium and sodium channels as templates. The energy of each model was Monte Carlo-minimized. The asparagines do not face the pore in the modeled states. In the open-channel models, the asparagine residue in a given repeat forms an inter-repeat H-bond with a polar residue, which is typically nine positions downstream from the conserved asparagine in the preceding repeat. The H-bonds, which are strengthened by surrounding hydrophobic residues, would stabilize the open channel and shape the open-pore geometry. According to our calculation, the latter is much more sensitive to mutations of the asparagines than the closed-pore geometry. Rearrangement of inter-repeat contacts may explain effects of these mutations on the voltage dependence of activation and inactivation and action of pore-targeting ligands.

Functional differences between neurotransmitter binding sites of muscle acetylcholine receptors.

A muscle acetylcholine receptor (AChR) has two neurotransmitter binding sites located in the extracellular domain, at αδ and either αε (adult) or αγ (fetal) subunit interfaces. We used single-channel electrophysiology to measure the effects of mutations of five conserved aromatic residues at each site with regard to their contribution to the difference in free energy of agonist binding to active versus resting receptors (ΔGB1). The two binding sites behave independently in both adult and fetal AChRs. For four different agonists, including ACh and choline, ΔGB1 is ∼-2 kcal/mol more favorable at αγ compared with at αε and αδ. Only three of the aromatics contribute significantly to ΔGB1 at the adult sites (αY190, αY198, and αW149), but all five do so at αγ (as well as αY93 and γW55). γW55 makes a particularly large contribution only at αγ that is coupled energetically to those contributions of some of the α-subunit aromatics. The hydroxyl and benzene groups of loop C residues αY190 and αY198 behave similarly with regard to ΔGB1 at all three kinds of site. ACh binding energies estimated from molecular dynamics simulations are consistent with experimental values from electrophysiology and suggest that the αγ site is more compact, better organized, and less dynamic than αε and αδ. We speculate that the different sensitivities of the fetal αγ site versus the adult αε and αδ sites to choline and ACh are important for the proper maturation and function of the neuromuscular synapse.

Catch-and-hold activation of muscle acetylcholine receptors having transmitter binding site mutations.

Agonists turn on receptors because their target sites have a higher affinity in the active versus resting conformation of the protein. We used single-channel electrophysiology to measure the lower-affinity (LA) and higher-affinity (HA) equilibrium dissociation constants for acetylcholine in adult-type muscle mouse nicotinic receptors (AChRs) having mutations of agonist binding site amino acids. For a series of agonists and for all mutations of αY93, αG147, αW149, αY190, αY198, εW55, and δW57, the change in LA binding energy was approximately half that in HA binding energy. The results were analyzed as a linear free energy relationship between LA and HA agonist binding, the slope of which (κ) gives the fraction of the overall binding chemical potential where the LA complex is established. The linear correlation between LA and HA binding energies suggests that the overall binding process is by an integrated mechanism (catch-and-hold). For the agonist and the above mutations, κ ∼ 0.5, but side-chain substitutions of two residues had a slope that was significantly higher (0.90; αG153) or lower (0.25; εP121). The results suggest that backbone rearrangements in loop B, loop C, and the non-α surface participate in both LA binding and the LA ↔ HA affinity switch. It appears that all of the intermediate steps in AChR activation comprise a single, energetically coupled process.

Fluoxetine blocks Nav1.5 channels via a mechanism similar to that of class 1 antiarrhythmics.

The voltage-gated Nav1.5 channel is essential for the propagation of action potentials in the heart. Malfunctions of this channel are known to cause hereditary diseases. It is a prime target for class 1 antiarrhythmic drugs and a number of antidepressants. Our study investigated the Nav1.5 blocking properties of fluoxetine, a selective serotonin reuptake inhibitor. Nav1.5 channels were expressed in HEK-293 cells, and Na(+) currents were recorded using the patch-clamp technique. Dose-response curves of racemic fluoxetine (IC50 = 39 µM) and its optical isomers had a similar IC50 [40 and 47 µM for the (+) and (-) isomers, respectively]. Norfluoxetine, a fluoxetine metabolite, had a higher affinity than fluoxetine, with an IC50 of 29 µM. Fluoxetine inhibited currents in a frequency-dependent manner, shifted steady-state inactivation to more hyperpolarized potentials, and slowed the recovery of Nav1.5 from inactivation. Mutating a phenylalanine (F1760) and a tyrosine (Y1767) in the S6 segment of domain (D) IV (DIVS6) significantly reduced the affinity of fluoxetine and its frequency-dependent inhibition. We used a noninactivating Nav1.5 mutant to show that fluoxetine displays open-channel block behavior. The molecular model of fluoxetine in Nav1.5 was in agreement with mutational experiments in which F1760 and Y1767 were found to be the key residues in binding fluoxetine. We concluded that fluoxetine blocks Nav1.5 by binding to the class 1 antiarrhythmic site. The blocking of cardiac Na(+) channels should be taken into consideration when prescribing fluoxetine alone or in association with other drugs that may be cardiotoxic or for patients with conduction disorders.

Energy for wild-type acetylcholine receptor channel gating from different choline derivatives.

Agonists, including the neurotransmitter acetylcholine (ACh), bind at two sites in the neuromuscular ACh receptor channel (AChR) to promote a reversible, global change in protein conformation that regulates the flow of ions across the muscle cell membrane. In the synaptic cleft, ACh is hydrolyzed to acetate and choline. Replacement of the transmitter's ester acetyl group with a hydroxyl (ACh→choline) results in a + 1.8 kcal/mol reduction in the energy for gating generated by each agonist molecule from a low- to high-affinity change of the transmitter binding site (ΔG(B)). To understand the distinct actions of structurally related agonist molecules, we measured ΔG(B) for 10 related choline derivatives. Replacing the hydroxyl group of choline with different substituents, such as hydrogen, chloride, methyl, or amine, increased the energy for gating (i.e., it made ΔG(B) more negative relative to choline). Extending the ethyl hydroxide tail of choline to propyl and butyl hydroxide also increased this energy. Our findings reveal the amount of energy that is available for the AChR conformational change provided by different, structurally related agonists. We speculate that a hydrogen bond between the choline hydroxyl and the backbone carbonyl of αW149 positions this agonist's quaternary ammonium group so as to reduce the cation-π interaction between this moiety and the aromatic groups at the binding site.

Sources of energy for gating by neurotransmitters in acetylcholine receptor channels.

Nicotinic acetylcholine receptors (AChRs) mediate signaling in the central and peripheral nervous systems. The AChR gating conformational change is powered by a low- to high-affinity change for neurotransmitters at two transmitter binding sites. We estimated (from single-channel currents) the components of energy for gating arising from binding site aromatic residues in the α-subunit. All mutations reduced the energy (TyrC1>>TrpB≈TyrC2>TyrA), with TyrC1 providing ~40% of the total. Considered one at a time, the fractional energy contributions from the aromatic rings were TrpB ~35%, TyrC1 ~28%, TyrC2 ~28%, and TyrA ~10%. Together, TrpB, TyrC1, and TyrC2 comprise an "aromatic triad" that provides much of the total energy from the transmitter for gating. Analysis of mutant pairs suggests that the energy contributions from some residues are nearly independent. Mutations of TyrC1 cause particularly large energy reductions because they remove two favorable and approximately equal interactions between the aromatic ring and the quaternary amine of the agonist and between the hydroxyl and αLysß7.

Design and control of acetylcholine receptor conformational change.

Allosteric proteins use energy derived from ligand binding to promote a global change in conformation. The "gating" equilibrium constant of acetylcholine receptor-channels (AChRs) is influenced by ligands, mutations, and membrane voltage. We engineered AChRs to have specific values of this constant by combining these perturbations, and then calculated the corresponding values for a reference condition. AChRs were designed to have specific rate and equilibrium constants simply by adding multiple, energetically independent mutations with known effects on gating. Mutations and depolarization (to remove channel block) changed the diliganded gating equilibrium constant only by changing the unliganded gating equilibrium constant (E(0)) and did not alter the energy from ligand binding. All of the tested perturbations were approximately energetically independent. We conclude that naturally occurring mutations mainly adjust E(0) and cause human disease because they generate AChRs that have physiologically inappropriate values of this constant. The results suggest that the energy associated with a structural change of a side chain in the gating isomerization is dissipated locally and is mainly independent of rigid body or normal mode motions of the protein. Gating rate and equilibrium constants are estimated for seven different AChR agonists using a stepwise engineering approach.

Subunit symmetry at the extracellular domain-transmembrane domain interface in acetylcholine receptor channel gating.

Transmitter molecules bind to synaptic acetylcholine receptor channels (AChRs) to promote a global channel-opening conformational change. Although the detailed mechanism that links ligand binding and channel gating is uncertain, the energy changes caused by mutations appear to be more symmetrical between subunits in the transmembrane domain compared with the extracellular domain. The only covalent connection between these domains is the pre-M1 linker, a stretch of five amino acids that joins strand ß10 with the M1 helix. In each subunit, this linker has a central Arg (Arg(3')), which only in the non-α-subunits is flanked by positively charged residues. Previous studies showed that mutations of Arg(3') in the α-subunit alter the gating equilibrium constant and reduce channel expression. We recorded single-channel currents and estimated the gating rate and equilibrium constants of adult mouse AChRs with mutations at the pre-M1 linker and the nearby residue Glu(45) in non-α-subunits. In all subunits, mutations of Arg(3') had similar effects as in the α-subunit. In the ε-subunit, mutations of the flanking residues and Glu(45) had only small effects, and there was no energy coupling between εGlu(45) and εArg(3'). The non-α-subunit Arg(3') residues had Φ-values that were similar to those for the α-subunit. The results suggest that there is a general symmetry between the AChR subunits during gating isomerization in this linker and that the central Arg is involved in expression more so than gating. The energy transfer through the AChR during gating appears to mainly involve Glu(45), but only in the α-subunits.

A homology model of the pore domain of a voltage-gated calcium channel is consistent with available SCAM data.

In the absence of x-ray structures of calcium channels, their homology models are used to rationalize experimental data and design new experiments. The modeling relies on sequence alignments between calcium and potassium channels. Zhen et al. (2005. J. Gen. Physiol. doi:10.1085/jgp.200509292) used the substituted cysteine accessibility method (SCAM) to identify pore-lining residues in the Ca(v)2.1 channel and concluded that their data are inconsistent with the symmetric architecture of the pore domain and published sequence alignments between calcium and potassium channels. Here, we have built K(v)1.2-based models of the Ca(v)2.1 channel with 2-(trimethylammonium)ethyl methanethiosulfonate (MTSET)-modified engineered cysteines and used Monte Carlo energy minimizations to predict their energetically optimal orientations. We found that depending on the position of an engineered cysteine in S6 and S5 helices, the ammonium group in the long flexible MTSET-modified side chain can orient into the inner pore, an interface between domains (repeats), or an interface between S5 and S6 helices. Different local environments of equivalent positions in the four repeats can lead to different SCAM results. The reported current inhibition by MTSET generally decreases with the predicted distances between the ammonium nitrogen and the pore axis. A possible explanation for outliers of this correlation is suggested. Our calculations rationalize the SCAM data, validate one of several published sequence alignments between calcium and potassium channels, and suggest similar spatial dispositions of S5 and S6 helices in voltage-gated potassium and calcium channels.

Access and binding of local anesthetics in the closed sodium channel.

Local anesthetics (LAs) are known to bind Na+ channels in the closed, open, and inactivated states and reach their binding sites via extracellular and intracellular access pathways. Despite intensive studies, no atomic-scale theory is available to explain the diverse experimental data on the LA actions. Here we attempt to contribute to this theory by simulating access and binding of LAs in the KcsA-based homology model of the closed Na+ channel. We used Monte Carlo minimizations to model the channel with representative local anesthetics N-(2,6-dimethylphenylcarbamoylmethyl)triethylammonium (QX-314), cocaine, and tetracaine. We found the nucleophilic central cavity to be a common binding region for the ammonium group of LAs, whose aromatic group can extend either along the pore axis (vertical binding mode) or to the III/IV domain interface (horizontal binding mode). The vertical mode was earlier predicted for the open channel, but only the horizontal mode is consistent with mutational data on the closed-channel block. To explore hypothetical access pathways of the permanently charged QX-314, we pulled the ligand via the selectivity filter, the closed activation gate, and the III/IV domain interface. Only the last pathway, which leads to the horizontal binding mode, did not impose steric obstacles. The LA ammonium group mobility within the central cavity was more restricted in the vertical mode than in the horizontal mode. Therefore, occupation of the selectivity-filter DEKA locus by a Na+ ion destabilizes the vertical mode, thus favoring the horizontal mode. LA binding in the closed channel requires the resident Na+ ion to leave the nucleophilic central cavity through the selectivity filter, whereas the LA egress should be coupled with reoccupation of the cavity by Na+. This hypothesis on the coupled movement of Na+ and LA in the closed channel explains seemingly contradictory data on how the outer-pore mutations as well as tetrodotoxin and micro-conotoxin binding affect the ingress and egress of LAs.

Chromanol 293B binding in KCNQ1 (Kv7.1) channels involves electrostatic interactions with a potassium ion in the selectivity filter.

The chromanol 293B (293B, trans-6-cyano-4-(N-ethylsulfonyl-N-methylamino)-3-hydroxy-2,2-dimethyl-chroman) is a lead compound of potential class III antiarrhythmics that inhibit cardiac I(Ks) potassium channels. These channels are formed by the coassembly of KCNQ1 (Kv7.1, KvLQT1) and KCNE1 subunits. Although homomeric KCNQ1 channels are the principal molecular targets, entry of KCNE1 to the channel complex enhances the chromanol block. Because closely related neuronal KCNQ2 potassium channels are insensitive to the drug, we used KCNQ1/KCNQ2 chimeras to identify the binding site of the inhibitor. We localized the putative drug receptor to the H5 selectivity filter and the S6 transmembrane segment. Single residues affecting 293B inhibition were subsequently identified through systematic exchange of amino acids that were either different in KCNQ1 and KCNQ2 or predicted by a docking model of 293B in the open and closed conformation of KCNQ1. Mutant channel proteins T312S, I337V, and F340Y displayed dramatically lowered sensitivity to chromanol block. The predicted drug binding receptor lies in the inner pore vestibule containing the lower part of the selectivity filter, and the S6 transmembrane domain also reported to be important for binding of benzodiazepines. We propose that the block of the ion permeation pathway involves hydrophobic interactions with the S6 transmembrane residues Ile337 and Phe340, and stabilization of chromanol 293B binding through electrostatic interactions of its oxygen atoms with the most internal potassium ion within the selectivity filter.

Monte Carlo-energy minimization of correolide in the Kv1.3 channel: possible role of potassium ion in ligand-receptor interactions.

BACKGROUND: Correolide, a nortriterpene isolated from the Costa Rican tree Spachea correa, is a novel immunosuppressant, which blocks Kv1.3 channels in human T lymphocytes. Earlier mutational studies suggest that correolide binds in the channel pore. Correolide has several nucleophilic groups, but the pore-lining helices in Kv1.3 are predominantly hydrophobic raising questions about the nature of correolide-channel interactions. RESULTS: We employed the method of Monte Carlo (MC) with energy minimization to search for optimal complexes of correolide in Kv1.2-based models of the open Kv1.3 with potassium binding sites 2/4 or 1/3/5 loaded with K+ ions. The energy was MC-minimized from many randomly generated starting positions and orientations of the ligand. In all the predicted low-energy complexes, oxygen atoms of correolide chelate a K+ ion. Correolide-sensing residues known from mutational analysis along with the ligand-bound K+ ion provide major contributions to the ligand-binding energy. Deficiency of K+ ions in the selectivity filter of C-type inactivated Kv1.3 would stabilize K+-bound correolide in the inner pore. CONCLUSION: Our study explains the paradox that cationic and nucleophilic ligands bind to the same region in the inner pore of K+ channels and suggests that a K+ ion is an important determinant of the correolide receptor and possibly receptors of other nucleophilic blockers of the inner pore of K+ channels.

Atomic determinants of state-dependent block of sodium channels by charged local anesthetics and benzocaine.

Molecular modeling predicts that a local anesthetic (LA) lidocaine binds to the resting and open Na(v)1.5 in different modes, interacting with LA-sensing residues known from experiments. Besides the major pathway via the open activation gate, LAs can reach the inner pore via a "sidewalk" between D3S6, D4S6, and D3P. The ammonium group of a cationic LA binds in the focus of the pore-helices macrodipoles, which also stabilize a Na(+) ion chelated by two benzocaine molecules. The LA's cationic group and a Na(+) ion in the selectivity filter repel each other suggesting that the Na(+) depletion upon slow inactivation would stabilize a LA, while a LA would stabilize slow-inactivated states.

KvAP-based model of the pore region of shaker potassium channel is consistent with cadmium- and ligand-binding experiments.

Potassium channels play fundamental roles in excitable cells. X-ray structures of bacterial potassium channels show that the pore-lining inner helices obstruct the cytoplasmic entrance to the closed channel KcsA, but diverge in widely open channels MthK and KvAP, suggesting a gating-hinge role for a conserved Gly in the inner helix. A different location of the gating hinge and a narrower open pore were proposed for voltage-gated Shaker potassium channels that have the Pro-473-Val-Pro motif. Two major observations back the proposal: cadmium ions lock mutant Val-476-Cys in the open state by bridging Cys-476 and His-486 in adjacent helices, and cadmium blocks the locked-open double mutant Val-474-Cys/Val-476-Cys by binding to Cys-474 residues. Here we used molecular modeling to show that the open Shaker should be as wide as KvAP to accommodate an open-channel blocker, correolide. We further built KvAP-, MthK-, and KcsA-based models of the Shaker mutants and Monte-Carlo-minimized them with constraints Cys-476-Cd(2+)-His-486. The latter were consistent with the KvAP-based model, causing a small-bend N-terminal to the Pro-473-Val-Pro motif. The constraints significantly distorted the MthK-based structure, making it similar to KvAP. The KcsA structure resisted the constraints. Two Cd(2+) ions easily block the locked-open KvAP-based model at Cys-474 residues, whereas constraining a single cadmium ion to four Cys-474 caused large conformational changes and electrostatic imbalance. Although mutual disposition of the voltage-sensor and pore domains in the KvAP x-ray structure is currently disputed, our results suggest that the pore-region domain retains a nativelike conformation in the crystal.