Physicochemical features of the HERG channel drug binding site

Blockade of hERG K(+) channels in the heart is an unintentional side effect of many drugs and can induce cardiac arrhythmia and sudden death. It has become common practice in the past few years to screen compounds for hERG channel activity early during the drug discovery process. Understanding the molecular basis of drug binding to hERG is crucial for the rational design of medications devoid of this activity. We previously identified 2 aromatic residues, Tyr-652 and Phe-656, located in the S6 domain of hERG, as critical sites of interaction with structurally diverse drugs. Here, Tyr-652 and Phe-656 were systematically mutated to different residues to determine how the physicochemical properties of the amino acid side group affected channel block by cisapride, terfenadine, and MK-499. The potency for block by all three drugs was well correlated with measures of hydrophobicity, especially the two-dimensional approximation of the van der Waals hydrophobic surface area of the side chain of residue 656. For residue 652, an aromatic side group was essential for high affinity block, suggesting the importance of a cation-pi interaction between Tyr-652 and the basic tertiary nitrogen of these drugs. hERG also lacks a Pro-Val-Pro motif common to the S6 domain of most other voltage-gated K(+) channels. Introduction of Pro-Val-Pro into hERG reduced sensitivity to drugs but also altered channel gating. Together, these findings assign specific residues to receptor fields predicted by pharmacophore models of hERG channel blockers and provide a refined molecular understanding of the drug binding site.     

High affinity HERG K(+) channel blockade by the antiarrhythmic agent dronedarone: resistance to mutations of the S6 residues Y652 and F656

Pharmacological inhibition of human-ether-a-go-go-related gene (HERG) K(+) channels by structurally and therapeutically diverse drugs is associated with the 'acquired' form of long QT syndrome and with potentially lethal cardiac arrhythmias. Two aromatic amino-acid residues (Y652 and F656) on the inner (S6) helices are considered to be key constituents of a high affinity drug binding site within the HERG channel pore cavity. Using wild-type (WT) and mutant HERG channels expressed in mammalian cell lines, we have investigated HERG channel current (I(HERG)) blockade at 37+/-1 degrees C by dronedarone (DRONED), a non-iodinated analogue of the Class III antiarrhythmic agent amiodarone (AMIOD). Under our conditions WT I(HERG) tails, measured at -40 mV following activating pulses to +30 mV, were blocked with IC(50) values of approximately 59 and 70 nM for DRONED and AMIOD, respectively. I(HERG) inhibition by DRONED was contingent upon channel gating, with block developing rapidly on membrane depolarization, but with no preference for activated over inactivated channels. High external [K(+)] (94 mM) reduced the potency of I(HERG) inhibition by both DRONED and AMIOD. Strikingly, mutagenesis to alanine of the S6 residue F656 (F656A) failed to eliminate blockade by both DRONED and AMIOD, whilst Y652A had comparatively little effect on DRONED but some effect on AMIOD. These findings demonstrate that high affinity drug blockade of I(HERG) can occur without a strong dependence on the Y652 and F656 aromatic amino-acid residues.     

Effect of sibutramine HCl on cardiac hERG K+ channel

Common clinically used drugs block the delayed rectifier K(+) channels and prolong the cardiac action potential duration associated with long QT syndrome. Here, we investigated the mechanism of hERG K(+) channel current (I(hERG)) blockade expressed in HEK-293 cells by sibutramine HCl, a serotonin-norepinephrine reuptake inhibitor. Sibutramine HCl inhibited I (hERG) in a concentration-dependent manner with the half-maximal inhibitory concentration (IC(50)) value of 2.5 microM at -40 mV. I(hERG) inhibition by sibutramine HCl showed weak voltage dependency, but the time-dependence of I(hERG) inhibition was developed relatively rapidly on membrane depolarization. On hERG channel gating for the S6 and pore regions, the S6 residue hERG mutant Y652A and F656A largely reduced the blocking potency of I(hERG), unlike the pore-region mutants T623A and S624A. These results indicate that sibutramine HCl preferentially inhibits the hERG potassium channel through the residue Y652 and F656, in a supratherapeutic concentration should be avoided by patients with high susceptibility for cardiac arrhythmia.     

The synergic modeling for the binding of fluoroquinolone antibiotics to the hERG potassium channel

The fluoroquinolone antibiotic binding site in the hERG potassium channel was examined for the residues involved and their position in the tetrameric channel. The blocking effect of the two fluoroquinolones levofloxacin and sparfloxacin to tandem dimers of the hERG mutants were evaluated electrophysiologically. The results indicated that two Tyr652s in the neighboring subunits and one or two Phe656s in the diagonal subunits contributed to the blockade in the case of both compounds, and Ser624 was also involved. The docking studies suggested that the protonated carboxyl group in the compounds strongly interacts with Phe656 as a π acceptor.     

Mutational analysis of block and facilitation of HERG current by a class III anti-arrhythmic agent, nifekalant

Chemicals and toxins are useful tools to elucidate the structure-function relationship of various proteins including ion channels. The HERG channel is blocked by many compounds and this may cause life-threatening cardiac arrhythmia. Besides block, some chemicals such as the class III anti-arrhythmic agent nifekalant stimulate HERG at low potentials by shifting its activation curve towards hyperpolarizing voltages. This is called "facilitation". Here, we report mutations and simulations analyzing the association between nifekalant and channel pore residues for block and facilitation. Alanine-scanning mutagenesis was performed in the pore region of HERG. The mutations at the base of the pore helix (T623A), the selectivity filter (V625A) and the S6 helix (G648A, Y652A and F656A) abolished and S624A attenuated both block and facilitation induced by the drug. On the other hand, the mutation of other residues caused either an increase or a decrease in nifekalant-induced facilitation without affecting block. An open-state homology model of the HERG pore suggested that T623, S624, Y652 and F656 faced the central cavity, and were positioned within geometrical range for the drug to be able to interact with all of them at the same time. Of these, S649 was the only polar residue located within possible interaction distance from the drug held in its blocking position. Further mutations and flexible-docking simulations suggest that the size, but not the polarity, of the side chain at S649 is critical for drug induced facilitation.     

Mutational Analysis of Block and Facilitation of HERG Current by A Class III Anti-Arrhythmic Agent,Nifekalant

Chemicals and toxins are useful tools to elucidate the structure-function relationship of various proteins including ion channels. The HERG channel is blocked by many compounds and this may cause life-threatening cardiac arrhythmia. Besides block, some chemicals such as the class III anti-arrhythmic agent nifekalant stimulate HERG at low potentials by shifting its activation curve towards hyperpolarizing voltages. This is called "facilitation". Here, we report mutations and simulations analyzing the association between nifekalant and channel pore residues for block and facilitation. Alanine-scanning mutagenesis was performed in the pore region of HERG. The mutations at the base of the pore helix (T623A), the selectivity filter (V625A) and the S6 helix (G648A, Y652A and F656A) abolished and S624A attenuated both block and facilitation induced by the drug. On the other hand, the mutation of other residues caused either an increase or a decrease in nifekalant-induced facilitation without affecting block. An open-state homology model of the HERG pore suggested that T623, S624, Y652 and F656 faced the central cavity, and were positioned within geometrical range for the drug to be able to interact with all of them at the same time. Of these, S649 was the only polar residue located within possible interaction distance from the drug held in its blocking position. Further mutations and flexible-docking simulations suggest that the size, but not the polarity, of the side chain at S649 is critical for drug induced facilitation.     

hERG K+ channel blockade by the antipsychotic drug thioridazine: An obligatory role for the S6 helix residue F656

The phenothiazine antipsychotic agent thioridazine has been linked with prolongation of the QT interval on the electrocardiogram, ventricular arrhythmias, and sudden death. Although thioridazine is known to inhibit cardiac hERG K(+) channels there is little mechanistic information on this action. We have investigated in detail hERG K(+) channel current (I(hERG)) blockade by thioridazine and identified a key molecular determinant of blockade. Whole-cell I(hERG) measurements were made at 37 degrees C from human embryonic kidney (HEK-293) cells expressing wild-type and mutant hERG channels. Thioridazine inhibited I(hERG) tails at -40mV following a 2s depolarization to +20mV with an IC(50) value of 80nM. Comparable levels of I(hERG) inhibition were seen with physiological command waveforms (ventricular and Purkinje fibre action potentials). Thioridazine block of I(hERG) was only weakly voltage-dependent, though the time dependence of I(hERG) inhibition indicated contingency of blockade upon channel gating. The S6 helix point mutation F656A almost completely abolished, and the Y652A mutation partially attenuated, I(hERG) inhibition by thioridazine. In summary, thioridazine is one of the most potent hERG K(+) channel blockers amongst antipsychotics, exhibiting characteristics of a preferential open/activated channel blocker and binding at a high affinity site in the hERG channel pore.     

Exploring blocker binding to a homology model of the open hERG K+ channel using docking and molecular dynamics methods

Binding of blockers to the human voltage-gated hERG potassium channel is studied using a combination of homology modelling, automated docking calculations and molecular dynamics simulations, where binding affinities are evaluated using the linear interaction energy method. A homology model was constructed based on the available crystal structure of the bacterial KvAP channel and the affinities of a series of sertindole analogues predicted using this model. The calculations reproduce the relative binding affinities of these compounds very well and indicate that both polar interactions near the intracellular opening of the selectivity filter as well as hydrophobic complementarity in the region around F656 are important for blocker binding. These results are consistent with recent alanine scanning mutation experiments on the blocking of the hERG channel by other compounds.     

A novel hypothesis for the binding mode of HERG channel blockers

We present a new docking model for HERG channel blockade. Our new model suggests three key interactions such that (1) a protonated nitrogen of the channel blocker forms a hydrogen bond with the carbonyl oxygen of HERG residue T623; (2) an aromatic moiety of the channel blocker makes a pi-pi interaction with the aromatic ring of HERG residue Y652; and (3) a hydrophobic group of the channel blocker forms a hydrophobic interaction with the benzene ring of HERG residue F656. The previous model assumes two interactions such that (1) a protonated nitrogen of the channel blocker forms a cation-pi interaction with the aromatic ring of HERG residue Y652; and (2) a hydrophobic group of the channel blocker forms a hydrophobic interaction with the benzene ring of HERG residue F656. To test these models, we classified 69 known HERG channel blockers into eight binding types based on their plausible binding modes, and further categorized them into two groups based on the number of interactions our model would predict with the HERG channel (two or three). We then compared the pIC(50) value distributions between these two groups. If the old hypothesis is correct, the distributions should not differ between the two groups (i.e., both groups show only two binding interactions). If our novel hypothesis is correct, the distributions should differ between Groups 1 and 2. Consistent with our hypothesis, the two groups differed with regard to pIC(50), and the group having more predicted interactions with the HERG channel had a higher mean pIC(50) value. Although additional work will be required to further validate our hypothesis, this improved understanding of the HERG channel blocker binding mode may help promote the development of in silico predictions methods for identifying potential HERG channel blockers.     

Erythromycin block of the HERG K+ channel: accessibility to F656 and Y652

The HERG potassium channel might have a non-canonical drug binding site, distinct from the channel's inner cavity, that could be responsible for elements of closed-state pharmacological inhibition of the channel. The macrolide antibiotic erythromycin is a drug that may block unconventionally because of its size. Here we used whole-cell patch-clamp recording at 37 degrees C from heterologously expressed HERG channels in a mammalian cell line to show that erythromycin either produces a rapid open-state-dependent HERG channel inhibition, or components of both open-state-dependent and closed-state-dependent inhibition. Alanine-substitution of HERG's canonical determinants of blockade revealed that Y652 was not important as a molecular determinant of blockade, and that mutation of F656 resulted in only weak attenuation of inhibition. In computer models of the channel, erythromycin could make several direct contacts with F656, but not with Y652, in the open-state model, and erythromycin was unable to fit into a closed-state channel model.     

Rutaecarpine targets hERG channels and participates in regulating electrophysiological properties leading to ventricular arrhythmia

Drug-mediated or medical condition-mediated disruption of hERG function accounts for the main cause of acquired long-QT syndrome (acLQTs), which predisposes affected individuals to ventricular arrhythmias (VA) and sudden death. Many Chinese herbal medicines, especially alkaloids, have risks of arrhythmia in clinical application. The characterized mechanisms behind this adverse effect are frequently associated with inhibition of cardiac hERG channels. The present study aimed to assess the potent effect of Rutaecarpine (Rut) on hERG channels. hERG-HEK293 cell was applied for evaluating the effect of Rut on hERG channels and the underlying mechanism. hERG current (I_hERG) was measured by patch-clamp technique. Protein levels were analysed by Western blot, and the phosphorylation of Sp1 was determined by immunoprecipitation. Optical mapping and programmed electrical stimulation were used to evaluate cardiac electrophysiological activities, such as APD, QT/QTc, occurrence of arrhythmia, phase singularities (PSs), and dominant frequency (DF). Our results demonstrated that Rut reduced the I_hERG by binding to F656 and Y652 amino acid residues of hERG channel instantaneously, subsequently accelerating the channel inactivation, and being trapped in the channel. The level of hERG channels was reduced by incubating with Rut for 24 hours, and Sp1 in nucleus was inhibited simultaneously. Mechanismly, Rut reduced threonine (Thr)/ tyrosine (Tyr) phosphorylation of Sp1 through PI3K/Akt pathway to regulate hERG channels expression. Cell-based model unables to fully reveal the pathological process of arrhythmia. In vivo study, we found that Rut prolonged QT/QTc intervals and increased induction rate of ventricular fibrillation (VF) in guinea pig heart after being dosed Rut for 2 weeks. The critical reasons led to increased incidence of arrhythmias eventually were prolonged APD₉₀ and APD₅₀ and the increase of DF, numbers of PSs, incidence of early after-depolarizations (EADs). Collectively, the results of this study suggest that Rut could reduce the I_hERG by binding to hERG channels through F656 and Y652 instantaneously. While, the PI3K/Akt/Sp1 axis may play an essential role in the regulation of hERG channels, from the perspective of the long-term effects of Rut (incubating for 24 hours). Importantly, the changes of electrophysiological properties by Rut were the main cause of VA.     

Flexibility of aromatic residues in the active-site gorge of acetylcholinesterase: X-ray versus molecular dynamics

The high aromatic content of the deep and narrow active-site gorge of acetylcholinesterase (AChE) is a remarkable feature of this enzyme. Here, we analyze conformational flexibility of the side chains of the 14 conserved aromatic residues in the active-site gorge of Torpedo californica AChE based on the 47 three-dimensional crystal structures available for the native enzyme, and for its complexes and conjugates, and on a 20-ns molecular dynamics (MD) trajectory of the native enzyme. The degree of flexibility of these 14 aromatic side chains is diverse. Although the side-chain conformations of F330 and W279 are both very flexible, the side-chain conformations of F120, W233, W432, Y70, Y121, F288, F290 and F331 appear to be fixed. Residues located on, or adjacent to, the Ω-loop (C67-C94), namely W84, Y130, Y442, and Y334, display different flexibilities in the MD simulations and in the crystal structures. An important outcome of our study is that the majority of the side-chain conformations observed in the 47 Torpedo californica AChE crystal structures are faithfully reproduced by the MD simulation on the native enzyme. Thus, the protein can assume these conformations even in the absence of the ligand that permitted their experimental detection. These observations are pertinent to structure-based drug design.     

A molecular switch driving inactivation in the cardiac k(+) channel HERG

K(+) channels control transmembrane action potentials by gating open or closed in response to external stimuli. Inactivation gating, involving a conformational change at the K(+) selectivity filter, has recently been recognized as a major K(+) channel regulatory mechanism. In the K(+) channel hERG, inactivation controls the length of the human cardiac action potential. Mutations impairing hERG inactivation cause life-threatening cardiac arrhythmia, which also occur as undesired side effects of drugs. In this paper, we report atomistic molecular dynamics simulations, complemented by mutational and electrophysiological studies, which suggest that the selectivity filter adopts a collapsed conformation in the inactivated state of hERG. The selectivity filter is gated by an intricate hydrogen bond network around residues S620 and N629. Mutations of this hydrogen bond network are shown to cause inactivation deficiency in electrophysiological measurements. In addition, drug-related conformational changes around the central cavity and pore helix provide a functional mechanism for newly discovered hERG activators.     

Molecular determinants of voltage-dependent human ether-a-go-go related gene (HERG) K+ channel block

The structural determinants for the voltage-dependent block of ion channels are poorly understood. Here we investigate the voltage-dependent block of wild-type and mutant human ether-a-go-go related gene (HERG) K(+) channels by the antimalarial compound chloroquine. The block of wild-type HERG channels expressed in Xenopus oocytes was enhanced as the membrane potential was progressively depolarized. The IC(50) was 8.4 +/- 0.9 microm when assessed during 4-s voltage clamp pulses to 0 mV. Chloroquine also slowed the apparent rate of HERG deactivation, reflecting the inability of drug-bound channels to close. Mutation to alanine of aromatic residues (Tyr-652 or Phe-656) located in the S6 domain of HERG greatly reduced the potency of channel block by chloroquine (IC(50) > 1 mm at 0 mV). However, mutation of Tyr-652 also altered the voltage dependence of the block. In contrast to wild-type HERG, block of Y652A HERG channels was diminished by progressive membrane depolarization, and complete relief from block was observed at +40 mV. HERG channel block was voltage-independent when the hydroxyl group of Tyr-652 was removed by mutating the residue to Phe. Together these findings indicate a critical role for Tyr-652 in voltage-dependent block of HERG channels. Molecular modeling was used to define energy-minimized dockings of chloroquine to the central cavity of HERG. Our experimental findings and modeling suggest that chloroquine preferentially blocks open HERG channels by cation-pi and pi-stacking interactions with Tyr-652 and Phe-656 of multiple subunits.     

Aromatic components of two ferric enterobactin binding sites in Escherichia coli FepA

Ferric enterobactin is a catecholate siderophore that binds with high affinity (Kd approximately 10-10 M) to the Escherichia coli outer membrane protein FepA. We studied the involvement of aromatic amino acids in its uptake by determining the binding affinities, kinetics and transport properties of site-directed mutants. We replaced seven aromatic residues (Y260, Y272, Y285, Y289, W297, Y309 and F329) in the central part of FepA primary structure with alanine, individually and in double combinations, and determined the ability of the mutant proteins to interact with ferric enterobactin and the protein toxins colicins B and D. All the constructs showed normal expression and localization. Among single mutants, Y260A and F329A were most detrimental, reducing the affinity between FepA and ferric enterobactin 100- and 10-fold respectively. Double substitutions involving Y260, Y272 and F329 impaired (100- to 2500-fold) adsorption of the iron chelate more strongly. For Y260A and Y272A, the drop in adsorption affinity caused commensurate decreases in transport efficiency, suggesting that the target residues primarily act in ligand binding. F329A, like R316A, showed greater impairment of transport than binding, intimating mechanistic involvement during ligand internalization. Furthermore, immunochemical studies localized F329 in the FepA ligand binding site. The mutagenesis results suggested the existence of dual ligand binding sites in the FepA vestibule, and measurements of the rate of ferric enterobactin adsorption to fluoresceinated FepA mutant proteins confirmed this conclusion. The initial, outermost site contains aromatic residues and probably functions through hydrophobic interactions, whereas the secondary site exists deeper in the vestibule, contains both charged and aromatic residues and probably acts through hydrophobic and electrostatic bonds.     

The flexibility of P-glycoprotein for its poly-specific drug binding from molecular dynamics simulations

The multidrug efflux pump P-glycoprotein (P-gp) contributes to multidrug resistance in about half of human cancers. Recently, high resolution X-ray crystal structures of mouse P-gp (inward-facing) were reported, which significantly facilitates the understanding of the function of P-gp and the structure-based design of inhibitors for P-gp. Here we perform 20?ns molecular dynamics simulations of inward-facing P-gp with/without ligand in explicit lipid and water to investigate the flexibility of P-gp for its poly-specific drug binding. By analyzing the interactions between P-gp and QZ59-RRR or QZ59-SSS, we summarize the important residues and the flexibility of different parts of P-gp. Particularly, the flexibility of the side chains of aromatic residues (Phe and Tyr) allows them to form rotamers with different orientations in the binding pocket, which plays a critical role for the poly-specificity of the drug-binding cavity of P-gp. MD simulations reveal that trans-membrane (TM) TM12 and TM6 are flexible and contribute to the poly-specific drug binding, while TM4 and TM5 are rigid and stabilize the whole structure. We also construct outward-facing P-gp based on the MsbA structure and perform 20?ns MD simulations. The comparison between the MD results for outward-facing P-gp and those for inward-facing P-gp shows that the TM parts in outward-facing P-gp undergo significant conformational change to facilitate the export of small molecules.     

Conformational dynamics of Tetracenomycin aromatase/cyclase regulate polyketide binding and enzyme aggregation propensity

BackgroundThe N-terminal domain of Tetracenomycin aromatase/cyclase (TcmN), an enzyme derived from Streptomyces glaucescens, is involved in polyketide cyclization, aromatization, and folding. Polyketides are a diverse class of secondary metabolites produced by certain groups of bacteria, fungi, and plants with various pharmaceutical applications. Examples include antibiotics, such as tetracycline, and anticancer drugs, such as doxorubicin. Because TcmN is a promising enzyme for in vitro production of polyketides, it is important to identify conditions that enhance its thermal resistance and optimize its function.MethodsTcmN unfolding, stability, and dynamics were evaluated by fluorescence spectroscopy, circular dichroism, nuclear magnetic resonance ¹⁵N relaxation experiments, and microsecond molecular dynamics (MD) simulations.ResultsTcmN thermal resistance was enhanced at low protein and high salt concentrations, was pH-dependent, and denaturation was irreversible. Conformational dynamics on the μs-ms timescale were detected for residues in the substrate-binding cavity, and two predominant conformers representing opened and closed cavity states were observed in the MD simulations.ConclusionBased on the results, a mechanism was proposed in which the thermodynamics and kinetics of the TcmN conformational equilibrium modulate enzyme function by favoring ligand binding and avoiding aggregation.General significanceUnderstanding the principles underlying TcmN stability and dynamics may help in designing mutants with optimal properties for biotechnological applications.     

Hydrophobic residues in small ankyrin 1 participate in binding to obscurin

Abstract Small ankyrin-1 is a splice variant of the ANK1 gene that binds to obscurin A. Previous studies have identified electrostatic interactions that contribute to this interaction. In addition, molecular dynamics (MD) simulations predict four hydrophobic residues in a 'hot spot' on the surface of the ankyrin-like repeats of sAnk1, near the charged residues involved in binding. We used site-directed mutagenesis, blot overlays and surface plasmon resonance assays to study the contribution of the hydrophobic residues, V70, F71, I102 and I103, to two different 30-mers of obscurin that bind sAnk1, Obsc????????? and Obsc?????????. Alanine mutations of each of the hydrophobic residues disrupted binding to the high affinity binding site, Obsc?????????. In contrast, V70A and I102A mutations had no effect on binding to the lower affinity site, Obsc?????????. Alanine mutagenesis of the five hydrophobic residues present in Obsc????????? showed that V6328, I6332, and V6334 were critical to sAnk1 binding. Individual alanine mutants of the six hydrophobic residues of Obsc????????? had no effect on binding to sAnk1, although a triple alanine mutant of residues V6233/I6234/I6235 decreased binding. We also examined a model of the Obsc?????????-sAnk1 complex in MD simulations and found I102 of sAnk1 to be within 2.2? of V6334 of Obsc?????????. In contrast to the I102A mutation, mutating I102 of sAnk1 to other hydrophobic amino acids such as phenylalanine or leucine did not disrupt binding to obscurin. Our results suggest that hydrophobic interactions contribute to the higher affinity of Obsc????????? for sAnk1 and to the dominant role exhibited by this sequence in binding.     

Electrostatic interaction-mediated conformational changes of adipocyte fatty acid binding protein probed by molecular dynamics simulation

Abstract

Adipocyte fatty acid binding protein (A-FABP) is a potential drug target for treatment of diabetes, obesity and atherosclerosis. Molecular dynamics (MD) simulations, principal component (PC) analysis and binding free energy calculations were combined to probe effect of electrostatic interactions of residues R78, R106 and R126 with inhibitors ZGB, ZGC and IBP on structural stability of three inhibitor/A-FABP complexes. The results indicate that mutation R126A produces significant influence on polar interactions of three inhibitors with A-FABP and these interactions are main force for driving the conformational change of A-FABP. Analyses on hydrogen bond interactions show that the decrease in hydrogen bonding interactions of residues R126 and Y128 with three inhibitors and the increase in that of K58 with inhibitors ZGC and IBP in the R126A mutated systems mostly regulate the conformational changes of A-FABP. This work shows that R126A can generate a significant perturbation on structural stability of A-FABP, which implies that R126 is of significance in inhibitor bindings. We expect that this study can provide a theoretical guidance for design of potent inhibitors targeting A-FABP.

Communicated by Ramaswamy H. Sarma     

Interactions of granisetron with an agonist-free 5-HT3A receptor model

A new homology model of type-3A serotonin receptors (5-HT(3A)Rs) was built on the basis of the electron microscopic structure of the nicotinic acetylcholine receptor and with an agonist-free binding cavity. The new model was used to re-evaluate the interactions of granisetron, a 5-HT(3A)R antagonist. Docking of granisetron identified two possible binding modes, including a newly identified region for antagonists formed by loop B, C, and E residues. Amino acid residues L184-D189 in loop B were mutated to alanine, while Y143 and Y153 in loop E were mutated to phenylalanine. Mutation H185A resulted in no detectable granisetron binding, while D189A resulted in a 22-fold reduction in affinity. Y143F and Y153F decreased granisetron affinity to the same extent as Y143A and Y153A mutations, supporting the role of the OH groups of these tyrosines in loop E. Modeling and mutation studies suggest that granisetron plays its antagonist role by hindering the closure of the back wall of the binding cavity.