Conformational dynamics of the nicotinic acetylcholine receptor channel: a 35-ns molecular dynamics simulation study

State-dependent binding of cholesterol and an anionic lipid to the muscle-type Torpedo nicotinic acetylcholine receptor

The ability of the Torpedo nicotinic acetylcholine receptor (nAChR) to undergo agonist-induced conformational transitions requires the presence of cholesterol and/or anionic lipids. Here we use recently solved structures along with multiscale molecular dynamics simulations to examine lipid binding to the nAChR in bilayers that have defined effects on nAChR function. We examine how phosphatidic acid and cholesterol, lipids that support conformational transitions, individually compete for binding with phosphatidylcholine, a lipid that does not. We also examine how the two lipids work synergistically to stabilize an agonist-responsive nAChR. We identify rapidly exchanging lipid binding sites, including both phospholipid sites with a high affinity for phosphatidic acid and promiscuous cholesterol binding sites in the grooves between adjacent transmembrane α-helices. A high affinity cholesterol site is confirmed in the inner leaflet framed by a key tryptophan residue on the MX α-helix. Our data provide insight into the dynamic nature of lipid-nAChR interactions and set the stage for a detailed understanding of the mechanisms by which lipids facilitate nAChR function at the neuromuscular junction.

Oleic Acid Could Act as a Channel Blocker in the Inhibition of nAChR: Insights from Molecular Dynamics Simulations

The activation of the muscular nicotinic acetylcholine receptor (nAChR) produces the opening of the channel, with the consequent increase in the permeability of cations, triggering an excitatory signal. Free fatty acids (FFA) are known to modulate the activity of the receptor as noncompetitive antagonists, acting at the membrane-AChR interface. We present molecular dynamics simulations of a model of nAChR in a desensitized closed state embedded in a lipid bilayer in which distinct membrane phospholipids were replaced by two different monounsaturated FFA that differ in the position of a double bond. This allowed us to detect and describe that the cis-18:1ω-9 FFA were located at the interface between the transmembrane segments of α2 and γ subunits diffused into the channel lumen with the consequent potential ability to block the channel to the passage of ions.

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'.

The functional role of the αM4 transmembrane helix in the muscle nicotinic acetylcholine receptor probed through mutagenesis and coevolutionary analyses

The activity of the muscle-type Torpedo nicotinic acetylcholine receptor (nAChR) is highly sensitive to lipids, but the underlying mechanisms remain poorly understood. The nAChR transmembrane α-helix, M4, is positioned at the perimeter of each subunit in direct contact with lipids and likely plays a central role in lipid sensing. To gain insight into the mechanisms underlying nAChR lipid sensing, we used homology modeling, coevolutionary analyses, site-directed mutagenesis, and electrophysiology to examine the role of the α-subunit M4 (αM4) in the function of the adult muscle nAChR. Ala substitutions for most αM4 residues, including those in clusters of polar residues at both the N and C termini, and deletion of up to 11 C-terminal residues had little impact on the agonist-induced response. Even Ala substitutions for coevolved pairs of residues at the interface between αM4 and the adjacent helices, αM1 and αM3, had little effect, although some impaired nAChR expression. On the other hand, Ala substitutions for Thr422 and Arg429 caused relatively large losses of function, suggesting functional roles for these specific residues. Ala substitutions for aromatic residues at the αM4-αM1/αM3 interface generally led to gains of function, as previously reported for the prokaryotic homolog, the Erwinia chrysanthemi ligand-gated ion channel (ELIC). The functional effects of individual Ala substitutions in αM4 were found to be additive, although not in a completely independent manner. Our results provide insight into the structural features of αM4 that are important. They also suggest how lipid-dependent changes in αM4 structure ultimately modify nAChR function.

A Review on the Receptor-ligand Molecular Interactions in the Nicotinic Receptor Signaling Systems

Nicotine is regarded as the main active addictive ingredient in tobacco products driving continued tobacco abuse behavior (smoking) to the addiction behavior, whereas nicotinic acetylcholine receptors (nAChR) is the crucial effective apparatus or molecular effector of nicotine and acetylcholine and other similar ligands. Many nAChR subunits have been revealed to bind to either neurotransmitters or exogenous ligands, such as nicotine and acetylcholine, being involved in the nicotinic receptor signal transduction. Therefore, the nicotinic receptor signalling molecules and the receptor-ligand molecular interactions between nAChRs and their ligands are universally regarded as crucial mediators of cellular functions and drug targets in medical treatment and clinical diagnosis. Given numerous endeavours have been made in defining the roles of nAChRs in response to nicotine and other addictive drugs, this review focuses on studies and reports in recent years on the receptor-ligand interactions between nAChR receptors and ligands, including lipid-nAChR and protein-nAChR molecular interactions, relevant signal transduction pathways and their molecular mechanisms in the nicotinic receptor signalling systems. All the references were carefully retrieved from the PubMed database by searching key words "nicotine", "acetylcholine", "nicotinic acetylcholine receptor(s)", "nAChR*", "protein and nAChR", "lipid and nAChR", "smok*" and "tobacco". All the relevant referred papers and reports retrieved were fully reviewed for manual inspection. This effort intend to get a quick insight and understanding of the nicotinic receptor signalling and their molecular interactions mechanisms. Understanding the cellular receptor-ligand interactions and molecular mechanisms between nAChRs and ligands will lead to a better translational and therapeutic operations and outcomes for the prevention and treatment of nicotine addiction and other chronic drug addictions in the brain's reward circuitry.

Towards the SMART workflow system for computational spectroscopy

Is it possible to convert highly specialized research in the field of computational spectroscopy into robust and user-friendly aids to experiments and industrial applications?

Raman scattering-based multiconformational analysis for probing the structural differences between acetylcholine and acetylthiocholine

Acetylcholine is the first discovered neurotransmitter that has received a great attention regarding its capability of binding to several cellular targets. The chemical composition of acetylcholine, including a positively charged trimethylammonium and a carbonyl group, as well as its conformational flexibility was pointed out as the key factors in the stabilization of its interactions. Here, the possibilities offered by a Raman scattering-based multiconformatioal analysis to access the most stable conformers of acetylcholine, is discussed. To control the validity of this protocol, acetylcholine and one of its closely structured analogues, acetylthiocholine, were simultaneously analyzed. Solution Raman spectra revealed distinct and well resolved strong markers for each molecule. Density functional theory calculations were consistent with the fact that the energy order of the low energy conformers is considerably affected by the acyloxy oxygen→sulfur atom substitution. Raman spectra were calculated on the basis of the thermal average of the spectra arising from the low energy conformers. It has been evidenced that the carbonyl and trimethylammonium groups are the most favorable hydration sites in aqueous environment. Taking into account the large gap between the carbonyl bond-stretch and aliphatic bending bands, Raman spectra also allowed separation of the HOH bending vibrations arising from the bound and bulk water molecules.

The molecular behavior of a single β-amyloid inside a dipalmitoylphosphatidylcholine bilayer at three different temperatures: An atomistic simulation study: Aβ interaction with DPPC: Atomistic simulation

The behavior of a single Aβ40 molecule within a dipalmitoylphosphatidylcholine (DPPC) bilayer was studied by all-atom molecular dynamics simulations. The effect of membrane structure was investigated on Aβ40 behavior, secondary structure, and insertion depth. Simulations were performed at three temperatures (323, 310, and 300 K) to probe three different bilayer fluidities. Results show that at all above temperatures, the peptide contains two short helices, coil, bend, and turn structures. At 300 K, the peptide contains a region with β structure in C-terminal region. Our results also show that Aβ decreases the bilayer thickness and the order of lipids in its vicinity which leads to water insertion into the bilayer and concomitant increase in the local fluidity. The peptide remains embedded in the bilayer at all temperatures, and become inserted into the bilayer up to several residues at 323 and 310 K. At 310 and 300 K, the dominant interaction energy between Aβ and bilayer changes from electrostatic to van der Waals. It can be proposed that at higher temperatures (e.g., 323 K), Lys28 and the C-terminal region of the peptide play the role of two anchors that keep Aβ inside the top leaflet. This study demonstrates that Aβ molecule can perturb the integrity of cellular membranes. Proteins 2017; 85:1298-1310. © 2017 Wiley Periodicals, Inc.

Fatty Acid Regulation of Voltage- and Ligand-Gated Ion Channel Function

Free fatty acids (FFA) are essential components of the cell, where they play a key role in lipid and carbohydrate metabolism, and most particularly in cell membranes, where they are central actors in shaping the physicochemical properties of the lipid bilayer and the cellular adaptation to the environment. FFA are continuously being produced and degraded, and a feedback regulatory function has been attributed to their turnover. The massive increase observed under some pathological conditions, especially in brain, has been interpreted as a protective mechanism possibly operative on ion channels, which in some cases is of stimulatory nature and in other cases inhibitory. Here we discuss the correlation between the structure of FFA and their ability to modulate protein function, evaluating the influence of saturation/unsaturation, number of double bonds, and cis vs. trans isomerism. We further focus on the mechanisms of FFA modulation operating on voltage-gated and ligand-gated ion channel function, contrasting the still conflicting evidence on direct vs. indirect mechanisms of action.

Probing the structure of the uncoupled nicotinic acetylcholine receptor

In the absence of activating anionic lipids and cholesterol, the nicotinic acetylcholine receptor (nAChR) from Torpedo adopts an uncoupled conformation that does not usually gate open in response to agonist. The uncoupled conformation binds both agonists and non-competitive channel blockers with a lower affinity than the desensitized state, consistent with both the extracellular agonist-binding and transmembrane channel-gating domains individually adopting resting-state like conformations. To test this hypothesis, we characterized the binding of the agonist, acetylcholine, and two fluorescent channel blockers, ethidium and crystal violet, to resting, desensitized and uncoupled nAChRs in reconstituted membranes. The measured K_d for acetylcholine binding to the uncoupled nAChR is similar to that for the resting state, confirming that the agonist binding site adopts a resting-state like conformation. Although both ethidium and crystal violet bind to the resting and desensitized channel pores with distinct affinities, no binding of either probe was detected to the uncoupled nAChR. Our data suggest that the transmembrane domain of the uncoupled nAChR adopts a conformation distinct from that of the resting and desensitized states. The lack of binding is consistent with a more constricted channel pore, possibly along the lines of what is observed in crystal structures of the prokaryotic homolog, ELIC.

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 nicotinic acetylcholine receptor reconstitution in Au(111)-supported thiolipid monolayers

The insertion and function of the muscle-type nicotinic acetylcholine receptor (nAChR) in Au(111)-supported thiolipid self-assembled monolayers have been studied by atomic force microscopy (AFM), surface plasmon resonance (SPR), and electrochemical techniques. It was possible for the first time to resolve the supramolecular arrangement of the protein spontaneously inserted in a thiolipid monolayer in an aqueous solution. Geometric supramolecular arrays of nAChRs were observed, most commonly in a triangular form compatible with three nAChR dimers of ∼20 nm each. Addition of the full agonist carbamoylcholine activated and opened the nAChR ion channel, as revealed by the increase in capacitance relative to that of the nAChR-thiolipid system under basal conditions. Thus, the self-assembled system appears to be a viable biomimetic model to measure ionic conductance mediated by ion-gated ion channels under different experimental conditions, with potential applications in biotechnology and pharmacology.

Role of the Fourth Transmembrane α Helix in the Allosteric Modulation of Pentameric Ligand-Gated Ion Channels

The gating of pentameric ligand-gated ion channels is sensitive to a variety of allosteric modulators that act on structures peripheral to those involved in the allosteric pathway leading from the agonist site to the channel gate. One such structure, the lipid-exposed transmembrane α helix, M4, is the target of lipids, neurosteroids, and disease-causing mutations. Here we show that M4 interactions with the adjacent transmembrane α helices, M1 and M3, modulate pLGIC function. Enhanced M4 interactions promote channel function while ineffective interactions reduce channel function. The interface chemistry governs the intrinsic strength of M4-M1/M3 inter-helical interactions, both influencing channel gating and imparting distinct susceptibilities to the potentiating effects of a lipid-facing M4 congenital myasthenic syndrome mutation. Through aromatic substitutions, functional studies, and molecular dynamics simulations, we elucidate a mechanism by which M4 modulates channel function.

The M4 Transmembrane α-Helix Contributes Differently to Both the Maturation and Function of Two Prokaryotic Pentameric Ligand-gated Ion Channels

The role of the outermost transmembrane α-helix in both the maturation and function of the prokaryotic pentameric ligand-gated ion channels, GLIC and ELIC, was examined by Ala scanning mutagenesis, deletion mutations, and mutant cycle analyses. Ala mutations at the M4-M1/M3 interface lead to loss-of-function phenotypes in GLIC, with the largest negative effects occurring near the M4 C terminus. In particular, two aromatic residues at the M4 C terminus form a network of π-π and/or cation-π interactions with residues on M3 and the β6-β7 loop that is essential for both maturation and function. M4-M1/M3 interactions appear to be optimized in GLIC with even subtle structural changes at this interface leading to detrimental effects. In contrast, mutations along the M4-M1/M3 interface of ELIC typically lead to gain-of-function phenotypes, suggesting that these interactions in ELIC are not optimized for channel function. In addition, no cluster of interacting residues involving the M4 C terminus, M3, and the β6-β7 loop was found, suggesting that the M4 C terminus plays little role in ELIC maturation or function. This study shows that M4 makes distinct contributions to the maturation and gating of these two closely related homologs, suggesting that GLIC and ELIC exhibit divergent features of channel function.

Phylogenetic conservation of protein-lipid motifs in pentameric ligand-gated ion channels

Using the crosstalk between the nicotinic acetylcholine receptor (nAChR) and its lipid microenvironment as a paradigm, this short overview analyzes the occurrence of structural motifs which appear not only to be conserved within the nAChR family and contemporary eukaryotic members of the pentameric ligand-gated ion channel (pLGIC) superfamily, but also extend to prokaryotic homologues found in bacteria. The evolutionarily conserved design is manifested in: 1) the concentric three-ring architecture of the transmembrane region, 2) the occurrence in this region of distinct lipid consensus motifs in prokaryotic and eukaryotic pLGIC and 3) the key participation of the outer TM4 ring in conveying the influence of the lipid membrane environment to the middle TM1-TM3 ring and this, in turn, to the inner TM2 channel-lining ring, which determines the ion selectivity of the channel. The preservation of these constant structural-functional features throughout such a long phylogenetic span likely points to the successful gain-of-function conferred by their early acquisition. This article is part of a Special Issue entitled: Lipid-protein interactions.

Nicotinic acetylcholine receptor-lipid interactions: Mechanistic insight and biological function

Membrane lipids are potent modulators of the nicotinic acetylcholine receptor (nAChR) from Torpedo. Lipids influence nAChR function by both conformational selection and kinetic mechanisms, stabilizing varying proportions of activatable versus non-activatable conformations, as well as influencing the transitions between these conformational states. Of note, some membranes stabilize an electrically silent uncoupled conformation that binds agonist but does not undergo agonist-induced conformational transitions. The uncoupled nAChR, however, does transition to activatable conformations in relatively thick lipid bilayers, such as those found in lipid rafts. In this review, we discuss current understanding of lipid-nAChR interactions in the context of increasingly available high resolution structural and functional data. These data highlight different sites of lipid action, including the lipid-exposed M4 transmembrane α-helix. Current evidence suggests that lipids alter nAChR function by modulating interactions between M4 and the adjacent transmembrane α-helices, M1 and M3. These interactions have also been implicated in both the folding and trafficking of nAChRs to the cell surface. We review current mechanistic understanding of lipid-nAChR interactions, and highlight potential biological roles for lipid-nAChR interactions in modulating the synaptic response. This article is part of a Special Issue entitled: Lipid-protein interactions.

Intramembrane aromatic interactions influence the lipid sensitivities of pentameric ligand-gated ion channels

Although the Torpedo nicotinic acetylcholine receptor (nAChR) reconstituted into phosphatidylcholine (PC) membranes lacking cholesterol and anionic lipids adopts a conformation where agonist binding is uncoupled from channel gating, the underlying mechanism remains to be defined. Here, we examine the mechanism behind lipid-dependent uncoupling by comparing the propensities of two prokaryotic homologs, Gloebacter and Erwinia ligand-gated ion channel (GLIC and ELIC, respectively), to adopt a similar uncoupled conformation. Membrane-reconstituted GLIC and ELIC both exhibit folded structures in the minimal PC membranes that stabilize an uncoupled nAChR. GLIC, with a large number of aromatic interactions at the interface between the outermost transmembrane α-helix, M4, and the adjacent transmembrane α-helices, M1 and M3, retains the ability to flux cations in this uncoupling PC membrane environment. In contrast, ELIC, with a level of aromatic interactions intermediate between that of the nAChR and GLIC, does not undergo agonist-induced channel gating, although it does not exhibit the expected biophysical characteristics of the uncoupled state. Engineering new aromatic interactions at the M4-M1/M3 interface to promote effective M4 interactions with M1/M3, however, increases the stability of the transmembrane domain to restore channel function. Our data provide direct evidence that M4 interactions with M1/M3 are modulated during lipid sensing. Aromatic residues strengthen M4 interactions with M1/M3 to reduce the sensitivities of pentameric ligand-gated ion channels to their surrounding membrane environment.

The role of the M4 lipid-sensor in the folding, trafficking, and allosteric modulation of nicotinic acetylcholine receptors

With the availability of high resolution structural data, increasing attention has focused on the mechanisms by which drugs and endogenous compounds allosterically modulate nicotinic acetylcholine receptor (nAChR) function. Lipids are potent modulators of the nAChR from Torpedo. Membrane lipids influence nAChR function by both conformational selection and kinetic mechanisms, stabilizing varying proportions of pre-existing resting, open, desensitized, and uncoupled conformations, as well as influencing the transitions between these conformational states. Structural and functional data highlight a role for the lipid-exposed M4 transmembrane α-helix of each subunit in lipid sensing, and suggest that lipids influence gating by altering the binding of M4 to the adjacent transmembrane α-helices, M1 and M3. M4 has also been implicated in both the folding and trafficking of nAChRs to the cell surface, as well as in the potentiation of nAChR gating by neurosteroids. Here, we discuss the roles of M4 in the folding, trafficking, and allosteric modulation of nAChRs. We also consider the hypothesis that variable chemistry at the M4-M1/M3 transmembrane α-helical interface in different nAChR subunits governs the capacity for potentiation by activating lipids. This article is part of the Special Issue entitled 'The Nicotinic Acetylcholine Receptor: From Molecular Biology to Cognition'.

Structural sensitivity of a prokaryotic pentameric ligand-gated ion channel to its membrane environment

Although the activity of the nicotinic acetylcholine receptor (nAChR) is exquisitely sensitive to its membrane environment, the underlying mechanisms remain poorly defined. The homologous prokaryotic pentameric ligand-gated ion channel, Gloebacter ligand-gated ion channel (GLIC), represents an excellent model for probing the molecular basis of nAChR sensitivity because of its high structural homology, relative ease of expression, and amenability to crystallographic analysis. We show here that membrane-reconstituted GLIC exhibits structural and biophysical properties similar to those of the membrane-reconstituted nAChR, although GLIC is substantially more thermally stable. GLIC, however, does not possess the same exquisite lipid sensitivity. In particular, GLIC does not exhibit the same propensity to adopt an uncoupled conformation where agonist binding is uncoupled from channel gating. Structural comparisons provide insight into the chemical features that may predispose the nAChR to the formation of an uncoupled state.

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.

The position of the double bond in monounsaturated free fatty acids is essential for the inhibition of the nicotinic acetylcholine receptor

Free fatty acids (FFAs) are non-competitive antagonists of the nicotinic acetylcholine receptor (AChR). Their site of action is supposedly located at the lipid-AChR interface. To elucidate the mechanism involved in this antagonism, we studied the effect that FFAs with a single double-bond at different positions (ω6, ω9, ω11 and ω13 cis-18:1) have on different AChR properties. Electrophysiological studies showed that only two FFAs (ω6 and ω9) reduced the duration of the channel open-state. The briefest component of the closed-time distribution remained unaltered, suggesting that ω6 and ω9 behave as allosteric blockers. Fluorescence resonance energy transfer studies indicated that all FFAs locate at the lipid-AChR interface, ω6 being restricted to annular sites and all others occupying non-annular sites. The perturbation of the native membrane order by FFAs was evaluated by DPH (1,6-diphenyl-1,3,5-hexatriene) and Laurdan fluorescence polarization studies, with the greatest decrease observed for ω9 and ω11. AChR conformational changes produced by FFAs present at the lipid bilayer were evaluated by fluorescence quenching studies of pyrene-labeled AChR and also using the AChR conformational-sensitive probe crystal violet. All cis-FFAs produced AChR conformational changes at the transmembrane level, but only ω9, ω11 and ω13 perturbed the resting state. Thus, the position and isomerism of the torsion angle of unsaturated FFAs are probably a key factor in terms of AChR blockage, suggesting that FFAs with a unique cis double bond at a superficial position inside the membrane directly inhibit AChR function by perturbing a potential conserved core structure for AChR gating at that level.

Molecular mechanisms of acetylcholine receptor-lipid interactions: from model membranes to human biology

Lipids are potent modulators of the Torpedo nicotinic acetylcholine receptor. Lipids influence nicotinic receptor function by allosteric mechanisms, stabilizing varying proportions of pre-existing resting, open, desensitized, and uncoupled conformations. Recent structures reveal that lipids could alter function by modulating transmembrane α-helix/α-helix packing, which in turn could alter the conformation of the allosteric interface that links the agonist-binding and transmembrane pore domains-this interface is essential in the coupling of agonist binding to channel gating. We discuss potential mechanisms by which lipids stabilize different conformational states in the context of the hypothesis that lipid-nicotinic receptor interactions modulate receptor function at biological synapses.

Tryptophan scanning mutagenesis reveals distortions in the helical structure of the δM4 transmembrane domain of the Torpedo californica nicotinic acetylcholine receptor

The lipid-protein interface is an important domain of the nicotinic acetylcholine receptor (nAChR) that has recently garnered increased relevance. Several studies have made significant advances toward determining the structure and dynamics of the lipid-exposed domains of the nAChR. However, there is still a need to gain insight into the mechanism by which lipid-protein interactions regulate the function and conformational transitions of the nAChR. In this study, we extended the tryptophan scanning mutagenesis (TrpScanM) approach to dissect secondary structure and monitor the conformational changes experienced by the δM4 transmembrane domain (TMD) of the Torpedo californica nAChR, and to identify which positions on this domain are potentially linked to the regulation of ion channel kinetics. The difference in oscillation patterns between the closed- and open-channel states suggests a substantial conformational change along this domain as a consequence of channel activation. Furthermore, TrpScanM revealed distortions along the helical structure of this TMD that are not present on current models of the nAChR. Our results show that a Thr-Pro motif at positions 462-463 markedly bends the helical structure of the TMD, consistent with the recent crystallographic structure of the GluCl Cys-loop receptor which reveals a highly bent TMD4 in each subunit. This Thr-Pro motif acts as a molecular hinge that delineates two gating blocks in the δM4 TMD. These results suggest a model in which a hinge-bending motion that tilts the helical structure is combined with a spring-like motion during transition between the closed- and open-channel states of the δM4 TMD.

Discrimination of agonists versus antagonists of nicotinic ligands based on docking onto AChBP structures

Numerous high-resolution crystallographic structures of the acetylcholine binding protein (AChBP), a molluscan cholinergic protein, homologous to the extracellular domain of nicotinic acetylcholine receptors, are available. This offers opportunities to model the interaction between various ligands and the acetylcholine binding site. Herein we present a study of the interplay between ligand binding and motions of the C-loop capping the binding site. Nicotinic agonists and antagonists were docked on AChBP X-ray structures. It is shown that the studied agonists and antagonists can be discriminated according to their higher affinities for structures respectively obtained in the presence of agonists or antagonists, highlighting the fact that AChBP structures retain a pharmacological footprint of the compound used in crystallography experiments. A detailed analysis of the binding site cavities suggests that this property is mainly related to the shape of the cavities.

NMR structure of the transmembrane domain of the n-acetylcholine receptor beta2 subunit

Nicotinic acetylcholine receptors (nAChRs) are involved in fast synaptic transmission in the central and peripheral nervous system. Among the many different types of subunits in nAChRs, the beta2 subunit often combines with the alpha4 subunit to form alpha4beta2 pentameric channels, the most abundant subtype of nAChRs in the brain. Besides computational predictions, there is limited experimental data available on the structure of the beta2 subunit. Using high-resolution NMR spectroscopy, we solved the structure of the entire transmembrane domain (TM1234) of the beta2 subunit. We found that TM1234 formed a four-helix bundle in the absence of the extracellular and intracellular domains. The structure exhibited many similarities to those previously determined for the Torpedo nAChR and the bacterial ion channel GLIC. We also assessed the influence of the fourth transmembrane helix (TM4) on the rest of the domain. Although secondary structures and tertiary arrangements were similar, the addition of TM4 caused dramatic changes in TM3 dynamics and subtle changes in TM1 and TM2. Taken together, this study suggests that the structures of the transmembrane domains of these proteins are largely shaped by determinants inherent in their sequence, but their dynamics may be sensitive to modulation by tertiary and quaternary contacts.

Complexes between photoactivated rhodopsin and transducin: progress and questions

Activation of GPCRs (G-protein-coupled receptors) leads to conformational changes that ultimately initiate signal transduction. Activated GPCRs transiently combine with and activate heterotrimeric G-proteins resulting in GTP replacement of GDP on the G-protein alpha subunit. Both the detailed structural changes essential for productive GDP/GTP exchange on the G-protein alpha subunit and the structure of the GPCR-G-protein complex itself have yet to be elucidated. Nevertheless, transient GPCR-G-protein complexes can be trapped by nucleotide depletion, yielding an empty-nucleotide G-protein-GPCR complex that can be isolated. Whereas early biochemical studies indicated formation of a complex between G-protein and activated receptor only, more recent results suggest that G-protein can bind to pre-activated states of receptor or even couple transiently to non-activated receptor to facilitate rapid responses to stimuli. Efficient and reproducible formation of physiologically relevant, conformationally homogenous GPCR-G-protein complexes is a prerequisite for structural studies designed to address these possibilities.

Structural model for the binding sites of allosterically potentiating ligands on nicotinic acetylcholine receptors

Current treatments of Alzheimer's disease include the allosteric potentiation of nicotinic acetylcholine receptor (nAChR) response. The location of the binding site for allosteric potentiating ligands (APLs) within the receptor is not yet fully understood. Based on homology models for the ligand binding domain of human alpha7, human alpha4beta2, and chicken alpha7 receptors, as well as blind docking experiments with galanthamine, physostigmine, codeine, and 5HT, we identified T197 as an essential element of the APL binding site at the outer surface of the ligand binding domain (LBD) of nAChR. We also found the previously known galanthamine binding site in the region of K123 at the inside of the receptor funnel, which, however, was shown to not be part of the APL site. Our results are verified by site-directed mutagenesis and electrophysiological experiments, and suggest that APL and ACh bind to different sites on nicotinic receptors and that allosteric potentiation may arise from a direct interplay between both these sites.

Boundary lipids in the nicotinic acetylcholine receptor microenvironment

The structural and functional properties of the nicotinic acetylcholine receptor (AChR), the archetype molecule in the superfamily of Cys-looped ligand-gated ion channels, are strongly dependent on the lipids in the vicinal microenvironment. The influence on receptor properties is mainly exerted by the AChR-vicinal ("shell" or "annular") lipids, which occur in the liquid-ordered phase as opposed to the more disordered and "fluid" bulk membrane lipids. Fluorescence studies from our laboratory have identified discrete sites for fatty acids, phospholipids, and cholesterol on the AChR protein, and electron-spin resonance spectroscopy has enabled the establishment of the stoichiometry and selectivity of the shell lipid for the AChR and the disclosure of lipid sites in the AChR transmembrane region. Experimental evidence supports the notion that the interface between the protein moiety and the adjacent lipid shell is the locus of a variety of pharmacologically relevant processes, including the action of steroids and other lipids. I surmise that the outermost ring of M4 helices constitutes the boundary interface, most suitable to convey the signals from the lipid microenvironment to the rest of the transmembrane region, and to the channel inner ring in particular.

Molecular-dynamics simulations of ELIC-a prokaryotic homologue of the nicotinic acetylcholine receptor

The ligand-gated ion channel from Erwinia chrysanthemi (ELIC) is a prokaryotic homolog of the eukaryotic nicotinic acetylcholine receptor (nAChR) that responds to the binding of neurotransmitter acetylcholine and mediates fast signal transmission. ELIC is similar to the nAChR in its primary sequence and overall subunit organization, but despite their structural similarity, it is not clear whether these two ligand-gated ion channels operate in a similar manner. Further, it is not known to what extent mechanistic insights gleaned from the ELIC structure translate to eukaryotic counterparts such as the nAChR. Here we use molecular-dynamics simulations to probe the conformational dynamics and hydration of the transmembrane pore of ELIC. The results are compared with those from our previous simulation of the human alpha7 nAChR. Overall, ELIC displays increased stability compared to the nAChR, whereas the two proteins exhibit remarkable similarity in their global motion and flexibility patterns. The majority of the increased stability of ELIC does not stem from the deficiency of the models used in the simulations, and but rather seems to have a structural basis. Slightly altered dynamical correlation features are also observed among several loops within the membrane region. In sharp contrast to the nAChR, ELIC is completely dehydrated from the pore center to the extracellular end throughout the simulation. Finally, the simulation of an ELIC mutant substantiates the important role of F246 on the stability, hydration and possibly function of the ELIC channel.

Fluorescence and molecular dynamics studies of the acetylcholine receptor γM4 transmembrane peptide in reconstituted systems

A combination of fluorescence spectroscopy and molecular dynamics (MD) is applied to assess the conformational dynamics of a peptide making up the outermost ring of the nicotinic acetylcholine receptor (AChR) transmembrane region and the effect of membrane thickness and cholesterol on the hydrophobic matching of this peptide. The fluorescence studies exploit the intrinsic fluorescence of the only tryptophan residue in a synthetic peptide corresponding to the fourth transmembrane domain of the AChR gamma subunit (gammaM4-Trp(6)) reconstituted in lipid bilayers of varying thickness, and combine this information with quenching studies using depth-sensitive phosphatidylcholine spin-labeled probes and acrylamide, polarization of fluorescence, and generalized polarization of Laurdan. A direct correlation was found between bilayer width and the depth of insertion of Trp(6). We further extend our recent MD study of the conformational dynamics of the AChR channel to focus on the crosstalk between M4 and the lipid-belt region. The isolated gammaM4 peptide is shown to possess considerable orientational flexibility while maintaining a linear alpha-helical structure, and to vary its tilt depending on bilayer width and cholesterol (Chol) content. MD studies also show that gammaM4 also establishes contacts with the other TM peptides on its inner face, stabilizing a shorter TM length that is still highly sensitive to the lipid environment. In the native membrane the topology of the M4 ring is likely to exhibit a similar behavior, dynamically modifying its tilt to match the hydrophobic thickness of the bilayer.

Structure and dynamics of the γM4 transmembrane domain of the acetylcholine receptor in lipid bilayers: insights into receptor assembly and function

A 28-mer peptide (gammaM4) corresponding to the fourth transmembrane segment of the nicotinic acetylcholine receptor (AChR) gamma-subunit, with a single tryptophan residue (Trp6), was reconstituted into lipid bilayers of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), loaded with either high or low amounts of cholesterol, i.e., in the conjugated liquid-ordered and liquid-disordered phases, respectively, at room temperature. By making use of the Trp intrinsic fluorescence, both steady-state and time-resolved fluorescence techniques were employed, namely, red-edge excitation shift effect, decay-associated spectra (DAS), and time-resolved anisotropy. The results obtained here, together with previous studies on the same reconstituted peptide, indicate that: (i) Trp6 is strongly anchored in the bilayer with a defined transverse location; (ii) the modifications in the measured DAS are related to the complex result of a self-quenching process on the decay parameters; (iii) the wobbling movement of the indole moiety of Trp6 is fast but severely restricted in amplitude; and, (iv) in the liquid-ordered phase, the bilayer properties and the tilt angle of the peptide enhance peptide-peptide interactions, with the formation of peptide rich patches and possibly some anti-parallel helix-helix aggregates, showing different dynamics from that of the peptide in the liquid-disordered phase where the peptide is randomly distributed.

A butterfly effect: highly insecticidal resistance caused by only a conservative residue mutated of drosophila melanogaster acetylcholinesterase

Acetylcholinesterase (AChE) and its mutation recently emerged as a significant research area, due to its resistance against organophosphate and carbamate insecticides. Residue G265, which is always a conservative residue, mutated to A265 is the most frequent mutant of AChE in Drosophila populations. However, only this mutation caused a 'butterfly effect' that gives high insecticidal resistance. Herein, the models of sensitive strain (Dm-S) and the resistance strain (Dm-R) were constructed, to give a total of 2000 ps molecular dynamics simulation and to reveal the insecticidal resistance mechanism, with implied, the active gorge of Dm-R was much less flexible than that of Dm-S. The "back door" channel was widened to accelerate the detoxication against insecticides by the conformation changing of W83 and I161. All the distances (S238-H480, S238-G150, S238-G151, Y71-M153) in Dm-R became smaller than those in Dm-S, which may deeply influence the binding between the insecticides and DmAChE.

In silico models for the human alpha4beta2 nicotinic acetylcholine receptor

The neuronal alpha4beta2 nicotinic acetylcholine receptor (nAChR) is one of the most widely expressed nAChR subtypes in the brain. Its subunits have high sequence identity (54 and 46% for alpha4 and beta2, respectively) with alpha and beta subunits in Torpedo nAChR. Using the known structure of the Torpedo nAChR as a template, the closed-channel structure of the alpha4beta2 nAChR was constructed through homology modeling. Normal-mode analysis was performed on this closed structure and the resulting lowest frequency mode was applied to it for a "twist-to-open" motion, which increased the minimum pore radius from 2.7 to 3.4 A and generated an open-channel model. Nicotine could bind to the predicted agonist binding sites in the open-channel model but not in the closed one. Both models were subsequently equilibrated in a ternary lipid mixture via extensive molecular dynamics (MD) simulations. Over the course of 11 ns MD simulations, the open channel remained open with filled water, but the closed channel showed a much lower water density at its hydrophobic gate comprised of residues alpha4-V259 and alpha4-L263 and their homologous residues in the beta2 subunits. Brownian dynamics simulations of Na+ permeation through the open channel demonstrated a current-voltage relationship that was consistent with experimental data on the conducting state of alpha4beta2 nAChR. Besides establishment of the well-equilibrated closed- and open-channel alpha4beta2 structural models, the MD simulations on these models provided valuable insights into critical factors that potentially modulate channel gating. Rotation and tilting of TM2 helices led to changes in orientations of pore-lining residue side chains. Without concerted movement, the reorientation of one or two hydrophobic side chains could be enough for channel opening. The closed- and open-channel structures exhibited distinct patterns of electrostatic interactions at the interface of extracellular and transmembrane domains that might regulate the signal propagation of agonist binding to channel opening. A potential prominent role of the beta2 subunit in channel gating was also elucidated in the study.

Behavioral differences between phosphatidic acid and phosphatidylcholine in the presence of the nicotinic acetylcholine receptor

It has been found experimentally that negatively charged phosphatidic acid (PA) lipids and cholesterol molecules stabilize the nicotinic acetylcholine receptor (nAChR) in a functional resting state that can participate in an agonist-induced conformational change. In this study, we compare phosphatidylcholine (PC) and PA lipid behavior in the presence of the nAChR to determine why PC lipids do not support a functional nAChR. For lipids that are located within 1.0 nm of the protein, both PC and PA lipids have very similar order-parameter and bilayer-thickness values, which indicate that the annular lipid properties are protein-dependent. The most significant difference between the PC and PA bilayers is the formation of a lipid domain around the protein, which is visible in the PA bilayer but not the PC bilayer. This suggests that the PA domain may help stabilize the nAChR resting state. The PA lipids in the microdomain have a decreased order compared to a homogeneous PA bilayer and the lipids near the protein attempt to increase the free space in their vicinity by residing in multiple lateral planes.

Molecular dynamics simulations of Na+/Cl(-)-dependent neurotransmitter transporters in a membrane-aqueous system

We have performed molecular dynamics simulations of a homology model of the human serotonin transporter (hSERT) in a membrane environment and in complex with either the natural substrate 5-HT or the selective serotonin reuptake inhibitor escitalopram. We have also included a transporter homologue, the Aquifex aeolicus leucine transporter (LeuT), in our study to evaluate the applicability of a simple and computationally attractive membrane system. Fluctuations in LeuT extracted from simulations are in good agreement with crystallographic B factors. Furthermore, key interactions identified in the X-ray structure of LeuT are maintained throughout the simulations indicating that our simple membrane system is suitable for studying the transmembrane protein hSERT in complex with 5-HT or escitalopram. For these transporter complexes, only relatively small fluctuations are observed in the ligand-binding cleft. Specific interactions responsible for ligand recognition, are identified in the hSERT-5HT and hSERT-escitalopram complexes. Our findings are in good agreement with predictions from mutagenesis studies.

Mechanics of channel gating of the nicotinic acetylcholine receptor

The nicotinic acetylcholine receptor (nAChR) is a key molecule involved in the propagation of signals in the central nervous system and peripheral synapses. Although numerous computational and experimental studies have been performed on this receptor, the structural dynamics of the receptor underlying the gating mechanism is still unclear. To address the mechanical fundamentals of nAChR gating, both conventional molecular dynamics (CMD) and steered rotation molecular dynamics (SRMD) simulations have been conducted on the cryo-electron microscopy (cryo-EM) structure of nAChR embedded in a dipalmitoylphosphatidylcholine (DPPC) bilayer and water molecules. A 30-ns CMD simulation revealed a collective motion amongst C-loops, M1, and M2 helices. The inward movement of C-loops accompanying the shrinking of acetylcholine (ACh) binding pockets induced an inward and upward motion of the outer beta-sheet composed of beta9 and beta10 strands, which in turn causes M1 and M2 to undergo anticlockwise motions around the pore axis. Rotational motion of the entire receptor around the pore axis and twisting motions among extracellular (EC), transmembrane (TM), and intracellular MA domains were also detected by the CMD simulation. Moreover, M2 helices undergo a local twisting motion synthesized by their bending vibration and rotation. The hinge of either twisting motion or bending vibration is located at the middle of M2, possibly the gate of the receptor. A complementary twisting-to-open motion throughout the receptor was detected by a normal mode analysis (NMA). To mimic the pulsive action of ACh binding, nonequilibrium MD simulations were performed by using the SRMD method developed in one of our laboratories. The result confirmed all the motions derived from the CMD simulation and NMA. In addition, the SRMD simulation indicated that the channel may undergo an open-close (O <--> C) motion. The present MD simulations explore the structural dynamics of the receptor under its gating process and provide a new insight into the gating mechanism of nAChR at the atomic level.

A comparative molecular dynamics analysis of the amyloid beta-peptide in a lipid bilayer

Because the amyloid beta-peptide (Abeta) functions as approximately half of the transmembrane domain of the amyloid precursor protein and interaction of Abeta with membranes is proposed to result in neurotoxicity, the association of Abeta with membranes likely is important in the etiology of Alzheimer's disease. Atomic details of the interaction of Abeta with membranes are not accessible with most experimental techniques, but computational methods can provide this information. Here, we present the results of ten 100-ns molecular dynamics (MD) simulations of the 40-residue amyloid beta-peptide (Abeta40) embedded in a dipalmitoylphosphatidylcholine (DPPC) bilayer. The present study examines the effects of insertion depth, protonation state of key residues, and ionic strength on Abeta40 in a DPPC bilayer. In all cases, a portion of the peptide remained embedded in the bilayer. In the case of deeper insertion depth, Abeta40 adopted a near-transmembrane orientation, drawing water molecules into the bilayer to associate with its charged amino acids. In the case of shallower insertion, the most widely-accepted construct, the peptide associated strongly with the membrane-water interface and the phosphatidylcholine headgroups of the bilayer. In most cases, significant disordering of the extracellular segment of the peptide was observed, and the brief appearance of a beta-strand was noted in one case. Our results compare well with a variety of experimental and computational findings. From this study, we conclude that Abeta associated with membranes is dynamic and capable of adopting a number of conformations, each of which may have significance in understanding the progression of Alzheimer's disease.

Dynamics of heteropentameric nicotinic acetylcholine receptor: implications of the gating mechanism

The dynamics characteristics of the currently available structure of Torpedo nicotinic acetylcholine receptor (nAChR), including the extracellular, transmembrane, and intracellular domains (ICDs), were analyzed using the Gaussian Network Model (GNM) and Anisotropic Network Model (ANM). We found that a symmetric quaternary twist motion, reported previously in the literature in a homopentameric receptor (Cheng et al. J Mol Biol 2006;355:310-324; Taly et al. Biophys J 2005;88:3954-3965), occurred also in the heteropentameric Torpedo nAChR. We believe, however, that the symmetric twist alone is not sufficient to explain a large body of experimental data indicating asymmetry and subunit nonequivalence during gating. Here we report our results supporting the hypothesis that a combination of symmetric and asymmetric motions opens the gate. We show that the asymmetric motion involves tilting of the TM2 helices. Furthermore, our study reveals three additional aspects of channel dynamics: (1) loop A serves as an allosteric mediator between the ligand binding loops and those at the domain interface, particularly the linker between TM2 and TM3; (2) the ICD can modulate the pore dynamics and thus should not be neglected in gating studies; and (3) the F loops, which are peculiarly longer and poorly-conserved in non-alpha-subunits, have important dynamical implications.

Nanosecond-timescale conformational dynamics of the human alpha7 nicotinic acetylcholine receptor

We explore the conformational dynamics of a homology model of the human alpha7 nicotinic acetylcholine receptor using molecular dynamics simulation and analyses of root mean-square fluctuations, block partitioning of segmental motion, and principal component analysis. The results reveal flexible regions and concerted global motions of the subunits encompassing extracellular and transmembrane domains of the subunits. The most relevant motions comprise a bending, hinged at the beta10-M1 region, accompanied by concerted tilting of the M2 helices that widens the intracellular end of the channel. Despite the nanosecond timescale, the observations suggest that tilting of the M2 helices may initiate opening of the pore. The results also reveal direct coupling between a twisting motion of the extracellular domain and dynamic changes of M2. Covariance analysis of interresidue motions shows that this coupling arises through a network of residues within the Cys and M2-M3 loops where Phe135 is stabilized within a hydrophobic pocket formed by Leu270 and Ile271. The resulting concerted motion causes a downward shift of the M2 helices that disrupts a hydrophobic girdle formed by 9' and 13' residues.

Pharmacophore-directed homology modeling and molecular dynamics simulation of G protein-coupled receptor: study of possible binding modes of 5-HT2C receptor agonists

A new pharmacophore-based modeling procedure, including homology modeling, pharmacophore study, flexible molecular docking, and long-time molecular dynamics (MD) simulations, was employed to construct the structure of the human 5-HT(2C) receptor and determine the characteristics of binding modes of 5-HT(2C) receptor agonists. An agonist-receptor complex has been constructed based on homology modeling and a pharmacophore hypothesis model based on some high active compounds. Then MD simulations of the ligand-receptor complex in an explicit membrane environment were carried out. The conformation of the 5-HT(2C) receptor during MD simulation was explored, and the stable binding modes of the studied agonist were determined. Flexible molecular docking of several structurally diverse agonists of the human 5-HT(2C) receptor was carried out, and the general binding modes of these agonists were investigated. According to the models presented in this work and the results of Flexi-Dock, the involvement of the amino acid residues Asp134, Ser138, Asn210, Asn331, Tyr358, Ile131, Ser132, Val135,Thr139, Ile189, Val202, Val208, Leu209, Phe214, Val215, Gly218, Ser219, Phe223, Trp324, Phe327, and Phe328 in agonist recognition was studied. The obtained binding modes of the human 5-HT(2C) receptor agonists have good agreement with the site-directed mutagenesis data and other studies.

Understanding the regulation mechanisms of PAF receptor by agonists and antagonists: Molecular modeling and molecular dynamics simulation studies

Platelet-activating factor receptor (PAFR) is a member of G-protein coupled receptor (GPCR) superfamily. Understanding the regulation mechanisms of PAFR by its agonists and antagonists at the atomic level is essential for designing PAFR antagonists as drug candidates for treating PAF-mediated diseases. In this study, a 3D model of PAFR was constructed by a hierarchical approach integrating homology modeling, molecular docking and molecular dynamics (MD) simulations. Based on the 3D model, regulation mechanisms of PAFR by agonists and antagonists were investigated via three 8-ns MD simulations on the systems of apo-PAFR, PAFR-PAF and PAFR-GB. The simulations revealed that binding of PAF to PAFR triggers the straightening process of the kinked helix VI, leading to its activated state. In contrast, binding of GB to PAFR locks PAFR in its inactive state.

Structure and dynamics of phospholipid bilayers using recently developed general all-atom force fields

Two fully hydrated pure-species phospholipids bilayers, 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and 1,2-dioleoyl-sn-glycero-3-phosphorylcholine (DOPC), in the fluid phase and explicit solvent have been studied using molecular dynamics simulation. Atom interactions were modeled using recently developed force fields based on AMBER with full atomistic details. Several representative liquid phase properties for the structure and dynamics of lipids with different length of hydrocarbon chains and different level of saturation have been reproduced without artificially biasing the system in order to match experimental data. In particular, as the new GAFF (General Amber Force Field) has not been explicitly developed to reproduce lipid characteristics and is naturally compatible with standard AMBER nucleic acids and proteins parameters, it is here proven a promising tool to study mixed lipid-protein processes as protein activity dependence on membrane composition, permeation of solute across membranes, and other cellular processes.

Blocking of the Nicotinic Acetylcholine Receptor Ion Channel by Chlorpromazine, a Noncompetitive Inhibitor: A Molecular Dynamics Simulation Study

A large series of pharmacological agents, distinct from the typical competitive antagonists, block in a noncompetitive manner the permeability response of the nicotinic acetylcholine receptor (nAChR) to the neurotransmitter acetylcholine. Taking the neuroleptic chlorpromazine (CPZ) as an example of such agents, the blocking mechanism of noncompetitive inhibitors to the ion channel pore of the nAChR has been explored at the atomic level using both conventional and steered molecular dynamics (MD) simulations. Repeated steered MD simulations have permitted calculation of the free energy (approximately 36 kJ/mol) of CPZ binding and identification of the optimal site in the region of the serine and leucine rings, at approximately 4 A from the pore entrance. Coulomb and the Lennard-Jones interactions between CPZ and the ion channel as well as the conformational fluctuations of CPZ were examined to assess the contribution of each to the binding of CPZ to the nAChR. The MD simulations disclose a dynamic interaction of CPZ binding to the nAChR ionic channel. The cationic ammonium head of CPZ forms strong hydrogen bonds with Glu262 (alpha), Asp268 (beta), Glu272 (beta), Ser276 (beta), Glu280 (delta), Gln271 (gamma), Glu275 (gamma), and Asn279 (gamma) nAChR residues. Finally, the conventional MD simulation of CPZ at its identified binding site demonstrates that the binding of CPZ not only blocks ion transport through the channel but also markedly inhibits the conformational transitions of the channel, necessary for nAChR to carry out its biological function.

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

The initial coupling between ligand binding and channel gating in the human α7 nicotinic acetylcholine receptor (nAChR) has been investigated with targeted molecular dynamics (TMD) simulation. During the simulation, eight residues at the tip of the C-loop in two alternating subunits were forced to move toward a ligand-bound conformation as captured in the crystallographic structure of acetylcholine binding protein (AChBP) in complex with carbamoylcholine. Comparison of apo- and ligand-bound AChBP structures shows only minor rearrangements distal from the ligand-binding site. In contrast, comparison of apo and TMD simulation structures of the nAChR reveals significant changes toward the bottom of the ligand-binding domain. These structural rearrangements are subsequently translated to the pore domain, leading to a partly open channel within 4 ns of TMD simulation. Furthermore, we confirmed that two highly conserved residue pairs, one located near the ligand-binding pocket (Lys145 and Tyr188), and the other located toward the bottom of the ligand-binding domain (Arg206 and Glu45), are likely to play important roles in coupling agonist binding to channel gating. Overall, our simulations suggest that gating movements of the α7 receptor may involve relatively small structural changes within the ligand-binding domain, implying that the gating transition is energy-efficient and can be easily modulated by agonist binding/unbinding.

Nicotinic acetylcholine receptors are ligand-gated ion channels responsible for neurotransmitter-mediated signal transduction at synapses throughout the central and peripheral nervous systems. Binding of neurotransmitter molecules to subunit interfaces in the N-terminal extracellular domain induces structural rearrangements of the membrane-spanning domain permitting the influx of cations. A full understanding of how the conformational changes propagate from the ligand-binding site to the pore domain is of great interest to biologists, yet remains to be established. Using a special simulation technique known as targeted molecular dynamics, Cheng and colleagues probed the early stages of ligand-induced conformational rearrangements that may lead to channel opening. During the simulation, Cheng et al. observed a sequence of conformational changes that stem from the ligand-binding site to the transmembrane domain resulting in a wider channel. From these results, they suggest that gating movements may entail only small structural changes in the ligand-binding domain, implying that channel gating is energy-efficient and can readily be modulated by the binding/unbinding of agonist molecules.