Molecular dynamics simulations of migration of ions and molecules through the acetylcholine receptor channel

A dynamic model of the channel of an acetylcholine receptor in a closed state has been proposed. The channel is formed by five a-helices of subunit M2 and stabilized by the cyclic hydrocarbon (CH2)105. The migration of charged and unchanged van der Waals particles with a diameter of 7.72 A equivalent to the diameter of a hydrated sodium ion has been studied. The migration occurred by the action of external force applied to the complex along the channel axis. In the closed state, the inhibition of ions is due to two components: electrostatic interaction and steric constraints. The van der Waals channel gate is formed by residues 13'-A-Val255, B-Val261, C-Val269, D-Val255, and E-Ile264, and the negatively changed residues occurring in the upper part of the channel have a great effect on ion selectivity.     

Acetylcholine receptor pore permeability studied by molecular dynamics simulation

A dynamic model of the closed-state pore of an acetylcholine receptor (five M2 α-helices stabilized with a (CH₂)₁₀₅ ring) is used to examine the migration of uncharged and charged probe particles equivalent to a hexahydrated sodium ion (van der Waals diameter 7.27 Å) propelled by varied external force along the channel axis. Ion movement through the pore is hindered by steric constraints and electrostatic interactions. The van der Waals gate is formed by helix residues 13′ (A-Val255, B-Val261, C-Val269, D-Val255, and E-Ile264), whereas the negatively charged residues in the upper part of the channel are important for ion selectivity.     

Regulation of channel function due to physical energetic coupling with a lipid bilayer

Regulation of membrane protein functions due to hydrophobic coupling with a lipid bilayer has been investigated. An energy formula describing interactions between lipid bilayer and integral ion channels with different structures, which is based on the screened Coulomb interaction approximation, has been developed. Here the interaction energy is represented as being due to charge-based interactions between channel and lipid bilayer. The hydrophobic bilayer thickness channel length mismatch is found to induce channel destabilization exponentially while negative lipid curvature linearly. Experimental parameters related to channel dynamics are consistent with theoretical predictions. To measure comparable energy parameters directly in the system and to elucidate the mechanism at an atomistic level we performed molecular dynamics (MD) simulations of the ion channel forming peptide–lipid complexes. MD simulations indicate that peptides and lipids experience electrostatic and van der Waals interactions for short period of time when found within each other’s proximity. The energies from these two interactions are found to be similar to the energies derived theoretically using the screened Coulomb and the van der Waals interactions between peptides (in ion channel) and lipids (in lipid bilayer) due to mainly their charge properties. The results of in silico MD studies taken together with experimental observable parameters and theoretical energetic predictions suggest that the peptides induce ion channels inside lipid membranes due to peptide–lipid physical interactions. This study provides a new insight helping better understand of the underlying mechanisms of membrane protein functions in cell membrane leading to important biological implications.     

Swelling of phospholipids by monovalent salt

Critical to biological processes such as membrane fusion and secretion, ion-lipid interactions at the membrane-water interface still raise many unanswered questions. Using reconstituted phosphatidylcholine membranes, we confirm here that multilamellar vesicles swell in salt solutions, a direct indication that salt modifies the interactions between neighboring membranes. By varying sample histories, and by comparing with data from ion carrier-containing bilayers, we eliminate the possibility that swelling is an equilibration artifact. Although both attractive and repulsive forces could be modified by salt, we show experimentally that swelling is driven primarily by weakening of the van der Waals attraction. To isolate the effect of salt on van der Waals interactions, we focus on high salt concentrations at which any possible electrostatic interactions are screened. By analysis of X-ray diffraction data, we show that salt does not alter membrane structure or bending rigidity, eliminating the possibility that repulsive fluctuation forces change with salt. By measuring changes in interbilayer separation with applied osmotic stress, we have determined, using the standard paradigm for bilayer interactions, that 1 M concentrations of KBr or KCl decrease the van der Waals strength by 50%. By weakening van der Waals attractions, salt increases energy barriers to membrane contact, possibly affecting cellular communication and biological signaling.     

Molecular mapping of general anesthetic sites in a voltage-gated ion channel

Several voltage-gated ion channels are modulated by clinically relevant doses of general anesthetics. However, the structural basis of this modulation is not well understood. Previous work suggested that n-alcohols and inhaled anesthetics stabilize the closed state of the Shaw2 voltage-gated (Kv) channel (K-Shaw2) by directly interacting with a discrete channel site. We hypothesize that the inhibition of K-Shaw2 channels by general anesthetics is governed by interactions between binding and effector sites involving components of the channel's activation gate. To investigate this hypothesis, we applied Ala/Val scanning mutagenesis to the S4-S5 linker and the post-PVP S6 segment, and conducted electrophysiological analysis to evaluate the energetic impact of the mutations on the inhibition of the K-Shaw2 channel by 1-butanol and halothane. These analyses identified residues that determine an apparent binding cooperativity and residue pairs that act in concert to modulate gating upon anesthetic binding. In some instances, due to their critical location, key residues also influence channel gating. Complementing these results, molecular dynamics simulations and in silico docking experiments helped us visualize possible anesthetic sites and interactions. We conclude that the inhibition of K-Shaw2 by general anesthetics results from allosteric interactions between distinct but contiguous binding and effector sites involving inter- and intrasubunit interfaces.     

Protein docking using continuum electrostatics and geometric fit

The computer program DOT quickly finds low-energy docked structures for two proteins by performing a systematic search over six degrees of freedom. A novel feature of DOT is its energy function, which is the sum of both a Poisson-Boltzmann electrostatic energy and a van der Waals energy, each represented as a grid-based correlation function. DOT evaluates the energy of interaction for many orientations of the moving molecule and maintains separate lists scored by either the electrostatic energy, the van der Waals energy or the composite sum of both. The free energy is obtained by summing the Boltzmann factor over all rotations at each grid point. Three important findings are presented. First, for a wide variety of protein-protein interactions, the composite-energy function is shown to produce larger clusters of correct answers than found by scoring with either van der Waals energy (geometric fit) or electrostatic energy alone. Second, free-energy clusters are demonstrated to be indicators of binding sites. Third, the contributions of electrostatic and attractive van der Waals energies to the total energy term appropriately reflect the nature of the various types of protein-protein interactions studied.     

Cross talk between activation and slow inactivation gates of Shaker potassium channels

This study addresses the energetic coupling between the activation and slow inactivation gates of Shaker potassium channels. To track the status of the activation gate in inactivated channels that are nonconducting, we used two functional assays: the accessibility of a cysteine residue engineered into the protein lining the pore cavity (V474C) and the liberation by depolarization of a Cs(+) ion trapped behind the closed activation gate. We determined that the rate of activation gate movement depends on the state of the inactivation gate. A closed inactivation gate favors faster opening and slower closing of the activation gate. We also show that hyperpolarization closes the activation gate long before a channel recovers from inactivation. Because activation and slow inactivation are ubiquitous gating processes in potassium channels, the cross talk between them is likely to be a fundamental factor in controlling ion flux across membranes.     

Computational studies of gramicidin permeation: An entryway sulfonate enhances cation occupancy at entry sites

The impact on the cation-transport free-energy profile of replacing the C-terminal ethanolamine in the gramicidin A channel with a taurine residue is studied using molecular dynamics simulations of gramicidin A (1JNO) embedded in a lipid bilayer (DMPC) with 1 mol/kg NaCl saline solution. The potential of mean force for ion transport is obtained by umbrella sampling. The presence of a negatively charged sulfonate group at the entrance of the gramicidin channel affects the depth and the location of the binding sites, producing a strong attraction for the cations in the bulk. The potential of mean force by the sulfonate acting directly through electrostatics and van der Waals interactions on the test ion is highly modulated by indirect effects (i.e., sulfonate effects on other components of the system that, in turn, affect the ion free-energy profile in the channel). Because the “entry” sites are located symmetrically at both entry and exit of the channel, the deeper free-energy wells should inhibit exit. Given that the channel has increased conductance experimentally, the simulation results suggest that the channel conductance is normally entry limited.     

Ion channel gating: insights via molecular simulations

Ion channels are gated, i.e. they can switch conformation between a closed and an open state. Molecular dynamics simulations may be used to study the conformational dynamics of ion channels and of simple channel models. Simulations on model nanopores reveal that a narrow (<4 A) hydrophobic region can form a functionally closed gate in the channel and can be opened by either a small (approximately 1 A) increase in pore radius or an increase in polarity. Modelling and simulation studies confirm the importance of hydrophobic gating in K channels, and support a model in which hinge-bending of the pore-lining M2 (or S6 in Kv channels) helices underlies channel gating. Simulations of a simple outer membrane protein, OmpA, indicate that a gate may also be formed by interactions of charged side chains within a pore, as is also the case in ClC channels.     

Molecular mechanical perspective on halogen bonding

The nature and strength of halogen bonding in halo molecule-Lewis base complexes were studied in terms of molecular mechanics using our recently developed positive extra-point (PEP) approach, in which the σ-hole on the halogen atom is represented by an extra point of positive charge. The contributions of the σ-hole (i.e., positively charged extra point) and the halogen atom to the strength of this noncovalent interaction were clarified using the atomic parameter contribution to the molecular interaction (APCtMI) approach. The molecular mechanical results revealed that the halogen bond is electrostatic and van der Waals in nature, and its strength depends on three types of interaction: (1) the attractive electrostatic interaction between the σ-hole and the Lewis base, (2) the repulsive electrostatic interaction between the negative halogen atom and the Lewis base, and (3) the repulsive/attractive van der Waals interactions between the halogen atom and the Lewis base. The strength of the halogen bond increases with increasing σ-hole size (i.e., magnitude of the extra-point charge) and increasing halogen atom size. The van der Waals interaction's contribution to the halogen bond strength is most favorable in chloro complexes, whereas the electrostatic interaction is dominant in iodo complexes. The idea that the chloromethane molecule can form a halogen bond with a Lewis base was revisited in terms of quantum mechanics and molecular mechanics. Although chloromethane does produce a positive region along the C-Cl axis, basis set superposition error corrected second-order M?ller-Plesset calculations showed that chloromethane-Lewis base complexes are unstable, producing halogen-Lewis base contacts longer than the sum of the van der Waals radii of the halogen and O/N atoms. Molecular mechanics using the APCtMI approach showed that electrostatic interactions between chloromethane and a Lewis base are unfavorable owing to the high negative charge on the chlorine atom, which overcomes the corresponding favorable van der Waals interactions.     

Energy of ligand-RNA complex formation

This work presents an energetic analysis of complex formation of eleven ligands with different structure and charge with RNA aptamers. The most entire set of the components of the total Gibbs energy of complex formation has been first calculated using different physical factors: van der Waals, electrostatic and hydrophobic interactions, hydrogen bonds and specific factors being predominantly of entropic character. The calculated Gibbs energy is found to be in a good agreement with experimental data. Different energy components which stabilize and destabilize complexes are lined up according to the degree of importance. The results obtained provide an understanding of the role of different physical interactions in a ligand-RNA complex formation.     

Molecular dynamics simulations of the Cx26 hemichannel: insights into voltage-dependent loop-gating

Loop-gating is one of two voltage-dependent mechanisms that regulate the open probability of connexin channels. The loop-gate permeability barrier is formed by a segment of the first extracellular loop (E1) (the parahelix) and appears to be accompanied by straightening of the bend angle between E1 and the first transmembrane domain (TM1). Here, all-atom molecular dynamics simulations are used to identify and characterize interacting van der Waals and electrostatic networks that stabilize the parahelices and TM1/E1 bend angles of the open Cx26 hemichannel. Dynamic fluctuations in an electrostatic network in each subunit are directly linked to the stability of parahelix structure and TM1/E1 bend angle in adjacent subunits. The electrostatic network includes charged residues that are pore-lining and thus positioned to be voltage sensors. We propose that the transition to the closed state is initiated by voltage-driven disruption of the networks that stabilize the open-state parahelix configuration, allowing the parahelix to protrude into the channel pore to form the loop-gate barrier. Straightening of the TM1/E1 bend appears to be a consequence of the reorganization of the interacting networks that accompany the conformational change of the parahelix. The electrostatic network extends across subunit boundaries, suggesting a concerted gating mechanism.     

Tricyclic antidepressants and mecamylamine bind to different sites in the human α4β2 nicotinic receptor ion channel

The interaction of tricyclic antidepressants with the human (h) α4β2 nicotinic acetylcholine receptor in different conformational states was compared with that for the noncompetitive antagonist mecamylamine by using functional and structural approaches. The results established that: (a) [³H]imipramine binds to hα4β2 receptors with relatively high affinity (K_d = 0.83 ± 0.08 μM), but imipramine does not differentiate between the desensitized and resting states, (b) although tricyclic antidepressants inhibit (±)-epibatidine-induced Ca²⁺ influx in HEK293-hα4β2 cells with potencies that are in the same concentration range as that for (±)-mecamylamine, tricyclic antidepressants inhibit [³H]imipramine binding to hα4β2 receptors with affinities >100-fold higher than that for (±)-mecamylamine. This can be explained by our docking results where imipramine interacts with the leucine (position 9′) and valine (position 13′) rings by van der Waals contacts, whereas mecamylamine interacts electrostatically with the outer ring (position 20′), (c) van der Waals interactions are in agreement with the thermodynamic results, indicating that imipramine interacts with the desensitized and resting receptors by a combination of enthalpic and entropic components. However, the entropic component is more important in the desensitized state, suggesting local conformational changes. In conclusion, our data indicate that tricyclic antidepressants and mecamylamine efficiently inhibit the ion channel by interacting at different luminal sites. The high proportion of protonated mecamylamine calculated at physiological pH suggests that this drug can be attracted to the channel mouth before binding deeper within the receptor ion channel finally blocking ion flux.     

A free energy analysis of nucleic acid base stacking in aqueous solution.

This paper reports a theoretical study of the free energy contributions to nucleic acid base stacking in aqueous solution. Electrostatic interactions are treated by using the finite difference Poisson-Boltzmann method and nonpolar effects are treated with explicit calculation of van der Waals interactions and/or free energy-surface area relationships. Although for some pairs of bases there is a favorable Coulombic interaction in the stacked conformation, generally the net effect of electrostatic interactions is to oppose stacking. This result is caused by the loss of favorable base-solvent electrostatic interactions, that accompany the partial removal of polar atoms from water in the stacked conformation. Nonpolar interactions, involving the hydrophobic effect and enhancement of van der Waals interactions caused by close-packing, drive stacking. The calculations qualitatively reproduce the experimental dependence of stacking free energy on purine-pyrimidine composition.     

Ion channels formed by amphipathic helical peptides

Channel forming peptides (CFPs) are amphipathic peptides, of length ca. 20 residues, which adopt an -helical conformation in the presence of lipid bilayers and form ion channels with electrophysiological properties comparable to those of ion channel proteins. We have modelled CFP channels as bundles of parallel trans-bilayer helices surrounding a central ion-permeable pore. Ion-channel interactions have been explored via accessible surface area calculations, and via evaluation of changes in van der Waals and electrostatic energies as a K⁺ ion is translated along the length of the pore. Two CFPs have been modelled: (a) zervamicin-A1-16, a synthetic apolar peptaibol related to alamethicin, and (b) -toxin from Staphylococcus aureus. Both of these CFPs have previously been shown to form ion channels in planar lipid bilayers, and have been shown to have predominantly helical conformations. Zervamicin-A1-16 channels were modelled as bundles of 4 to 8 parallel helices. Two related helix bundle geometries were explored. K⁺channel interactions have been shown to involve exposed backbone carbonyl oxygen atoms. -Toxin channels were modelled as bundles of 6 parallel helices. Residues Q3, D11 and D18 generate favourable K⁺-channel interactions. Rotation of W15 about its C-C bond has been shown to be capable of occluding the central pore, and is discussed as a possible model for sidechain conformational changes in relation to ion channel gating.     

Thermodynamic coupling between activation and inactivation gating in potassium channels revealed by free energy molecular dynamics simulations

The amount of ionic current flowing through K(+) channels is determined by the interplay between two separate time-dependent processes: activation and inactivation gating. Activation is concerned with the stimulus-dependent opening of the main intracellular gate, whereas inactivation is a spontaneous conformational transition of the selectivity filter toward a nonconductive state occurring on a variety of timescales. A recent analysis of multiple x-ray structures of open and partially open KcsA channels revealed the mechanism by which movements of the inner activation gate, formed by the inner helices from the four subunits of the pore domain, bias the conformational changes at the selectivity filter toward a nonconductive inactivated state. This analysis highlighted the important role of Phe103, a residue located along the inner helix, near the hinge position associated with the opening of the intracellular gate. In the present study, we use free energy perturbation molecular dynamics simulations (FEP/MD) to quantitatively elucidate the thermodynamic basis for the coupling between the intracellular gate and the selectivity filter. The results of the FEP/MD calculations are in good agreement with experiments, and further analysis of the repulsive, van der Waals dispersive, and electrostatic free energy contributions reveals that the energetic basis underlying the absence of inactivation in the F103A mutation in KcsA is the absence of the unfavorable steric interaction occurring with the large Ile100 side chain in a neighboring subunit when the intracellular gate is open and the selectivity filter is in a conductive conformation. Macroscopic current analysis shows that the I100A mutant indeed relieves inactivation in KcsA, but to a lesser extent than the F103A mutant.     

NMR investigations of protein-carbohydrate interactions: refined three- dimensional structure of the complex between hevein and methyl beta- chitobioside

The specific interaction of hevein with GlcNAc-containing oligosaccharides has been analyzed by1H-NMR spectroscopy. The association constants for the binding of hevein to a variety of ligands have been estimated from1H-NMR titration experiments. The association constants increase in the order GlcNAc-alpha(1-->6)-Man < GlcNAc < benzyl-beta-GlcNAc < p-nitrophenyl-beta-GlcNAc < chitobiose < p- nitrophenyl-beta-chitobioside < methyl-beta-chitobioside < chitotriose. Entropy and enthalpy of binding for different complexes have been obtained from van't Hoff analysis. The driving force for the binding process is provided by a negative DeltaH0which is partially compensated by negative DeltaS0. These negative signs indicate that hydrogen bonding and van der Waals forces are the major interactions stabilizing the complex. NOESY NMR experiments in water solution provided 475 accurate protein proton-proton distance constraints after employing the MARDIGRAS program. In addition, 15 unambiguous protein/carbohydrate NOEs were detected. All the experimental constraints were used in a refinement protocol including restrained molecular dynamics in order to determine the highly refined solution conformation of this protein- carbohydrate complex. With regard to the NMR structure of the free protein, no important changes in the protein nOe's were observed, indicating that carbohydrate-induced conformational changes are small. The average backbone rmsd of the 20 refined structures was 0.055 nm, while the heavy atom rmsd was 0.116 nm. It can be deduced that both hydrogen bonds and van der Waals contacts confer stability to the complex. A comparison of the three-dimensional structure of hevein in solution to those reported for wheat germ agglutinin (WGA) and hevein itself in the solid state has also been performed. The polypeptide conformation has also been compared to the NMR-derived structure of a smaller antifungical peptide, Ac-AMP2.      

Probing structure-function relationships and gating mechanisms in the CorA Mg2+ transport system

Recent crystal structures of the CorA Mg(2+) transport protein from Thermotoga maritima (TmCorA) revealed an unusually long ion pore putatively gated by hydrophobic residues near the intracellular end and by universally conserved asparagine residues at the periplasmic entrance. A conformational change observed in an isolated funnel domain structure also led to a proposal for the structural basis of gating. Because understanding the molecular mechanisms underlying ion channel and transporter gating remains an important challenge, we have undertaken a structure-guided engineering approach to probe structure-function relationships in TmCorA. The intracellular funnel domain is shown to constitute an allosteric regulatory module that can be engineered to promote an activated or closed state. A periplasmic gate centered about a proline-induced kink of the pore-lining helix is described where "helix-straightening" mutations produce a dramatic gain-of-function. Mutation to the narrowest constriction along the pore demonstrates that a hydrophobic gate is operational within this Mg(2+)-selective transport protein and likely forms an energetic barrier to ion flux. We also provide evidence that highly conserved acidic residues found in the short periplasmic loop are not essential for TmCorA function or Mg(2+) selectivity but may be required for proper protein folding and stability. This work extends our gating model for the CorA-Alr1-Mrs2 superfamily and reveals features that are characteristic of an ion channel. Aspects of these results that have broader implications for a range of channel and transporter families are highlighted.     

Investigation of Atomic Level Patterns in Protein—Small Ligand Interactions

Background

Shape complementarity and non-covalent interactions are believed to drive protein-ligand interaction. To date protein-protein, protein-DNA, and protein-RNA interactions were systematically investigated, which is in contrast to interactions with small ligands. We investigate the role of covalent and non-covalent bonds in protein-small ligand interactions using a comprehensive dataset of 2,320 complexes.

Methodology and Principal Findings

We show that protein-ligand interactions are governed by different forces for different ligand types, i.e., protein-organic compound interactions are governed by hydrogen bonds, van der Waals contacts, and covalent bonds; protein-metal ion interactions are dominated by electrostatic force and coordination bonds; protein-anion interactions are established with electrostatic force, hydrogen bonds, and van der Waals contacts; and protein-inorganic cluster interactions are driven by coordination bonds. We extracted several frequently occurring atomic-level patterns concerning these interactions. For instance, 73% of investigated covalent bonds were summarized with just three patterns in which bonds are formed between thiol of Cys and carbon or sulfur atoms of ligands, and nitrogen of Lys and carbon of ligands. Similar patterns were found for the coordination bonds. Hydrogen bonds occur in 67% of protein-organic compound complexes and 66% of them are formed between NH- group of protein residues and oxygen atom of ligands. We quantify relative abundance of specific interaction types and discuss their characteristic features. The extracted protein-organic compound patterns are shown to complement and improve a geometric approach for prediction of binding sites.

Conclusions and Significance

We show that for a given type (group) of ligands and type of the interaction force, majority of protein-ligand interactions are repetitive and could be summarized with several simple atomic-level patterns. We summarize and analyze 10 frequently occurring interaction patterns that cover 56% of all considered complexes and we show a practical application for the patterns that concerns interactions with organic compounds.