The surface accessibility of the glycine receptor M2-M3 loop is increased in the channel open state

Mutations in the extracellular M2-M3 loop of the glycine receptor (GlyR) alpha1 subunit have been shown previously to affect channel gating. In this study, the substituted cysteine accessibility method was used to investigate whether a structural rearrangement of the M2-M3 loop accompanies GlyR activation. All residues from R271C to V277C were covalently modified by both positively charged methanethiosulfonate ethyltrimethylammonium (MTSET) and negatively charged methanethiosulfonate ethylsulfonate (MTSES), implying that these residues form an irregular surface loop. The MTSET modification rate of all residues from R271C to K276C was faster in the glycine-bound state than in the unliganded state. MTSES modification of A272C, L274C, and V277C was also faster in the glycine-bound state. These results demonstrate that the surface accessibility of the M2-M3 loop is increased as the channel transitions from the closed to the open state, implying that either the loop itself or an overlying domain moves during channel activation.

Nuclear hourglass technique: An approach that detects electrically open nuclear pores in Xenopus laevis oocyte

Nuclear pore complexes (NPCs) mediate both active transport and passive diffusion across the nuclear envelope (NE). Determination of NE electrical conductance, however, has been confounded by the lack of an appropriate technical approach. The nuclear patch clamp technique is restricted to preparations with electrically closed NPCs, and microelectrode techniques fail to resolve the extremely low input resistance of large oocyte nuclei. To address the problem, we have developed an approach for measuring the NE electrical conductance of Xenopus laevis oocyte nuclei. The method uses a tapered glass tube, which narrows in its middle part to 2/3 of the diameter of the nucleus. The isolated nucleus is sucked into the narrow part of the capillary by gentle fluid movement, while the resulting change in electrical resistance is monitored. NE electrical conductance was unexpectedly large (7.9 ± 0.34 S/cm2). Evaluation of NPC density by atomic force microscopy showed that this conductance corresponded to 3.7 × 106 NPCs. In contrast to earlier conclusions drawn from nuclear patch clamp experiments, NPCs were in an electrically “open” state with a mean single NPC electrical conductance of 1.7 ± 0.07 nS. Enabling or blocking of active NPC transport (accomplished by the addition of cytosolic extracts or gp62-directed antibodies) revealed this large NPC conductance to be independent of the activation state of the transport machinery located in the center of NPCs. We conclude that peripheral channels, which are presumed to reside in the NPC subunits, establish a high ionic permeability that is virtually independent of the active protein transport mechanism.

Comparative surface accessibility of a pore-lining threonine residue (T6') in the glycine and GABA(A) receptors

The substituted cysteine accessibility method was used to probe the surface exposure of a pore-lining threonine residue (T6') common to both the glycine receptor (GlyR) and gamma-aminobutyric acid, type A receptor (GABAAR) chloride channels. This residue lies close to the channel activation gate, the ionic selectivity filter, and the main pore blocker binding site. Despite their high amino acid sequence homologies and common role in conducting chloride ions, recent studies have suggested that the GlyRs and GABA(A)Rs have divergent open state pore structures at the 6' position. When both the human alpha1(T6'C) homomeric GlyR and the rat alpha1(T6'C)beta1(T6'C) heteromeric GABA(A)R were expressed in human embryonic kidney 293 cells, their 6' residue surface accessibilities differed significantly in the closed state. However, when a soluble cysteine-modifying compound was applied in the presence of saturating agonist concentrations, both receptors were locked into the open state. This action was not induced by oxidizing agents in either receptor. These results provide evidence for a conserved pore opening mechanism in anion-selective members of the ligand-gated ion channel family. The results also indicate that the GABA(A)R pore structure at the 6' level may vary between different expression systems.

Comparative surface accessibility of a pore-lining threonine residue (T6') in the glycine and GABA(A) receptors

The substituted cysteine accessibility method was used to probe the surface exposure of a pore-lining threonine residue (T6’) common to both the glycine receptor (GlyR) and GABAA receptor (GABAAR) chloride channels. This residue lies close to the channel activation gate, the ionic selectivity filter and the main pore blocker binding site. Recent studies have suggested that the GlyRs and GABAARs have divergent open state pore structures at the 6’ position. When both the human a1T6’C homomeric GlyR and the rat a1T6’Cb1T6’C heteromeric GABAAR were expressed in HEK293 cells, their 6’ residue surface accessibilities differed significantly in the closed state. However, when a soluble cysteine-modifying compound was applied in the presence of saturating agonist concentrations, both receptors were locked into the open state. This action was not induced by oxidising agents in either receptor. These results provide evidence for a conserved pore opening mechanism in anion-selective members of the ligand-gated ion channel family. The results also indicate that the GABAAR pore structure at the 6’ level may vary between different expression systems.

A single beta subunit M2 domain residue controls the picrotoxin sensitivity of alpha beta heteromeric glycine receptor chloride channels

This study investigated the residues responsible for the reduced picrotoxin sensitivity of the alpha beta heteromeric glycine receptor relative to the alpha homomeric receptor. By analogy with structurally related receptors, the beta subunit M2 domain residues P278 and F282 were considered the most likely candidates for mediating this effect. These residues align with G254 and T258 of the alpha subunit. The T258A, T258C and T258F mutations dramatically reduced the picrotoxin sensitivity of the alpha homomeric receptor. Furthermore, the converse F282T mutation in the beta subunit increased the picrotoxin sensitivity of the alpha beta heteromeric receptor. The P278G mutation in the beta subunit did not affect the picrotoxin sensitivity of the alpha beta heteromer. Thus, a ring of five threonines at the M2 domain depth corresponding to alpha subunit T258 is specifically required for picrotoxin sensitivity. Mutations to alpha subunit T258 also profoundly influenced the apparent glycine affinity. A substituted cysteine accessibility analysis revealed that the T258C sidechain increases its pore exposure in the channel open state. This provides further evidence for an allosteric mechanism of picrotoxin inhibition, but renders it unlikely that picrotoxin las an allosterically acting 'competitive' antagonist) binds to this residue.

Asymmetric contribution of α and β subunits to the activation of αβ heteromeric glycine receptors

This study investigated the role of beta subunits in the activation of alphabeta heteromeric glycine receptor (GlyR) chloride channels recombinantly expressed in HEK293 cells. The approach involved incorporating mutations into corresponding positions in alpha and beta subunits and comparing their effects on receptor function. Although cysteine-substitution mutations to residues in the N-terminal half of the alpha subunit M2-M3 loop dramatically impaired the gating efficacy, the same mutations exerted little effect when incorporated into corresponding positions of the beta subunit. Furthermore, although the alpha subunit M2-M3 loop cysteines were modified by a cysteine-specific reagent, the corresponding beta subunit cysteines showed no evidence of reactivity. These observations suggest structural or functional differences between alpha and beta subunit M2-M3 loops. In addition, a threonine-->leucine mutation at the 9' position in the beta subunit M2 pore-lining domain dramatically increased the glycine sensitivity. By analogy with the effects of the same mutation in other ligand-gated ion channels, it was concluded that the mutation affected the GlyR activation mechanism. This supports the idea that the GlyR beta subunit is involved in receptor gating. In conclusion, this study demonstrates that beta subunits contribute to the activation of the GlyR, but that their involvement in this process is significantly different to that of the alpha subunit.

Large scale assessment of regulatory evolution and transcriptome complexity in mammals

AbstractIn addition to genetic changes affecting the function of gene products, changes in gene expression have been suggested to underlie many or even most of the phenotypic differences among mammals. However, detailed gene expression comparisons were, until recently, restricted to closely related species, owing to technological limitations. Thus, we took advantage of the latest technologies (RNA-Seq) to generate extensive qualitative and quantitative transcriptome data for a unique collection of somatic and germline tissues from representatives of all major mammalian lineages (placental mammals, marsupials and monotremes) and birds, the evolutionary outgroup.In the first major project of my thesis, we performed global comparative analyses of gene expression levels based on these data. Our analyses provided fundamental insights into the dynamics of transcriptome change during mammalian evolution (e.g., the rate of expression change across species, tissues and chromosomes) and allowed the exploration of the functional relevance and phenotypic implications of transcription changes at a genome-wide scale (e.g., we identified numerous potentially selectively driven expression switches).In a second project of my thesis, which was also based on the unique transcriptome data generated in the context of the first project we focused on the evolution of alternative splicing in mammals. Alternative splicing contributes to transcriptome complexity by generating several transcript isoforms from a single gene, which can, thus, perform various functions. To complete the global comparative analysis of gene expression changes, we explored patterns of alternative splicing evolution. This work uncovered several general and unexpected patterns of alternative splicing evolution (e.g., we found that alternative splicing evolves extremely rapidly) as well as a large number of conserved alternative isoforms that may be crucial for the functioning of mammalian organs.Finally, the third and final project of my PhD consisted in analyzing in detail the unique functional and evolutionary properties of the testis by exploring the extent of its transcriptome complexity. This organ was previously shown to evolve rapidly both at the phenotypic and molecular level, apparently because of the specific pressures that act on this organ and are associated with its reproductive function. Moreover, my analyses of the amniote tissue transcriptome data described above, revealed strikingly widespread transcriptional activity of both functional and nonfunctional genomic elements in the testis compared to the other organs. To elucidate the cellular source and mechanisms underlying this promiscuous transcription in the testis, we generated deep coverage RNA-Seq data for all major testis cell types as well as epigenetic data (DNA and histone methylation) using the mouse as model system. The integration of these complete dataset revealed that meiotic and especially post-meiotic germ cells are the major contributors to the widespread functional and nonfunctional transcriptome complexity of the testis, and that this "promiscuous" spermatogenic transcription is resulting, at least partially, from an overall transcriptionally permissive chromatin state. We hypothesize that this particular open state of the chromatin results from the extensive chromatin remodeling that occurs during spermatogenesis which ultimately leads to the replacement of histones by protamines in the mature spermatozoa. Our results have important functional and evolutionary implications (e.g., regarding new gene birth and testicular gene expression evolution).Generally, these three large-scale projects of my thesis provide complete and massive datasets that constitute valuables resources for further functional and evolutionary analyses of mammalian genomes.

The epithelial sodium channel mediates the directionality of galvanotaxis in human keratinocytes.

Cellular directional migration in an electric field (galvanotaxis) is one of the mechanisms guiding cell movement in embryogenesis and in skin epidermal repair. The epithelial sodium channel (ENaC), in addition to its function of regulating sodium transport in kidney, has recently been found to modulate cell locomotory speed. Here we tested whether ENaC has an additional function of mediating the directional migration of galvanotaxis in keratinocytes. Genetic depletion of ENaC completely blocks only galvanotaxis and does not decrease migration speed. Overexpression of ENaC is sufficient to drive galvanotaxis in otherwise unresponsive cells. Pharmacologic blockade or maintenance of the open state of ENaC also decreases or increases, respectively, galvanotaxis, suggesting that the channel open state is responsible for the response. Stable lamellipodial extensions formed at the cathodal sides of wild-type cells at the start of galvanotaxis; these were absent in the ENaC knockout keratinocytes, suggesting that ENaC mediates galvanotaxis by generating stable lamellipodia that steer cell migration. We provide evidence that ENaC is required for directional migration of keratinocytes in an electric field, supporting a role for ENaC in skin wound healing.

Modulation of the c-Jun N-terminal kinase activity in the embryonic heart in response to anoxia-reoxygenation: involvement of the Ca2+ and mitoKATP channels.

Whether the response of the fetal heart to ischemia-reperfusion is associated with activation of the c-Jun N-terminal kinase (JNK) pathway is not known. In contrast, involvement of the sarcolemmal L-type Ca2+ channel (LCC) and the mitochondrial KATP (mitoKATP) channel has been established. This work aimed at investigating the profile of JNK activity during anoxia-reoxygenation and its modulation by LCC and mitoK(ATP) channel. Hearts isolated from 4-day-old chick embryos were submitted to anoxia (30 min) and reoxygenation (60 min). Using the kinase assay method, the profile of JNK activity in the ventricle was determined every 10 min throughout anoxia-reoxygenation. Effects on JNK activity of the LCC blocker verapamil (10 nM), the mitoK(ATP) channel opener diazoxide (50 microM) and the blocker 5-hydroxydecanoate (5-HD, 500 microM), the mitochondrial Ca2+ uniporter (MCU) inhibitor Ru360 (10 microM), and the antioxidant N-(2-mercaptopropionyl) glycine (MPG, 1 mM) were determined. In untreated hearts, JNK activity was increased by 40% during anoxia and peaked fivefold relative to basal level after 30-40 min reoxygenation. This peak value was reduced by half by diazoxide and was tripled by 5-HD. Furthermore, the 5-HD-mediated stimulation of JNK activity during reoxygenation was abolished by diazoxide, verapamil or Ru360. MPG had no effect on JNK activity, whatever the conditions. None of the tested pharmacological agents altered JNK activity under basal normoxic conditions. Thus, in the embryonic heart, JNK activity exhibits a characteristic pattern during anoxia and reoxygenation and the respective open-state of LCC, MCU and mitoKATP channel can be a major determinant of JNK activity in a ROS-independent manner.

Extracellular Subunit Interactions Control Transitions between Functional States of Acid-sensing Ion Channel 1a.

Acid-sensing ion channels (ASICs) are neuronal, voltage-independent Na(+) channels that are transiently activated by extracellular acidification. They are involved in pain sensation, the expression of fear, and in neurodegeneration after ischemic stroke. Our study investigates the role of extracellular subunit interactions in ASIC1a function. We identified two regions involved in critical intersubunit interactions. First, formation of an engineered disulfide bond between the palm and thumb domains leads to partial channel closure. Second, linking Glu-235 of a finger loop to either one of two different residues of the knuckle of a neighboring subunit opens the channel at physiological pH or disrupts its activity. This suggests that one finger-knuckle disulfide bond (E235C/K393C) sets the channel in an open state, whereas the other (E235C/Y389C) switches the channel to a non-conducting state. Voltage-clamp fluorometry experiments indicate that both the finger loop and the knuckle move away from the β-ball residue Trp-233 during acidification and subsequent desensitization. Together, these observations reveal that ASIC1a opening is accompanied by a distance increase between adjacent thumb and palm domains as well as a movement of Glu-235 relative to the knuckle helix. Our study identifies subunit interactions in the extracellular loop and shows that dynamic changes of these interactions are critical for normal ASIC function.

Études de type structure fonction du couplage électromécanique et de la coopérativité sous-unitaire chez les canaux potassiques dépendants du voltage

Les canaux potassiques voltage-dépendants forment des tétramères dont chaque sous-unité comporte six segments transmembranaires (S1 à S6). Le pore, formé des segments S5-S6 de chaque sous-unité, est entouré de quatre domaines responsables de la sensibilité au potentiel membranaire, les senseurs de voltage (VS; S1-S4). Lors d’une dépolarisation membranaire, le mouvement des résidus chargés situés dans le VS entraine un mouvement de charges détectable en électrophysiologie, le courant de « gating ». L’activation du VS conduit à l'ouverture du pore, qui se traduit par un changement de conformation en C-terminal du segment S6. Pour élucider les principes qui sous-tendent le couplage électromécanique entre ces deux domaines, nous avons étudié deux régions présumées responsables du couplage chez les canaux de type Shaker K+, soit la région carboxy-terminale du segment S6 et le lien peptidique reliant les segments transmembranaire S4-S5 (S4-5L). Avec la technique du « cut-open voltage clamp fluorometry » (COVCF), nous avons pu déterminer que l’interaction inter-sous-unitaire RELY, formée par des acides aminés situés sur le lien S4-5L et S6 de deux sous-unités voisines, est impliquée dans le développement de la composante lente observée lors du retour des charges de « gating » vers leur état de repos, le « OFF-gating ». Nous avons observé que l’introduction de mutations dans la région RELY module la force de ces interactions moléculaires et élimine l’asymétrie observée dans les courants de « gating » de type sauvage. D’ailleurs, nous démontrons que ce couplage inter-sous-unitaire est responsable de la stabilisation du pore dans l’état ouvert. Nous avons également identifié une interaction intra-sous-unitaire entre les résidus I384 situé sur le lien S4-5L et F484 sur le segment S6 d’une même sous-unité. La déstabilisation de cette interaction hydrophobique découple complètement le mouvement des senseurs de voltage et l'ouverture du pore. Sans cette interaction, l’énergie nécessaire pour activer les VS est moindre en raison de l’absence du poids mécanique appliqué par le pore. De plus, l’abolition du couplage électromécanique élimine également le « mode shift », soit le déplacement de la dépendance au voltage des charges de transfert (QV) vers des potentiels hyperpolarisants. Ceci indique que le poids mécanique du pore imposé au VS entraine le « mode shift », en modulant la conformation intrinsèque du VS par un processus allostérique.

Étude structurale et fonctionnelle du canal potassium dépendant du voltage KvAP

Les canaux ioniques dépendants du voltage sont responsables de l'initiation et de la propagation des potentiels d'action dans les cellules excitables. De nombreuses maladies héréditaires (channelopathies) sont associées à un contrôle défectueux du voltage par ces canaux (arythmies, épilepsie, etc.). L’établissement de la relation structure-fonction exacte de ces canaux est donc crucial pour le développement de nouveaux agents thérapeutiques spécifiques. Dans ce contexte, le canal procaryote dépendant du voltage et sélectif au potassium KvAP a servi de modèle d’étude afin d’approfondir i) le processus du couplage électromécanique, ii) l’influence des lipides sur l’activité voltage-dépendante et iii) l’inactivation de type closed-state. Afin de pallier à l’absence de données structurales dynamiques du côté cytosolique ainsi que de structure cristalline dans l’état fermé, nous avons mesuré le mouvement du linker S4-S5 durant le gating par spectroscopie de fluorescence (LRET). Pour ce faire, nous avons utilisé une technique novatrice du contrôle de l’état conformationnel du canal en utilisant les lipides (phospholipides et non phospholipides) au lieu du voltage. Un modèle dans l’état fermé a ainsi été produit et a démontré qu’un mouvement latéral modeste de 4 Å du linker S4-S5 est suffisant pour mener à la fermeture du pore de conduction. Les interactions lipides - canaux jouent un rôle déterminant dans la régulation de la fonction des canaux ioniques mais ne sont pas encore bien caractérisées. Nous avons donc également étudié l’influence de différents lipides sur l’activation voltage - dépendante de KvAP et mis en évidence deux sites distincts d’interactions menant à des effets différents : au niveau du senseur de voltage, menant au déplacement de la courbe conductance-voltage, et du côté intracellulaire, influençant le degré de la pente de cette même courbe. Nous avons également démontré que l’échange de lipides autour de KvAP est extrêmement limité et affiche une dépendance à l’état conformationnel du canal, ne se produisant que dans l’état ouvert. KvAP possède une inactivation lente particulière, accessible depuis l'état ouvert. Nous avons étudié les effets de la composition lipidique et de la température sur l'entrée dans l'état inactivé et le temps de récupération. Nous avons également utilisé la spectroscopie de fluorescence (quenching) en voltage imposé afin d'élucider les bases moléculaires de l’inactivation de type closed-state. Nous avons identifié une position à la base de l’hélice S4 qui semble impliquée à la fois dans le mécanisme responsable de ce type d'inactivation et dans la récupération particulièrement lente qui est typique du canal KvAP.

Estudo das estruturas terciária e quaternária de proteínas por espalhamento de raios X a baixos ângulos integrado a técnicas de modelagem

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Na +-coupled betaine symporters involved in osmotic stress response in prokaryotes and eukaryotes

The betaine/GABA transporter BGT1 is one of the most important osmolyte transporters in the kidney. BGT1 is a member of the neurotransmitter sodium symporter (NSS) family, facilitates Na+/Cl--coupled betaine uptake to cope with hyperosmotic stress. Betaine transport in kidney cells is upregulated under hypertonic conditions by a yet unknown mechanism when increasing amounts of intracellular BGT1 are inserted into the plasma membrane. Re-establishing isotonicity results in ensuing depletion of BGT1 from the membrane. BGT1 phosphorylation on serines and threonines might be a regulation mechanism. In the present study, four potential PKC phosphorylation sites were mutated to alanines and the responses to PKC activators, phorbol 12-myristate acetate (PMA) and dioctanoyl-sn-glycerol (DOG) were determined. GABA-sensitive currents were diminished after 30 min preincubation with these PKC activators. Staurosporine blocked the response to DOG. Three mutants evoked normal GABA-sensitive currents but currents in oocytes expressing the mutant T40A were greatly diminished. [3H]GABA uptake was also determined in HEK-293 cells expressing EGFP-tagged BGT1 with the same mutations. Three mutants showed normal upregulation of GABA uptake after hypertonic stress, and downregulation by PMA was normal compared to EGFP-BGT1. In contrast, GABA uptake by the T40A mutant showed no response to hypertonicity or PMA. Confocal microscopy of the EGFP-BGT1 mutants expressed in MDCK cells, grown on glass or filters, revealed that T40A was present in the cytoplasm after 24 h hypertonic stress while the other mutants and EGFP-BGT1 were predominantely present in the plasma membrane. All four mutants co-migrated with EGFP-BGT1 on Western blots suggesting they are full-length proteins. In conclusion, T235, S428, and S564 are not involved in downregulation of BGT1 due to phosphorylation by PKC. However, T40 near the N-terminus may be part of a hot spot important for normal trafficking or insertion of BGT1 into the plasma membrane. Additionally, a link between substrate transport regulation, insertion of BGT1 into the plasma membrane and N-glycosylation in the extracellular loop 2 (EL2) could be revealed. The functional importance of two predicted N-glycosylation sites, which are conserved in EL2 within the NSS family were investigated for trafficking, transport and regulated plasma membrane insertion by immunogold-labelling, electron microscopy, mutagenesis, two-electrode voltage clamp measurements in Xenopus laevis oocytes and uptake of radioactive-labelled substrate into MDCK cells. Trafficking and plasma membrane insertion of BGT1 was clearly promoted by proper N-glycosylation in both, oocytes and MDCK cells. De-glycosylation with PNGase F or tunicamycin led to a decrease in substrate affinity and transport rate. Mutagenesis studies revealed that in BGT1 N183 is the major N-glycosylation site responsible for full protein activity. Replacement of N183 with aspartate resulted in a mutant, which was not able to bind N-glycans suggesting that N171 is a non-glycosylated site in BGT1. N183D exhibited close to WT transport properties in oocytes. Surprisingly, in MDCK cells plasma membrane insertion of the N183D mutant was no longer regulated by osmotic stress indicating unambiguously that association with N-glycans at this position is linked to osmotic stress-induced transport regulation in BGT1. The molecular transport mechanism of BGT1 remains largely unknown in the absence of a crystal structure. Therefore investigating the structure-function relationship of BGT1 by a combination of structural biology (2D and 3D crystallization) and membrane protein biochemistry (cell culture, substrate transport by radioactive labeled GABA uptake into cells and proteoliposomes) was the aim of this work. While the functional assays are well established, structure determination of eukaryotic membrane transporters is still a challenge. Therefore, a suitable heterologous expression system could be defined, starting with cloning and overexpression of an optimized gene. The achieved expression levels in P. pastoris were high enough to proceed with isolation of BGT1. Furthermore, purification protocols could be established and resulted in pure protein, which could even be reconstituted in an active form. The quality and homogeneity of the protein allowed already 2D and 3D crystallization, in which initial crystals could be obtained. Interestingly, the striking structural similarity of BGT1 to the bacterial betaine transporter BetP, which became a paradigm for osmoregulated betaine transport, provided information on substrate coordination in BGT1. The structure of a BetP mutant that showed activity for GABA was solved to 3.2Å in complex with GABA in an inward facing open state. This structure shed some light into the molecular transport mechanisms in BGT1 and might help in future to design conformationally locked BGT1 to enforce the on-going structure determination.

Block of the hERG channel by bupivacaine: Electrophysiological and modeling insights towards stereochemical optimization

The hERG voltage-gated potassium channel mediates the cardiac I(Kr) current, which is crucial for the duration of the cardiac action potential. Undesired block of the channel by certain drugs may prolong the QT interval and increase the risk of malignant ventricular arrhythmias. Although the molecular determinants of hERG block have been intensively studied, not much is known about its stereoselectivity. Levo-(S)-bupivacaine was the first drug reported to have a higher affinity to block hERG than its enantiomer. This study strives to understand the principles underlying the stereoselectivity of bupivacaine block with the help of mutagenesis analyses and molecular modeling simulations. Electrophysiological measurements of mutated hERG channels allowed for the identification of residues involved in bupivacaine binding and stereoselectivity. Docking and molecular mechanics simulations for both enantiomers of bupivacaine and terfenadine (a non-stereoselective blocker) were performed inside an open-state model of the hERG channel. The predicted binding modes enabled a clear depiction of ligand-protein interactions. Estimated binding affinities for both enantiomers were consistent with electrophysiological measurements. A similar computational procedure was applied to bupivacaine enantiomers towards two mutated hERG channels (Tyr652Ala and Phe656Ala). This study confirmed, at the molecular level, that bupivacaine stereoselectively binds the hERG channel. These results help to lay the foundation for structural guidelines to optimize the cardiotoxic profile of drug candidates in silico.