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

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.

Modulation of a rapid neurotransmitter receptor-ion channel by membrane lipids

Membrane lipids modulate the proteins embedded in the bilayer matrix by two non-exclusive mechanisms: direct or indirect. The latter comprise those effects mediated by the physicochemical state of the membrane bilayer, whereas direct modulation entails the more specific regulatory effects transduced via recognition sites on the target membrane protein. The nicotinic acetylcholine receptor (nAChR), the paradigm member of the pentameric ligand-gated ion channel (pLGIC) superfamily of rapid neurotransmitter receptors, is modulated by both mechanisms. Reciprocally, the nAChR protein exerts influence on its surrounding interstitial lipids. Folding, conformational equilibria, ligand binding, ion permeation, topography, and diffusion of the nAChR are modulated by membrane lipids. The knowledge gained from biophysical studies of this prototypic membrane protein can be applied to other neurotransmitter receptors and most other integral membrane proteins.

Simple Does Not Mean Trivial: Behavior of Phosphatidic Acid in Lipid Mono- and Bilayers

Phosphatidic acid (PA) is one of the simplest membrane phospholipids, yet it plays a crucial role in various biologically relevant processes that take place in cells. Since PA generation may be triggered by a variety of factors, very often of antagonistic character, the specific nature of physiological responses driven by PA is not clear. In order to shed more light on these issues, we carried out a systematic characterization of membranes containing one of the three biologically significant PA molecular species. The effect of these molecules on the properties of membranes composed of phosphatidylcholine and/or cholesterol was assessed in a multidisciplinary approach, including molecular dynamic simulations, flicker noise spectroscopy, and Langmuir monolayer isotherms. The first enables the determination of various macroscopic and microscopic parameters such as lateral diffusion, membrane thickness, and defect analysis. The obtained data revealed a strong interaction between unsaturated PA species and phosphatidylcholine. On the other hand, the behavior of saturated PA was greatly influenced by cholesterol. Additionally, a strong effect on mechanical properties was observed in the case of three-component systems, which could not be explained by the simple extrapolation of parameters of the corresponding two-component systems. Our data show that various PA species are not equivalent in terms of their influence on lipid mono- and bilayers and that membrane composition/properties, particularly those related to the presence of cholesterol, may strongly modulate PA behavior.

Lipid-modulation of membrane insertion and refolding of the apoptotic inhibitor Bcl-xL

Bcl-xL is a member of the Bcl-2 family of apoptotic regulators, responsible for inhibiting the permeabilization of the mitochondrial outer membrane, and a promising anti-cancer target. Bcl-xL exists in the following conformations, each believed to play a role in the inhibition of apoptosis: (a) a soluble folded conformation, (b) a membrane-anchored (by its C-terminal α8 helix) form, which retains the same fold as in solution and (c) refolded membrane-inserted conformations, for which no structural data are available. Previous studies established that in the cell Bcl-xL exists in a dynamic equilibrium between soluble and membranous states, however, no direct evidence exists in support of either anchored or inserted conformation of the membranous state in vivo. In this in vitro study, we employed a combination of fluorescence and EPR spectroscopy to characterize structural features of the bilayer-inserted conformation of Bcl-xL and the lipid modulation of its membrane insertion transition. Our results indicate that the core hydrophobic helix α6 inserts into the bilayer without adopting a transmembrane orientation. This insertion disrupts the packing of Bcl-xL and releases the regulatory N-terminal BH4 domain (α1) from the rest of the protein structure. Our data demonstrate that both insertion and refolding of Bcl-xL are modulated by lipid composition, which brings the apparent pK_a of insertion to the threshold of physiological pH. We hypothesize that conformational rearrangements associated with the bilayer insertion of Bcl-xL result in its switching to a so-called non-canonical mode of apoptotic inhibition. Presented results suggest that the alteration in lipid composition before and during apoptosis can serve as an additional factor regulating the permeabilization of the mitochondrial outer membrane.

Boundary lipids of the nicotinic acetylcholine receptor: Spontaneous partitioning via coarse-grained molecular dynamics simulation

Reconstituted nicotinic acetylcholine receptors (nAChRs) exhibit significant gain-of-function upon addition of cholesterol to reconstitution mixtures, and cholesterol affects the organization of nAChRs within domain-forming membranes, but whether nAChR partitions to cholesterol-rich liquid-ordered ("raft" or l_o) domains or cholesterol-poor liquid-disordered (l_do) domains is unknown. We use coarse-grained molecular dynamics simulations to observe spontaneous interactions of cholesterol, saturated lipids, and polyunsaturated (PUFA) lipids with nAChRs. In binary Dipalmitoylphosphatidylcholine:Cholesterol (DPPC:CHOL) mixtures, both CHOL and DPPC acyl chains were observed spontaneously entering deep "non-annular" cavities in the nAChR TMD, particularly at the subunit interface and the β subunit center, facilitated by the low amino acid density in the cryo-EM structure of nAChR in a native membrane. Cholesterol was highly enriched in the annulus around the TMD, but this effect extended over (at most) 5-10 Å. In domain-forming ternary mixtures containing PUFAs, the presence of a single receptor did not significantly affect the likelihood of domain formation. nAChR partitioned to any cholesterol-poor l_do domain that was present, regardless of whether the l_do or l_o domain lipids had PC or PE headgroups. Enrichment of PUFAs among boundary lipids was positively correlated with their propensity for demixing from cholesterol-rich phases. Long n-3 chains (tested here with Docosahexaenoic Acid, DHA) were highly enriched in annular and non-annular embedded sites, partially displacing cholesterol and completely displacing DPPC, and occupying sites even deeper within the bundle. Shorter n-6 chains were far less effective at displacing cholesterol from non-annular sites.

A novel identification approach for discovery of 5-HydroxyTriptamine 2A antagonists: combination of 2D/3D similarity screening, molecular docking and molecular dynamics

5-HydroxyTriptamine 2A antagonists are potential targets for treatment of various cerebrovascular and cardiovascular disorders. In this study, we have developed and performed a unique screening pipeline for filtering ZINC database compounds on the basis of similarities to known antagonists to determine novel small molecule antagonists of 5-HydroxyTriptamine 2A. The screening pipeline is based on 2D similarity, 3D dissimilarity and a combination of 2D/3D similarity. The shortlisted compounds were docked to a 5-HydroxyTriptamine 2A homology-based model, and complexes with low binding energies (287 complexes) were selected for molecular dynamics (MD) simulations in a lipid bilayer. The MD simulations of the shortlisted compounds in complex with 5-HydroxyTriptamine 2A confirmed the stability of the complexes and revealed novel interaction insights. The receptor residues S239, N343, S242, S159, Y370 and D155 predominantly participate in hydrogen bonding. π-π stacking is observed in F339, F340, F234, W151 and W336, whereas hydrophobic interactions are observed amongst V156, F339, F234, V362, V366, F340, V235, I152 and W151. The known and potential antagonists shortlisted by us have similar overlapping molecular interaction patterns. The 287 potential 5-HydroxyTriptamine 2A antagonists may be experimentally verified.

Cholesterol affects C₆₀ translocation across lipid bilayers

Cholesterol plays an important role in regulating the structural properties of phospholipid membranes and further influences the permeability of molecules and nanoparticles. However, nanoparticles' translocation across phospholipid membranes in the presence of cholesterol on the molecular scale is rarely studied. Here, we performed coarse-grained molecular dynamics simulations to probe the translocation of C60, one of the most popular nanoparticles, across dipalmitoylphosphatidylcholine bilayers with different concentrations of cholesterol molecules (0-50 mol%). The results reveal that the presence of cholesterol molecules induces lower area per lipid, larger bilayer thickness, and more ordered orientation of lipid tails. The higher the concentration of cholesterol molecules, the more significant is the condensing effect of lipid bilayer as just mentioned. Besides, dynamic processes, free energy profiles and permeability coefficients further indicate that the permeability of C60 decreases with increasing cholesterol concentration, which can be explained by the condensation effect and reduced free volume. Our researches provide an explicit description of the impact of cholesterol on C60 translocation across lipid bilayers.

Binding of β-amyloid (1-42) peptide to negatively charged phospholipid membranes in the liquid-ordered state: modeling and experimental studies

To explore the initial stages of amyloid β peptide (Aβ42) deposition on membranes, we have studied the interaction of Aβ42 in the monomeric form with lipid monolayers and with bilayers in either the liquid-disordered or the liquid-ordered (L(o)) state, containing negatively charged phospholipids. Molecular dynamics (MD) simulations of the system have been performed, as well as experimental measurements. For bilayers in the L(o) state, in the absence of the negatively charged lipids, interaction is weak and it cannot be detected by isothermal calorimetry. However, in the presence of phosphatidic acid, or of cardiolipin, interaction is detected by different methods and in all cases interaction is strongest with lower (2.5-5 mol%) than higher (10-20 mol%) proportions of negatively charged phospholipids. Liquid-disordered bilayers consistently allowed a higher Aβ42 binding than L(o) ones. Thioflavin T assays and infrared spectroscopy confirmed a higher proportion of β-sheet formation under conditions when higher peptide binding was measured. The experimental results were supported by MD simulations. We used 100 ns MD to examine interactions between Aβ42 and three different 512 lipid bilayers consisting of palmitoylsphingomyelin, dimyristoyl phosphatidic acid, and cholesterol in three different proportions. MD pictures are different for the low- and high-charge bilayers, in the former case the peptide is bound through many contact points to the bilayer, whereas for the bilayer containing 20 mol% anionic phospholipid only a small fragment of the peptide appears to be bound. The MD results indicate that the binding and fibril formation on the membrane surface depends on the composition of the bilayer, and is the result of a subtle balance of many inter- and intramolecular interactions between the Aβ42 and membrane.

Cannabinoids suppress inflammatory and neuropathic pain by targeting α3 glycine receptors

Systemic and intrathecal administration of derivatives of a nonpsychoactive component of marijuana significantly suppresses chronic inflammatory and neuropathic pain, without causing analgesic tolerance, in several rodent models.

Certain types of nonpsychoactive cannabinoids can potentiate glycine receptors (GlyRs), an important target for nociceptive regulation at the spinal level. However, little is known about the potential and mechanism of glycinergic cannabinoids for chronic pain treatment. We report that systemic and intrathecal administration of cannabidiol (CBD), a major nonpsychoactive component of marijuana, and its modified derivatives significantly suppress chronic inflammatory and neuropathic pain without causing apparent analgesic tolerance in rodents. The cannabinoids significantly potentiate glycine currents in dorsal horn neurons in rat spinal cord slices. The analgesic potency of 11 structurally similar cannabinoids is positively correlated with cannabinoid potentiation of the α3 GlyRs. In contrast, the cannabinoid analgesia is neither correlated with their binding affinity for CB1 and CB2 receptors nor with their psychoactive side effects. NMR analysis reveals a direct interaction between CBD and S296 in the third transmembrane domain of purified α3 GlyR. The cannabinoid-induced analgesic effect is absent in mice lacking the α3 GlyRs. Our findings suggest that the α3 GlyRs mediate glycinergic cannabinoid-induced suppression of chronic pain. These cannabinoids may represent a novel class of therapeutic agents for the treatment of chronic pain and other diseases involving GlyR dysfunction.

Molecular interactions of the physiological anti-hypertensive peptide catestatin with the neuronal nicotinic acetylcholine receptor

Catestatin (CST), a chromogranin-A-derived peptide, is a potent endogenous inhibitor of the neuronal nicotinic acetylcholine receptor (nAChR). It exerts an anti-hypertensive effect by acting as a 'physiological brake' on transmitter release into the circulation. However, the mechanism of interaction of CST with nAChR is only partially understood. To unravel molecular interactions of the wild-type human CST (CST-WT) as well as its naturally occurring variants (CST-364S and CST-370L, which have Gly→Ser and Pro→Leu substitutions, respectively) with the human α3β4 nAChR, we generated a homology-modeled human α3β4 nAChR structure and solution structures of CST peptides. Docking and molecular dynamics simulations showed that ~90% of interacting residues were within 15 N-terminal residues of CST peptides. The rank order of binding affinity of these peptides with nAChR was: CST-370L>CST-WT>CST-364S; the extent of occlusion of the receptor pore by these peptides was also in the same order. In corroboration with computational predictions, circular dichroism analysis revealed significant differences in global structures of CST peptides (e.g. the order of α-helical content was: CST-370L>CST-WT>CST-364S). Consistently, CST peptides blocked various stages of nAChR signal transduction, such as nicotine- or acetylcholine-evoked inward current, rise in intracellular Ca(2+) and catecholamine secretion in or from neuron-differentiated PC12 cells, in the same rank order. Taken together, this study shows molecular interactions between human CST peptides and human α3β4 nAChR, and demonstrates that alterations in the CST secondary structure lead to the gain of potency for CST-370L and loss of potency for CST-364S. These findings have implications for understanding the nicotinic cholinergic signaling in humans.

The role of structured water in mediating general anesthetic action on alpha4beta2 nAChR

Water is an essential component for many biological processes. Pauling proposed that water might play a critical role in general anesthesia by forming water clathrates around anesthetic molecules. To examine potential involvement of water in general anesthesia, we analyzed water within alpha4beta2 nAChR, a putative protein target hypersensitive to volatile anesthetics. Experimental structure-derived closed- and open-channel nAChR systems in a fully hydrated lipid bilayer were examined using all-atom molecular dynamics simulations. At the majority of binding sites in alpha4beta2 nAChR, halothane replaced the slow-exchanging water molecules and caused a regional water population decrease. Only two binding sites had an increased quantity of water in the presence of halothane, where water arrangements resemble clathrate-like structures. The small number of such clathrate-like water clusters suggests that the formation of water clathrates is unlikely to be a primary cause for anesthesia. Despite the decrease in water population at most of the halothane binding sites, the number of sites that can be occupied transiently by water is actually increased in the presence of halothane. Many of these water sites were located between two subunits or in regions containing agonist binding sites or critical structural elements for transducing agonist binding to channel gating. Changes in water sites in the presence of halothane affected water-mediated protein-protein interactions and the protein dynamics, which can have direct impact on protein function. Collectively, water contributes to the action of anesthetics in proteins by mediating interactions between protein subunits and altering protein dynamics, instead of forming water clathrates around 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.

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.

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.

Anionic lipid and cholesterol interactions with alpha4beta2 nAChR: insights from MD simulations

Anionic lipids and cholesterols (CHOL) are critical to the function of nicotinic acetylcholine receptors (nAChR). We investigated their interactions with an open- and closed-channel alpha4beta2 nAChR by over 10 ns molecular dynamics simulations in a ternary lipid mixture of 1-palmitoyl-2-oleoyl phosphatidylcholine (POPC), 1-palmitoyl-2-oleoyl phosphatidic acid (POPA), and CHOL with a ratio of 3:1:1 (Haddadian et al., J. Phys. Chem. B 2008, 112, 13981). On average there were 65 and 74 interfacial lipids around the closed- and open-channel alpha4beta2 nAChR, respectively, in the equilibrated simulation systems. In the open-channel system, 42% of the interfacial POPA had acyl chains partially inserted into intra- or intersubunit cavities, as compared to only 7% in the closed-channel alpha4beta2. No CHOL was found in cavities within single subunits, though some CHOL infiltrated into the gaps between subunits. Because of its smaller headgroup, POPA could access some nonannular sites where POPC could not easily reach due to steric exclusion. Furthermore, POPA acted not only as an acceptor for hydrogen bonding (H bonding) as POPC did, but also as a donor through its hydroxyl group for H bonding with the backbone of the protein. The charged headgroup of POPA allowed the lipid to form stable salt bridges with conserved Arg and Lys residues at the interfaces of the transmembrane (TM) and extracellular (EC) or intracellular (IC) domains of the alpha4beta2. A higher number of salt bridges and hydrogen bonds (H bonds) between POPA and the alpha4beta2 nAChR were found in the open system than in the closed system, suggesting a potential role of POPA in the equilibrium between different channel states. Most interfacial POPA molecules showed lower order parameters than the bulk POPA due to the mixed effect of gauche defects, hydrophobic mismatch, and the lipid orientations near the magic angle. These unique properties enable the interfacial POPA to achieve what POPC cannot with regard to specific interactions with the protein, thereby making POPA essential for the function of nAChR.

Release of 11-cis-retinal from cellular retinaldehyde-binding protein by acidic lipids

PURPOSE

To determine molecular mechanisms for the release of 11-cis-retinal from the binding pocket of cellular retinaldehyde-binding protein (CRALBP).

METHODS

Binding of CRALBP to lipid surfaces was assessed with a lipid-immunoblot assay. Lipids were presented to CRALBP as small unilamellar vesicles (SUVs) consisting of phosphatidylcholine (PC) plus other lipids. Release of 9-cis-retinal or 11-cis-retinal from CRALBP was measured with spectral and high performance liquid chromatography (HPLC) assays based on the protection of the protein-bound retinal carbonyl group from reaction with NH(2)OH. The electrostatic surface potential of CRALBP was calculated from a model of its structure using the program CCP4mg.

RESULTS

Incubation of CRALBP.11-cis-retinal with lipids absorbed on nitrocellulose revealed binding to the acidic lipids, phosphatidic acid (PA)>phosphatidylinositol 3,4,5-trisphosphate [PI(3,4,5)P(3)]>phosphatidylserine (PS)> PI(4,5)P(2) and little or no binding to PC, phosphatidylethanolamine (PE), or PI(4)P. 11-cis-retinal was released during incubation of CRALBP with SUVs consisting of PC plus 50 mol% PA but not during incubation with those composed of 100 mol% PC. The efficacy of release of 9-cis-retinal or 11-cis-retinal from CRALBP by phospholipid-containing SUVs generally paralleled that of the binding of CRALBP to the lipids (PA>PS>PI>>PC). Examination of the electrostatic surface potential of the protein structure revealed a basic recess on one face of the protein, which may bind acidic lipids.

CONCLUSIONS

Our results identify the first physiologic substances that release 11-cis-retinal from CRALBP. PA and PS are relatively minor membrane lipids that can be generated in the cytoplasmic leaflet of the plasma membrane in response to various signal transduction pathways, where they could interact with cytosolic CRALBP. The mechanism for release of retinal from CRALBP by acidic lipids remains to be determined but could involve binding of the acidic lipid in the 11-cis-retinal binding site or to the positive basic recess on the protein surface. These results open a new facet in our understanding of how CRALBP functions in the regeneration of visual pigments.

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.

Examining the contributions of lipid shape and headgroup charge on bilayer behavior

To better understand bilayer property dependency on lipid electrostatics and headgroup size, we use atomistic molecular dynamics simulations to study negatively charged and neutral lipid membranes. We compare the negatively charged phosphatidic acid (PA), which at physiological pH and salt concentration has a negative spontaneous curvature, with the negatively charged phosphatidylglycerol (PG) and neutrally charged phosphatidylcholine (PC), both of which have zero spontaneous curvature. The PA lipids are simulated using two different sets of partial charges for the headgroup and the varied charge distribution between the two PA systems results in significantly different locations for the Na(+) ions relative to the water/membrane interface. For one PA system, the Na(+) ions are localized around the phosphate group. In the second PA system, the Na(+) ions are located near the ester carbonyl atoms, which coincides with the preferred location site for the PG Na(+) ions. We find that the Na(+) ion location has a larger effect on bilayer fluidity properties than lipid headgroup size, where the A(lipid) and acyl chain order parameter values are more similar between the PA and PG bilayers that have Na(+) ions located near the ester groups than between the two PA bilayers.