NMR resolved multiple anesthetic binding sites in the TM domains of the α4ß2 nAChR.

The α4ß2 nicotinic acetylcholine receptor (nAChR) has significant roles in nervous system function and disease. It is also a molecular target of general anesthetics. Anesthetics inhibit the α4ß2 nAChR at clinically relevant concentrations, but their binding sites in α4ß2 remain unclear. The recently determined NMR structures of the α4ß2 nAChR transmembrane (TM) domains provide valuable frameworks for identifying the binding sites. In this study, we performed solution NMR experiments on the α4ß2 TM domains in the absence and presence of halothane and ketamine. Both anesthetics were found in an intra-subunit cavity near the extracellular end of the ß2 transmembrane helices, homologous to a common anesthetic binding site observed in X-ray structures of anesthetic-bound GLIC (Nury et al., [32]). Halothane, but not ketamine, was also found in cavities adjacent to the common anesthetic site at the interface of α4 and ß2. In addition, both anesthetics bound to cavities near the ion selectivity filter at the intracellular end of the TM domains. Anesthetic binding induced profound changes in protein conformational exchanges. A number of residues, close to or remote from the binding sites, showed resonance signal splitting from single to double peaks, signifying that anesthetics decreased conformation exchange rates. It was also evident that anesthetics shifted population of two conformations. Altogether, the study comprehensively resolved anesthetic binding sites in the α4ß2 nAChR. Furthermore, the study provided compelling experimental evidence of anesthetic-induced changes in protein dynamics, especially near regions of the hydrophobic gate and ion selectivity filter that directly regulate channel functions.

NMR structures of the transmembrane domains of the α4ß2 nAChR.

The α4ß2 nicotinic acetylcholine receptor (nAChR) is the predominant heteromeric subtype of nAChRs in the brain, which has been implicated in numerous neurological conditions. The structural information specifically for the α4ß2 and other neuronal nAChRs is presently limited. In this study, we determined structures of the transmembrane (TM) domains of the α4 and ß2 subunits in lauryldimethylamine-oxide (LDAO) micelles using solution NMR spectroscopy. NMR experiments and size exclusion chromatography-multi-angle light scattering (SEC-MALS) analysis demonstrated that the TM domains of α4 and ß2 interacted with each other and spontaneously formed pentameric assemblies in the LDAO micelles. The Na(+) flux assay revealed that α4ß2 formed Na(+) permeable channels in lipid vesicles. Efflux of Na(+) through the α4ß2 channels reduced intra-vesicle Sodium Green™ fluorescence in a time-dependent manner that was not observed in vesicles without incorporating α4ß2. The study provides structural insight into the TM domains of the α4ß2 nAChR. It offers a valuable structural framework for rationalizing extensive biochemical data collected previously on the α4ß2 nAChR and for designing new therapeutic modulators.

Asymmetric ligand binding facilitates conformational transitions in pentameric ligand-gated ion channels.

The anesthetic propofol inhibits the currents of the homopentameric ligand-gated ion channel GLIC, yet the crystal structure of GLIC with five propofol molecules bound symmetrically shows an open-channel conformation. To address this dilemma and determine if the symmetry of propofol binding sites affects the channel conformational transition, we performed a total of 1.5 µs of molecular dynamics simulations for different GLIC systems with propofol occupancies of 0, 1, 2, 3, and 5. GLIC without propofol binding or with five propofol molecules bound symmetrically, showed similar channel conformation and hydration status over multiple replicates of 100-ns simulations. In contrast, asymmetric binding to one, two or three equivalent sites in different subunits accelerated the channel dehydration, increased the conformational heterogeneity of the pore-lining TM2 helices, and shifted the lateral and radial tilting angles of TM2 toward a closed-channel conformation. The results differentiate two groups of systems based on the propofol binding symmetry. The difference between symmetric and asymmetric groups is correlated with the variance in the propofol-binding cavity adjacent to the hydrophobic gate and the force imposed by the bound propofol. Asymmetrically bound propofol produced greater variance in the cavity size that could further elevate the conformation heterogeneity. The force trajectory generated by propofol in each subunit over the course of a simulation exhibits an ellipsoidal shape, which has the larger component tangential to the pore. Asymmetric propofol binding creates an unbalanced force that expedites the channel conformation transitions. The findings from this study not only suggest that asymmetric binding underlies the propofol functional inhibition of GLIC, but also advocate for the role of symmetry breaking in facilitating channel conformational transitions.

Isoflurane alters the structure and dynamics of GLIC.

Pentameric ligand-gated ion channels are targets of general anesthetics. Although the search for discrete anesthetic binding sites has achieved some degree of success, little is known regarding how anesthetics work after the events of binding. Using the crystal structures of the bacterial Gloeobacter violaceus pentameric ligand-gated ion channel (GLIC), which is sensitive to a variety of general anesthetics, we performed multiple molecular dynamics simulations in the presence and absence of the general anesthetic isoflurane. Isoflurane bound to several locations within GLIC, including the transmembrane pocket identified crystallographically, the extracellular (EC) domain, and the interface of the EC and transmembrane domains. Isoflurane also entered the channel after the pore was dehydrated in one of the simulations. Isoflurane disrupted the quaternary structure of GLIC, as evidenced in a striking association between the binding and breakage of intersubunit salt bridges in the EC domain. The pore-lining helix experienced lateral and inward radial tilting motion that contributed to the channel closure. Isoflurane binding introduced strong anticorrelated motions between different subunits of GLIC. The demonstrated structural and dynamical modulations by isoflurane aid in the understanding of the underlying mechanism of anesthetic inhibition of GLIC and possibly other homologous pentameric ligand-gated ion channels.

General anesthetic binding to neuronal alpha4beta2 nicotinic acetylcholine receptor and its effects on global dynamics.

The neuronal alpha4beta2 nicotinic acetylcholine receptor (nAChR) is a target for general anesthetics. Currently available experimental structural information is inadequate to understand where anesthetics bind and how they modulate the receptor motions essential to function. Using our published open-channel structure model of alpha4beta2 nAChR, we identified and evaluated six amphiphilic interaction sites for the volatile anesthetic halothane via flexible ligand docking and subsequent 20-ns molecular dynamics simulations. Halothane binding energies ranged from -6.8 to -2.4 kcal/mol. The primary binding sites were located at the interface of extracellular and transmembrane domains, where halothane perturbed conformations of, and widened the gap among, the Cys loop, the beta1-beta2 loop, and the TM2-TM3 linker. The halothane with the highest binding affinity at the interface between the alpha4 and beta2 subunits altered interactions between the protein and nearby lipids by competing for hydrogen bonds. Gaussian network model analyses of the alpha4beta2 nAChR structures at the end of 20-ns simulations in the absence or presence of halothane revealed profound changes in protein residue mobility. The concerted motions critical to protein function were also perturbed considerably. Halothane's effect on protein dynamics was not confined to the residues adjacent to the binding sites, indicating an action on a more global scale.

Molecular dynamics simulations of ternary membrane mixture: phosphatidylcholine, phosphatidic acid, and cholesterol.

A ternary mixture of 1-palmitoyl-2-oleoyl phosphatidylcholine (POPC), 1-palmitoyl-2-oleoyl phosphatidic acid (POPA), and cholesterol (CHOL) works effectively for a functional conformation of nicotinic acetylcholine receptor (nAChR) that can undergo agonist-induced conformation changes, but POPC alone can stabilize only a desensitized state of nAChR. To gain insights into the lipid mixture that has strong impact to nAChR functions, we performed more than 50 ns all atom molecular dynamic (MD) simulations at 303 K on a fully hydrated bilayer consisting of 240 POPC, 80 POPA, and 80 CHOL (3:1:1). The MD simulation revealed various interactions between different types of molecular pairs that ultimately regulated lipid organization. The heterogeneous interactions among three different constituents resulted in a broad spectrum of lipid properties, including extensive distributions of average area per lipid and varied lipid ordering as a function of lipid closeness to CHOL. Higher percentage of POPA than POPC had close association with CHOL, which coincided with relatively higher ordering of POPA molecules in their acyl chains near lipid head groups. Lower fraction of gauche dihedrals was also found in the same region of POPA. Although the CHOL molecules had the effects on the enhancement of surrounding lipid order, relatively low lipid order parameters and high fraction of gauche bonds were observed in the ternary mixture. Collectively, these results suggest that the dynamical structure of the ternary system could be determinant for a functional nAChR.

Modeling the effects of torsional disorder on the spectra of poly- and oligo-(p-phenyleneethynylenes).

The absorption spectra of phenyleneethynylene oligomers show an unusual change in shape with oligomer length. The unusual aspects of the spectra arise from rotation of the phenylene rings about the long axis of the oligomer. In the ground electronic state, the barrier to this rotation is low and the spectra in room temperature come from an ensemble of different structures. In the excited electronic state, the barrier to rotation is substantially higher, giving rise to strong nonlinear electron-phonon coupling. A multidimensional semiempirical model that includes these effects is developed for the photophysics of phenyleneethynylene oligomers. The ground-state energy is modeled with a molecular mechanics expression, and the excitation energy is modeled with an exciton model. Intermediate Neglect of Differential Overlap/Singles Configuration Interaction (INDO/SCI) calculations verify the exciton model and provide initial estimates of the model parameters. These parameters generate the qualitative features seen in experimental spectra. Inclusion of entropy effects from the multiple torsional coordinates is essential. Refinement of the parameters yields quantitative agreement with experiment. This agreement shows that coupling to torsional motion is a major factor in the spectroscopy and photophysics of these conjugated polymers.

Mechanistic insights into xenon inhibition of NMDA receptors from MD simulations.

Inhibition of N-methyl-D-aspartate (NMDA) receptors has been viewed as a primary cause of xenon anesthesia, yet the mechanism is unclear. Here, we investigated interactions between xenon and the ligand-binding domain (LBD) of a NMDA receptor and examined xenon-induced structural and dynamical changes that are relevant to functional changes of the NMDA receptor. Several comparative molecular dynamics simulations were performed on two X-ray structures representing the open- and closed-cleft LBD of the NMDA receptor. We identified plausible xenon action sites in the LBD, including those nearby agonist sites, in the hinge region, and at the interface between two subunits. The xenon-binding energy varies from -5.3 to -0.7 kcal/mol. Xenon's effect on the NMDA receptor is conformation-dependent and is produced through both competitive and noncompetitive mechanisms. Xenon can promote cleft opening in the absence of agonists and consequently stabilizes the closed channel. Xenon can also bind at the interface of two subunits, alter the intersubunit interaction, and lead to a reduction of the distance between two GT linkers. This reduction corresponds to a rearrangement of the channel toward a direction of pore size decreasing, implying a closed or desensitized channel. In addition to these noncompetitive actions, xenon was found to weaken the glutamate binding, which could lead to low agonist efficacy and appear as competitive inhibition.

Higher susceptibility to halothane modulation in open- than in closed-channel alpha4beta2 nAChR revealed by molecular dynamics simulations.

The neuronal alpha4beta2 nicotinic acetylcholine receptor (nAChR) is a potential molecular target for general anesthetics. It is unclear, however, whether anesthetic action produces the same effect on the open and closed channels. Computations parallel to our previous open channel study (J. Phys. Chem. B 2009, 113, 12581) were performed on the closed-channel alpha4beta2 nAChR to investigate the conformation-dependent anesthetic effects on channel structures and dynamics. Flexible ligand docking and over 20 ns molecular dynamics simulations revealed similar halothane-binding sites in the closed and open channels. The sites with relatively high binding affinities (approximately -6.0 kcal/mol) were identified at the interface of extracellular (EC) and transmembrane (TM) domains or at the interface between alpha4 and beta2 subunits. Despite similar sites for halothane binding, the closed-channel conformation showed much less sensitivity than the open channel to the structural and dynamical perturbations from halothane. Compared to the systems without anesthetics, the amount of water inside the pore decreased by 22% in the presence of halothane in the open channel but only by 6% in the closed channel. Comparison of the nonbonded interactions at the EC/TM interfaces suggested that the beta2 subunits were more prone than the alpha4 subunits to halothane binding. In addition, our data support the notion that halothane exerts its effect by disturbing the quaternary structure and dynamics of the channel. The study concludes that sensitivity and global dynamics responsiveness of alpha4beta2 nAChR to halothane are conformation dependent. The effect of halothane on the global dynamics of the open-channel conformation might also account for the action of other inhaled general anesthetics.

Unresponsive correlated motion in alpha7 nAChR to halothane binding explains its functional insensitivity to volatile anesthetics.

Neuronal nicotinic acetylcholine receptors (nAChRs) have been implicated as targets for general anesthetics, but the functional responses to anesthetic modulation vary considerably among different subtypes of nAChRs. Inhaled general anesthetics, such as halothane, could effectively inhibit the channel activity of the alpha4beta2 nAChR but not the homologous alpha7 nAChR. To understand why alpha7 is insensitive to inhaled general anesthetics, we performed multiple sets of 20 ns molecular dynamics (MD) simulations on the closed- and open-channel alpha7 in the absence and presence of halothane and critically compared the results with those from our studies on the alpha4beta2 nAChR (Liu et al. J. Phys. Chem. B 2009, 113, 12581 and Liu et al. J. Phys. Chem. B 2010, 114, 626). Several halothane binding sites with fairly high binding affinities were identified in both closed- and open-channel alpha7, suggesting that a lack of sensitive functional responses of the alpha7 nAChR to halothane in the previous experiments was unlikely due to a lack of halothane interaction with alpha7. The binding affinities of halothane in alpha7 seemed to be protein conformation-dependent. Overall, halothane affinity was higher in the closed-channel alpha7. Halothane binding to alpha7 did not induce profound changes in alpha7 structure and dynamics that could be related to the channel function. In contrast, correlated motion of the open-channel alpha4beta2 was reduced substantially in the presence of halothane, primarily due to the more susceptible nature of beta2 to anesthetic modulation. The amphiphilic extracellular and transmembrane domain interface of the beta2 subunit is attractive to halothane and is susceptible to halothane perturbation, which may be why alpha4beta2 is functionally more sensitive to halothane than alpha7.