Modeling negative ion defect migration through the gramicidin A channel

The results of potential of mean force (PMF) calculations for the distinct stages of proton conduction through the gramicidin A channel, including proton migration, reorientation of the water file and negative ion defect migration, are presented. The negative ion defect migration mechanism was hypothesized in experimental studies but was not considered previously in molecular dynamics simulations. The model system consisted of the peptide chains constructed on the base of the structure PDBID:1JNO, the inner file of nine water molecules and external clusters of water molecules placed at both ends of the channel. Potential energy functions were computed with the CHARMM/PM6/TIP3P parameters. The results obtained for proton migration and water file reorientation are basically consistent with those reported previously by Pómès and Roux (Biophys J 82:2304, 2002) within the similar approach. For the newly considered mechanism of negative ion defect migration from the channel center to the end of the water file we obtain the energy 3.8 kcal mol⁻¹ which is not considerably different from the activation energy of water reorientation, 5.4 kcal mol⁻¹. Therefore this mechanism may principally compete for the rate-limiting step in proton conduction in gramicidin.     

Bioactivity of topologically confined gramicidin A dimers

The d-/l-peptide gramicidin A (gA) is well known as a pivotal ion channel model and shows a broad spectrum of bioactivities such as antibiosis, antimalarial activity, as well as hemolysis. We applied inter-chain disulfide bonds to constrain the conformational freedom of gA into parallel and antiparallel dimeric topologies. Albeit the constructs were not found to be monoconformational, CD- and IR-spectroscopic studies suggested that this strategy indeed restricted the conformational space of the d-/l-peptide construct, and that β-helical secondary structures prevail. Correlative testing of gA dimers in antimicrobial, antimalarial, and ion conduction assays suggested that the tail-to-tail antiparallel single stranded β^6.3 helix dominantly mediates the bioactivity of gA. Other conformers are unlikely to contribute to these activities. From these investigations, only weakly ion conducting gA dimers were identified that retained nM antimalarial activity.     

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.     

Photolysis of gramicidin A channels in lipid bilayers

We have tested the hypothesis that peptide tryptophan groups can control the ionic conductance of transmembrane channels. We report here that single gramicidin A channels change conductance state when the peptide tryptophans are flash photolyzed with ultraviolet light. The current flow through planar lipid bilayers containing multiple gramicidin A channels decreases irreversibly when exposed to ultraviolet light. The current-loss action spectrum peaks sharply at the 280 nm absorption maximum of the gramicidin A tryptophans. Gramicidin channel sensitivity to ultraviolet light is found to be about 20-fold higher than that of frog node sodium channels which is even more than expected based on the high tryptophan content of gramicidin. Channels which survive an ultraviolet light exposure exist in a wide variety of different low-conductance forms. The broad distribution of the single channel conductance of these partially photolyzed channels is attributable to the loss of different combinations of the dimer's normal complement of eight tryptophans per channel. Flash photolysis of single channels results in discrete conductance state changes. Partially photolyzed single channels manifest a further conductance cascade when exposed to a second flash of ultraviolet light. Analysis of the photolysis conductance turn-off process indicates that gramicidin A is a multistate electrochemical unit where the peptide tryptophan groups can modulate the flow of ions through the transmembrane channel.     

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.     

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.     

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.     

Three-dimensional Poisson-Nernst-Planck theory studies: influence of membrane electrostatics on gramicidin A channel conductance

A recently introduced real-space lattice methodology for solving the three-dimensional Poisson-Nernst-Planck equations is used to compute current-voltage relations for ion permeation through the gramicidin A ion channel embedded in membranes characterized by surface dipoles and/or surface charge. Comparisons to a variety of experimental results, presented herein, have proven largely successful. Strengths and weaknesses of the method are discussed.     

Single channel properties of d-Leu2-gramicidin A side chain modulation of channel lifetime

The channel forming properties of synthetic gramicidin A and dLeu²-gramicidin A were compared in black lipid membranes. The most probable single channel conductance was identical for both derivatives but in each case a distribution of smaller channel sizes was observed. However, the lifetime of the channel formed by dLeu²-gramicidin A was considerably shorter than for gramicidin A. The dLeu² substitution is considered to interfere with the head to head hydrogen bonding which forms the conducting dimer, thus destabilizing the dimeric structure of the channel and reducing the lifetime. This represents the first demonstration of side-chain modulation of channel lifetime.     

Effects of gramicidin and tryptophan-N-formylated gramicidin on the sodium and potassium content of human erythrocytes

Summary

In order to get a better understanding in the mechanism by which tryptophan-N-formylated gramicidin (NFG) and gramicidin kill the malaria parasite Plasmodium falciparum in vitro, we studied the capacity of these peptides to change the potassium, as well as the sodium, composition of normal human erythrocytes, and their ability to cause cell lysis. It is shown that both peptides are able to induce potassium leakage from, and sodium flux into, erythrocytes in such a manner that it is most likely that they are able to form cation channels in the membrane of these cells. For both peptides, potassium efflux proceeds at a faster rate than sodium influx, but this difference is greater for NFG than for gramicidin. This explains the observation that gramicidin is more lytic than NFG is, even when comparing concentrations that show the same antimalarial activity. The finding that gramicidin is approximately 10 times more active than NFG in causing potassium efflux from normal erythrocytes, as well as in killing the malaria parasite, supports the hypothesis that peptideinduced parasite death is related to their capacity to induce potassium leakage from infected erythrocytes. Finally, the observation that erythrocytes are able to restore their normal ion contents after losing more than 50% of their potassium content by incubation with NFG or gramicidin, suggests that, in vivo, and upon treatment with drug concentrations that cause full inhibition of parasite growth, these cells would not be irreversibly damaged by action of the drugs.     

The structure of the gramicidin A transmembrane channel

Gramicidin A is a linear polypeptide antibiotic that facilitates the diffusion of monovalent cations across lipid bilayer membranes by forming channels. It has been proposed that the conducting channel is a dimer which is in equilibrium with nonconducting monomers in the membrane. To directly test this model in several independent ways, we have prepared and purified a series of gramicidin C derivatives. All of these derivatives are fully active analogs of gramicidin A, and each derivative has a useful chromophore esterified to the phenolic hydroxyl of tyrosine #11. Simultaneous conductance and fluorescence measurements on planar lipid bi-layer membranes containing dansyl gramicidin C yielded four conclusions: (1) A plot of the logarithm of the membrane conductance versus the logarithm of the membrane fluorescence had a slope of 2.0 ± 0.3, over a concentration range for which nearly all the gramicidin was monomeric. Hence, the active channel is a dimer of the nonconducting species. (2) In a membrane in which nearly all of the gramicidin was dimeric, the number of channels was approximately equal to the number of dimers. Thus, most dimers are active channels and so it should be feasible to carry out spectroscopic studies of the conformation of the transmembrane channel. (3) The association constant for dimerization is more than 1,000-fold larger in a glycerolester membrane with 26 Å-hydrocarbon thickness than in a 47 Å-glycerolester membrane. The dimerization constant in a 48 Å-phosphatidyl choline membrane was 200 times larger than in a 47 Å-glycerolester membrane, showing that it depends on the type of lipid as well as on the thickness of the hydrocarbon core. (4) We were readily able to detect 10^?14 mole cm^?2 of dansyl gramicidin C in a bilayer membrane, which corresponds to 60 fluorescent molecules per square μm. The fluorescent techniques described here should be sufficiently sensitive for fluorescence studies of reconstituted gates and receptors in planar bilayer membranes. An alternative method of determining the number of molecules of gramicidin in the channel is to measure the fraction of hybrid channels present in a mixture of 2 chemically different gramicidins. The single-channel conductance of p-phenylazo-benzene-sulfonyl ester gramicidin C (PABS gramicidin C) was found to be 0.68 that of gramicidin A. In membranes containing a mixture of these 2 gramicidins, a hybrid channel was evident in addition to 2 pure channels. The hybrid channel conductance was 0.82 that of gramicidin A. Fluorescence energy transfer from dansyl gramicidin C to diethylamino-phenylazobenzene-sulfonyl ester gramicidin C (DPBS gramicidin C), provided an independent way to measure the fraction of hybrid channels on liposomes. For both techniques the fraction of hybrid channels was found to be 2ad where a² and d² were the fractions of the 2 kinds of pure channels. This result strongly supports a dimer channel and the hybrid data excludes the possibility of a tetramer channel. The study of hybrid species by conductance and fluorescence techniques should be generally useful in elucidating the subunit structure of oligomeric assemblies in membranes. The various models which have been proposed for the conformation of the gramicidin transmembrane channel are briefly discussed.     

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.     

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.     

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.     

Computed energy profiles for ion transport in the gramicidin A channel

Energy profiles have been established for the transport of Na⁺, K⁺ and Cs⁺ through the gramicidin A channel. In Urry’s head-to-headβ _3.3 ^6.3 left-handed helical dimer structure, using a refined methodology for the calculation of intra and intermolecular interactions. The computations show the important role, for the energy profile and the position of the possible binding site, of the flexible ethanolamine chain, whose significance was till now overlooked. The calculations indicate that the barriers at the entrance and at the center of the channel should be in the order Na⁺ XXX K⁺ XXX Cs⁺ but predict also that the energies of the binding sites should be the greatest for Na⁺ and, then, comparable for K⁺ and Cs⁺. The indications concerning the barriers are confirmed by experiment.     

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 states 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.     

Asymmetry of the gramicidin channel in bilayers of asymmetric lipid composition: I. Single channel conductance

Summary Gramicidin-doped asymmetric bilayers made by the Montal-Mueller method exhibited an asymmetric current-voltage relationship. The asymmetric conductance was shown to be the product of two components, a rectifying single-channel conductance and an asymmetric voltage dependence of the reaction which leads to the conducting channel. The single-channel conductance was asymmetric in both asymmetric bilayers made of charged lipids and asymmetric bilayers made only of neutral lipids. The single-channel asymmetry decreased with increasing ion concentration. From the comparison of the singlechannel conductance in symmetric and asymmetric bilayers and the dependence of the asymmetry on the solution ion concentrations, it was concluded that (1) the rate of ion entry into the channel is dependent on the lipid composition of the membrane and is asymmetric in asymmetric bilayers; (2) the entry step is rate determining at low ion concentrations; and (3) at higher ion concentrations the rate-determining step is the translocation across the main barrier in the membrane; and this translocation appears insensitive to lipid asymmetry.     

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.     

Gravity dependency of the gramicidin A channel conductivity

Theoretical investigations involving the membrane-solution interface have revealed that the density of the solution varies appreciably within interfacial layers adjacent to charged membrane surfaces. The hypothesis that gravity interacts with this configuration and modifies transport rates across horizontal and vertical membranes differently was supported by initial experiments with gramicidin A channels in phosphatidylserine (PS) membranes in 0.1 M KCl. Channel conductivity was found to be about 1.6 times higher in horizontal membranes than in vertical membranes. Here we present the results of further experiments with gramicidin A channels (incorporated into charged PS- and uncharged phosphatidylcholine (PC) membranes in KCl- and CsCl-solutions) to demonstrate that the hypothesis is more generally applicable. Again, channel conductivity was found to be higher in horizontal PS membranes by a factor of between 1.20 and 1.75 in 0.1 M CsCl. No difference in channel conductivity was found for uncharged PC membranes in 0.1 M KCl and in 0.1 M CsCl. However, for PC membranes in 0.05 M KCl the channel conductivity was significantly higher in horizontal membranes by a factor of between 1.07 and 1.14. These results are consistent with the results of our model calculations of layer density and extension, which showed that the layer formation is enhanced by increasing membrane surface charge and decreasing electrolyte ion concentration. The mechanism of gravity interaction with membrane transport processes via interface reactions might be utilized by biological systems for orientational behaviour in the gravity field, which has been observed even for cellular systems. Received: 16 October 1995 / Accepted: 23 April 1996     

The dimeric nature of the gramicidin A transmembrane channel: Conductance and fluorescence energy transfer studies of hybrid channels

Gramicidin A, a linear peptide antibiotic, makes membranes permeable to alkali cations and hydrogen ions by forming transmembrane channels. We report here conductance and fluorescence energy transfer studies of channels containing two kinds of gramicidin. These studies of hybrid channels were designed to determine the number of molecules in a channel. The gramicidins studied were gramicidin A, dansyl gramicidin C, the p-phenylazobenzene sulfonyl derivative of gramicidin C (PABS⁴ gramicidin C), and the 4-(diethylamino)-phenylazobenzene-4-sulfonyl chloride derivative of gramicidin C (DPBS gramicidin C). The dansyl, PABS and DPBS groups were linked to the hydroxyl group of tyrosine 11 in gramicidin C. The single-channel conductance of PABS gramicidin C in planar bilayer membranes is 0.68 that of gramicidin A. Membranes containing both PABS gramicidin C and gramicidin A exhibit three kinds of channels: a pure gramicidin A, a pure PABS gramicidin C channel, and a hybrid channel with an intermediate conductance (0.82 that of gramicidin A). The dependence of the frequencies of these three kinds of channels on the mole fractions of gramicidin A and PABS gramicidin C in the membrane-forming solution fits a dimer model. Fluorescence energy transfer was used as a complementary means of ascertaining the frequency of hybrid channels. Dansyl gramicidin C was the fluorescent energy donor and DPBS gramicidin C was the energy acceptor. The efficiency of energy transfer between these chromophores in hybrid channels in liposomes was 75%. The relative quantum yield of the dansyl fluorescence was measured as a function of the mole fraction of DPBS gramicidin C. These fluorescence studies, like the single-channel conductance measurements, showed that there are two molecules of gramicidin in a channel. The study of hybrid species by conductance and fluorescence techniques should be generally useful in elucidating the subunit structure of oligomeric assemblies in membranes.