Determination of the structure of a membrane-incorporated ion channel. Solid-state nuclear magnetic resonance studies of gramicidin A

NMR techniques for investigating antimicrobial peptides in model membranes and bacterial cells

AMPs are short, mainly cationic membrane-active peptides found in all living organism. They perform diverse roles including signaling and acting as a line of defense against bacterial infections. AMPs have been extensively investigated as templates to facilitate the development of novel antimicrobial therapeutics. Understanding the interplay between these membrane-active peptides and the lipid membranes is considered to be a significant step in elucidating the specific mechanism of action of AMPs against prokaryotic and eukaryotic cells to aid the development of new therapeutics. In this review, we have provided a brief overview of various NMR techniques commonly used for studying AMP structure and AMP-membrane interactions in model membranes and whole cells.

Determining the helical tilt angle and dynamic properties of the transmembrane domains of pinholin S2168 using mechanical alignment EPR spectroscopy

The lytic cycle of bacteriophage φ21 for the infected E. coli is initiated by pinholin S²¹, which determines the timing of host cell lysis through the function of pinholin (S²¹68) and antipinholin (S²¹71). The activity of pinholin or antipinholin directly depends on the function of two transmembrane domains (TMDs) within the membrane. For active pinholin, TMD1 externalizes and lies on the surface while TMD2 remains incorporated inside the membrane forming the lining of the small pinhole. In this study, spin labeled pinholin TMDs were incorporated separately into mechanically aligned POPC (1-palmitoyl-2-oleoyl-glycero-3-phosphocholine) lipid bilayers and investigated with electron paramagnetic resonance (EPR) spectroscopy to determine the topology of both TMD1 and TMD2 with respect to the lipid bilayer; the TOAC (2,2,6,6-tetramethyl-N-oxyl-4-amino-4-carboxylic acid) spin label was used here because it attaches to the backbone of a peptide and is very rigid. TMD2 was found to be nearly colinear with the bilayer normal (n) with a helical tilt angle of 16 ± 4° while TMD1 lies on or near the surface with a helical tilt angle of 84 ± 4°. The order parameters (~0.6 for both TMDs) obtained from our alignment study were reasonable, which indicates the samples incorporated inside the membrane were well aligned with respect to the magnetic field (B₀). The data obtained from this study supports previous findings on pinholin: TMD1 partially externalizes from the lipid bilayer and interacts with the membrane surface, whereas TMD2 remains buried in the lipid bilayer in the active conformation of pinholin S²¹68. In this study, the helical tilt angle of TMD1 was measured for the first time. For TMD2 our experimental data corroborates the findings of the previously reported helical tilt angle by the Ulrich group.

Chemically modified and conjugated antimicrobial peptides against superbugs

Antimicrobial resistance (AMR) is one of the greatest threats to human health that, by 2050, will lead to more deaths from bacterial infections than cancer. New antimicrobial agents, both broad-spectrum and selective, that do not induce AMR are urgently required. Antimicrobial peptides (AMPs) are a novel class of alternatives that possess potent activity against a wide range of Gram-negative and positive bacteria with little or no capacity to induce AMR. This has stimulated substantial chemical development of novel peptide-based antibiotics possessing improved therapeutic index. This review summarises recent synthetic efforts and their impact on analogue design as well as their various applications in AMP development. It includes modifications that have been reported to enhance antimicrobial activity including lipidation, glycosylation and multimerization through to the broad application of novel bio-orthogonal chemistry, as well as perspectives on the direction of future research. The subject area is primarily the development of next-generation antimicrobial agents through selective, rational chemical modification of AMPs. The review further serves as a guide toward the most promising directions in this field to stimulate broad scientific attention, and will lead to new, effective and selective solutions for the several biomedical challenges to which antimicrobial peptidomimetics are being applied.

A theoretical assessment of structure determination of multi-span membrane proteins by oriented sample solid-state NMR spectroscopy

Oriented sample solid state NMR (OS-ssNMR) spectroscopy allows direct determination of the structure and topology of membrane proteins reconstituted into aligned lipid bilayers. While OS-ssNMR theoretically has no upper size limit, its application to multi-span membrane proteins has not been established since most studies have been restricted to single or dual span proteins and peptides. Here, we present a critical assessment of the application of this method to multi-span membrane proteins. We used molecular dynamics simulations to back-calculate [¹⁵N-¹H] separated local field (SLF) spectra from a G protein-coupled receptor (GPCR) and show that fully resolved spectra can be obtained theoretically for a multi-span membrane protein with currently achievable resonance linewidths.

Peptides derived from the C-terminal domain of HIV-1 Viral Protein R in lipid bilayers: Structure, membrane positioning and gene delivery

Viral protein R (Vpr) is a small accessory protein of 96 amino acids that is present in Human and simian immunodeficiency viruses. Among the very different properties that Vpr possesses we can find cell penetrating capabilities. Based on this and on its capacity to interact with nucleic acids we previously investigated the DNA transfection properties of Vpr and subfragments thereof. We found that fragments of the C-terminal helical domain of Vpr are able to deliver efficiently plasmid DNA into different cell lines. As the amphipathic helix may play a role in the interactions with membranes, we investigated whether insertion of a proline residue in the α-helix modifies the transfection properties of Vpr. Unexpectedly, we found that the resulting Vpr55-82 Pro70 peptide was even more efficient than wild type Vpr55-82 in the gene delivery assays. Using circular dichroism, light scattering and solid-state NMR techniques, we characterized the secondary structure and interactions of Vpr and several mutants with model membranes. A model is proposed where the proline shifts the dissociation equilibrium of the peptide-cargo complex and thereby its endosomal release.

Mechanistic Landscape of Membrane-Permeabilizing Peptides

Membrane permeabilizing peptides (MPPs) are as ubiquitous as the lipid bilayer membranes they act upon. Produced by all forms of life, most membrane permeabilizing peptides are used offensively or defensively against the membranes of other organisms. Just as nature has found many uses for them, translational scientists have worked for decades to design or optimize membrane permeabilizing peptides for applications in the laboratory and in the clinic ranging from antibacterial and antiviral therapy and prophylaxis to anticancer therapeutics and drug delivery. Here, we review the field of membrane permeabilizing peptides. We discuss the diversity of their sources and structures, the systems and methods used to measure their activities, and the behaviors that are observed. We discuss the fact that "mechanism" is not a discrete or a static entity for an MPP but rather the result of a heterogeneous and dynamic ensemble of structural states that vary in response to many different experimental conditions. This has led to an almost complete lack of discrete three-dimensional active structures among the thousands of known MPPs and a lack of useful or predictive sequence-structure-function relationship rules. Ultimately, we discuss how it may be more useful to think of membrane permeabilizing peptides mechanisms as broad regions of a mechanistic landscape rather than discrete molecular processes.

Dynamic membrane interactions of antibacterial and antifungal biomolecules, and amyloid peptides, revealed by solid-state NMR spectroscopy

A variety of biomolecules acting on the cell membrane folds into a biologically active structure in the membrane environment. It is, therefore, important to determine the structures and dynamics of such biomolecules in a membrane environment. While several biophysical techniques are used to obtain low-resolution information, solid-state NMR spectroscopy is one of the most powerful means for determining the structure and dynamics of membrane bound biomolecules such as antibacterial biomolecules and amyloidogenic proteins; unlike X-ray crystallography and solution NMR spectroscopy, applications of solid-state NMR spectroscopy are not limited by non-crystalline, non-soluble nature or molecular size of membrane-associated biomolecules. This review article focuses on the applications of solid-state NMR techniques to study a few selected antibacterial and amyloid peptides. Solid-state NMR studies revealing the membrane inserted bent α-helical structure associated with the hemolytic activity of bee venom melittin and the chemical shift oscillation analysis used to determine the transmembrane structure (with α-helix and 310-helix in the N- and C-termini, respectively) of antibiotic peptide alamethicin are discussed in detail. Oligomerization of an amyloidogenic islet amyloid polypeptide (IAPP, or also known as amylin) resulting from its aggregation in a membrane environment, molecular interactions of the antifungal natural product amphotericin B with ergosterol in lipid bilayers, and the mechanism of lipid raft formation by sphingomyelin studied using solid state NMR methods are also discussed in this review article. This article is part of a Special Issue entitled "Biophysical Exploration of Dynamical Ordering of Biomolecular Systems" edited by Dr. Koichi Kato.

Antimicrobial Peptide Structures: From Model Membranes to Live Cells

The rise in antibiotic resistance has led to a renewed interest in antimicrobial peptides (AMPs) that target membranes. The mode of action of AMPs involves the disruption of the lipid bilayer and leads to growth inhibition and death of the bacteria. However, details at the molecular level of how these peptides kill bacteria and the reasons for the observed differences in selectivity remain unclear. Structural information is crucial for defining the molecular mechanism by which these peptides recognize, self-assemble and interact with a particular lipid membrane. Solid-state NMR is a non-invasive technique that allows the study of the structural details of lipid-peptide and peptide-peptide interactions. Following on from studies of antibiotic and lytic peptides, gramicidin A and melittin, respectively, we investigated maculatin 1.1, an AMP from the skin of Australian tree frogs that acts against Gram-positive bacteria. By using perdeuterated phospholipids and specifically labelled peptides, ² H, ³¹ P and {³¹ P}¹⁵ N REDOR solid-state NMR experiments have been used to localize, maculatin 1.1 in neutral and anionic model membranes. However, the structure, location and activity depend on the composition of the model membrane and current advances in solid-state NMR spectroscopy now allow structure determination of AMPs in live bacteria.

Concepts and Methods of Solid-State NMR Spectroscopy Applied to Biomembranes

Concepts of solid-state NMR spectroscopy and applications to fluid membranes are reviewed in this paper. Membrane lipids with ²H-labeled acyl chains or polar head groups are studied using ²H NMR to yield knowledge of their atomistic structures in relation to equilibrium properties. This review demonstrates the principles and applications of solid-state NMR by unifying dipolar and quadrupolar interactions and highlights the unique features offered by solid-state ²H NMR with experimental illustrations. For randomly oriented multilamellar lipids or aligned membranes, solid-state ²H NMR enables direct measurement of residual quadrupolar couplings (RQCs) due to individual C-²H-labeled segments. The distribution of RQC values gives nearly complete profiles of the segmental order parameters S_CD⁽ⁱ⁾ as a function of acyl segment position (i). Alternatively, one can measure residual dipolar couplings (RDCs) for natural abundance lipid samples to obtain segmental S_CH order parameters. A theoretical mean-torque model provides acyl-packing profiles representing the cumulative chain extension along the normal to the aqueous interface. Equilibrium structural properties of fluid bilayers and various thermodynamic quantities can then be calculated, which describe the interactions with cholesterol, detergents, peptides, and integral membrane proteins and formation of lipid rafts. One can also obtain direct information for membrane-bound peptides or proteins by measuring RDCs using magic-angle spinning (MAS) in combination with dipolar recoupling methods. Solid-state NMR methods have been extensively applied to characterize model membranes and membrane-bound peptides and proteins, giving unique information on their conformations, orientations, and interactions in the natural liquid-crystalline state.

Gramicidin A disassembles large conductive clusters of its lysine-substituted derivatives in lipid membranes

N-terminally substituted lysine derivatives of gramicidin A (gA), [Lys1]gA and [Lys3]gA, but not glutamate- or aspartate-substituted peptides have been previously shown to cause the leakage of carboxyfluorescein from liposomes. Here, the leakage induction was also observed for [Arg1]gA and [Arg3]gA, while [His1]gA and [His3]gA were inactive at neutral pH. The Lys3-containing analogue with all tryptophans replaced by isoleucines did not induce liposome leakage, similar to gA. This suggests that the presence of both tryptophans and N-terminal cationic residues is critical for pore formation. Remarkably, the addition of gA blocked the leakage induced by [Lys3]gA. By examining with fluorescence correlation spectroscopy the peptide-induced leakage of fluorescent markers from liposomes, we estimated the diameter of pores responsible for the leakage to be about 1.6 nm. Transmission electron cryo-microscopy imaging of liposomes with [Lys3]gA showed that the liposomal membranes contained high electron density particles with a size of about 40 Å, suggesting the formation of peptide clusters. No such clusterization was observed in liposomes incorporating gA or a mixture of gA with [Lys3]gA. Three-dimensional reconstruction of the clusters was compatible with their pentameric arrangement. Based on experimental data and computational modeling, we suggest that the large pore formed by [Lys3]gA represents a barrel-stave oligomeric cluster formed by antiparallel double-stranded helical dimers (DH). In a tentative model, the pentamer of dimers may be stabilized by aromatic Trp-Trp and cation-π Trp-Lys interactions between the neighboring DHs. The inhibiting effect of gA on the [Lys3]gA-induced leakage can be attributed to breaking of cation-π interactions, which prevents peptide clusterization and pore formation.

An ensemble dynamics approach to decipher solid-state NMR observables of membrane proteins

Solid-state NMR (SSNMR) is an invaluable tool for determining orientations of membrane proteins and peptides in lipid bilayers. Such orientational descriptions provide essential information about membrane protein functions. However, when a semi-static single conformer model is used to interpret various SSNMR observables, important dynamics information can be missing, and, sometimes, even orientational information can be misinterpreted. In addition, over the last decade, molecular dynamics (MD) simulation and semi-static SSNMR interpretation have shown certain levels of discrepancies in terms of transmembrane helix orientation and dynamics. Dynamic fitting models have recently been proposed to resolve these discrepancies by taking into account transmembrane helix whole body motions using additional parameters. As an alternative approach, we have developed SSNMR ensemble dynamics (SSNMR-ED) using multiple conformer models, which generates an ensemble of structures that satisfies the experimental observables without any fitting parameters. In this review, various computational methods for determining transmembrane helix orientations are discussed, and the distributions of VpuTM (from HIV-1) and WALP23 (a synthetic peptide) orientations from SSNMR-ED simulations are compared with those from MD simulations and semi-static/dynamic fitting models. Such comparisons illustrate that SSNMR-ED can be used as a general means to extract both membrane protein structure and dynamics from the SSNMR measurements. This article is part of a Special Issue entitled: Membrane protein structure and function.

Solid-state C NMR spectroscopy of a C carbonyl-labeled polypeptide

High resolution structural elucidation of macromolecular structure by solid-state nuclear magnetic resonance requires the preparation of uniformly aligned samples that are isotopically labeled. In addition, to use the chemical shift interaction as a high resolution constraint requires an in situ tensor characterization for each site of interest. For (13)C in the peptide backbone, this characterization is complicated by the presence of dipolar coupled (14)N from the peptide bond. Here the (13)C(1)-Gly(2) site in gramicidin A is studied both as a dry powder and in a fully hydrated lipid bilayer environment. Linewidths reported for the oriented samples are a factor of five narrower than those reported elsewhere, and previous misinterpretations of the linewidths are corrected. The observed frequency from oriented samples is shown to be consistent with the recently determined structure for this site in the gramicidin backbone. It is also shown that, whereas a dipolar coupling between (13)C and (14)N is apparent in dry preparations of the polypeptide, in a hydrated bilayer the dipolar coupling is absent, presumably due to a ;self-decoupling' mechanism.

Application of the fuzzy-oil-drop model to membrane protein simulation

The analysis of structural properties and biological activity of membrane proteins requires long lasting simulation of molecular dynamics. The large number of atoms present in protein molecule, membrane (phospholipids), and water environment makes the simulation of large scale. The implementation of simplified model representing the natural environment for membrane proteins is presented and compared with the vacuum simulation and simulation in the presence of water molecules and membrane phospholipids presented explicite. The comparative structural analysis and computational times for these three models makes the simplified model promising.

Structural and dynamic studies of the γ-M4 trans-membrane domain of the nicotinic acetylcholine receptor

A structural characterization of a synthetic peptide corresponding to the fourth transmembrane domain (M4-TMD) of the gamma-subunit of the nicotinic acetylcholine receptor from Torpedo californica has been undertaken. Solid-state NMR and CD spectroscopy studies indicate that upon reconstitution into lipid vesicles or magnetically aligned lipid bilayers, the synthetic M4-TMD adopts a linear alpha-helical conformation with the helix aligned within 15 degrees of the membrane normal. Furthermore, analysis of the motional averaging of anisotropic interactions present in the solid-state NMR spectra of the reconstituted peptide, indicate that the dynamics of the peptide within the bilayer are highly sensitive to the phase adopted by the lipid bilayer, providing an insight into how the interaction of lipids with this domain may play a important role in the modulation of this receptor by its lipid environment.

Resolution enhancement in solid-state NMR of oriented membrane proteins by anisotropic differential linebroadening

We demonstrate that a significant improvement in the spectral resolution may be achieved in solid-state NMR experiments of proteins in inhomogeneously disordered oriented lipid bilayers. Using 1H homonuclear decoupling instead of standard 1H heteronuclear decoupling, the 15N line widths may be reduced by up to seven times for such samples. For large oriented membrane proteins, such resolution enhancements may be crucial for assignment and structural interpretation.

Sensitized photoinactivation of minigramicidin channels in bilayer lipid membranes

The method of sensitized photoinactivation based on the photosensitized damage of gramicidin A (gA) molecules was applied here to study ionic channels formed by minigramicidin (the 11-residue analogue of gramicidin A) in a planar bilayer lipid membrane (BLM) of different thickness. Irradiation of BLM with a single flash of visible light in the presence of a photosensitizer (aluminum phthalocyanine or Rose Bengal) generating singlet oxygen provoked a decrease in the minigramicidin-induced electric current across BLM, the kinetics of which had the characteristic time of several seconds, as observed with gA. For gA, there is good correlation between the characteristic time of photoinactivation and the single-channel lifetime. In contrast to the covalent dimer of gA characterized by extremely long single-channel lifetime and the absence of current relaxation upon flash excitation, the covalent head-to-head dimer of minigramicidin displayed the flash-induced current decrease with the kinetics being strongly dependent on the membrane thickness. The current decrease became slower both upon increasing the concentration of the minigramicidin covalent dimer and upon including cholesterol in the membrane composition. These data in combination with the quadratic dependence of the current on the peptide concentration can be rationalized by hypothesizing that the macroscopic current across BLM measured at high concentrations of the peptide is provided by dimers of minigramicidin covalent dimers in the double beta(5.7)-helical conformation having the lifetime of about 0.4 s, while single channels with the lifetime of 0.01 s, observed at a very low peptide concentration, correspond to the single-stranded beta(6.3)-helical conformation. Alternatively the results can be explained by clustering of channels at high concentrations of the minigramicidin covalent dimer.

The orientation effect of gramicidin A on bicelles and Eu3+-doped bicelles as studied by solid-state NMR and FT-IR spectroscopy

We have explored the effect of gramicidin A (gA) on bicelle (Bic) orientation in the absence and presence of Eu(3+) by (31)P and (2)H NMR at different DMPC/gA ratios. FT-IR spectroscopy was used to assess the lipid chain ordering and verify the transmembrane peptide conformation. Our results show a time-dependent flipping of the bilayer normal alignment at high temperatures and high proportion of gA. The results are explained by both the diamagnetic susceptibility anisotropy of the beta(6.3) helical peptides and viscosity of the lipid mixture. The concentration effect of gramicidin on Bic/Eu(3+) is compared to that on Eu(3+)-doped DMPC liposomes. The Bic/Eu(3+) system is no longer oriented in the presence of gA and adopts a vesicular morphology while the peptide incorporation induces the formation of ellipsoidal DMPC/Eu(3+) assemblies aligned with their normal parallel to the magnetic field. The difference is explained in terms of lipid chain disorder and size of the bilayers.

Structure and function of "metalloantibiotics"

Although most antibiotics do not need metal ions for their biological activities, there are a number of antibiotics that require metal ions to function properly, such as bleomycin (BLM), streptonigrin (SN), and bacitracin. The coordinated metal ions in these antibiotics play an important role in maintaining proper structure and/or function of these antibiotics. Removal of the metal ions from these antibiotics can cause changes in structure and/or function of these antibiotics. Similar to the case of "metalloproteins," these antibiotics are dubbed "metalloantibiotics" which are the title subjects of this review. Metalloantibiotics can interact with several different kinds of biomolecules, including DNA, RNA, proteins, receptors, and lipids, rendering their unique and specific bioactivities. In addition to the microbial-originated metalloantibiotics, many metalloantibiotic derivatives and metal complexes of synthetic ligands also show antibacterial, antiviral, and anti-neoplastic activities which are also briefly discussed to provide a broad sense of the term "metalloantibiotics."

Gramicidin derivatives as membrane-based pH sensors

Ion channels provide a means for sensitive pH measurement at membrane interfaces. Detailed knowledge of the structure and function of gramicidin channels permits the engineering of pH-sensitive derivatives. Two derivatives, gramicidin-ethylenediamine and gramicidin-histamine, are shown to exhibit pH-dependent single-channel behaviour over the pH ranges 9-11 and 6.5-8.5, respectively. Thermal isomerization of a carbamate group at the entrance of the channels leads to a pattern of steps in single-channel recordings. The size of the steps depends on the time-averaged degree of protonation of the appended group (ethylenediamine or histamine). Measurement of the size of the steps thus permits single-molecule pH sensing under symmetrical pH conditions or in the presence of a pH gradient.

Design and characterization of gramicidin channels

This article summarizes methods for the chemical synthesis and biophysical characterization of gramicidins with varying sequences and labels. The family of gramicidin channels has developed into a powerful model system for understanding fundamental properties, interactions, and dynamics of proteins and lipids generally, and ion channels specifically, in biological membranes.

The structure, dynamics and orientation of antimicrobial peptides in membranes by multidimensional solid-state NMR spectroscopy

Linear peptide antibiotics have been isolated from amphibians, insects and humans and used as templates to design cheaper and more potent analogues for medical applications. Peptides such as cecropins or magainins are < or = 40 amino acids in length. Many of them have been prepared by solid-phase peptide synthesis with isotopic labels incorporated at selected sites. Structural analysis by solid-state NMR spectroscopy and other biophysical techniques indicates that these peptide antibiotics strongly interact with lipid membranes. In bilayer environments they exhibit amphipathic alpha-helical conformations and alignments of the helix axis parallel to the membrane surface. This contrasts the transmembrane orientations observed for alamethicin or gramicidin A. Models that have been proposed to explain the antibiotic and pore-forming activities of membrane-associated peptides, as well as other experimental results, include transmembrane helical bundles, wormholes, carpets, detergent-like effects or the in-plane diffusion of peptide-induced bilayer instabilities.

The binding site of sodium in the gramicidin A channel

The available information concerning the structure and location of the main binding site for sodium in the gramicidin A channel is reviewed and discussed. Results from molecular dynamics simulations using an atomic model of the channel embedded in a lipid bilayer are compared with experimental observations. The combined information from experiment and simulation suggests that the main binding sites for sodium are near the channel's mouth, approximately 9.2 A from the centre of the dimer channel, although the motion along the axis could be as large as 1 to 2 A. In the binding site, the sodium ion is lying off axis, making contact with two carbonyl oxygens and two single-file water molecules. The main channel ligand is provided by the carbonyl group of the Leu10-Trp11 peptide linkage, which exhibits the largest deflection from the ion-free channel structure.

13C solid-state NMR of gramicidin A in a lipid membrane

The natural-abundance 13C NMR spectrum of gramicidin A in a lipid membrane was acquired under magic-angle spinning conditions. With fast sample spinning (15 kHz) at approximately 65 degrees C the peaks from several of the aliphatic, beta-, alpha-, aromatic, and carbonyl carbons in the peptide could be resolved. The resolution in the 13C spectrum was superior that observed with 1H NMR under similar conditions. The 13C linewidths were in the range 30-100 Hz, except for the alpha- and beta-carbons, the widths of which were approximately 350 Hz. The beta-sheet-like local structure of gramicidin A was observed as an upfield shift of the gramicidin alpha and carbonyl resonances. Under slow sample spinning (500 Hz), the intensity of the spinning sidebands from 13C in the backbone carbonyls was used to determine the residual chemical shift tensor. As expected, the elements of the residual chemical shift tensor were consistent with the single-stranded, right-handed beta6.3 helix structure proposed for gramicidin A in lipid membranes.

Characterization of the secondary structure and membrane interaction of the putative membrane anchor domains of prostaglandin I2 synthase and cytochrome P450 2C1

Prostaglandin I2 synthase (PGIS) produces prostaglandin I2 (PGI2) which has opposite actions on platelet aggregatory and vasoconstrictive properties compared to thromboxane A2 (TXA2) produced from the same substrate by another P450 enzyme, thromboxane A2 synthase (TXAS). PGIS and TXAS have only 16% amino acid sequence identity. Hydropathy analysis suggests that the putative NH2-terminal membrane anchor domain of PGIS is similar to many other membrane-bound microsomal P450s, which are believed to be anchored by a single transmembrane segment, and thus different from the TXAS anchor, which appears to have two transmembrane segments. To characterize the membrane anchor function of the PGIS NH2-terminal region, we have used the peptidoliposome reconstitution assay to identify the membrane anchor segment in the PGIS NH2-terminal domain and compared it with the anchor segment of P450 2C1. Four peptides, mimicking putative NH2-terminal membrane anchor segments of PGIS and P450 2C1, containing residues 1-28 (PGIS-LP1 and P450 2C1-LP1) or residues 25-54 (PGIS-LP2 and P450 2C1-LP2), were synthesized and their ability to insert in a lipid bilayer was evaluated. The results indicated that both LP1 peptides of PGIS and P450 2C1 became bound to the lipid bilayer, whereas both LP2 peptides did not bind the lipid. The two LP1 peptides were further characterized as to their conformation using CD spectroscopy. Helical structure induced in these peptides by addition of trifluoroethanol, dodecylphosphocholine, or incorporation into liposomes indicated that these segments tend to adopt a helical structure in a hydrophobic environment and thus could function as membrane anchor segments. These results support the hypothesis that PGIS and TXAS interact with the endoplasmic reticulum membrane in different ways, in which the NH2-terminal anchor domain of PGIS, as with P450 2C1, appears to have a single transmembrane segment.

The binding site of sodium in the gramicidin A channel: comparison of molecular dynamics with solid-state NMR data

The location of the main binding site for sodium in the gramicidin A (GA) channel was investigated with molecular dynamics simulations, using an atomic model of the channel embedded in a fully hydrated dimyristoyl phosphatidycholine (DMPC) bilayer. Twenty-four separate simulations in which a sodium was restrained at different locations along the channel axis were generated. The results are compared with carbonyl 13C chemical shift anisotropy solid-state NMR experimental data previously obtained with oriented GA:DMPC samples. Predictions are made for other solid-state NMR properties that could be observed experimentally. The combined information from experiment and simulation strongly suggests that the main binding sites for sodium are near the channel's mouth, approximately 9.2 A from the center of the dimer channel. The 13C chemical shift anisotropy of Leu10 is the most affected by the presence of a sodium ion in the binding site. In the binding site, the sodium ion is lying off-axis, making contact with two carbonyl oxygens and two single-file water molecules. The main channel ligand is provided by the carbonyl group of the Leu10-Trp11 peptide linkage, which exhibits the largest deviation from the ion-free channel structure. Transient contacts with the carbonyl group of Val8 and Trp15 are also present. The influence of the tryptophan side chains on the channel conductance is examined based on the current information about the binding site.

Structure, energetics, and dynamics of lipid-protein interactions: A molecular dynamics study of the gramicidin A channel in a DMPC bilayer

The microscopic details of lipid-protein interactions are examined using molecular dynamics simulations of the gramicidin A channel embedded in a fully hydrated dimyristoyl phosphatidylcholine (DMPC) bilayer. A novel construction protocol was used to assemble the initial configurations of the membrane protein complex for the simulations. Three hundred systems were constructed with different initial lipid placement and conformations. Seven systems were simulated with molecular dynamics. One system was simulated for a total of 600 psec, four were simulated for 300 psec, and two for 100 psec. Analysis of the resulting trajectories shows that the bulk solvent-membrane interface region is much broader than traditionally pictured in simplified continuum theories: its width is almost 15 angstroms. In addition, lipid-protein interactions are far more varied, both structurally and energetically, than is usually assumed: the total interaction energy between the gramicidin A and the individual lipids varies from 0 to -50 kcal/mol. The deuterium quadrupolar splittings of the lipid acyl chains calculated from the trajectories are in good agreement with experimental data. The lipid chains in direct contact with the GA are ordered but the effect is not uniform due to the irregular surface of the protein. Energy decompositions shows that the most energetically favorable interactions between lipid and protein involve nearly equal contributions from van der Waals and electrostatic interactions. The tryptophans, located near the bulk-membrane interface, appear to be particularly important in mediating both hydrogen bonding interactions with the lipid glycerol backbone and water and also in forming favorable van der Waals contacts with the hydrocarbon chains. In contrast, the interactions of the leucine residues with the lipids, also located near the interface, are dominated by van der Waals interactions with the hydrocarbon lipid chains.

High resolution 1H nuclear magnetic resonance of a transmembrane peptide

Although the strong 1H-1H dipolar interaction is known to result in severe homogeneous broadening of the 1H nuclear magnetic resonance (NMR) spectra of ordered systems, in the fluid phase of biological and model membranes the rapid, axially symmetric reorientation of the molecules about the local bilayer normal projects the dipolar interaction onto the motional symmetry axis. Because the linewidth then scales as (3 cos2 theta-1)/2, where theta is the angle between the local bilayer normal and the magnetic field, the dipolar broadening has been reduced to an "inhomogeneous" broadening by the rapid axial reorientation. It is then possible to obtain high resolution 1H-NMR spectra of membrane components by using magic angle spinning (MAS). Although the rapid axial reorientation effectively eliminates the homogeneous dipolar broadening, including that due to n = 0 rotational resonances, the linewidths observed in both lipids and peptides are dominated by low frequency motions. For small peptides the most likely slow motions are either a "wobble" or reorientation of the molecular diffusion axis relative to the local bilayer normal, or the reorientation of the local bilayer normal itself through surface undulations or lateral diffusion over the curved surface. These motions render the peptide 1H-NMR lines too broad to be observed at low spinning speeds. However, the linewidths due to these slow motions are very sensitive to spinning rate, so that at higher speeds the lines become readily visible. The synthetic amphiphilic peptide K2GL20K2A-amide (peptide-20) has been incorporated into bilayers of 1,2-di-d 27-myristoyl-sn-glycero-3-phosphocholine (DMPC-d54) and studied by high speed 1H-MAS-NMR. The linewidths observed for this transbilayer peptide, although too broad to be observable at spinning rates below -5 kHz, are reduced to 68 Hz at a spinning speed of 14 kHz (at 500C). Further improvements in spinning speed and modifications in sample composition designed to reduce the effectiveness of the slow motions responsible for the linewidth should result in significant further reduction in peptide linewidths. With this technique, there is now the potential for the use of 1H-MAS-NMR for the study of conformation, folding, and dynamics of small membrane peptides and protein fragments.

Species heterogeneity of Gly-11 gramicidin A incorporated into sodium dodecyl sulfate micelles

Evidence is presented for species heterogeneity of the gly-11 analog of gramicidin A incorporated into sodium dodecyl sulfate (SDS) micelles. The evidence for species heterogeneity has been obtained using one-dimensional (1D) 1H NMR spectroscopy. The 1D spectra of the indole NH moiety of tryptophans 9, 13, and 15 show the presence of more than one species. It has been found that the heterogeneity is dependent upon the gly-11/SDS molar ratio. At high SDS concentration (i.e., gly-11/SDS of 3 mM/700 mM) the heterogeneity almost completely disappears. The temperature dependence of these 1H NMR signals suggests that the two species do not interconvert. The results of nuclear Overhauser effect spectroscopy NMR experiments indicate that one species is embedded within the micelle, while the other is nearer the aqueous interface. The importance of side chain interactions with the membrane environment in producing stable, solubilized species of small peptides in SDS micelles is illustrated.

Ion transport in the gramicidin channel: molecular dynamics study of single and double occupancy

The structural and thermodynamic factors responsible for the singly and doubly occupied saturation states of the gramicidin channel are investigated with molecular dynamics simulations and free energy perturbation methods. The relative free energy of binding of all of the five common cations Li+, Na+, K+, Rb+, and Cs+ is calculated in the singly and doubly occupied channel and in bulk water. The atomic system, which includes the gramicidin channel, a model membrane made of neutral Lennard-Jones particles and 190 explicit water molecules to form the bulk region, is similar to the one used in previous work to calculate the free energy profile of a Na+ ion along the axis of the channel. In all of the calculations, the ions are positioned in the main binding sites located near the entrances of the channel. The calculations reveal that the doubly occupied state is relatively more favorable for the larger ions. Thermodynamic decomposition is used to show that the origin of the trend observed in the calculations is due to the loss of favorable interactions between the ion and the single file water molecules inside the channel. Small ions are better solvated by the internal water molecules in the singly occupied state than in the doubly occupied state; bigger ions are solvated almost as well in both occupation states. Water-channel interactions play a role in the channel response. The observed trends are related to general thermodynamical properties of electrolyte solutions.

Molecular dynamics simulation of the gramicidin channel in a phospholipid bilayer

A molecular dynamics simulation of the gramicidin A channel in an explicit dimyristoyl phosphatidylcholine bilayer was generated to study the details of lipid-protein interactions at the microscopic level. Solid-state NMR properties of the channel averaged over the 500-psec trajectory are in excellent agreement with available experimental data. In contrast with the assumptions of macroscopic models, the membrane/solution interface region is found to be at least 12 A thick. The tryptophan side chains, located within the interface, are found to form hydrogen bonds with the ester carbonyl groups of the lipids and with water, suggesting their important contribution to the stability of membrane proteins. Individual lipid-protein interactions are seen to vary from near 0 to -50 kcal/mol. The most strongly interacting conformations are short-lived and have a nearly equal contribution from both van der Waals and electrostatic energies. This approach for performing molecular dynamics simulations of membrane proteins in explicit phospholipid bilayers should help in studying the structure, dynamics, and energetics of lipid-protein interactions.

Sodium ion binding in the gramicidin A channel. Solid-state NMR studies of the tryptophan residues

Gramicidin A analogs, labeled with 13C in the backbone carbonyl groups and the C-2 indole carbons of the tryptophan-11 and tryptophan-13 residues, were synthesized using t-Boc-protected amino acids. The purified analogs were incorporated into phosphatidylcholine bilayers at a 1:15 molar ratio and macroscopically aligned between glass coverslips. The orientations of the labeled groups within the channel were investigated using solid-state NMR and the effect of a monovalent ion (Na+) on the orientation of these groups determined. The presence of sodium ions did not perturb the 13C spectra of the tryptophan carbonyl groups. These results contrast with earlier results in which the Leu-10, Leu-12, and Leu-14 carbonyl groups were found to be significantly affected by the presence of sodium ions and imply that the tryptophan carbonyl groups are not directly involved in ion binding. The channel form of gramicidin A has been demonstrated to be the right-handed form of the beta 6.3 helix: consequently, the tryptophan carbonyls would be directed away from the entrance to the channel and take part in internal hydrogen bonding, so that the presence of cations in the channel would have less effect than on the outer leucine residues. Sodium ions also had no effect on the C-2 indole resonance of the tryptophan side chains. However, a small change was observed in Trp-11 when the ether lipid, ditetradecylphosphatidylcholine, was substituted for the ester lipid, dimyristoylphosphatidylcholine, indicating some sensitivity of the gramicidin side chains to the surrounding lipid.

Conformation states of gramicidin A along the pathway to the formation of channels in model membranes determined by 2D NMR and circular dichroism spectroscopy

Gramicidin A incorporated into SDS (sodium dodecyl sulfate) micelles exists as a right-handed, N-to-N-terminal beta 6.3 helical dimer [Lomize, A. L., Orechov, V. Yu., & Arseniev, A.S. (1992) Bioorg. Khim. 18, 182-189]. In the incorporation procedure to achieve the ion channel state of gramicidin A in SDS micelles, trifluoroethanol (TFE) is used to solubilize the hydrophobic peptide before addition to the aqueous/micelle solution. The conformational transition of gramicidin A to form ion channels in SDS micelles, i.e., in TFE and 10% TFE/water, has been investigated using 2D NMR and CD spectroscopy. In neat TFE, gramicidin A was found to be monomeric and may possibly exist in an equilibrium of rapidly interconverting conformers of at least three different forms believed to be left- and/or right-handed alpha and beta 4.4 helices. It was found that the interconversion between these conformers was slowed down in 55% TFE as evident by the observation of at least three different sets of d alpha N COSY peaks although CD gave a net spectrum similar to that in neat TFE. In 10% TFE gramicidin A spontaneously forms a precipitate. The precipitated species were isolated and solubilized in dioxane where gramicidin conformers undergo very slow interconversion and could be characterized by NMR. At least seven different gramicidin A conformations were found in 10% TFE. Four of thes are the same types of double helices as previously found in ethanol (i.e., a symmetric left-handed parallel beta 5.6 double helix, an unsymmetric left-handed parallel beta 5.6 double helix, a symmetric left-handed antiparallel beta 5.6 double helix, a symmetric right-handed parallel beta 5.6 double helix); the fifth is possibly a symmetric right-handed antiparallel beta 5.6 double helix. There is also evidence for the presence of at least one form of monomeric species. Previous observation on the solvent history dependence in the ease of channel incorporation may be explained by the presence of several different folding pathways to channel formation. To test this proposal, the conformation of gramicidin A in 10% DMSO and 10% methanol was studied. In the former environment, the major form was a random coil with a minor population of double-stranded helices, while in the latter, NMR spectra indicate the presence of the same double-helical conformers as found in neat methanol.

Alpha-helical hydrophobic polypeptides form proton-selective channels in lipid bilayers

Proton translocation is important in membrane-mediated processes such as ATP-dependent proton pumps, ATP synthesis, bacteriorhodopsin, and cytochrome oxidase function. The fundamental mechanism, however, is poorly understood. To test the theoretical possibility that bundles of hydrophobic alpha-helices could provide a low energy pathway for ion translocation through the lipid bilayer, polyamino acids were incorporated into extruded liposomes and planar lipid membranes, and proton translocation was measured. Liposomes with incorporated long-chain poly-L-alanine or poly-L-leucine were found to have proton permeability coefficients 5 to 7 times greater than control liposomes, whereas short-chain polyamino acids had relatively little effect. Potassium permeability was not increased markedly by any of the polyamino acids tested. Analytical thin layer chromatography measurements of lipid content and a fluorescamine assay for amino acids showed that there were approximately 135 polyleucine or 65 polyalanine molecules associated with each liposome. Fourier transform infrared spectroscopy indicated that a major fraction of the long-chain hydrophobic peptides existed in an alpha-helical conformation. Single-channel recording in both 0.1 N HCl and 0.1 M KCl was also used to determine whether proton-conducting channels formed in planar lipid membranes (phosphatidylcholine/phosphatidylethanolamine, 1:1). Poly-L-leucine and poly-L-alanine in HCl caused a 10- to 30-fold increase in frequency of conductive events compared to that seen in KCl or by the other polyamino acids in either solution. This finding correlates well with the liposome observations in which these two polyamino acids caused the largest increase in membrane proton permeability but had little effect on potassium permeability. Poly-L-leucine was considerably more conductive than poly-L-alanine due primarily to larger event amplitudes and, to a lesser extent, a higher event frequency. Poly-L-leucine caused two populations of conductive events, one in the 0.1-0.5 pA range, and one in the 1.0-5.0 pA range, whereas nearly all events caused by poly-L-alanine were in the 0.1-0.5 pA range at an applied voltage of +60 mV. The channel-like activity appeared to switch between conductive and nonconductive states, with most open-times in the range of 50-200 ms. We conclude that hydrophobic polyamino acids produce proton-conducting defects in lipid bilayers that may be used to model functional proton channels in biological membranes.

The structure of an integral membrane peptide: a deuterium NMR study of gramicidin

Solid state deuterium NMR was employed on oriented multilamellar dispersions consisting of 1,2-dilauryl-sn-glycero-3-phosphatidylcholine and deuterium (2H) exchange-labeled gramicidin D, at a lipid to protein molar ratio (L/P) of 15:1, in order to study the dynamic structure of the channel conformation of gramicidin in a liquid crystalline phase. The corresponding spectra were used to discriminate between several structural models for the channel structure of gramicidin (based on the left- and right-handed beta 6.3 LD helix) and other models based on a structure obtained from high resolution NMR. The oriented spectrum is complicated by the fact that many of the doublets, corresponding to the 20 exchangeable sites, partially overlap. Furthermore, the asymmetry parameter, eta, of the electric field gradient tensor of the amide deuterons is large (approximately 0.2) and many of the amide groups are involved in hydrogen bonding, which is known to affect the quadrupole coupling constant. In order to account for these complications in simulating the spectra in the fast motional regime, an ab initio program called Gaussian 90 was employed, which permitted us to calculate, by quantum mechanical means, the complete electric field gradient tensor for each residue in gramicidin (using two structural models). Our results indicated that the left-handed helical models were inconsistent with our observed spectra, whereas a model based on the high-resolution structure derived by Arseniev and coworkers, but relaxed by a simple energy minimization procedure, was consistent with our observed spectra. The molecular order parameter was then estimated from the motional narrowing assuming the relaxed (right-handed) Arseniev structure. Our resultant order parameter of SZZ = 0.91 translates into an rms angle of 14 degrees, formed by the helix axis and the local bilayer normal. The strong resemblance between our spectra (and also those reported for gramicidin in 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC) multilayers) and the spectra of the same peptide incorporated in a lyotropic nematic phase, suggests that the lyotropic nematic phase simulates the local environment of the lipid bilayer.

Two-dimensional rotational-echo double resonance of Val1-[1-13C]Gly2-[15N]Ala3-gramicidin A in multilamellar dimyristoylphosphatidylcholine dispersions

The dipolar coupling between the Gly2 13C-1 carbon and Ala3 15N-amide nitrogen was used to investigate the conformation and dynamics of the Gly2-Ala3 13C-15N peptide bond in Val1-[1-13C]Gly2-[15N]Ala3-gramicidin A incorporated into multilamellar dispersions of dimyristoylphosphatidylcholine. Measurement of the 13C-15N dipolar coupling constant D of the labeled gramicidin in a powder and the effective dipolar coupling constant De in a multilamellar dispersion was accomplished by two-dimensional rotational-echo double-resonance (2D REDOR) NMR, a magic-angle spinning experiment designed to measure weak dipolar coupling constants. The magnitudes of D and De were measured by the mirror-symmetric form of 2D REDOR, and the signs of D and De were determined relative to the sign of the isotropic indirect spin-spin coupling constant J by the mirror-asymmetric form of 2D REDOR. From knowledge of the magnitudes of D and De, four possible values were calculated for the angle between the Gly2-Ala3 13C-15N peptide bond and the gramicidin helical axis. Additional knowledge of the signs of D and De permitted the set of possible values for the peptide bond angle to be reduced to a single angle and its supplement (64 degrees, 116 degrees). This information about the Gly2-Ala3 13C-15N peptide bond angle eliminates the double-stranded, helical dimers and the left-handed, single-stranded, beta 6.3 helical dimer but supports the right-handed, single-stranded, beta 6.3 helical dimer as the structural model for gramicidin in multilamellar dispersions.

2H NMR lineshapes of immobilized uniaxially oriented membrane proteins

As a method for the structure determination of integral membrane proteins or other large macromolecular complexes, a solid state 2H NMR approach is presented, capable of measuring the orientations of individual chemical bond vectors. In an immobilized uniaxially oriented sample, the bond angle of a deuterium-labelled methyl group relative to the axis of ordering can be calculated from the quadrupole splitting in the "zero-tilt" spectrum where the sample normal is aligned parallel to the spectrometer field direction. However, since positive and negative values of this splitting cannot be distinguished, there may appear to be two solutions, of which only one describes the correct molecular geometry. We show that it is possible to determine the bond angle uniquely between 0 degree and 90 degrees, by analysing the lineshapes of a tilt series of spectra acquired over different sample inclinations. The lineshape equation describing such oriented 2H NMR spectra will be derived (for asymmetry parameter eta = 0) and discussed, with an illustration of the various linebroadening effects from which the orientational distribution function in the macroscopically ordered system can be determined. This strategy is then applied to specifically deuterium-labelled retinal in dark-adapted bacteriorhodopsin, prepared in a uniaxially oriented sample from purple membrane fragments. From the quadrupole splitting in the zero-tilt spectrum and by lineshape simulations, the deuteromethyl group at C20 on retinal is found to make an angle of 32 degrees +/- 1 degree with the membrane normal, and the sample mosaic spread to be around +/- 8 degrees. The resulting orientation of retinal is in excellent agreement with its known structure in bacteriorhodopsin, and together with the results on other methyl groups it will be possible to construct a detailed picture of the chromophore in the protein binding pocket.

Solid state 13C NMR of unlabeled phosphatidylcholine bilayers: spectral assignments and measurement of carbon-phosphorus dipolar couplings and 13C chemical shift anisotropies

The direct measurement of 13C chemical shift anisotropies (CSA) and 31P-13C dipolar splitting in random dispersions of unlabeled L alpha-phase phosphatidylcholine (PC) has traditionally been difficult because of extreme spectral boradening due to anisotropy. In this study, mixtures of dimyristoyl phosphatidylcholine (DMPC) with three different detergents known to promote the magnetic orientation of DMPC were employed to eliminate the powder-pattern nature of signals without totally averaging out spectral anisotropy. The detergents utilized were CHAPSO, Triton X-100, and dihexanoylphosphatidylcholine (DHPC). Using such mixtures, many of the individual 13C resonances from DMPC were resolved and a number of 13C-31P dipolar couplings were evident. In addition, differing line widths were observed for the components of some dipolar doublets, suggestive of dipolar/chemical shift anisotropy (CSA) relaxation interference effects. Oriented sample resonance assignments were made by varying the CHAPSO or DHPC to DMPC ratio to systematically scale overall bilayer order towards the isotropic limit. In this manner, peaks could be identified based upon extrapolation to their isotropic positions, for which assignments have previously been made (Lee, C.W.B., and R.G. Griffin. 1989. Biophys. J. 55:355-358; Forbes, J., J. Bowers, X. Shan, L. Moran, E. Oldfield, and M.A. Moscarello. 1988. J. Chem. Soc., Faraday, Trans. 1 84:3821-3849). It was observed that the plots of CSA or dipolar coupling versus overall bilayer order obtained from DHPC and CHAPSO titrations were linear. Estimates of the intrinsic dipolar couplings and chemical shift anisotropies for pure DMPC bilayers were made by extrapolating shifts and couplings from the detergent titrations to zero detergent. Both detergent titrations led to similar "intrinsic" CSAs and dipolar couplings. Results extracted from an oriented Triton-DMPC mixture also led to similar estimates for the detergent-free DMPC shifts and couplings. The results from these experiments were found to compare favorably with limited measurements made from pure L alpha PC. This detergent-based method for assigning spectra and for determining dipolar couplings and CSA in detergent-free systems should be extendable to other lipid systems. The resulting data set from this study may prove useful in future modeling of the structure and dynamics of DMPC bilayers. In addition, the fact that experiments utilizing each of the three detergents led to similar estimates for the spectral parameters of pure DMPC, and the fact that spectral parameter versus bilayer order plots were linear, indicate that the averaged conformation and dynamics of DMPC in the presence of the three detergents are very similar to those of pure L alpha DMPC.

Energy-minimized conformation of gramicidin-like channels. I. Infinitely long poly-(L,D)-alanine beta 6.3-helix

The energy-minimized conformation of an infinitely long poly-(L,D)-alanine in single-stranded beta 6.3-helix was calculated by the molecular mechanics method. When energy minimization was started from a wide range of initial geometries, six optimized conformations were obtained and identified as the right- and left-handed counterparts of the beta 4.5-, beta 6.3-, and beta 8.2-helices. It was found that their conformation energies increase in this order, the beta 4.5-helix having the lowest energy. The backbone dihedral angles of the energy-minimized beta 6.3-helix were: phi L = -116 degrees (or -131 degrees), psi L = 122 degrees (or 111 degrees), phi D = 131 degrees (or 116 degrees), psi D = -111 degrees (or -122 degrees), omega L = 173 degrees (or 173 degrees), and omega D = -173 degrees (or -173 degrees) for the right-handed (or left-handed) helix. This helix was composed of 6.30 residues/turn with a pitch of 4.97 A. All the alpha-carbons of L- and D-configurations appeared on one common circular helix. Interestingly, small deviations (approximately 7 degrees) of the peptide bonds from the planar structure caused a considerable lowering of the conformation energy, and, at the same time, they produced more favorable fitting of the hydrogen bonds; the carbonyl oxygens and the nearest-neighbor alpha-hydrogens also took more favorable relative positions.

Model ion channels: gramicidin and alamethicin

We have discussed in some detail a variety of experimental studies which were designed to elucidate the conformational and dynamic properties of gramicidin and alamethicin. Although the behavior of these peptides is by no means fully characterized, these studies have already permitted aspects of ion channel activity to be understood in molecular terms. Studies with gramicidin in a variety of organic solutions have revealed conformational heterogeneity of this peptide; at least five major isomers exist, several of which have been characterized in detail using NMR spectroscopy and X-ray crystallography. When added to lipid membranes gramicidin undergoes a further conformational conversion. Although the conformation of gramicidin in membranes is not as well characterized as the solution conformation(s) and an X-ray structure is not yet available, detailed data, particularly from solid-state NMR studies, continue to become available and a right-handed beta 6.3 helical conformation of the peptide backbone is now generally accepted. Two of these beta 6.3 helices joined at their N-termini are believed to form the conducting channel. The conformational behavior of the side-chains of gramicidin in the membrane-bound form is not well established and several NMR, CD, fluorescence and theoretical studies are now focussed on this. Although the side-chains do not directly contact the permeating ions, they can have distinct effects on conductance and selectivity by altering the electrostatic environment sensed by the ion. The dynamics of both side-chain and backbone conformations of gramicidin appear critical to a detailed understanding of the ion transport process in this channel. As the description of the membrane-bound conformation of gramicidin becomes more detailed, simulations of ion transport using computational methods are likely to improve and will further our understanding of the processes of ion transport. As well as internal motion of the backbone and side-chains, gramicidin undergoes rotational and translational motion in the plane of the membrane. These motions do not appear to be essential for the process of ion transport but can affect channel lifetime since lifetime is determined by the rate of association and dissociation of gramicidin monomers. Gramicidin-membrane interactions are also likely to be involved in the frequency of occurrence of channel subconductance states, the frequency of channel flickering and fundamentally in the stability of the membrane-bound gramicidin conformation. Alamethicin forms channels in membranes which are strongly voltage-dependent. The molecular origin of voltage-dependent conductances has been a fundamental problem in biophysics for many years.(ABSTRACT TRUNCATED AT 400 WORDS)