Targeted molecular dynamics study of C-loop closure and channel gating in nicotinic receptors

Molecular Dynamics Refinement of Open State Serotonin 5-HT3A Receptor Structures

Pentameric ligand-gated ion channels play an important role in mediating fast neurotransmissions. As a member of this receptor family, cation-selective 5-HT₃ receptors are a clinical target for treating nausea and vomiting associated with chemotherapy and radiation therapy (Thompson and Lummis, 2006). Multiple cryo-electron microscopy (cryo-EM) structures of 5-HT₃ receptors have been determined in distinct functional states (e.g., open, closed, etc.) (Basak et al., 2018; Basak et al., 2018; Polovinkin et al., 2018; Zhang et al., 2015). However, recent work has shown that the transmembrane pores of the open 5-HT₃ receptor structures rapidly collapse and become artificially asymmetric in molecular dynamics (MD) simulations. To avoid this hydrophobic collapse, Dämgen and Biggin developed an equilibration protocol that led to a stable open state structure of the glycine receptor in MD simulations (Dämgen and Biggin, 2020). However, the protocol failed to yield open-like structures of the 5-HT₃ receptor in our simulations. Here, we present a refined equilibration protocol that involves the rearrangement of the transmembrane helices to achieve stable open state structures of the 5-HT₃ receptor that allow both water and ion permeation through the channel. Notably, channel gating is mediated through collective movement of the transmembrane helices, involving not only pore lining M2 helices but also their cross-talk with the adjacent M1 and M3 helices. Thus, the successful application of our refined equilibration protocol underscores the importance of the conformational coupling between the transmembrane helices in stabilizing open-like structures of the 5-HT₃ receptor.

Predicting the structural basis of targeted protein degradation by integrating molecular dynamics simulations with structural mass spectrometry

Targeted protein degradation (TPD) is a promising approach in drug discovery for degrading proteins implicated in diseases. A key step in this process is the formation of a ternary complex where a heterobifunctional molecule induces proximity of an E3 ligase to a protein of interest (POI), thus facilitating ubiquitin transfer to the POI. In this work, we characterize 3 steps in the TPD process. (1) We simulate the ternary complex formation of SMARCA2 bromodomain and VHL E3 ligase by combining hydrogen-deuterium exchange mass spectrometry with weighted ensemble molecular dynamics (MD). (2) We characterize the conformational heterogeneity of the ternary complex using Hamiltonian replica exchange simulations and small-angle X-ray scattering. (3) We assess the ubiquitination of the POI in the context of the full Cullin-RING Ligase, confirming experimental ubiquitinomics results. Differences in degradation efficiency can be explained by the proximity of lysine residues on the POI relative to ubiquitin.

Elephants in the Dark: Insights and Incongruities in Pentameric Ligand-gated Ion Channel Models

The superfamily of pentameric ligand-gated ion channels (pLGICs) comprises key players in electrochemical signal transduction across evolution, including historic model systems for receptor allostery and targets for drug development. Accordingly, structural studies of these channels have steadily increased, and now approach 250 depositions in the protein data bank. This review contextualizes currently available structures in the pLGIC family, focusing on morphology, ligand binding, and gating in three model subfamilies: the prokaryotic channel GLIC, the cation-selective nicotinic acetylcholine receptor, and the anion-selective glycine receptor. Common themes include the challenging process of capturing and annotating channels in distinct functional states; partially conserved gating mechanisms, including remodeling at the extracellular/transmembrane-domain interface; and diversity beyond the protein level, arising from posttranslational modifications, ligands, lipids, and signaling partners. Interpreting pLGIC structures can be compared to describing an elephant in the dark, relying on touch alone to comprehend the many parts of a monumental beast: each structure represents a snapshot in time under specific experimental conditions, which must be integrated with further structure, function, and simulations data to build a comprehensive model, and understand how one channel may fundamentally differ from another.

Atomic-Resolution 1.3 Å Crystal Structure, Inhibition by Sulfate, and Molecular Dynamics of the Bacterial Enzyme DapE

We report the atomic-resolution (1.3 Å) X-ray crystal structure of an open conformation of the dapE-encoded N-succinyl-l,l-diaminopimelic acid desuccinylase (DapE, EC 3.5.1.18) from Neisseria meningitidis. This structure [Protein Data Bank (PDB) entry 5UEJ] contains two bound sulfate ions in the active site that mimic the binding of the terminal carboxylates of the N-succinyl-l,l-diaminopimelic acid (l,l-SDAP) substrate. We demonstrated inhibition of DapE by sulfate (IC₅₀ = 13.8 ± 2.8 mM). Comparison with other DapE structures in the PDB demonstrates the flexibility of the interdomain connections of this protein. This high-resolution structure was then utilized as the starting point for targeted molecular dynamics experiments revealing the conformational change from the open form to the closed form that occurs when DapE binds l,l-SDAP and cleaves the amide bond. These simulations demonstrated closure from the open to the closed conformation, the change in RMS throughout the closure, and the independence in the movement of the two DapE subunits. This conformational change occurred in two phases with the catalytic domains moving toward the dimerization domains first, followed by a rotation of catalytic domains relative to the dimerization domains. Although there were no targeting forces, the substrate moved closer to the active site and bound more tightly during the closure event.

Nicotinic receptor pharmacology in silico: Insights and challenges

Nicotinic acetylcholine receptors (nAChR) are homo- or hetero-pentameric ligand-gated ion channels of the Cys-loop superfamily and play important roles in the nervous system and muscles. Studies on nAChR benefit from in silico modeling due to the lack of high-resolution structures for most receptor subtypes and challenges in experiments addressing the complex mechanism of activation involving allosteric sites. Although there is myriad of computational modeling studies on nAChR, the multitude of the methods and parameters used in these studies makes modeling nAChR a daunting task, particularly for the non-experts in the field. To address this problem, the modeling literature on Torpedo nAChR and α7 nAChR were focused on as examples of heteromeric and homomeric nAChR, and the key in silico modeling studies between the years 1995-2019 were concisely reviewed. This was followed by a critical analysis of these studies by comparing the findings with each other and with the emerging experimental and computational data on nAChR. Based on these critical analyses, suggestions were made to guide the future researchers in the field of in silico modeling of nAChR. This article is part of the special issue on 'Contemporary Advances in Nicotine Neuropharmacology'.

Dynamics and Interactions of GPI-Linked lynx1 Protein with/without Nicotinic Acetylcholine Receptor in Membrane Bilayers

Nicotinic acetylcholine receptors (nAChRs) participate in diverse biological processes, such as mood, learning, and addiction. Glycosylphosphatidylinositol-linked lynx1 is an allosteric modulator of nAChR function, including shifts in agonist sensitivity, reduced desensitization, and slower recovery from desensitization. This modulation is thought to be achieved by lynx1's interaction with nAChR subunits, particularly at the α4:α4 interface. In this study, we used molecular modeling and simulation to study the structure, dynamics, and interactions of lynx1 when bound to nAChRs, as well as unbound, monomeric lynx1 in membranes. Though lynx1 structures are similar in both states, its dynamics is more restricted in the bound state than in the unbound one. When bound, interactions between lynx1 and nAChR are observed to be maintained throughout the simulations. Of particular note, lynx1 demonstrates prolonged interactions with the receptor C-loop in one of the nAChR α4 subunits, a region important for agonist binding and possibly the transition between open/closed states. During interactions with lynx1, an α4 C-loop tends to be restricted in either a closed or open state, whereas the C-loop state transitions are more evident when lynx1 is unbound. Interestingly, the conformational change of the C-loop is stochastic, suggesting that lynx1 can influence nAChR (critical for its multimodal action), for instance, by shifting its agonist sensitivity and recovery from desensitization.

In Silico Modeling of the α7 Nicotinic Acetylcholine Receptor: New Pharmacological Challenges Associated with Multiple Modes of Signaling

The α7 nicotinic acetylcholine receptor is a homopentameric ion-channel of the Cys-loop superfamily characterized by its low probability of opening, high calcium permeability, and rapid desensitization. The α7 receptor has been targeted for the treatment of the cognitive symptoms of schizophrenia, depression, and Alzheimer's disease, but it is also involved in inflammatory modulation as a part of the cholinergic anti-inflammatory pathway. Despite its functional importance, in silico studies of the α7 receptor cannot produce a general model explaining the structural features of receptor activation, nor predict the mode of action for various ligand classes. Two particular problems in modeling the α7 nAChR are the absence of a high-resolution structure and the presence of five potentially nonequivalent orthosteric ligand binding sites. There is wide variability regarding the templates used for homology modeling, types of ligands investigated, simulation methods, and simulation times. However, a systematic survey focusing on the methodological similarities and differences in modeling α7 has not been done. In this work, we make a critical analysis of the modeling literature of α7 nAChR by comparing the findings of computational studies with each other and with experimental studies under the main topics of structural studies, ligand binding studies, and comparisons with other nAChR. In light of our findings, we also summarize current problems in the field and make suggestions for future studies concerning modeling of the α7 receptor.

The C loop at the orthosteric binding site is critically involved in GABAA receptor gating

GABA_A receptors (GABA_ARs) play a crucial role in mammalian adult brain inhibition. The dysfunction of GABAergic drive is related to such disorders as epilepsy, schizophrenia, and depression. Substantial progress has recently been made in describing the static structure of GABA_ARs, but the molecular mechanisms that underlie the activation process remain elusive. The C loop of the GABA_AR structure shows the largest movement upon ligand binding to the orthosteric binding site, a phenomenon that is referred to as "capping." The C loop is known to be involved in agonist binding, but its role in the gating of Cys-loop receptors is still debated. Herein, we investigated this issue by analyzing the impact of a β₂F200 residue mutation of the C loop on gating properties of α₁β₂γ₂ GABA_ARs. Extensive analyses and the modeling of current responses to saturating agonist application demonstrated that this mutation strongly affected preactivation, opening, closing and desensitization, i.e. all considered gating steps. Single-channel analysis revealed that the β₂F200 mutation slowed all shut time components, and open times were shortened. Model fitting of these single-channel data further confirmed that the β₂F200 mutation strongly affected all of the gating characteristics. We also found that this mutation altered receptor sensitivity to the benzodiazepine flurazepam, which was attributable to a change in preactivation kinetics. In silico analysis indicated that the β₂F200 mutation resulted in distortion of the C loop structure, causing the movement of its tip from the binding site. Altogether, we provide the first evidence that C loop critically controls GABA_AR gating.

Blockade of Human α7 Nicotinic Acetylcholine Receptor by α-Conotoxin ImI Dendrimer: Insight from Computational Simulations

Nicotinic acetylcholine receptors (nAChRs) are ligand-gated ion channels that are involved in fast synaptic transmission and mediated physiological activities in the nervous system. α-Conotoxin ImI exhibits subtype-specific blockade towards homomeric α7 and α9 receptors. In this study, we established a method to build a 2×ImI-dendrimer/h (human) α7 nAChR model, and based on this model, we systematically investigated the molecular interactions between the 2×ImI-dendrimer and hα7 nAChR. Our results suggest that the 2×ImI-dendrimer possessed much stronger potency towards hα7 nAChR than the α-ImI monomer and demonstrated that the linker between α-ImI contributed to the potency of the 2×ImI-dendrimer by forming a stable hydrogen-bond network with hα7 nAChR. Overall, this study provides novel insights into the binding mechanism of α-ImI dendrimer to hα7 nAChR, and the methodology reported here opens an avenue for the design of more selective dendrimers with potential usage as drug/gene carriers, macromolecular drugs, and molecular probes.

A single molecular distance predicts agonist binding energy in nicotinic receptors

Agonists turn on receptors because they bind more strongly to active (R*) versus resting (R) conformations of their target sites. Here, to explore how agonists activate neuromuscular acetylcholine receptors, we built homology models of R and R* neurotransmitter binding sites, docked ligands to those sites, ran molecular dynamics simulations to relax ("equilibrate") the structures, measured binding site structural parameters, and correlated them with experimental agonist binding energies. Each binding pocket is a pyramid formed by five aromatic amino acids and covered partially by loop C. We found that in R* versus R, loop C is displaced outward, the pocket is smaller and skewed, the agonist orientation is reversed, and a key nitrogen atom in the agonist is closer to the pocket center (distance dx) and a tryptophan pair but farther from αY190. Of these differences, the change in dx shows the largest correlation with experimental binding energy and provides a good estimate of agonist affinity, efficacy, and efficiency. Indeed, concentration-response curves can be calculated from just dx values. The contraction and twist of the binding pocket upon activation resemble gating rearrangements of the extracellular domain of related receptors at a smaller scale.

Closed-Locked and Apo-Resting State Structures of the Human α7 Nicotinic Receptor: A Computational Study

Nicotinic acetylcholine receptors, belonging to the Cys-loop superfamily of ligand-gated ion channels (LGICs), are membrane proteins present in neurons and at neuromuscular junctions. They are responsible for signal transmission, and their function is regulated by neurotransmitters, agonists, and antagonists drugs. A detailed knowledge of their conformational transition in response to ligand binding is critical to understanding the basis of ligand-receptor interaction, in view of new pharmacological approaches to control receptor activity. However, the scarcity of experimentally derived structures of human channels makes this perspective extremely challenging. To contribute overcoming this issue, we have recently reported structural models for the open and the desensitized states of the human α7 nicotinic receptor. Here, we provide all-atom structural models of the same receptor in two different nonconductive states. The first structure, built via homology modeling and relaxed with extensive Molecular Dynamics simulations, represents the receptor bound to the natural antagonist α-conotoxin ImI. After comparison with available experimental data and computational models of other eukaryotic LGICs, we deem it consistent with the "closed-locked" state. The second model, obtained with simulations from the spontaneous relaxation of the open, agonist-bound α7 structure after ligand removal, recapitulates the characteristics of the apo-resting state of the receptor. These results add to our previous work on the active and desensitized state conformations, contributing to the structural characterization of the conformational landscape of the human α7 receptor and suggesting benchmarks to discriminate among conformations found in experiments or in simulations of LGICs. In particular key interactions at the interface between the extracellular domain and the transmembrane domain are identified, that could be critical to the α7 receptor function.

The neural γ2α1β2α1β2 gamma amino butyric acid ion channel receptor: structural analysis of the effects of the ivermectin molecule and disulfide bridges

While ~30% of the human genome encodes membrane proteins, only a handful of structures of membrane proteins have been resolved to high resolution. Here, we studied the structure of a member of the Cys-loop ligand gated ion channel protein superfamily of receptors, human type A γ₂α₁β₂α₁β₂ gamma amino butyric acid receptor complex in a lipid bilayer environment. Studying the correlation between the structure and function of the gamma amino butyric acid receptor may enhance our understanding of the molecular basis of ion channel dysfunctions linked with epilepsy, ataxia, migraine, schizophrenia and other neurodegenerative diseases. The structure of human γ₂α₁β₂α₁β₂ has been modeled based on the X-ray structure of the Caenorhabditis elegans glutamate-gated chloride channel via homology modeling. The template provided the first inhibitory channel structure for the Cys-loop superfamily of ligand-gated ion channels. The only available template structure before this glutamate-gated chloride channel was a cation selective channel which had very low sequence identity with gamma aminobutyric acid receptor. Here, our aim was to study the effect of structural corrections originating from modeling on a more reliable template structure. The homology model was analyzed for structural properties via a 100 ns molecular dynamics (MD) study. Due to the structural shifts and the removal of an open channel potentiator molecule, ivermectin, from the template structure, helical packing changes were observed in the transmembrane segment. Namely removal of ivermectin molecule caused a closure around the Leu 9 position along the ion channel. In terms of the structural shifts, there are three potential disulfide bridges between the M1 and M3 helices of the γ₂ and 2 α₁ subunits in the model. The effect of these disulfide bridges was investigated via monitoring the differences in root mean square fluctuations (RMSF) of individual amino acids and principal component analysis of the MD trajectory of the two homology models-one with the disulfide bridge and one with protonated Cys residues. In all subunit types, RMSF of the transmembrane domain helices are reduced in the presence of disulfide bridges. Additionally, loop A, loop F and loop C fluctuations were affected in the extracellular domain. In cross-correlation analysis of the trajectory, the two model structures displayed different coupling in between the M2-M3 linker region, protruding from the membrane, and the β1-β2/D loop and cys-loop regions in the extracellular domain. Correlations of the C loop, which collapses directly over the bound ligand molecule, were also affected by differences in the packing of transmembrane helices. Finally, more localized correlations were observed in the transmembrane helices when disulfide bridges were present in the model. The differences observed in this study suggest that dynamic coupling at the interface of extracellular and ion channel domains differs from the coupling introduced by disulfide bridges in the transmembrane region. We hope that this hypothesis will be tested experimentally in the near future.

Molecular mechanism of GPCR-mediated arrestin activation

Despite intense interest in discovering drugs that cause G-protein-coupled receptors (GPCRs) to selectively stimulate or block arrestin signalling, the structural mechanism of receptor-mediated arrestin activation remains unclear^1,2. Here we reveal this mechanism through extensive atomic-level simulations of arrestin. We find that the receptor's transmembrane core and cytoplasmic tail-which bind distinct surfaces on arrestin-can each independently stimulate arrestin activation. We confirm this unanticipated role of the receptor core, and the allosteric coupling between these distant surfaces of arrestin, using site-directed fluorescence spectroscopy. The effect of the receptor core on arrestin conformation is mediated primarily by interactions of the intracellular loops of the receptor with the arrestin body, rather than the marked finger-loop rearrangement that is observed upon receptor binding. In the absence of a receptor, arrestin frequently adopts active conformations when its own C-terminal tail is disengaged, which may explain why certain arrestins remain active long after receptor dissociation. Our results, which suggest that diverse receptor binding modes can activate arrestin, provide a structural foundation for the design of functionally selective ('biased') GPCR-targeted ligands with desired effects on arrestin signalling.

SPIDR: small-molecule peptide-influenced drug repurposing

Background

Conventional de novo drug design is costly and time consuming, making it accessible to only the best resourced research organizations. An emergent approach to new drug development is drug repurposing, in which compounds that have already gone through some level of clinical testing are examined for efficacy against diseases divergent than their original application. Repurposing of existing drugs circumvents the time and considerable cost of early stages of drug development, and can be accelerated by using software to screen existing chemical databases to identify suitable drug candidates.

Results

Small-molecule Peptide-Influenced Drug Repurposing (SPIDR) was developed to identify small molecule drugs that target a specific receptor by exploring the conformational binding space of peptide ligands. SPIDR was tested using the potent and selective 16-amino acid peptide α-conotoxin MII ligand and the α₃β₂-nicotinic acetylcholine receptor (nAChR) isoform. SPIDR incorporates a genetic algorithm-based, heuristic search procedure, which was used to explore the ligand binding domain of the α₃β₂-nAChR isoform using a library consisting of 640,000 α-conotoxin MII peptide analogs. The peptides that exhibited the highest affinity for α₃β₂-nAChR were used as models for a small-molecule structure similarity search of the PubChem Compound database. SPIDR incorporates the SimSearcher utility, which generates shape distribution signatures of molecules and employs multi-level K-means clustering to insure fast database queries. SPIDR identified non-peptide drugs with estimated binding affinities nearly double that of the native α-conotoxin MII peptide.

Conclusions

SPIDR has been generalized and integrated into DockoMatic v 2.1. This software contains an intuitive graphical interface for peptide mutant screening workflow and facilitates mapping, clustering, and searching of local molecular databases, making DockoMatic a valuable tool for researchers in drug design and repurposing.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2153-y) contains supplementary material, which is available to authorized users.

Competitive docking model for prediction of the human nicotinic acetylcholine receptor α7 binding of tobacco constituents

The detrimental health effects associated with tobacco use constitute a major public health concern. The addiction associated with nicotine found in tobacco products has led to difficulty in quitting among users. Nicotinic acetylcholine receptors (nAChRs) are the targets of nicotine and are responsible for addiction to tobacco products. However, it is unknown if the other >8000 tobacco constituents are addictive. Since it is time-consuming and costly to experimentally assess addictive potential of such larger number of chemicals, computationally predicting human nAChRs binding is important for in silico evaluation of addiction potential of tobacco constituents and needs structures of human nAChRs. Therefore, we constructed three-dimensional structures of the ligand binding domain of human nAChR α7 subtype and then developed a predictive model based on the constructed structures to predict human nAChR α7 binding activity of tobacco constituents. The predictive model correctly predicted 11 out of 12 test compounds to be binders of nAChR α7. The model is a useful tool for high-throughput screening of potential addictive tobacco constituents. These results could inform regulatory science research by providing a new validated predictive tool using cutting-edge computational methodology to high-throughput screen tobacco additives and constituents for their binding interaction with the human α7 nicotinic receptor. The tool represents a prediction model capable of screening thousands of chemicals found in tobacco products for addiction potential, which improves the understanding of the potential effects of additives.

Single-molecule techniques in biophysics: a review of the progress in methods and applications

Single-molecule biophysics has transformed our understanding of biology, but also of the physics of life. More exotic than simple soft matter, biomatter lives far from thermal equilibrium, covering multiple lengths from the nanoscale of single molecules to up to several orders of magnitude higher in cells, tissues and organisms. Biomolecules are often characterized by underlying instability: multiple metastable free energy states exist, separated by levels of just a few multiples of the thermal energy scale k _B T, where k _B is the Boltzmann constant and T absolute temperature, implying complex inter-conversion kinetics in the relatively hot, wet environment of active biological matter. A key benefit of single-molecule biophysics techniques is their ability to probe heterogeneity of free energy states across a molecular population, too challenging in general for conventional ensemble average approaches. Parallel developments in experimental and computational techniques have catalysed the birth of multiplexed, correlative techniques to tackle previously intractable biological questions. Experimentally, progress has been driven by improvements in sensitivity and speed of detectors, and the stability and efficiency of light sources, probes and microfluidics. We discuss the motivation and requirements for these recent experiments, including the underpinning mathematics. These methods are broadly divided into tools which detect molecules and those which manipulate them. For the former we discuss the progress of super-resolution microscopy, transformative for addressing many longstanding questions in the life sciences, and for the latter we include progress in 'force spectroscopy' techniques that mechanically perturb molecules. We also consider in silico progress of single-molecule computational physics, and how simulation and experimentation may be drawn together to give a more complete understanding. Increasingly, combinatorial techniques are now used, including correlative atomic force microscopy and fluorescence imaging, to probe questions closer to native physiological behaviour. We identify the trade-offs, limitations and applications of these techniques, and discuss exciting new directions.

Deciphering Intrinsic Inter-subunit Couplings that Lead to Sequential Hydrolysis of F1-ATPase Ring

Rotary sequential hydrolysis of the metabolic machine F₁-ATPase is a prominent manifestation of high coordination among multiple chemical sites in ring-shaped molecular machines, and it is also functionally essential for F₁ to tightly couple chemical reactions and central γ-shaft rotation. High-speed AFM experiments have identified that sequential hydrolysis is maintained in the F₁ stator ring even in the absence of the γ-rotor. To explore the origins of intrinsic sequential performance, we computationally investigated essential inter-subunit couplings on the hexameric ring of mitochondrial and bacterial F₁. We first reproduced in stochastic Monte Carlo simulations the experimentally determined sequential hydrolysis schemes by kinetically imposing inter-subunit couplings and following subsequent tri-site ATP hydrolysis cycles on the F₁ ring. We found that the key couplings to support the sequential hydrolysis are those that accelerate neighbor-site ADP and Pi release upon a certain ATP binding or hydrolysis reaction. The kinetically identified couplings were then examined in atomistic molecular dynamics simulations at a coarse-grained level to reveal the underlying structural mechanisms. To do that, we enforced targeted conformational changes of ATP binding or hydrolysis to one chemical site on the F₁ ring and monitored the ensuing conformational responses of the neighboring sites using structure-based simulations. Notably, we found asymmetrical neighbor-site opening that facilitates ADP release upon enforced ATP binding. We also captured a complete charge-hopping process of the Pi release subsequent to enforced ATP hydrolysis in the neighbor site, confirming recent single-molecule analyses with regard to the role of ATP hydrolysis in F₁. Our studies therefore elucidate both the coordinated chemical kinetics and structural dynamics mechanisms underpinning the sequential operation of the F₁ ring.

Exploring the Relationship between Nicotinic Acetylcholine Receptor Ligand Size, Efficiency, Efficacy, and C-Loop Opening

Nicotinic acetylcholine receptors (nAChRs) are ligand-gated ion channels mediating fundamental physiological activities in the nervous system and have become important targets for drug design. For a long time, the acetylcholine binding protein (AChBP) has been used as a surrogate to study the nAChR structure-function. Taking advantage of more than 100 AChBP crystal structures in the Protein DataBank (PDB), we explored the relationship between the size, efficiency, and efficacy of nAChR ligands and the C-loop movement. We found that the size of the ligand is correlated with the opening of the C-loop, which can be used in selecting AChBP crystal structures with appropriate C-loop opening to be used for nAChR ligand docking. Ligand size and C-loop opening are reversely correlated with the ligand efficiency rather than the binding affinity. Ligand efficiency could be accurately predicted using simple computational docking, giving a correlation coefficients (R²) up to 0.73. The efficacy of nAChR ligands might be related to ligand size, C-loop opening, and ligand efficiency. Results from this study are useful for engineering the binding affinity and efficacy of nAChR ligands.

Exploring the binding energy profiles of full agonists, partial agonists, and antagonists of the α7 nicotinic acetylcholine receptor

Nicotinic acetylcholine receptors (nAChRs) belong to the Cys-loop receptor family and are important drug targets for the treatment of neurological diseases. However, the precise determinants of the binding efficacies of ligands for these receptors are unclear. Therefore, in this study, the binding energy profiles of various ligands (full agonists, partial agonists, and antagonists) were quantified by docking those ligands with structural ensembles of the α7 nAChR exhibiting different degrees of C-loop closure. This approximate treatment of interactions suggested that full agonists, partial agonists, and antagonists of the α7 nAChR possess distinctive binding energy profiles. Results from docking revealed that ligand binding efficacy may be related to the capacity of the ligand to stabilize conformational states with a closed C loop.

Genetic Algorithm Managed Peptide Mutant Screening: Optimizing Peptide Ligands for Targeted Receptor Binding

This study demonstrates the utility of genetic algorithms to search exceptionally large and otherwise intractable mutant libraries for sequences with optimal binding affinities for target receptors. The Genetic Algorithm Managed Peptide Mutant Screening (GAMPMS) program was used to search an α-conotoxin (α-CTx) MII mutant library of approximately 41 billion possible peptide sequences for those exhibiting the greatest binding affinity for the α₃β₂-nicotinic acetylcholine receptor (nAChR) isoform. A series of top resulting peptide ligands with high sequence homology was obtained, with each mutant having an estimated ΔG_bind approximately double that of the potent native α-CTx MII ligand. A consensus sequence from the top GAMPMS results was subjected to more rigorous binding free energy calculations by molecular dynamics and compared to α-CTx MII and other related variants for binding with α₃β₂-nAChR. In this study, the efficiency of GAMPMS to substantially reduce the sample population size through evolutionary selection criteria to produce ligands with higher predicted binding affinity is demonstrated.

Loop-driven conformational transition between the alternative and collapsed form of prethrombin-2: targeted molecular dynamics study

Two distinct crystal structures of prethrombin-2, the alternative and collapsed forms, are elucidated by X-ray crystallogrphy. We analyzed the conformational transition from the alternative to the collapsed form employing targeted molecular dynamics (TMD) simulation. Despite small RMSD difference in the two X-ray crystal structures, some hydrophobic residues (W60d, W148, W215, and F227) show a significant difference between the two conformations. TMD simulation shows that the four hydrophobic residues undergo concerted movement from dimer to trimer transition via tetramer state in the conformational change from the alternative to the collapsed form. We reveal that the concerted movement of the four hydrophobic residues is controlled by movement of specific loop regions behind. In this paper, we propose a sequential scenario for the conformational transition from the alternative form to the collapsed form, which is partially supported by the mutant W148A simulation.

The Free Energy Landscape of GABA Binding to a Pentameric Ligand-Gated Ion Channel and Its Disruption by Mutations

Pentameric ligand-gated ion channels (pLGICs) of the Cys-loop superfamily are important neuroreceptors that mediate fast synaptic transmission. They are activated by the binding of a neurotransmitter, but the details of this process are still not fully understood. As a prototypical pLGIC, here we choose the insect resistance to dieldrin (RDL) receptor involved in resistance to insecticides and investigate the binding of the neurotransmitter GABA to its extracellular domain at the atomistic level. We achieve this by means of μ-sec funnel-metadynamics simulations, which efficiently enhance the sampling of bound and unbound states by using a funnel-shaped restraining potential to limit the exploration in the solvent. We reveal the sequence of events in the binding process from the capture of GABA from the solvent to its pinning between the charged residues Arg111 and Glu204 in the binding pocket. We characterize the associated free energy landscapes in the wild-type RDL receptor and in two mutant forms, where the key residues Arg111 and Glu204 are mutated to Ala. Experimentally these mutations produce nonfunctional channels, which is reflected in the reduced ligand binding affinities due to the loss of essential interactions. We also analyze the dynamical behavior of the crucial loop C, whose opening allows the access of GABA to the binding site and closure locks the ligand into the protein. The RDL receptor shares structural and functional features with other pLGICs; hence, our work outlines a valuable protocol to study the binding of ligands to pLGICs beyond conventional docking and molecular dynamics techniques.

Study on the Application of the Combination of TMD Simulation and Umbrella Sampling in PMF Calculation for Molecular Conformational Transitions

Free energy calculations of the potential of mean force (PMF) based on the combination of targeted molecular dynamics (TMD) simulations and umbrella samplings as a function of physical coordinates have been applied to explore the detailed pathways and the corresponding free energy profiles for the conformational transition processes of the butane molecule and the 35-residue villin headpiece subdomain (HP35). The accurate PMF profiles for describing the dihedral rotation of butane under both coordinates of dihedral rotation and root mean square deviation (RMSD) variation were obtained based on the different umbrella samplings from the same TMD simulations. The initial structures for the umbrella samplings can be conveniently selected from the TMD trajectories. For the application of this computational method in the unfolding process of the HP35 protein, the PMF calculation along with the coordinate of the radius of gyration (R_g) presents the gradual increase of free energies by about 1 kcal/mol with the energy fluctuations. The feature of conformational transition for the unfolding process of the HP35 protein shows that the spherical structure extends and the middle α-helix unfolds firstly, followed by the unfolding of other α-helices. The computational method for the PMF calculations based on the combination of TMD simulations and umbrella samplings provided a valuable strategy in investigating detailed conformational transition pathways for other allosteric processes.

Secondary Interaction Interfaces with PCNA Control Conformational Switching of DNA Polymerase PolB from Polymerization to Editing

Replicative DNA polymerases (Pols) frequently possess two distinct DNA processing activities: DNA synthesis (polymerization) and proofreading (3'-5' exonuclease activity). The polymerase and exonuclease reactions are performed alternately and are spatially separated in different protein domains. Thus, the growing DNA primer terminus has to undergo dynamic conformational switching between two distinct functional sites on the polymerase. Furthermore, the transition from polymerization (pol) mode to exonuclease (exo) mode must occur in the context of a DNA Pol holoenzyme, wherein the polymerase is physically associated with processivity factor proliferating cell nuclear antigen (PCNA) and primer-template DNA. The mechanism of this conformational switching and the role that PCNA plays in it have remained obscure, largely due to the dynamic nature of ternary Pol/PCNA/DNA assemblies. Here, we present computational models of ternary assemblies for archaeal polymerase PolB. We have combined all available structural information for the binary complexes with electron microscopy data and have refined atomistic models for ternary PolB/PCNA/DNA assemblies in pol and exo modes using molecular dynamics simulations. In addition to the canonical PIP-box/interdomain connector loop (IDCL) interface of PolB with PCNA, contact analysis of the simulation trajectories revealed new secondary binding interfaces, distinct between the pol and exo states. Using targeted molecular dynamics, we explored the conformational transition from pol to exo mode. We identified a hinge region between the thumb and palm domains of PolB that is critical for conformational switching. With the thumb domain anchored onto the PCNA surface, the neighboring palm domain executed rotational motion around the hinge, bringing the core of PolB down toward PCNA to form a new interface with the clamp. A helix from PolB containing a patch of arginine residues was involved in the binding, locking the complex in the exo mode conformation. Together, these results provide a structural view of how the transition between the pol and exo states of PolB is coordinated through PCNA to achieve efficient proofreading.

Allosteric Modulation of Alpha7 Nicotinic Receptors: Mechanistic Insight through Metadynamics and Essential Dynamics

Increasing attention has recently been devoted to allosteric modulators, as they can provide inherent advantages over classic receptor agonists. In the field of nicotinic receptors (nAChRs), the main advantage is that allosteric modulators can trigger pharmacological responses, limiting receptor desensitization. Most of the known allosteric ligands are "positive allosteric modulators" (PAMs), which increase both sensitivity to receptor agonists and current amplitude. Intriguingly, some allosteric modulators are also able to activate the α7 receptor (α7-nAChR) even in the absence of orthosteric agonists. These compounds have been named "ago-allosteric modulators" and GAT107 has been studied in depth because of its unique mechanism of action. We here investigate by molecular dynamics simulations, metadynamics, and essential dynamics the activation mechanism of α7-nAChR, in the presence of different nicotinic modulators. We determine the free energy profiles associated with the closed-to-open motion of the loop C, and we highlight mechanistic differences observed in the presence of different modulators. In particular, we demonstrate that GAT107 triggers conformational motions and cross-talk similar to those observed when the α7-nACh receptor is in complex with both an agonist and an allosteric modulator.

Pathways and Barriers for Ion Translocation through the 5-HT3A Receptor Channel

Pentameric ligand gated ion channels (pLGICs) are ionotropic receptors that mediate fast intercellular communications at synaptic level and include either cation selective (e.g., nAChR and 5-HT₃) or anion selective (e.g., GlyR, GABA_A and GluCl) membrane channels. Among others, 5-HT₃ is one of the most studied members, since its first cloning back in 1991, and a large number of studies have successfully pinpointed protein residues critical for its activation and channel gating. In addition, 5-HT₃ is also the target of a few pharmacological treatments due to the demonstrated benefits of its modulation in clinical trials. Nonetheless, a detailed molecular analysis of important protein features, such as the origin of its ion selectivity and the rather low conductance as compared to other channel homologues, has been unfeasible until the recent crystallization of the mouse 5-HT₃A receptor. Here, we present extended molecular dynamics simulations and free energy calculations of the whole 5-HT₃A protein with the aim of better understanding its ion transport properties, such as the pathways for ion permeation into the receptor body and the complex nature of the selectivity filter. Our investigation unravels previously unpredicted structural features of the 5-HT₃A receptor, such as the existence of alternative intersubunit pathways for ion translocation at the interface between the extracellular and the transmembrane domains, in addition to the one along the channel main axis. Moreover, our study offers a molecular interpretation of the role played by an arginine triplet located in the intracellular domain on determining the characteristic low conductance of the 5-HT₃A receptor, as evidenced in previous experiments. In view of these results, possible implications on other members of the superfamily are suggested.

Functional Impact of 14 Single Nucleotide Polymorphisms Causing Missense Mutations of Human α7 Nicotinic Receptor

The α7nicotinic receptor (nAChR) is a major subtype of the nAChRs in the central nervous system, and the receptor plays an important role in brain function. In the dbSNP database, there are 55 single nucleotide polymorphisms (SNPs) that cause missense mutations of the human α7nAChR in the coding region. In this study, we tested the impact of 14 SNPs that cause missense mutations in the agonist binding site or the coupling region between binding site and channel gate on the receptor function. The wild type or mutant receptors were expressed or co-expressed in Xenopus oocytes, and the agonist-induced currents were tested using two-electrode voltage clamp. Our results demonstrated that 6 mutants were nonfunctional, 4 mutants had reduced current expression, and 1 mutants altered ACh and nicotine efficacy in the opposite direction, and one additional mutant had slightly reduced agonist sensitivity. Interestingly, the function of most of these nonfunctional mutants could be rescued by α7nAChR positive allosteric modulator PNU-120596 and agonist-PAM 4BP-TQS. Finally, when coexpressed with the wild type, the nonfunctional mutants could also influence the receptor function. These changes of the receptor properties by the mutations could potentially have an impact on the physiological function of the α7nAChR-mediated cholinergic synaptic transmission and anti-inflammatory effects in the human SNP carriers. Rescuing the nonfunctional mutants could provide a novel way to treat the related disorders.

Conformational gating dynamics in the GluCl anion-selective chloride channel

Cys-loop receptors are central to propagation of signals in the nervous system. The gating of the membrane-spanning pore is triggered by structural rearrangements in the agonist-binding site, located some 50 Å away from the pore. A sequential conformational change, propagating from the ligand-binding site to the pore, has been proposed to govern gating in all Cys-loop receptors. Here, we identify structural and dynamic components of the conformational gating in the eukaryotic glutamate-gated chloride channel (GluCl) by means of molecular dynamics (MD) simulations with and without the l-glutamate agonist bound. A significant increase in pore opening and accompanying hydration is observed in the presence of glutamate. Potential of mean force calculations reveal that the barrier for ion passage drops from 15 kcal/mol to 5-10 kcal/mol with the agonist bound. This appears to be explained by agonist binding that leads to significant changes in the intersubunit hydrogen-bonding pattern, which induce a slight tilt of the extracellular domain relative to the transmembrane domain in the simulations. This rearrangement is subtle, but correspond to the direction of the quaternary twist observed as a key difference between open and closed X-ray structures. While the full reversible gating is still a much slower process, the observed structural dynamics sheds new light on the early stages of how the agonist influences the extracellular domain, how the extracellular domain interacts with the transmembrane domain, and how changes in the transmembrane domain alter the free energy of ion passage.

A Structural Model of the Human α7 Nicotinic Receptor in an Open Conformation

Nicotinic acetylcholine receptors (nAchRs) are ligand-gated ion channels that regulate chemical transmission at the neuromuscular junction. Structural information is available at low resolution from open and closed forms of an eukaryotic receptor, and at high resolution from other members of the same structural family, two prokaryotic orthologs and an eukaryotic GluCl channel. Structures of human channels however are still lacking. Homology modeling and Molecular Dynamics simulations are valuable tools to predict structures of unknown proteins, however, for the case of human nAchRs, they have been unsuccessful in providing a stable open structure so far. This is due to different problems with the template structures: on one side the homology with prokaryotic species is too low, while on the other the open eukaryotic GluCl proved itself unstable in several MD studies and collapsed to a dehydrated, non-conductive conformation, even when bound to an agonist. Aim of this work is to obtain, by a mixing of state-of-the-art homology and simulation techniques, a plausible prediction of the structure (still unknown) of the open state of human α7 nAChR complexed with epibatidine, from which it is possible to start structural and functional test studies. To prevent channel closure we employ a restraint that keeps the transmembrane pore open, and obtain in this way a stable, hydrated conformation. To further validate this conformation, we run four long, unbiased simulations starting from configurations chosen at random along the restrained trajectory. The channel remains stable and hydrated over the whole runs. This allows to assess the stability of the putative open conformation over a cumulative time of 1 μs, 800 ns of which are of unbiased simulation. Mostly based on the analysis of pore hydration and size, we suggest that the obtained structure has reasonable chances to be (at least one of the possible) structures of the channel in the open conformation.

Identification of the Conformational transition pathway in PIP2 Opening Kir Channels

The gating of Kir channels depends critically on phosphatidylinositol 4,5-bisphosphate (PIP2), but the detailed mechanism by which PIP2 regulates Kir channels remains obscure. Here, we performed a series of Targeted molecular dynamics simulations on the full-length Kir2.1 channel and, for the first time, were able to achieve the transition from the closed to the open state. Our data show that with the upward motion of the cytoplasmic domain (CTD) the structure of the C-Linker changes from a loop to a helix. The twisting of the C-linker triggers the rotation of the CTD, which induces a small downward movement of the CTD and an upward motion of the slide helix toward the membrane that pulls the inner helix gate open. At the same time, the rotation of the CTD breaks the interaction between the CD- and G-loops thus releasing the G-loop. The G-loop then bounces away from the CD-loop, which leads to the opening of the G-loop gate and the full opening of the pore. We identified a series of interaction networks, between the N-terminus, CD loop, C linker and G loop one by one, which exquisitely regulates the global conformational changes during the opening of Kir channels by PIP2.

Substrate binding accelerates the conformational transitions and substrate dissociation in multidrug efflux transporter AcrB

The tripartite efflux pump assembly AcrAB-TolC is the major multidrug resistance transporter in E. coli. The inner membrane transporter AcrB is a homotrimer, energized by the proton movement down the transmembrane electrochemical gradient. The asymmetric crystal structures of AcrB with three monomers in distinct conformational states [access (A), binding (B) and extrusion (E)] support a functional rotating mechanism, in which each monomer of AcrB cycles among the three states in a concerted way. However, the relationship between the conformational changes during functional rotation and drug translocation has not been totally understood. Here, we explored the conformational changes of the AcrB homotrimer during the ABE to BEA transition in different substrate-binding states using targeted MD simulations. It was found that the dissociation of substrate from the distal binding pocket of B monomer is closely related to the concerted conformational changes in the translocation pathway, especially the side chain reorientation of Phe628 and Tyr327. A second substrate binding at the proximal binding pocket of A monomer evidently accelerates the conformational transitions as well as substrate dissociation in B monomer. The acceleration effect of the multi-substrate binding mode provides a molecular explanation for the positive cooperativity observed in the kinetic studies of substrate efflux and deepens our understanding of the functional rotating mechanism of AcrB.

Constrained geometric simulation of the nicotinic acetylcholine receptor

Constrained geometric simulations have been performed for the recently published closed-channel state of the nicotinic acetylcholine receptor. These simulations support the theory that correlated motion in the flexible β-sheet structure of the extracellular domain helps to communicate a "conformational wave", spreading from the acetylcholine binding pocket. Furthermore, we have identified key residues that act at the interface between subunits and between domains that could potentially facilitate rapid communication between the binding site and the transmembrane gate.

Conformational changes in acetylcholine binding protein investigated by temperature accelerated molecular dynamics

Despite the large number of studies available on nicotinic acetylcholine receptors, a complete account of the mechanistic aspects of their gating transition in response to ligand binding still remains elusive. As a first step toward dissecting the transition mechanism by accelerated sampling techniques, we study the ligand-induced conformational changes of the acetylcholine binding protein (AChBP), a widely accepted model for the full receptor extracellular domain. Using unbiased Molecular Dynamics (MD) and Temperature Accelerated Molecular Dynamics (TAMD) simulations we investigate the AChBP transition between the apo and the agonist-bound state. In long standard MD simulations, both conformations of the native protein are stable, while the agonist-bound structure evolves toward the apo one if the orientation of few key sidechains in the orthosteric cavity is modified. Conversely, TAMD simulations initiated from the native conformations are able to produce the spontaneous transition. With respect to the modified conformations, TAMD accelerates the transition by at least a factor 10. The analysis of some specific residue-residue interactions points out that the transition mechanism is based on the disruption/formation of few key hydrogen bonds. Finally, while early events of ligand dissociation are observed already in standard MD, TAMD accelerates the ligand detachment and, at the highest TAMD effective temperature, it is able to produce a complete dissociation path in one AChBP subunit.

Computing the relative stabilities and the per-residue components in protein conformational changes

Protein molecules often undergo conformational changes. In order to gain insights into the forces that drive such changes, it would be useful to have a method that computes the per-residue contributions to the conversion free energy. Here, we describe the "confine-convert-release" (CCR) method, which is applicable to large conformational changes. We show that CCR correctly predicts the stable states of several "chameleon" sequences that have previously been challenging for molecular simulations. CCR can often discriminate better from worse predictions of native protein models in critical assessment of protein structure prediction (CASP). We show how the total conversion free energies can be parsed into per-residue free-energy components. Such parsing gives insights into which amino acids are most responsible for given transformations. For example, here we are able to "reverse-engineer" the known design principles of the chameleon proteins. This opens up the possibility for systematic improvements in structure-prediction scoring functions, in the design of protein conformational switches, and in interpreting protein mechanisms at the amino-acid level.

Surface-Tension Replica-Exchange Molecular Dynamics Method for Enhanced Sampling of Biological Membrane Systems

Conformational sampling is fundamentally important for simulating complex biomolecular systems. The generalized-ensemble algorithm, especially the temperature replica-exchange molecular dynamics method (T-REMD), is one of the most powerful methods to explore structures of biomolecules such as proteins, nucleic acids, carbohydrates, and also of lipid membranes. T-REMD simulations have focused on soluble proteins rather than membrane proteins or lipid bilayers, because explicit membranes do not keep their structural integrity at high temperature. Here, we propose a new generalized-ensemble algorithm for membrane systems, which we call the surface-tension REMD method. Each replica is simulated in the NPγT ensemble, and surface tensions in a pair of replicas are exchanged at certain intervals to enhance conformational sampling of the target membrane system. We test the method on two biological membrane systems: a fully hydrated DPPC (1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine) lipid bilayer and a WALP23-POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) membrane system. During these simulations, a random walk in surface tension space is realized. Large-scale lateral deformation (shrinking and stretching) of the membranes takes place in all of the replicas without collapse of the lipid bilayer structure. There is accelerated lateral diffusion of DPPC lipid molecules compared with conventional MD simulation, and a much wider range of tilt angle of the WALP23 peptide is sampled due to large deformation of the POPC lipid bilayer and through peptide-lipid interactions. Our method could be applicable to a wide variety of biological membrane systems.

Interactions of acetylcholine binding site residues contributing to nicotinic acetylcholine receptor gating: role of residues Y93, Y190, K145 and D200

The nicotinic acetylcholine receptor exhibits multiple conformational states, resting (channel closed), active (channel open) and desensitized (channel closed). The resting state may be distinguished from the active and desensitized states by the orientation of loop C in the extracellular ligand binding domain (LBD). Homology modeling was used to generate structures of the Torpedo californica α2βδγ nAChR that initially represent the resting state (loop C open) and the desensitized state (loop C closed). Molecular dynamics (MD) simulations were performed on the extracellular LBD on each nAChR conformational state, with and without the agonist anabaseine present in each binding site (the αγ and the αδ sites). Three MD simulations of 10ns each were performed for each of the four conditions. Comparison of dynamics revealed that in the presence of agonist, loop C was drawn inward and attains a more stable conformation. Examination of side-chain interactions revealed that residue αY190 exhibited hydrogen-bonding interactions either with residue αY93 in the ligand binding site or with residue αK145 proximal to the binding site. αK145 also exhibited side chain (salt bridge) interactions with αD200 and main chain interactions with αY93. Residues αW149, αY198, γY116/δT119, γL118/δL121 and γL108/δL111 appear to play the role of stabilizing ligand in the binding site. In MD simulations for the desensitized state, the effect of ligand upon the interactions among αK145, αY190, and αY93 as well as ligand-hydrogen-bonding to αW149 were more pronounced at the αγ interface than at the αδ interface. Differences in affinity for the desensitized state were determined experimentally to be 10-fold. The changes in side chain interactions observed for the two conformations and induced by ligand support a model wherein hydrogen bond interactions between αD200 and αY93 are broken and rearrange to form a salt-bridge between αK145 and αD200 and hydrogen bond interactions between αY93 and αY190 and between αK145 and αY190.

Computational gibberellin-binding channel discovery unraveling the unexpected perception mechanism of hormone signal by gibberellin receptor

Gibberellins (GAs) are phytohormones essential for many developmental processes in plants. In this work, fundamental mechanism of hormone perception by receptor GID1 has been studied by performing computational simulations, revealing a new GA-binding channel of GID1 and a novel hormone perception mechanism involving only one conformational state of GID1. The novel hormone perception mechanism demonstrated here is remarkably different from the previously proposed/speculated mechanism [Murase et al., Nature 2008, 456, 459] involving two conformational states ("OPEN" and "CLOSED") of GID1. According to the new perception mechanism, GA acts as a "conformational stabilizer," rather than the previously speculated "allosteric inducer," to induce the recognition of protein DELLA by GID1. The novel mechanistic insights obtained in this study provide a new starting point for further studies on the detailed molecular mechanisms of GID1 interacting with DELLA and various hormones and for mechanism-based rational design of novel, potent growth regulators that target crops and ornamental plants.

New horizons for congenital myasthenic syndromes

During the past five years an increasing number of patients have been diagnosed with congenital myasthenic syndromes (CMS) and a number of novel syndromes have been recognized and investigated. This presentation focuses on the CMS caused by defects in choline acetyltransferase, novel fast-channel syndromes that hinder isomerization of the acetylcholine receptor from the closed to the open state, the consequences of deleterious mutations in the intermediate filament linker plectin, altered neuromuscular transmission in a centronuclear myopathy, and two recently identified CMS caused by congenital defects in glycosylation.

GABA binding to an insect GABA receptor: a molecular dynamics and mutagenesis study

RDL receptors are GABA-activated inhibitory Cys-loop receptors found throughout the insect CNS. They are a key target for insecticides. Here, we characterize the GABA binding site in RDL receptors using computational and electrophysiological techniques. A homology model of the extracellular domain of RDL was generated and GABA docked into the binding site. Molecular dynamics simulations predicted critical GABA binding interactions with aromatic residues F206, Y254, and Y109 and hydrophilic residues E204, S176, R111, R166, S176, and T251. These residues were mutated, expressed in Xenopus oocytes, and their functions assessed using electrophysiology. The data support the binding mechanism provided by the simulations, which predict that GABA forms many interactions with binding site residues, the most significant of which are cation-π interactions with F206 and Y254, H-bonds with E204, S205, R111, S176, T251, and ionic interactions with R111 and E204. These findings clarify the roles of a range of residues in binding GABA in the RDL receptor, and also show that molecular dynamics simulations are a useful tool to identify specific interactions in Cys-loop receptors.

Size matters in activation/inhibition of ligand-gated ion channels

Cys loop, glutamate, and P2X receptors are ligand-gated ion channels (LGICs) with 5, 4, and 3 protomers, respectively. There is now growing atomic level understanding of their gating mechanisms. Although each family is unique in the architecture of the ligand-binding pocket, the pathway for motions to propagate from ligand-binding domain to transmembrane domain, and the gating motions of the transmembrane domain, there are common features among the LGICs, which are the focus of the present review. In particular, agonists and competitive antagonists apparently induce opposite motions of the binding pocket. A simple way to control the motional direction is ligand size. Agonists, usually small, induce closure of the binding pocket, leading to opening of the channel pore, whereas antagonists, usually large, induce opening of the binding pocket, thereby stabilizing the closed pore. A cross-family comparison of the gating mechanisms of the LGICs, focusing in particular on the role played by ligand size, provides new insight on channel activation/inhibition and design of pharmacological compounds.

α6 nAChR subunit residues that confer α-conotoxin BuIA selectivity

Nicotinic acetylcholine receptors (nAChRs) containing α6 and/or α4 subunits modulate the release of dopamine. However, few compounds can effectively discriminate between ligand-binding sites that contain α6 vs. α4 nAChR subunits. Using a chimeric (α6/α4) subunit, we showed that α-conotoxin BuIA binds the extracellular rat α6β2 vs. α4β2 interface with ∼60,000-fold selectivity. Chimeras containing residues from the α6 subunit were inserted into the homologous position of the α4 subunit to identify critical sequence segments. The region between residues 184 and 207 in the α6 subunit accounted for the potency difference. Chimeras within this region followed by point mutations were constructed for further definition. α6 Lys185, Thr187, and Ile188 form a triad of key residues that influence BuIA binding; when these 3 α6 residues were inserted into the α4 subunit, there was an ∼2000-fold increase in toxin potency. We used a crystal structure of BuIA bound to the acetylcholine-binding protein together with the structure of the Torepedo marmorata nAChR to build a homology model of BuIA bound to the interface between α6 and β2 subunits. The results indicate that the triad of α6 residues lies outside the C loop and is distantly located from bound BuIA (>10 Å). This suggests that alterations in potency are not caused by the direct interaction between the triad and BuIA. Instead, alterations in C-loop 3-dimensional structure and/or flexibility may account for differential potency. Thr198 and Tyr205 also contributed to BuIA potency. In addition, Thr198 caused BuIA potency differences between the closely related α6 and α3 subunits. Together, the findings provide insight into differences between the α6 and other α subunits that may be exploited by α-conotoxins to achieve binding selectivity.

End-plate acetylcholine receptor: structure, mechanism, pharmacology, and disease

The synapse is a localized neurohumoral contact between a neuron and an effector cell and may be considered the quantum of fast intercellular communication. Analogously, the postsynaptic neurotransmitter receptor may be considered the quantum of fast chemical to electrical transduction. Our understanding of postsynaptic receptors began to develop about a hundred years ago with the demonstration that electrical stimulation of the vagus nerve released acetylcholine and slowed the heart beat. During the past 50 years, advances in understanding postsynaptic receptors increased at a rapid pace, owing largely to studies of the acetylcholine receptor (AChR) at the motor endplate. The endplate AChR belongs to a large superfamily of neurotransmitter receptors, called Cys-loop receptors, and has served as an exemplar receptor for probing fundamental structures and mechanisms that underlie fast synaptic transmission in the central and peripheral nervous systems. Recent studies provide an increasingly detailed picture of the structure of the AChR and the symphony of molecular motions that underpin its remarkably fast and efficient chemoelectrical transduction.

Myasthenic syndrome AChRα C-loop mutant disrupts initiation of channel gating

Congenital myasthenic syndromes (CMSs) are neuromuscular disorders that can be caused by defects in ace-tylcholine receptor (AChR) function. Disease-associated point mutants can reveal the unsuspected functional significance of mutated residues. We identified two pathogenic mutations in the extracellular domain of the AChR α subunit (AChRα) in a patient with myasthenic symptoms since birth: a V188M mutation in the C-loop and a heteroallelic G74C mutation in the main immunogenic region. The G74C mutation markedly reduced surface AChR expression in cultured cells, whereas the V188M mutant was expressed robustly but had severely impaired kinetics. Single-channel patch-clamp analysis indicated that V188M markedly decreased the apparent AChR channel opening rate and gating efficiency. Mutant cycle analysis of energetic coupling among conserved residues within or dispersed around the AChRα C-loop revealed that V188 is functionally linked to Y190 in the C-loop and to D200 in β-strand 10, which connects to the M1 transmembrane domain. Furthermore, V188M weakens inter-residue coupling of K145 in β-strand 7 with Y190 and with D200. Cumulatively, these results indicate that V188 of AChRα is part of an interdependent tetrad that contributes to rearrangement of the C-loop during the initial coupling of agonist binding to channel gating.