Water-mediated conformational transitions in nicotinic receptor M2 helix bundles: a molecular dynamics study

Melittin: a membrane-active peptide with diverse functions

Melittin is the principal toxic component in the venom of the European honey bee Apis mellifera and is a cationic, hemolytic peptide. It is a small linear peptide composed of 26 amino acid residues in which the amino-terminal region is predominantly hydrophobic whereas the carboxy-terminal region is hydrophilic due to the presence of a stretch of positively charged amino acids. This amphiphilic property of melittin has resulted in melittin being used as a suitable model peptide for monitoring lipid-protein interactions in membranes. In this review, the solution and membrane properties of melittin are highlighted, with an emphasis on melittin-membrane interaction using biophysical approaches. The recent applications of melittin in various cellular processes are discussed.

Dynamics of a membrane-bound tryptophan analog in environments of varying hydration: a fluorescence approach

Tryptophan octyl ester (TOE) represents an important model for membrane-bound tryptophan residues. In this article, we have employed a combination of wavelength-selective fluorescence and time-resolved fluorescence spectroscopies to monitor the effect of varying degrees of hydration on the dynamics of TOE in reverse micellar environments formed by sodium bis(2-ethylhexyl) sulfosuccinate (AOT) in isooctane. Our results show that TOE exhibits red edge excitation shift (REES) and other wavelength-selective fluorescence effects when bound to reverse micelles of AOT. Fluorescence parameters such as intensity, emission maximum, anisotropy, and lifetime of TOE in reverse micelles of AOT depend on [water]/[surfactant] molar ratio (w (o)). These results are relevant and potentially useful for analyzing dynamics of proteins or peptides bound to membranes or membrane-mimetic media under conditions of changing hydration.

Effect of graded hydration on the dynamics of an ion channel peptide: a fluorescence approach

Water plays an important role in determining the folding, structure, dynamics, and, in turn, the function of proteins. We have utilized a combination of fluorescence approaches such as the wavelength-selective fluorescence approach to monitor the effect of varying degrees of hydration on the organization and dynamics of the functionally important tryptophan residues of gramicidin in reverse micelles formed by sodium bis(2-ethylhexyl) sulfosuccinate. Our results show that tryptophans in gramicidin, present in the single-stranded beta6.3 conformation, experience slow solvent relaxation giving rise to red-edge excitation shift (REES). In addition, changes in fluorescence polarization with increasing excitation or emission wavelength reinforce that the gramicidin tryptophans are localized in motionally restricted regions of the reverse micelle. Interestingly, the extent of REES is found to be independent of the [water]/[surfactant] molar ratio (w(o)). We attribute this to heterogeneity in gramicidin tryptophan localization. Fluorescence intensity and mean fluorescence lifetime of the gramicidin tryptophans show significant reductions with increasing w(o) indicating sensitivity to increased polarity. Since the dynamics of hydration is related to folding, structure, and eventually function of proteins, we conclude that REES could prove to be a potentially sensitive tool to explore the dynamics of proteins under conditions of changing hydration.

Organization and dynamics of tryptophan residues in erythroid spectrin: novel structural features of denatured spectrin revealed by the wavelength-selective fluorescence approach

We have investigated the organization and dynamics of the functionally important tryptophan residues of erythroid spectrin in native and denatured conditions utilizing the wavelength-selective fluorescence approach. We observed a red edge excitation shift (REES) of 4 nm for the tryptophans in the case of spectrin in its native state. This indicates that tryptophans in spectrin are localized in a microenvironment of restricted mobility, and that the regions surrounding the spectrin tryptophans offer considerable restriction to the reorientational motion of the water dipoles around the excited state tryptophans. Interestingly, spectrin exhibits a REES of 3 nm even when denatured in 8 M urea. This represents the first report of a denatured protein displaying REES. Observation of REES in the denatured state implies that some of the structural and dynamic features of this microenvironment around the spectrin tryptophans are retained even when the protein is denatured. Fluorescence quenching data of denatured spectrin support this conclusion. In addition, we have deduced the organization and dynamics of the hydrophobic binding site of the polarity-sensitive fluorescent probe PRODAN that binds erythroid spectrin with high affinity. When bound to spectrin, PRODAN exhibits a REES of 9 nm. Because PRODAN binds to a hydrophobic site in spectrin, such a result would directly imply that this region of spectrin offers considerable restriction to the reorientational motion of the solvent dipoles around the excited state fluorophore. The results of our study could provide vital insight into the role of tryptophans in the stability and folding of spectrin.

Exploring membrane organization and dynamics by the wavelength-selective fluorescence approach

Wavelength-selective fluorescence comprises a set of approaches based on the red edge effect in fluorescence spectroscopy which can be used to directly monitor the environment and dynamics around a fluorophore in a complex biological system. A shift in the wavelength of maximum fluorescence emission toward higher wavelengths, caused by a shift in the excitation wavelength toward the red edge of absorption band, is termed red edge excitation shift (REES). This effect is mostly observed with polar fluorophores in motionally restricted media such as very viscous solutions or condensed phases where the dipolar relaxation time for the solvent shell around a fluorophore is comparable to or longer than its fluorescence lifetime. REES arises from slow rates of solvent relaxation (reorientation) around an excited state fluorophore which is a function of the motional restriction imposed on the solvent molecules in the immediate vicinity of the fluorophore. Utilizing this approach, it becomes possible to probe the mobility parameters of the environment itself (which is represented by the relaxing solvent molecules) using the fluorophore merely as a reporter group. Further, since the ubiquitous solvent for biological systems is water, the information obtained in such cases will come from the otherwise 'optically silent' water molecules. This makes REES and related techniques extremely useful since hydration plays a crucial modulatory role in a large number of important cellular events, including lipid-protein interactions and ion transport. The interfacial region in membranes, characterized by unique motional and dielectric characteristics, represents an appropriate environment for displaying wavelength-selective fluorescence effects. The application of REES and related techniques (wavelength-selective fluorescence approach) as a powerful tool to monitor the organization and dynamics of probes and peptides bound to membranes, micelles, and reverse micelles is discussed.

Passive water and urea permeability of a human Na(+)-glutamate cotransporter expressed in Xenopus oocytes

The human Na(+)-glutamate transporter (EAAT1) was expressed in Xenopus laevis oocytes. The passive water permeability, L(p), was derived from volume changes of the oocyte induced by changes in the external osmolarity. Oocytes were subjected to two-electrode voltage clamp. In the presence of Na(+), the EAAT1-specific (defined in Discussion) L(p) increased linearly with positive clamp potentials, the L(p) being around 23 % larger at +50 mV than at -50 mV. L-Glutamate increased the EAAT1-specific L(p) by up to 40 %. The K(0.5) for the glutamate-dependent increase was 20 +/- 6 microM, which is similar to the K(0.5) value for glutamate activation of transport. The specific inhibitor DL-threo-beta-benzyloxyaspartate (TBOA) reduced the EAAT1-specific L(p) to 72 %. EAAT1 supported passive fluxes of [(14)C]urea and [(14)C]glycerol. The [(14)C]urea flux was increased in the presence of glutamate. The data suggest that the permeability depends on the conformational equilibrium of the EAAT1. At positive potentials and in the presence of Na(+) and glutamate, the pore is enlarged and water and urea penetrate more readily. The L(p) was larger when measured with urea or glycerol as osmolytes as compared with mannitol. Apparently, the properties of the pore are not uniform along its length. The outer section may accommodate urea and glycerol in an osmotically active form, giving rise to larger water fluxes. The physiological role of EAAT1 for water homeostasis in the central nervous system is discussed.

Molecular modeling of ligand-gated ion channels: progress and challenges

There has been rapid progress in molecular modeling of LGICs in recent years. The convergence of improved software for molecular mechanics/dynamics, techniques of chimeric substitution and site-directed mutations, and the first X-ray structures of transmembrane ion channels will make it possible to build reasonable models of neuronal ion channels well in advance of publication of their crystal structures. These models will not only serve as guides for future site-directed mutagenesis, but they will also be a starting point for understanding the dynamics of ion channel gating.

Anesthetics and ion channels: molecular models and sites of action

The mechanisms of general anesthesia in the central nervous system are finally yielding to molecular examination. As a result of research during the past several decades, a group of ligand-gated ion channels have emerged as plausible targets for general anesthetics. Molecular biology techniques have greatly accelerated attempts to classify ligand-gated ion channel sensitivity to general anesthetics, and have identified the sites of receptor subunits critical for anesthetic modulation using chimeric and mutated receptors. The experimental data have facilitated the construction of tenable molecular models for anesthetic binding sites, which in turn allows structural predictions to be tested. In vivo significance of a putative anesthetic target can now be examined by targeted gene manipulations in mice. In this review, we summarize from a molecular perspective recent advances in our understanding of mechanisms of action of general anesthetics on ligand-gated ion channels.

Specific binding sites for alcohols and anesthetics on ligand-gated ion channels

Ligand-gated ion channels are a target for inhaled anesthetics and alcohols in the central nervous system. The inhibitory strychnine-sensitive glycine and gamma-aminobutyric acid type A receptors are positively modulated by anesthetics and alcohols, and site-directed mutagenesis techniques have identified amino acid residues important for the action of volatile anesthetics and alcohols in these receptors. A key question is whether these amino acids are part of an alcohol/anesthetic-binding site. In the present study, we used an alkanethiol anesthetic to covalently label its binding site by mutating selected amino acids to cysteine. We demonstrated that the anesthetic propanethiol, or alternatively, propyl methanethiosulfonate, covalently binds to cysteine residues introduced into a specific second transmembrane site in glycine receptor and gamma-aminobutyric acid type A receptor subunits and irreversibly enhances receptor function. Moreover, upon permanent occupation of the site by propyl disulfide, the usual ability of octanol, enflurane, and isoflurane to potentiate the function of the ion channels was lost. This approach provides strong evidence that the actions of anesthetics in these receptors are due to binding at a single site.

Intersegment hydrogen bonds as possible structural determinants of the N/Q/R site in glutamate receptors

Specific electrophysiological and pharmacological properties of ionic channels in NMDA, AMPA, and kainate subtypes of ionotropic glutamate receptors (GluRs) are determined by the Asn (N), Gln (Q), and Arg (R) residues located at homologous positions of the pore-lining M2 segments (the N/Q/R site). Presumably, the N/Q/R site is located at the apex of the reentrant membrane loop and forms the narrowest constriction of the pore. Although the shorter Asn residues are expected to protrude in the pore to a lesser extent than the longer Gln residues, the effective dimension of the NMDA channel (corresponding to the size of the largest permeant organic cation) is, surprisingly, smaller than that of the AMPA channel. To explain this paradox, we propose that the N/Q/R residues form macrocyclic structures (rings) stabilized by H-bonds between a NH(2) group in the side chain of a given M2 segment and a C==O group of the main chain in the adjacent M2 segment. Using Monte Carlo minimization, we have explored conformational properties of the rings. In the Asn, but not in the Gln ring, the side-chain oxygens protruding into the pore may facilitate ion permeation and accept H-bonds from the blocking drugs. In this way, the model explains different electrophysiological and pharmacological properties of NMDA and non-NMDA GluR channels. The ring of H-bonded polar residues at the pore narrowing resembles the ring of four Thr(75) residues observed in the crystallographic structure of the KcsA K(+) channel.

Channels formed by the transmembrane helix of phospholamban: a simulation study

Phospholamban is a small membrane protein which can form cation selective ion channels in lipid bilayers. Each subunit contains a single, largely hydrophobic transmembrane helix. The helices are thought to assemble as a pentameric and approximately parallel bundle surrounding a central pore. A model of this assembly (PDB code IPSL) has been used as the starting point for molecular dynamics (MD) simulations of a system consisting of the pentameric helix bundle, plus 217 water molecules located within and at either mouth of the pore. Interhelix distance restraints were employed to maintain the integrity of the helix bundle during a 500 ps MD simulation. Water molecules within the pore exhibited reduced diffusional and rotational mobility. Interactions between the alpha-helix dipoles and the water dipoles, the latter aligned anti-parallel to the former, contribute to the stability of the system. Analysis of the potential energy of interaction of a K+ ion as it was moved through the pore suggested that unfavourable interactions of the cation with the aligned helix dipoles at the N-terminal mouth were overcome by favourable ion-water interactions. Comparable analysis for a Cl ion revealed that the ion-(pore + water) interactions were unfavourable along the whole of the pore, increasingly so from the N- to the C-terminal mouth. Overall, the interaction energy profiles were consistent with a pore selective for cations over anions. Pore radius profiles were used to predict a channel conductance of 50 to 70 ps in 0.2 M KCl, which compares well with an experimental value of 100 ps.

Molecular dynamics study of water and Na+ ions in models of the pore region of the nicotinic acetylcholine receptor

The nicotinic acetylcholine receptor (nAChR) is an integral membrane protein that forms ligand-gated and cation-selective channels. The central pore is lined by a bundle of five approximately parallel M2 helices, one from each subunit. Candidate model structures of the solvated pore region of a homopentameric (alpha7)5 nAChR channel in the open state, and in two possible forms of the closed state, have been studied using molecular dynamics simulations with restraining potentials. It is found that the mobility of the water is substantially lower within the pore than in bulk, and the water molecules become aligned with the M2 helix dipoles. Hydrogen-bonding patterns in the pore, especially around pore-lining charged and hydrophilic residues, and around exposed regions of the helix backbone, have been determined. Initial studies of systems containing both water and sodium ions together within the pore region have also been conducted. A sodium ion has been introduced into the solvated models at various points along the pore axis and its energy profile evaluated. It is found that the ion causes only a local perturbation of the water structure. The results of these calculations have been used to examine the effectiveness of the central ring of leucines as a component of a gate in the closed-channel model.

Alamethicin channels - modelling via restrained molecular dynamics simulations

Alamethicin channels have been modelled as approximately parallel bundles of transbilayer helices containing between N = 4 and 8 helices per bundle. Initial models were generated by in vacuo restrained molecular dynamics (MD) simulations, and were refined by 60 ps MD simulations with water molecules present within and at the mouths of the central pore. The helix bundles were stabilized by networks of H-bonds between intra-pore water molecules and Gln-7 side-chains. Channel conductances were predicted on the basis of pore radius profiles, and suggested that the N = 4 bundle formed an occluded pore, whereas pores with N > or = 5 helices per bundle were open. Continuum electrostatics calculations suggested that the N = 6 pore is cation-selective, whereas pores with N > or = 7 helices per bundle were predicted to be somewhat less ion-selective.

Molecular dynamics simulations of water within models of ion channels

The transbilayer pores formed by ion channel proteins contain extended columns of water molecules. The dynamic properties of such waters have been suggested to differ from those of water in its bulk state. Molecular dynamics simulations of ion channel models solvated within and at the mouths of their pores are used to investigate the dynamics and structure of intra-pore water. Three classes of channel model are investigated: a) parallel bundles of hydrophobic (Ala20) alpha-helices; b) eight-stranded hydrophobic (Ala10) antiparallel beta-barrels; and c) parallel bundles of amphipathic alpha-helices (namely, delta-toxin, alamethicin, and nicotinic acetylcholine receptor M2 helix). The self-diffusion coefficients of water molecules within the pores are reduced significantly relative to bulk water in all of the models. Water rotational reorientation rates are also reduced within the pores, particularly in those pores formed by alpha-helix bundles. In the narrowest pore (that of the Ala20 pentameric helix bundle) self-diffusion coefficients and reorientation rates of intra-pore waters are reduced by approximately an order of magnitude relative to bulk solvent. In Ala20 helix bundles the water dipoles orient antiparallel to the helix dipoles. Such dipole/dipole interaction between water and pore may explain how water-filled ion channels may be formed by hydrophobic helices. In the bundles of amphipathic helices the orientation of water dipoles is modulated by the presence of charged side chains. No preferential orientation of water dipoles relative to the pore axis is observed in the hydrophobic beta-barrel models.