A lipid dependence in the formation of twin ion channels

A gramicidin A derivative with a polyether linkage between both ethanolamine termini was synthesized and its ion channel properties were studied. The compound showed a duplication in the state of conductance for alkali cations in thick DOPC bilayer membranes, which is interpreted as the occurrence of twin-channels. In thinner DMPC membranes mono-channels were dominant. The influence of hydrophobic coupling on the mono channel/twin channel equilibrium is discussed.     

Proton relaxation mechanisms and the measurements of r phi, r psi and transannular interproton distances in gramicidin S

Monoselective, Rio(SE), biselective, Rio(i,j), and nonselective proton spin-lattice relaxation rates have been measured for dilute solutions of gramicidin S in dimethyl sulfoxide and used to evaluate cross-relaxation rates (sigma ij = Rio(i,j)-Rio(SE)) and Fi ratios (Fi = Ri(NS)/Rio(SE)). The cross-relaxation parameters, sigma, and Fi ratios measured for backbone gramicidin S protons predict that the same correlation time, tau c = 1.2 X 10(-9)s, modulates all the dipolar proton-proton interactions and that these interactions represent the main source for the proton spin-lattice relaxation process. The larger relaxation rates for amide versus alpha-protons of the backbone are attributed to dipolar relaxation between 14N and its directly bonded protons and is an approximate measure of the extent of this. The intrabackbone proton-proton distances, evaluated from sigma values, were consistent with the antiparallel beta-plated sheet/beta II'-turn conformation previously proposed for gramicidin S in solution.     

Comparison of the gramicidin a potassium/sodium permeability and single channel conductance ratio

If the ion concentration is low enough that most channels are unoccupied, then the ‘independence relations’ should be satisfied and the permeability ratio should equal the conductance ratio. It has been previously reported that for the gramicidin A channel these ratios for Na⁺ and K⁺ were not equal at concentrations as low as 10 mM. However, these ratios were not measured at the same applied potential, as is required by the theory. Instead, the conductance ratio was measured at 100 mV and corrected using calculated current-voltage relations. In this report the comparison between permeability and conductance ratios is reexamined using data obtained at the correct potential. There is no significant difference in the ratios at 10 mM when they are measured at the same voltage. This implies that most channels are not occupied by sodium or potassium ions at 10 mM.     

Constructing the suitable initial configuration of the membrane-protein system in molecular dynamics simulations

A method for constructing the suitable initial configuration of the membrane-protein system for molecular dynamics (MD) simulations is presented. This method could provide some hydrated initial configurations and help us to determine the best surface area of the system by contracting the surface area and comparing the optimized lowest energy of the system by energy minimization. The gramicidin A (GA) channel in;the fully hydrated dimyristoylphosphatidylcholine (DMPC) bilayer was used as our model. Three configurations with different surface areas were selected and applied for one 400 ps and two 300 ps MD simulations at constant pressure and temperature. All simulations were fairly stable without any constraints. Through analysis of the MD trajectories we found that the system with the best surface area was more stable than the other two systems, whose sizes were changed in the simulations. Further analysis of the bilayer normal length and the order parameters of the lipid alkyl tails indicates that the system with the best surface area shows some characteristics of the L_α phase, while both the smaller and the larger size systems have distinct deviations from the L_α phase that we expect. This illustrates that the correct surface area and the suitable initial configuration have an important influence on the phase of the membrane in the MD simulation. In addition, by comparing the root mean square differences of GA relative to the initial structure and interaction energy between different components of the system for all three systems, we find that the state of the DMPC bilayer has exerted a significant influence on the structure of GA. All these results demonstrate the validity of our method for constructing the initial configuration of the membrane-protein system for MD simulations. Received: 10 September 1998 / Revised version: 19 March 1999 / Accepted: 19 March 1999     

Test of molecular dynamics force fields in gramicidin A

The force fields commonly used in molecular dynamics simulations of proteins are optimized under bulk conditions. Whether the same force fields can be used in simulations of membrane proteins is not well established, although they are increasingly being used for such purposes. Here we consider ion permeation in the gramicidin A channel as a test of the AMBER force field in a membrane environment. The potentials of mean force for potassium ions are calculated along the channel axis and compared with the one deduced from the experimental conductance data. The calculated result indicates a rather large central barrier similar to those obtained from other force fields, which are incompatible with the conductance data. We suggest that lack of polarizability is the most likely cause of this problem, and, therefore, urge development of polarizable force fields for simulations of membrane proteins.     

Chemical and Photochemical Modification of Colicin E1 and Gramicidin A in Bilayer Lipid Membranes

Chemical modification and photodynamic treatment of the colicin E1 channel-forming domain (P178) in vesicular and planar bilayer lipid membranes (BLMs) was used to elucidate the role of tryptophan residues in colicin E1 channel activity. Modification of colicin tryptophan residues by N-bromosuccinimide (NBS), as judged by the loss of tryptophan fluorescence, resulted in complete suppression of wild-type P178 channel activity in BLMs formed from fully saturated (diphytanoyl) phospholipids, both at the macroscopic-current and single-channel levels. The similar effect on both the tryptophan fluorescence and the electric current across BLM was observed also after NBS treatment of gramicidin channels. Of the single-tryptophan P178 mutants studied, W460 showed the highest sensitivity to NBS treatment, pointing to the importance of the water-exposed Trp460 in colicin channel activity. In line with previous work, the photodynamic treatment (illumination with visible light in the presence of a photosensitizer) led to suppression of P178 channel activity in diphytanoyl-phospholipid membranes concomitant with the damage to tryptophan residues detected here by a decrease in tryptophan fluorescence. The present work revealed novel effects: activation of P178 channels as a result of both NBS and photodynamic treatments was observed with BLMs formed from unsaturated (dioleoyl) phospholipids. These phenomena are ascribed to the effect of oxidative modification of double-bond-containing lipids on P178 channel formation. The pronounced stimulation of the colicin-mediated ionic current observed after both pretreatment with NBS and sensitized photomodification of the BLMs support the idea that distortion of membrane structure can facilitate channel formation.Abbreviations: AlP_cS₃, almininum trisulfophthalocyanine; BLM, bilayer lipid membrane; DOPC, dioleoylphosphatidylcholine; DOPG, dioleoylphosphatidyl-glycerol; DPhPG, diphytanoylphos-phatidylglycerol; DPhPg, diphytanoylphosphatidylcholine; gA, gramicidin A; NBS, N-bromosuccinimideThis revised version was published online in August 2005 with a corrected cover date.     

Difference in association of Tl(I) with gramicidin A and gramicidin B in trifluoroethanol determined by Tl-205 NMR

The thallium-205 chemical shift was determined as a function of temperature for the thallium(I) complexes of gramicidin A and gramicidin B in 2,2,2-trifluoroethanol. From the difference in magnitude of the induced chemical shift it was determined that gramicidin B does not bind the Tl(I) ion as well as does gramicidin A. This result may explain the lower single-channel conductance of gramicidin B relative to gramicidin A. Cabon-13 NMR studies strongly indicate that the binding site for gramicidin A and B is at teh tryptophan end of the molecule and that replacement of tryptophan residue at position 11 in gramicidin A with a phenylalanine to form gramicidin B produces a significant structural change at the tryptophan end of the molecule, but has little effect on the N-terminus.     

Macromolecular structural elucidation with solid-state NMR-derived orientational constraints

Summary The complete structure determination of a polypeptide in a lipid bilayer environment is demonstrated built solely upon orientational constraints derived from solid-state NMR observations. Such constraints are obtained from isotopically labeled samples uniformly aligned with respect to the B₀ field. Each observation constrains the molecular frame with respect to B₀ and the bilayer normal, which are arranged to be parallel. These constraints are not only very precise (a few tenths of a degree), but also very accurate. This is clearly demonstrated as the back bone structure is assembled sequentially and the i to i+6 hydrogen bonds in this structure of the gramicidin channel are shown on average to be within 0.5 Å of ideal geometry. Similarly, the side chains are assembled independently and in a radial direction from the backbone. The lack of considerable atomic overlap between side chains also demonstrates the accuracy of the constraints. Through this complete structure, solid-state NMR is demonstrated as an approach for determining three-dimensional macromolecular structure.     

Conformational studies on the gramicidin A transmembrane channel in lipid micelles and liposomes

The interaction of gramicidin A with lysolecithin micelles and with lecithin liposomes is demonstrated by circular dichroism to result in several metastable conformational states. A stable state can be obtained after extensive heating when the gramicidin A was added dry or in ethanol solution to the phospholipid dispersion but the stable state is readily obtained when gramicidin A is added in a trifluoroethanol solution. The circular dichroism of the stable conformational state is characterized by negative ellipticity below 205 nm and principally by a positive 220 nm band on which is superposed a weak 230 nm band (the latter likely arising from tryptophan side chains). The stable conformational state is considered to be that of the functional transmembrane channel primarily on the basis of extensive studies on its interaction with sodium ions.     

Influence of protein flexibility on the electrostatic energy landscape in gramicidin A

We describe an electrostatic model of the gramicidin A channel that allows protein atoms to move in response to the presence of a permeating ion. To do this, molecular dynamics simulations are carried out with a permeating ion at various positions within the channel. Then an ensemble of atomic coordinates taken from the simulations are used to construct energy profiles using macroscopic electrostatic calculations. The energy profiles constructed are compared to experimentally-determined conductance data by inserting them into Brownian dynamics simulations. We find that the energy landscape seen by a permeating ion changes significantly when we allow the protein atoms to move rather than using a rigid protein structure. However, the model developed cannot satisfactorily reproduce all of the experimental data. Thus, even when protein atoms are allowed to move, the dielectric model used in our electrostatic calculations breaks down when modeling the gramicidin channel.     

A genomic reservoir for Tnfrsf genes is developmentally regulated and imprinted in the mouse

Tumor necrosis factor receptor superfamily is composed of at least 26 members in the mouse, three of which exist as a cluster within the imprinted Kcnq1 domain on chromosome 7. Tnfrsf22, 23 and 26 contain typical cystein-rich domains and Tnfrsf22 and 23 can bind ligands but have no signaling capacity. Thus, they are assumed to be decoy receptors. The developmental expression profile of these genes is unknown and knowledge of their imprinting patterns is incomplete and controversial. We found that all three genes are expressed during mouse embryonic development, and that they have a strong maternal bias, indicating that they may be affected by the KvDMR, the Kcnq1 imprinting control region. We found expression of an antisense non-coding RNA, {"type":"entrez-nucleotide","attrs":{"text":"AK155734","term_id":"74212930","term_text":"AK155734"}}AK155734, in embryos and some neonatal tissues. This RNA overlaps the Tnfrsf22 and possibly the Tnfrsf23 coding regions and is also expressed with a maternal bias. We were interested in exploring the evolutionary origins of the three Tnfrsf genes, because they are absent in the orthologous human Kcnq1 domain. To determine whether the genes were deleted from humans or acquired in the rodent lineage, we performed phylogenetic analyses. Our data suggest that TNFRSF sequences were duplicated and/or degenerated or eliminated from the KCNQ1 region several times during the evolution of mammals. In humans, multiple mutations (point mutations and/or deletions) have accumulated on the ancestral TNFRSF, leaving a single short non-functional sequence.     

A genomic reservoir for Tnfrsf genes is developmentally regulated and imprinted in the mouse

Tumor necrosis factor receptor superfamily is composed of at least 26 members in the mouse, three of which exist as a cluster within the imprinted Kcnq1 domain on chromosome 7. Tnfrsf22, 23 and 26 contain typical cystein-rich domains and Tnfrsf22 and 23 can bind ligands but have no signaling capacity. Thus, they are assumed to be decoy receptors. The developmental expression profile of these genes is unknown and knowledge of their imprinting patterns is incomplete and controversial. We found that all three genes are expressed during mouse embryonic development, and that they have a strong maternal bias, indicating that they may be affected by the KvDMR, the Kcnq1 imprinting control region. We found expression of an antisense non-coding RNA, AK155734, in embryos and some neonatal tissues. This RNA overlaps the Tnfrsf22 and possibly the Tnfrsf23 coding regions and is also expressed with a maternal bias. We were interested in exploring the evolutionary origins of the three Tnfrsf genes, because they are absent in the orthologous human Kcnq1 domain. To determine whether the genes were deleted from humans or acquired in the rodent lineage, we performed phylogenetic analyses. Our data suggest that TNFRSF sequences were duplicated and/or degenerated or eliminated from the KCNQ1 region several times during the evolution of mammals. In humans, multiple mutations (point mutations and/or deletions) have accumulated on the ancestral TNFRSF, leaving a single short non-functional sequence.     

Chemical characterization of cyanogen bromide fragments from the beta-chain of human complement factor C3

The isolated beta-chain of human complement factor C3 (C3 beta) was fragmented by cyanogen bromide. Nine fragments were defined by gel filtration and high-pressure liquid chromatography, and characterized with respect to their Mr, amino acid composition and N-terminal amino acid sequence. Approx. 30% of the primary structure of C3 beta was determined. Alignment of the 3 N-terminal fragments allowed determination of 61 of the amino terminal residues of C3 beta. This region demonstrated 40% homology with the sequence in the N-terminal segment of the alpha-chain of the cobra venom factor.     

Molecular dynamics - potential of mean force calculations as a tool for understanding ion permeation and selectivity in narrow channels

Ion channels catalyze the permeation of charged molecules across cell membranes and are essential for many vital physiological functions, including nerve and muscle activity. To understand better the mechanisms underlying ion conduction and valence selectivity of narrow ion channels, we have employed free energy techniques to calculate the potential of mean force (PMF) for ion movement through the prototypical gramicidin A channel. Employing modern all-atom molecular dynamics (MD) force fields with umbrella sampling methods that incorporate one hundred 1-2 ns trajectories, we find that it is possible to achieve semi-quantitative agreement with experimental binding and conductance measurements. We also examine the sensitivity of the MD-PMF results to the choice of MD force field and compare PMFs for potassium, calcium and chloride ions to explore the basis for the valence selectivity of this narrow and uncharged ion channel. A large central barrier is observed for both anions and divalent ions, consistent with lack of experimental conductance. Neither anion or divalent cation is seen to be stabilized inside the channel relative to the bulk electrolyte and each leads to large disruptions to the protein and membrane structure when held deep inside the channel. Weak binding of calcium ions outside the channel corresponds to a free energy well that is too shallow to demonstrate channel blocking. Our findings emphasize the success of the MD-PMF approach and the sensitivity of ion energetics to the choice of biomolecular force field.     

Structural comparison in solution of a native and retro peptide derived from the third helix of Staphylococcus aureus protein A,domain B: retro peptides,a useful tool for the discrimination of helix stabilization factors dependent on the peptide chain orientation

A peptide fragment corresponding to the third helix of Staphylococcus Aureus protein A, domain B, was chosen to study the effect of the main‒chain direction upon secondary structure formation and stability, applying the retro‒enantio concept. For this purpose, two peptides consisting of the native (Ln) and reversed (Lr) sequences were synthesized and their conformational preferences analysed by CD and NMR spectroscopy. A combination of CD and NMR data, such as molar ellipcitity, NOE connectivities, Hα and NH chemical shifts, ³J_αN coupling constants and amide temperature coefficients indicated the presence of nascent helices for both Ln and Lr in water, stabilized upon addition of the fluorinated solvents TFE and HFIP. Helix formation and stabilization appeared to be very similar in both normal and retro peptides, despite the unfavourable charge–macrodipole interactions and bad N-capping in the retro peptide. Thus, these helix stabilization factors are not a secondary structure as determined for this specific peptide. In general, the synthesis and confirmational analysis of peptide pairs with opposite main‒chain directions, normal and retro peptides, could be useful in the determination of secondary structure stabilization factors dependent on the direction. © 1997 European Peptide Society and John Wiley & Sons, Ltd.     

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.     

A high resolution structure of the putative hinge region in M2 channel-lining segments of the nicotinic acetylcholine receptor

Transmembrane peptide helices play key roles in signal transduction across cell membranes, yet little is known about their high-resolution structure or the role membrane composition plays in their association, structure, dynamics and ultimately their performance. Using magic angle spinning (MAS) homonuclear dipolar recoupling experiments, the backbone structure at positions L10, L11, and A12 of the M2 ion channel peptide was determined in two lipid systems. Their measurements are in agreement with M2 forming transmembrane helices, but the torsion angles vary considerably from common α−helical values. These measurements show remarkable agreement with a previous computational model of M2 peptides forming a pore domain in which their helices are kinked near the central leucine, L11 [R. Sankararamakrishnan, C. Adcock, M.S.P. Sansom, The pore domain of the nicotinic acetylcholine receptor: Molecular modeling, pore dimensions, and electrostatics, Biophys. J. 71 (1996) 1659-1671]. The generation of high resolution data for transmembrane helices is of critical importance in refining structures for membrane protein and developing models of helix packing interactions.     

Use of reverse micelles in membrane protein structural biology

Membrane protein structural biology is a rapidly developing field with fundamental importance for elucidating key biological and biophysical processes including signal transduction, intercellular communication, and cellular transport. In addition to the intrinsic interest in this area of research, structural studies of membrane proteins have direct significance on the development of therapeutics that impact human health in diverse and important ways. In this article we demonstrate the potential of investigating the structure of membrane proteins using the reverse micelle forming surfactant dioctyl sulfosuccinate (AOT) in application to the prototypical model ion channel gramicidin A. Reverse micelles are surfactant based nanoparticles which have been employed to investigate fundamental physical properties of biomolecules. The results of this solution NMR based study indicate that the AOT reverse micelle system is capable of refolding and stabilizing relatively high concentrations of the native conformation of gramicidin A. Importantly, pulsed-field-gradient NMR diffusion and NOESY experiments reveal stable gramicidin A homodimer interactions that bridge reverse micelle particles. The spectroscopic benefit of reverse micelle-membrane protein solubilization is also explored, and significant enhancement over commonly used micelle based mimetic systems is demonstrated. These results establish the effectiveness of reverse micelle based studies of membrane proteins, and illustrate that membrane proteins solubilized by reverse micelles are compatible with high resolution solution NMR techniques. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Single channel and monolayer studies of acylated gramicidin A: influence of the length of the alkyl group

Three different gramicidin A analogues bearing acyl chains of various length on the ethanolamine moiety have been studied by investigating their single channel behaviour and their monolayer properties. It is shown that the single channel conductance does not depend on the substitution of the ethanolamine OH group and that the channel lifetime is roughly proportional to the length of the alkyl chain. The monolayer study indicates that acylation of gramicidin A produces compounds which have medium-dependent conformations. These acylated compounds are miscible with lipids, while GA is not, and the surface potential is not modified by the esterification of the alcohol group. Offprint requests to: F. Heitz     

A method for the facile solid-phase synthesis of gramicidin S and its analogs

Summary A simple method is described for the facile synthesis of gramicidin S and six other analogs, using standard solidphase synthetic technology and a single solution-phase cyclization step. The peptides were purified to homogeneity and characterized by plasma desorption time-of-flight mass spectrometry and NMR spectroscopy. Complete ¹H NMR assignments for all seven peptides (in aqueous solution) are presented. Unlike previous approaches, the presented method is simple, automatable, rapid (less than three days), high-yielding, requires no side-chain protection during cyclization, and appears to be generally applicable to the preparation of a variety of related head-to-tail cyclic peptides.Abbreviations Boc t-butyloxycarbonyl - BOP benzotriazoyl N-oxytris(dimethylamino)phosphonium hexafluorophosphate - Bzl benzyl - DCC N,N-dicyclohexylcarbodiimide - DCM dichloromethane - DIEA N,N-diisopropylethylamine - DMF N,N-dimethylformamide - DQF-COSY double-quantum-filtered correlation spectroscopy - DSS 2,2-dimethyl-2-silapentane-5-sulfonate, sodium salt - EDAC 1-ethyl-3-(3-dimethylaminopropyl)carbodiimide - HBTU 2-(1H-benzotriazol-1-yl)-1,1,3,3-tetramethyl uronium hexafluorophosphate - HOBt 1-hydroxybenzotriazole - 4-MeBzl 4-methylbenzyl - NHS N-hydroxysuccinimide - NOESY nuclear Overhauser effect spectroscopy - PAM phenylacetamidomethyl (resin) - RP-HPLC reversed-phase high-performance liquid chromatography - TFA trifluoroacetic acid - TOCSY total correlation spectroscopy - Tos tosyl - Troc 2,2,2-trichloroethylcarbamate.