International Union of Basic and Clinical Pharmacology. CX. Classification of Receptors for 5-hydroxytryptamine; Pharmacology and Function.

5-HT receptors expressed throughout the human body are targets for established therapeutics and various drugs in development. Their diversity of structure and function reflects the important role 5-HT receptors play in physiologic and pathophysiological processes. The present review offers a framework for the official receptor nomenclature and a detailed understanding of each of the 14 5-HT receptor subtypes, their roles in the systems of the body, and, where appropriate, the (potential) utility of therapeutics targeting these receptors. SIGNIFICANCE STATEMENT: This review provides a comprehensive account of the classification and function of 5-hydroxytryptamine receptors, including how they are targeted for therapeutic benefit.

Regulation of a pentameric ligand-gated ion channel by a semiconserved cationic lipid-binding site.

Pentameric ligand-gated ion channels (pLGICs) are crucial mediators of electrochemical signal transduction in various organisms from bacteria to humans. Lipids play an important role in regulating pLGIC function, yet the structural bases for specific pLGIC-lipid interactions remain poorly understood. The bacterial channel ELIC recapitulates several properties of eukaryotic pLGICs, including activation by the neurotransmitter GABA and binding and modulation by lipids, offering a simplified model system for structure-function relationship studies. In this study, functional effects of noncanonical amino acid substitution of a potential lipid-interacting residue (W206) at the top of the M1-helix, combined with detergent interactions observed in recent X-ray structures, are consistent with this region being the location of a lipid-binding site on the outward face of the ELIC transmembrane domain. Coarse-grained and atomistic molecular dynamics simulations revealed preferential binding of lipids containing a positive charge, particularly involving interactions with residue W206, consistent with cation-π binding. Polar contacts from other regions of the protein, particularly M3 residue Q264, further support lipid binding via headgroup ester linkages. Aromatic residues were identified at analogous sites in a handful of eukaryotic family members, including the human GABAA receptor ε subunit, suggesting conservation of relevant interactions in other evolutionary branches. Further mutagenesis experiments indicated that mutations at this site in ε-containing GABAA receptors can change the apparent affinity of the agonist response to GABA, suggesting a potential role of this site in channel gating. In conclusion, this work details type-specific lipid interactions, which adds to our growing understanding of how lipids modulate pLGICs.

Snake venom phospholipase A2s exhibit strong virucidal activity against SARS-CoV-2 and inhibit the viral spike glycoprotein interaction with ACE2.

The COVID-19 pandemic caused by SARS-CoV-2 requires new treatments both to alleviate the symptoms and to prevent the spread of this disease. Previous studies demonstrated good antiviral and virucidal activity of phospholipase A2s (PLA2s) from snake venoms against viruses from different families but there was no data for coronaviruses. Here we show that PLA2s from snake venoms protect Vero E6 cells against SARS-CoV-2 cytopathic effects. PLA2s showed low cytotoxicity to Vero E6 cells with some activity at micromolar concentrations, but strong antiviral activity at nanomolar concentrations. Dimeric PLA2 from the viper Vipera nikolskii and its subunits manifested especially potent virucidal effects, which were related to their phospholipolytic activity, and inhibited cell-cell fusion mediated by the SARS-CoV-2 spike glycoprotein. Moreover, PLA2s interfered with binding both of an antibody against ACE2 and of the receptor-binding domain of the glycoprotein S to 293T/ACE2 cells. This is the first demonstration of a detrimental effect of PLA2s on ß-coronaviruses. Thus, snake PLA2s are promising for the development of antiviral drugs that target the viral envelope, and could also prove to be useful tools to study the interaction of viruses with host cells.

Ecological change in dynamic environments: Accounting for temporal environmental variability in studies of ocean change biology.

The environmental conditions in the ocean have long been considered relatively more stable through time compared to the conditions on land. Advances in sensing technologies, however, are increasingly revealing substantial fluctuations in abiotic factors over ecologically and evolutionarily relevant timescales in the ocean, leading to a growing recognition of the dynamism of the marine environment as well as new questions about how this dynamism may influence species' vulnerability to global environmental change. In some instances, the diurnal or seasonal variability in major environmental change drivers, such as temperature, pH and seawater carbonate chemistry, and dissolved oxygen, can exceed the changes expected with continued anthropogenic global change. While ocean global change biologists have begun to experimentally test how variability in environmental conditions mediates species' responses to changes in the mean, the extensive literature on species' adaptations to temporal variability in their environment and the implications of this variability for their evolutionary responses has not been well integrated into the field. Here, we review the physiological mechanisms underlying species' responses to changes in temperature, pCO2 /pH (and other carbonate parameters), and dissolved oxygen, and discuss what is known about behavioral, plastic, and evolutionary strategies for dealing with variable environments. In addition, we discuss how exposure to variability may influence species' responses to changes in the mean conditions and highlight key research needs for ocean global change biology.

Multiple regions in the extracellular domain of the glycine receptor determine receptor activity.

Glycine receptors (GlyRs) are Cys-loop receptors that mediate fast synaptic inhibition in the brain stem and spinal cord. They are involved in the generation of motor rhythm, reflex circuit coordination, and sensory signal processing and therefore represent targets for therapeutic interventions. The extracellular domains (ECDs) of Cys-loop receptors typically contain many aromatic amino acids, but only those in the receptor binding pocket have been extensively studied. Here, we show that many Phe residues in the ECD that are not located in the binding pocket are also involved in GlyR function. We examined these Phe residues by creating several GlyR variants, characterizing these variants with the two-electrode voltage clamp technique in Xenopus oocytes, and interpreting changes in receptor parameters by using currently available structural information on the open and closed states of the GlyR. Substitution of six of the eight Phe residues in the ECD with Ala resulted in loss of function or significantly increased the EC50 and also altered the maximal response to the partial GlyR agonist taurine compared with glycine in those receptor variants that were functional. Substitutions with other amino acids, combined with examination of nearby residues that could potentially interact with these Phe residues, suggested interactions that could be important for GlyR function, and possibly similar interactions could contribute to the function of other members of the Cys-loop receptor family. Overall, our results suggest that many ECD regions are important for GlyR function and that these regions could inform the design of therapeutic agents targeting GlyR activity.

Probing Proline Residues in the Prokaryotic Ligand-Gated Ion Channel, ELIC.

Erwinia ligand-gated ion channel (ELIC) is a bacterial homologue of vertebrate pentameric ligand-gated ion channels (pLGICs) and has proven to be a valuable model for understanding the structure and function of this important protein family. There is nevertheless still a question about whether molecular details can be accurately extrapolated from this protein to those found in eukaryotes. Here we explore the role of proline residues (Pros) in ELIC by creating mutant receptors, expressing them in Xenopus laevis oocytes, and using whole-cell voltage-clamp electrophysiology to monitor channel activity. In contrast to eukaryotic pLGICs, proline-to-alanine (Pro-to-Ala) substitution in ELIC mostly resulted in gain of function, and even altering highly conserved Pro residues in M1 and the M2-M3 loop did not ablate function. These substitutions also mostly resulted in ablation of the modulation by Ca2+ observed in wild-type receptors. Substitution of the Pro in the "Cys loop", however, did result in nonfunctional receptors. Probing this residue with noncanonical amino acids revealed a requirement for a substituted amine at this position, as well as a general preference for Pro analogues with greater intrinsic cis biases. We propose there is likely a cis bond at the apex of the Cys loop in this protein, which is consistent with some, but not all, findings from other pLGICs. Overall, the data show that the roles of proline residues are less critical in ELIC than in other pLGICs, supporting other studies that suggest caution must be applied in using data from this prokaryotic receptor to understand molecular details of eukaryotic pLGIC receptor function.

Identification of Novel Functionally Important Aromatic Residue Interactions in the Extracellular Domain of the Glycine Receptor.

The extracellular domains (ECDs) of Cys-loop receptors contain many aromatic amino acids, but only relatively few have been well studied. Here we explore the roles of Tyr and Trp residues in the ECD of the glycine receptor and show that four such residues that have not been previously studied (Y24, Y58, W170, and Y197) contribute significantly to the function of the protein. The residues were studied by creating mutant receptors, characterizing them using two-electrode voltage clamp in Xenopus oocytes, and interpreting changes in receptor parameters using structural information about the open and closed states of the receptor. Alanine substitution of all these residues ablates function or increases the glycine EC50. There are also a number of changes in the relative maximal responses to taurine, a partial agonist, compared to glycine. Further mutations, in combination with structural information, suggest Y24 contributes to an anion-π interaction with a binding loop D residue, Y58 to an S-π interaction stabilizing the Cys loop, W170 to hydrophobic interactions stabilizing the hydrophobic interior of the subunit, and Y197 to a hydrogen bond linking binding loops B and C. These interactions appear to be broadly conserved in other Cys-loop receptors. Thus, we have identified new regions of the glycine receptor that are important contributors to receptor activation and are likely also to contribute to function in other members of this important protein family.

Modeling and Mutational Analysis of the Binding Mode for the Multimodal Antidepressant Drug Vortioxetine to the Human 5-HT3A Receptor.

5-Hydroxytryptamine3 (5-HT3) receptors are ligand-gated ion channels that mediate neurotransmission by serotonin in the central nervous system. Pharmacological inhibition of 5-HT3 receptor activity has therapeutic potential in several psychiatric diseases, including depression and anxiety. The recently approved multimodal antidepressant vortioxetine has potent inhibitory activity at 5-HT3 receptors. Vortioxetine has an inhibitory mechanism that differs from classic 5-HT3 receptor competitive antagonists despite being believed to bind in the same binding site. Specifically, vortioxetine shows partial agonist activity followed by persistent and insurmountable inhibition. We have investigated the binding mode of vortioxetine at the human 5-HT3A receptor through computational and in vitro experiments to provide insight into the molecular mechanisms behind the unique pharmacological profile of the drug. We find that vortioxetine binds in a manner different from currently known 5-HT3A orthosteric ligands. Specifically, while the binding pattern of vortioxetine mimics some aspects of both the setron class of competitive antagonists and 5-hydroxytryptamine (5-HT) with regards to interactions with residues of the aromatic box motif in the orthosteric binding site, vortioxetine also forms interactions with residues not previously described to be important for the binding of either setrons or 5-HT such as Val202 on Loop F. Our results expand the framework for understanding how orthosteric ligands drive 5-HT3 receptor function, which is of importance for the potential future development of novel classes of 5-HT3 receptor antagonists.

The roles of aromatic residues in the glycine receptor transmembrane domain.

BACKGROUND: Cys-loop receptors play important roles in fast neuronal signal transmission. Functional receptors are pentamers, with each subunit having an extracellular, transmembrane (TM) and intracellular domain. Each TM domain contains 4 α-helices (M1-M4) joined by loops of varying lengths. Many of the amino acid residues that constitute these α-helices are hydrophobic, and there has been particular interest in aromatic residues, especially those in M4, which have the potential to contribute to the assembly and function of the receptor via a range of interactions with nearby residues. RESULTS: Here we show that many aromatic residues in the M1, M3 and M4 α-helices of the glycine receptor are involved in the function of the receptor. The residues were explored by creating a range of mutant receptors, characterising them using two electrode voltage clamp in Xenopus oocytes, and interpreting changes in receptor parameters using currently available structural information on the open and closed states of the receptor. For 7 residues function was ablated with an Ala substitution: 3 Tyr residues at the extracellular end of M1, 2 Trp residues located towards the centers of M1 and M3, and a Phe and a Tyr residue in M4. For many of these an alternative aromatic residue restored wild-type-like function indicating the importance of the π ring. EC50s were increased with Ala substitution of 8 other aromatic residues, with those in M1 and M4 also having reduced currents, indicating a role in receptor assembly. The structure shows many potential interactions with nearby residues, especially between those that form the M1/M3/M4 interface, and we identify those that are supported by the functional data. CONCLUSION: The data reveal the importance and interactions of aromatic residues in the GlyR M1, M3 and M4 α-helices, many of which are essential for receptor function.

Unexpected resilience of a seagrass system exposed to global stressors.

Despite a growing interest in identifying tipping points in response to environmental change, our understanding of the ecological mechanisms underlying nonlinear ecosystem dynamics is limited. Ecosystems governed by strong species interactions can provide important insight into how nonlinear relationships between organisms and their environment propagate through ecosystems, and the potential for environmentally mediated species interactions to drive or protect against sudden ecosystem shifts. Here, we experimentally determine the functional relationships (i.e., the shapes of the relationships between predictor and response variables) of a seagrass assemblage with well-defined species interactions to ocean acidification (enrichment of CO2 ) in isolation and in combination with nutrient loading. We demonstrate that the effect of ocean acidification on grazer biomass (Phyllaplysia taylori and Idotea resecata) was quadratic, with the peak of grazer biomass at mid-pH levels. Algal grazing was negatively affected by nutrients, potentially due to low grazer affinity for macroalgae (Ulva intestinalis), as recruitment of both macroalgae and diatoms were favored in elevated nutrient conditions. This led to an exponential increase in macroalgal and epiphyte biomass with ocean acidification, regardless of nutrient concentration. When left unchecked, algae can cause declines in seagrass productivity and persistence through shading and competition. Despite quadratic and exponential functional relationships to stressors that could cause a nonlinear decrease in seagrass biomass, productivity of our model seagrass-the eelgrass (Zostera marina)- remained highly resilient to increasing acidification. These results suggest that important species interactions governing ecosystem dynamics may shift with environmental change, and ecosystem state may be decoupled from ecological responses at lower levels of organization.

The Proton Responsiveness in the Extracellular Domain of GLIC Differs in the Presence of the ELIC Transmembrane Domain.

Prokaryotic homologues of Cys-loop receptors have proven to be useful in understanding their eukaryotic counterparts, but even the best studied of these, Gloeobacter ligand-gated ion channel (GLIC), is still not yet fully understood. GLIC is activated by protons with a pH50 between 5 and 6, implicating a histidine residue in its activation, but although a histidine residue (His11') in the pore-forming α-helix (M2) is known to be involved in gating, the His in the extracellular domain (ECD), His127, is not. Nevertheless, there is evidence from a GLIC-glycine chimera for a proton sensitive residue or region in the GLIC extracellular domain. Here we create a novel chimeric receptor with the ECD of GLIC and the transmembrane domain of ELIC (GELIC). Expression of this receptor in oocytes reveals proton activation, although the pH50 (6.7) differs from that of GLIC (5.4). Exploration of protonatable residues in the ECD reveals that the pKas of five Asp residues (31, 49, 91, 136, and 178) differ between the open and closed states of GLIC. Substitution of these residues with Ala or Asn shows somewhat similar effects for GLIC and GELIC in Asp91 mutants, but different effects for the others. Overall, the data suggest that protonation of residues in the ECD is a requirement for channel opening in GELIC but plays only a minor role in GLIC, where gating may be largely driven via protonation of the His residue in its pore.

Perturbation of Critical Prolines in Gloeobacter violaceus Ligand-gated Ion Channel (GLIC) Supports Conserved Gating Motions among Cys-loop Receptors.

Gloeobacter violaceus ligand-gated ion channel (GLIC) has served as a valuable structural and functional model for the eukaryotic Cys-loop receptor superfamily. In Cys-loop and other receptors, we have previously demonstrated the crucial roles played by several conserved prolines. Here we explore the role of prolines in the gating transitions of GLIC. As conventional substitutions at some positions resulted in nonfunctional proteins, we used in vivo non-canonical amino acid mutagenesis to determine the specific structural requirements at these sites. Receptors were expressed heterologously in Xenopus laevis oocytes, and whole-cell electrophysiology was used to monitor channel activity. Pro-119 in the Cys-loop, Pro-198 and Pro-203 in the M1 helix, and Pro-299 in the M4 helix were sensitive to substitution, and distinct roles in receptor activity were revealed for each. In the context of the available structural data for GLIC, the behaviors of Pro-119, Pro-203, and Pro-299 mutants are consistent with earlier proline mutagenesis work. However, the Pro-198 site displays a unique phenotype that gives evidence of the importance of the region surrounding this residue for the correct functioning of GLIC.

Disturbed neuronal ER-Golgi sorting of unassembled glycine receptors suggests altered subcellular processing is a cause of human hyperekplexia.

Recent studies on the pathogenic mechanisms of recessive hyperekplexia indicate disturbances in glycine receptor (GlyR) α1 biogenesis. Here, we examine the properties of a range of novel glycine receptor mutants identified in human hyperekplexia patients using expression in transfected cell lines and primary neurons. All of the novel mutants localized in the large extracellular domain of the GlyR α1 have reduced cell surface expression with a high proportion of receptors being retained in the ER, although there is forward trafficking of glycosylated subpopulations into the ER-Golgi intermediate compartment and cis-Golgi compartment. CD spectroscopy revealed that the mutant receptors have proportions of secondary structural elements similar to wild-type receptors. Two mutants in loop B (G160R, T162M) were functional, but none of those in loop D/ß2-3 were. One nonfunctional truncated mutant (R316X) could be rescued by coexpression with the lacking C-terminal domain. We conclude that a proportion of GlyR α1 mutants can be transported to the plasma membrane but do not necessarily form functional ion channels. We suggest that loop D/ß2-3 is an important determinant for GlyR trafficking and functionality, whereas alterations to loop B alter agonist potencies, indicating that residues here are critical elements in ligand binding.

Crotonic Acid Blocks the Gloeobacter Ligand-Gated Ion Channel (GLIC) via the Extracellular Domain.

Cys-loop receptors play important roles in signal transduction in multicellular organisms, but similar proteins exist in prokaryotes, the best studied of which is the Gloeobacter ligand-gated ion channel (GLIC). GLIC is activated by protons with 50% activation (pH50) at pH 5.5, and while a histidine residue in its pore-forming α-helix (M2) is known to be involved in gating, there is also evidence of a proton-sensitive region in the extracellular domain. However, this proton-sensitive region does not appear to be located in the region of GLIC equivalent to the agonist binding site in related proteins. Here we explore functional effects of a range of compounds that could bind to this site and show that some GABA analogues, the most potent of which is crotonic acid, inhibit GLIC function. Mutagenesis and docking studies suggest crotonic acid can bind to this region of the protein and, when bound, can allosterically inhibit GLIC function. These data therefore suggest that there is a transduction pathway from the orthosteric binding site to the pore in GLIC, as exists in related eukaryotic ligand-gated ion channels, and thus provide further evidence that this prokaryotic receptor is a good model for understanding this family of proteins.

Probing residues in the pore-forming (M2) domain of the Cys-loop receptor homologue GLIC reveals some unusual features.

Cys-loop receptors play important roles in signal transduction. The Gloeobacter ligand-gated ion channel (GLIC) pore binds similar compounds to Cys-loop receptor pores, but has the advantage of known structures in open and closed states. GLIC is activated by protons with a pEC50 of 5.4, and has a histidine residue (His 11') in its pore-forming α-helix (M2) which is involved in gating. Here we explore the role of this His and other M2 residues using two-electrode voltage clamp of mutant receptors expressed in oocytes. We show that 11'His is very sensitive to substitution; replacement with a range of amino acids ablates function. Similarly altering its location in M2 to the 8', 9', 10', 12', 13' or 14' positions ablated function. Most substitutions of Ser6' or Ile9' were also non-functional, although not Ile9'Leu and Ile9'Val. Unexpectedly, an Ile9'His substitution was constitutively active at pH 7, but closed as [H+] increased, with a pIC50 of 5.8. Substitution at 2', 5' and 7' had little effect on pEC50. Overall the data show Ser6' and His11' are critical for the function of the receptor, and thus distinguish the roles of these M2 residues from those of Cys-loop receptors, where substitutions are mostly well tolerated. These data suggest modellers should be aware of these atypical features when using the GLIC pore as a model for Cys-loop receptor pores.

Structural basis of ligand recognition in 5-HT3 receptors.

The 5-HT(3) receptor is a pentameric serotonin-gated ion channel, which mediates rapid excitatory neurotransmission and is the target of a therapeutically important class of anti-emetic drugs, such as granisetron. We report crystal structures of a binding protein engineered to recognize the agonist serotonin and the antagonist granisetron with affinities comparable to the 5-HT(3) receptor. In the serotonin-bound structure, we observe hydrophilic interactions with loop E-binding site residues, which might enable transitions to channel opening. In the granisetron-bound structure, we observe a critical cation-π interaction between the indazole moiety of the ligand and a cationic centre in loop D, which is uniquely present in the 5-HT(3) receptor. We use a series of chemically tuned granisetron analogues to demonstrate the energetic contribution of this electrostatic interaction to high-affinity ligand binding in the human 5-HT(3) receptor. Our study offers the first structural perspective on recognition of serotonin and antagonism by anti-emetics in the 5-HT(3) receptor.

Pentameric ligand-gated ion channel ELIC is activated by GABA and modulated by benzodiazepines.

GABA(A) receptors are pentameric ligand-gated ion channels involved in fast inhibitory neurotransmission and are allosterically modulated by the anxiolytic, anticonvulsant, and sedative-hypnotic benzodiazepines. Here we show that the prokaryotic homolog ELIC also is activated by GABA and is modulated by benzodiazepines with effects comparable to those at GABA(A) receptors. Crystal structures reveal important features of GABA recognition and indicate that benzodiazepines, depending on their concentration, occupy two possible sites in ELIC. An intrasubunit site is adjacent to the GABA-recognition site but faces the channel vestibule. A second intersubunit site partially overlaps with the GABA site and likely corresponds to a low-affinity benzodiazepine-binding site in GABA(A) receptors that mediates inhibitory effects of the benzodiazepine flurazepam. Our study offers a structural view how GABA and benzodiazepines are recognized at a GABA-activated ion channel.

Phenylalanine in the pore of the Erwinia ligand-gated ion channel modulates picrotoxinin potency but not receptor function.

The Erwinia ligand-gated ion channel (ELIC) is a bacterial homologue of eukaryotic Cys-loop ligand-gated ion channels. This protein has the potential to be a useful model for Cys-loop receptors but is unusual in that it has an aromatic residue (Phe) facing into the pore, leading to some predictions that this protein is incapable of ion flux. Subsequent studies have shown this is not the case, so here we probe the role of this residue by examining the function of the ELIC in cases in which the Phe has been substituted with a range of alternative amino acids, expressed in Xenopus oocytes and functionally examined. Most of the mutations have little effect on the GABA EC50, but the potency of the weak pore-blocking antagonist picrotoxinin at F16'A-, F16'D-, F16'S-, and F16'T-containing receptors was increased to levels comparable with those of Cys-loop receptors, suggesting that this antagonist can enter the pore only when residue 16' is small. T6'S has no effect on picrotoxinin potency when expressed alone but abolishes the increased potency when combined with F16'S, indicating that the inhibitor binds at position 6', as in Cys-loop receptors, if it can enter the pore. Overall, the data support the proposal that the ELIC pore is a good model for Cys-loop receptor pores if the role of F16' is taken into consideration.

Multisite binding of a general anesthetic to the prokaryotic pentameric Erwinia chrysanthemi ligand-gated ion channel (ELIC).

Pentameric ligand-gated ion channels (pLGICs), such as nicotinic acetylcholine, glycine, γ-aminobutyric acid GABA(A/C) receptors, and the Gloeobacter violaceus ligand-gated ion channel (GLIC), are receptors that contain multiple allosteric binding sites for a variety of therapeutics, including general anesthetics. Here, we report the x-ray crystal structure of the Erwinia chrysanthemi ligand-gated ion channel (ELIC) in complex with a derivative of chloroform, which reveals important features of anesthetic recognition, involving multiple binding at three different sites. One site is located in the channel pore and equates with a noncompetitive inhibitor site found in many pLGICs. A second transmembrane site is novel and is located in the lower part of the transmembrane domain, at an interface formed between adjacent subunits. A third site is also novel and is located in the extracellular domain in a hydrophobic pocket between the ß7-ß10 strands. Together, these results extend our understanding of pLGIC modulation and reveal several specific binding interactions that may contribute to modulator recognition, further substantiating a multisite model of allosteric modulation in this family of ion channels.

Insights into the binding of GABA to the insect RDL receptor from atomistic simulations: a comparison of models.

The resistance to dieldrin (RDL) receptor is an insect pentameric ligand-gated ion channel (pLGIC). It is activated by the neurotransmitter γ-aminobutyric acid (GABA) binding to its extracellular domain; hence elucidating the atomistic details of this interaction is important for understanding how the RDL receptor functions. As no high resolution structures are currently available, we built homology models of the extracellular domain of the RDL receptor using different templates, including the widely used acetylcholine binding protein and two pLGICs, the Erwinia Chrysanthemi ligand-gated ion channel (ELIC) and the more recently resolved GluCl. We then docked GABA into the selected three dimensional structures, which we used as starting points for classical molecular dynamics simulations. This allowed us to analyze in detail the behavior of GABA in the binding sites, including the hydrogen bond and cation-π interaction networks it formed, the conformers it visited and the possible role of water molecules in mediating the interactions; we also estimated the binding free energies. The models were all stable and showed common features, including interactions consistent with experimental data and similar to other pLGICs; differences could be attributed to the quality of the models, which increases with increasing sequence identity, and the use of a pLGIC template. We supplemented the molecular dynamics information with metadynamics, a rare event method, by exploring the free energy landscape of GABA binding to the RDL receptor. Overall, we show that the GluCl template provided the best models. GABA forming direct salt-bridges with Arg211 and Glu204, and cation-π interactions with an aromatic cage including Tyr109, Phe206 and Tyr254, represents a favorable binding arrangement, and the interaction with Glu204 can also be mediated by a water molecule.