Water permeation through gramicidin A: desformylation and the double helix: a molecular dynamics study

Multinanosecond molecular dynamics simulations of gramicidin A embedded in a dimyristoylphosphatidylcholine bilayer show a remarkable structural stability for both experimentally determined conformations: the head-to-head helical dimer and the double helix. Water permeability was found to be much higher in the double helical conformation, which is explained by lower hydrogen bond-mediated enthalpic barriers at the channel entrance and its larger pore size. Free-energy perturbation calculations show that the double helical structure is stabilized by the positive charges at the N termini introduced by the desformylation, whereas the helical dimer is destabilized. Together with the recent experimental observation that desformyl gramicidin conducts water hundredfold better than gramicidin, this suggests that desformyl gramicidin A predominantly occurs in the double helical conformation.     

Crystals of a trypsin-modified alkaline phosphatase. Preliminary crystallographic characterization

Trypsin-modified alkaline phosphatase from Escherichia coli has been crystallized in a form distinct from the two known crystal forms of the native enzyme. The large well diffracting crystals belong to the orthorhombic space group P2(1)2(1)2(1), possess unit cell dimensions a = 56.0 A, b = 136.0 A, c = 283.9 A with 2 dimers per asymmetric unit, and are suitable for high resolution x-ray crystallographic studies. The observed structural and functional differences between the native and modified molecules are a result of peptide bond cleavage at Arg10-Ala11 with loss of the NH2-terminal decapeptide in both subunits of the dimer.     

Purification and characterization of the trifunctional beta-subunit of anthranilate synthase from Neurospora crassa

The trifunctional beta-subunit of anthranilate synthase complex of Neurospora crassa has been purified from a mutant which produces no detectable alpha-subunit. The isolated beta-subunit appeared to be a highly asymmetric dimer with a s20,w of 7.35 and an apparent molecular weight of 200,000 as determined by gel filtration on Sephacryl S-300 compared with a monomer molecular weight of approximately 84,000 Da as determined by sodium dodecyl sulfate-gel electrophoresis. The purified subunit was cleaved by elastase, trypsin, or chymotrypsin into fragments which retained the three enzyme activities. After elastase digestion, two active fragments were separated by gel filtration and ion exchange chromatography. A 30,000-Da fragment, which behaved as a monomer on gel filtration, interacted with free alpha-subunit to produce glutamine-dependent anthranilate synthase activity. A second 56,000-Da fragment, which behaved as an asymmetric dimer (apparent molecular weight 140,000) on gel filtration, retained both N-(5'-phosphoribosyl)anthranilate isomerase and indole-3-glycerol phosphate synthase activity. The failure to detect an NH2-terminal amino acid residue on either the intact beta-subunit or the 30,000-Da complementing fragment, while the 56,000-Da fragment possessed an NH2-terminal histidine residue, indicated that the complementing fragment was derived from the NH2-terminal sequence of the beta-subunit.     

Electric field effects on membranes: gramicidin A as a test ground

Electric fields due to transmembrane potential differences or ionic gradients across the membrane are presumably crucial for many reactions across membranes or close to membranes like signal transduction, control of ion channels or the generation of neural impulses. Molecular dynamics simulations have been used to study the influence of external electric fields on a mixed gramicidin/phospholipid bilayer system. At high field strengths, formation of membrane electropores occurred both close and distal to the gramicidin. Gramicidin was found to stabilize the membrane adjacent to the protein but also at larger distances of up to 2-3 nm. As a result, membrane pore formation was found to be significantly suppressed for the mixed gramicidin/DMPC system. Moderate field strengths only weakly affected the structure and dynamics of the gramicidin. Spontaneous potassium passage events in external electric fields were observed for both the head-to-head helical conformation as well as for the double helical conformation of gramicidin A. The double-helical conformation was found to facilitate ion passage compared to the head-to-head helical dimer.     

Energetics of Shaker K channels block by inactivation peptides

A synthetic peptide of the NH2-terminal inactivation domain of the ShB channel blocks Shaker channels which have an NH2-terminal deletion and mimics many of the characteristics of the intramolecular inactivation reaction. To investigate the role of electrostatic interactions in both peptide block and the inactivation process we measured the kinetics of block of macroscopic currents recorded from the intact ShB channel, and from ShB delta 6-46 channels in the presence of peptides, at different ionic strengths. The rate of inactivation and the association rate constants (k(on)) for the ShB peptides decreased with increasing ionic strength. k(on) for a more positively charged peptide was more steeply dependent on ionic strength consistent with a simple electrostatic mechanism of enhanced diffusion. This suggests that a rate limiting step in the inactivation process is the diffusion of the NH2-terminal domain towards the pore. The dissociation rates (k(off)) were insensitive to ionic strength. The temperature dependence of k(on) for the ShB peptide was very high, (Q10 = 5.0 +/- 0.58), whereas k(off) was relatively temperature insensitive (Q10 approximately 1.1). The results suggest that at higher temperatures the proportion of time either the peptide or channel spends in the correct conformation for binding is increased. There were two components to the time course of recovery from block by the ShB peptide, indicating two distinct blocked states, one of which has similar kinetics and dependence on external K+ concentration as the inactivated state of ShB. The other is voltage- dependent and at -120 mV is very unstable. Increasing the net charge on the peptide did not increase sensitivity to knock-off by external K+. We propose that the free peptide, having fewer constraints than the tethered NH2-terminal domain binds to a similar site on the channel in at least two different conformations.     

Peptide backbone chemistry and membrane channel function: effects of a single amide-to-ester replacement on gramicidin channel structure and function

To examine the structural and functional importance of backbone amide groups in ion channels for subunit folding, hydrogen bonding, ion solvation, and ion permeation, we replaced the peptide bond between Val(1) and Gly(2) in gramicidin A by an ester bond. The substitution is at the junction between the two channel subunits, where it removes an intramolecular hydrogen bond between the NH of Gly(2) and the C==O of Val(7) and perturbs an intermolecular hydrogen bond between the C==O of Val(1) in one subunit and the NH of Ala(5) in the other subunit. The substitution thus perturbs not only subunit folding but also dimer assembly, in addition to any effects on ion permeation. This backbone modification has large effects on channel function: It alters channel stability, as monitored by the channel forming ability and channel lifetime, and ion permeability, as monitored by changes in single-channel conductance and cation permeability ratios. In fact, the homodimeric channels, with two ester-containing subunits, have lifetimes so short that it becomes impossible to characterize them in any detail. The peptide --> ester substitution, however, does not affect the basic subunit fold because heterodimeric channels can form between a subunit with an ester bond and a native subunit. These heterodimeric channels, with only a single ester bond, are more easily characterized; the lone ester reduces the single-channel conductance about 4-fold and the lifetime about 200-fold as compared to the native homodimeric channels. The altered channel function results from a perturbation/disruption of the hydrogen bond network that stabilizes the backbone, as well as the membrane-spanning dimer, and that forms the lining of the ion-conducting pore. Molecular dynamics simulations show the expected destabilization of the modified heterodimeric or homodimeric channels, but the changes in backbone structure and dynamics are remarkably small. The ester bond is somewhat unstable, which precluded further structural characterization. The lability also led to a hydrolysis product that terminates with an alcohol and lacks formyl-Val. Symmetric channels formed by the hydrolyzed product again have short lifetimes, but the channels are distinctly different from those formed by the ester gramicidin A. Furthermore, well-behaved asymmetric channels form between the hydrolysis product and reference subunits that have either an L- or a D-residue at the formyl-NH-terminus.     

Binding of alkaline cations to the double-helical form of gramicidin.

Gramicidin is a polypeptide antibiotic that forms monovalent cation-specific channels in membrane environments. In organic solvents and in lipids containing unsaturated fatty acid chains, it forms a double-helical "pore" structure, in which two monomers are intertwined. This form of gramicidin can bind two cations inside its lumen, and the crystal structures of both an ion complex and an ion-free form have been determined. In this study, we have used circular dichroism (CD) spectroscopy to examine the binding mechanism and the binding constants (K1 and K2) of cations to gramicidin in the double helical form in methanol solution. The dramatic change in optical rotation in the far-ultraviolet CD spectrum of gramicidin provides a useful tool for monitoring the binding. The binding mechanism appears to involve a large conformation change associated with the binding of ions to the first of the two sites. The calculated values for the K1 binding constants for alkaline cations are considerably smaller than the K2 binding constants. The order of binding affinity for alkaline cations is similar to that for the helical dimer "channel" form of gramicidin, i.e., Cs+ approximately Rb+ > > K+ > Li+, but in comparison to the helical dimer form, the binding to double-helical dimers is dominated by a cation size-dependent conformational change in the gramicidin structure.     

Modeling the ion channel structure of cecropin.

Atomic-scale computer models were developed for how cecropin peptides may assemble in membranes to form two types of ion channels. The models are based on experimental data and physiochemical principles. Initially, cecropin peptides, in a helix-bend-helix motif, were arranged as antiparallel dimers to position conserved residues of adjacent monomers in contact. The dimers were postulated to bind to the membrane with the NH2-terminal helices sunken into the head-group layer and the COOH-terminal helices spanning the hydrophobic core. This causes a thinning of the top lipid layer of the membrane. A collection of the membrane bound dimers were then used to form the type I channel structure, with the pore formed by the transmembrane COOH-terminal helices. Type I channels were then assembled into a hexagonal lattice to explain the large number of peptides that bind to the bacterium. A concerted conformational change of a type I channel leads to the larger type II channel, in which the pore is formed by the NH2-terminal helices. By having the dimers move together, the NH2-terminal helices are inserted into the hydrophobic core without having to desolvate the charged residues. It is also shown how this could bring lipid head-groups into the pore lining.     

A lattice relaxation algorithm for three-dimensional Poisson-Nernst-Planck theory with application to ion transport through the gramicidin A channel.

A lattice relaxation algorithm is developed to solve the Poisson-Nernst-Planck (PNP) equations for ion transport through arbitrary three-dimensional volumes. Calculations of systems characterized by simple parallel plate and cylindrical pore geometries are presented in order to calibrate the accuracy of the method. A study of ion transport through gramicidin A dimer is carried out within this PNP framework. Good agreement with experimental measurements is obtained. Strengths and weaknesses of the PNP approach are discussed.     

Distinction between dipolar and inductive effects in modulating the conductance of gramicidin channels.

The ion permeability of transmembrane channels formed by the linear gramicidins is altered by amino acid sequence substitutions. We have previously shown that the polarity of the side chain at position one is important in modulating a channel's conductance and ion selectivity [Russel et al. (1986) Biophys. J. 49, 673-686]. Changes in polarity could alter ion permeability by (through-space) ion-dipole interactions or by (through-bond) inductive electron shifts. We have addressed this question by investigating the permeability characteristics of channels formed by gramicidins where the NH2-terminal amino acid is either phenylalanine or one of a series of substituted phenylalanines: p-hydroxy-, p-methoxy-, o-fluoro-, m-fluoro-, or p-fluorophenylalanine. The electron-donating or -withdrawing nature, as quantified by the Hammett constant, ranges from -0.37 to +0.34 for these side chains. Channels formed by these gramicidins show a more than 2.5-fold variation in their Na+ conductance, but the conductance variations do not rank in the order of the Hammett constants of the side chains. Inductive effects cannot therefore be of primary importance in the modulation of the gramicidin single-channel conductance by these side chains. The results support previous suggestions that electrostatic interactions between side chain dipoles and permeating ions can modify the energy profile for ion movement through the gramicidin channel and thus alter the conductance.     

Multiple solution conformations and internal rotations of the decapeptide gramicidin S.

The conformations of every C alpha H-C beta H2 moiety of the peptide gramicidin S are reported. Internal rotation occurs, but distinct preferences for one side chain rotamer, greater than 80%, are found for the D-phenylalanine and ornithine residues. Leucine and valine exhibit more extensive averaging while proline is shown to be at least 90% in the Ramachandran B conformation. The data are consistent with the coexistence of many tertiary conformations of gramicidin S; the statistical weights of the twelve major tertiary conformations consistent with the rotamer populations are reported. The relative statistical weights of the tertiary conformers depend upon temperature and solvent. A comparison of the conclusions from this publication and conformations derived by energy minimization procedures is made. Partial agreement was found, but the calculations have not yet predicted the wealth of coexisting tertiary conformations nor accounted for the subtle effects of solvent. It is proposed that a more complete picture of the conformational dynamics of gramicidin S and other peptides will result from calculations which use as a basis the extensive data reported here.     

Modulation of troponin T molecular conformation and flexibility by metal ion binding to the NH2-terminal variable region

Troponin T (TnT) plays an allosteric signal transduction role in the actin thin-filament-based Ca(2+)-regulation of striated muscle contraction. Developmentally regulated alternative RNA splicing produces TnT isoforms differing in their NH(2)-terminal structure. Physical property variations of the NH(2)-terminal hypervariable region of TnT may have a role in tuning the Ca(2+)-sensitivity and overall cooperativity of the muscle. We have previously demonstrated that metal ion or monoclonal antibody binding to the NH(2)-terminal region can modulate the epitopic conformation and troponin I and tropomyosin binding affinity of TnT. To further establish the molecular basis of this conformational and functional modulation, we have characterized the NH(2)-terminal variable region-originated secondary conformational effect in TnT using fluorescence spectral analysis. The chicken fast skeletal muscle TnT isoform, TnT8e16, containing a cluster of transition-metal ion binding sites (Tx) in the NH(2)-terminal variable region was used in this study. TnT8e16 was titrated for Cu(II) binding-induced changes in fluorescence intensity and anisotropy of the COOH-domain Trp residues (W234, W236, and W285), which demonstrated considerable environmental sensitivity in TnT denaturation studies. Nonlinear Stern-Volmer plots of Trp quenching indicated a metal ion binding-induced conformational change in TnT. Fluorescence anisotropy changes upon metal ion binding indicated a decrease in the mobility of the Trp residues and an increase in the flexibility of fluorescein-labeled Cys263 in the COOH domain. These data support a model that the alternatively spliced NH(2)-terminal variable region of TnT modulates conformation and flexibility of other domains of the protein.     

Gramicidin Channels Are Internally Gated

Gramicidin channels are archetypal molecular subjects for solid-state NMR studies and investigations of single-channel or cation conductance. Until now, the transitions between on and off conductance states have been thought, based on multichannel studies, to represent monomer ↔ dimer reactions. Here we use a single-molecule deposition method (vesicle fusion to a planar bilayer) to show that gramicidin dimer channels do not normally dissociate when conductance terminates. Furthermore, the observation of two ¹³C peaks in solid-state NMR indicates very stable dichotomous conformations for both the first and second peptide bonds in the monomers, and a two-dimensional chemical exchange spectrum with a 12-s mixing time demonstrates that the Val₁ carbonyl conformations exchange slowly, with lifetimes of several seconds. It is proposed that gramicidin channels are gated by small conformational changes in the channel near the permeation pathway. These studies demonstrate how regulation of conformations governing closed ↔ open transitions may be achieved and studied at the molecular level.     

Proton conduction in gramicidin A and in its dioxolane-linked dimer in different lipid bilayers.

Gramicidin A (gA) molecules were covalently linked with a dioxolane ring. Dioxolane-linked gA dimers formed ion channels, selective for monovalent cations, in planar lipid bilayers. The main goal of this study was to compare the functional single ion channel properties of natural gA and its covalently linked dimer in two different lipid bilayers and HCl concentrations (10-8000 mM). Two ion channels with different gating and conductance properties were identified in bilayers from the product of dimerization reaction. The most commonly observed and most stable gramicidin A dimer is the main object of this study. This gramicidin dimer remained in the open state most of the time, with brief closing flickers (tau(closed) approximately 30 micros). The frequency of closing flickers increased with transmembrane potential, making the mean open time moderately voltage dependent (tau(open) changed approximately 1.43-fold/100 mV). Such gating behavior is markedly different from what is seen in natural gA channels. In PEPC (phosphatidylethanolamine-phosphatidylcholine) bilayers, single-channel current-voltage relationships had an ohmic behavior at low voltages, and a marked sublinearity at relatively higher voltages. This behavior contrasts with what was previously described in GMO (glycerylmonooleate) bilayers. In PEPC bilayers, the linear conductance of single-channel proton currents at different proton concentrations was essentially the same for both natural and gA dimers. g(max) and K(D), obtained from fitting experimental points to a Langmuir adsorption isotherm, were approximately 1500 pS and 300 mM, respectively, for both the natural gA and its dimer. In GMO bilayers, however, proton affinities of gA and the dioxolane-dimer were significantly lower (K(D) of approximately 1 and 1.5 M, respectively), and the g(max) higher (approximately 1750 and 2150 pS, respectively) than in PEPC bilayers. Furthermore, the relationship between single-channel conductance and proton concentration was linear at low bulk concentrations of H+ (0.01-2 M) and saturated at concentrations of more than 3 M. It is concluded that 1) The mobility of protons in gramicidin A channels in different lipid bilayers is remarkably similar to proton mobilities in aqueous solutions. In particular, at high concentrations of HCl, proton mobilities in gramicidin A channel and in solution differ by only 25%. 2) Differences between proton conductances in gramicidin A channels in GMO and PEPC cannot be explained by surface charge effects on PEPC membranes. It is proposed that protonated phospholipids adjacent to the mouth of the pore act as an additional source of protons for conduction through gA channels in relation to GMO bilayers. 3) Some experimental results cannot be reconciled with simple alterations in access resistance to proton flow in gA channels. Said differences could be explained if the structure and/or dynamics of water molecules inside gramicidin A channels is modulated by the lipid environment and by modifications in the structure of gA channels. 4) The dioxolane ring is probably responsible for the closing flickers seen in the dimer channel. However, other factors can also influence closing flickers.     

The preference of tryptophan for membrane interfaces: insights from N-methylation of tryptophans in gramicidin channels

To better understand the structural and functional roles of tryptophan at the membrane/water interface in membrane proteins, we examined the structural and functional consequences of Trp --> 1-methyl-tryptophan substitutions in membrane-spanning gramicidin A channels. Gramicidin A channels are miniproteins that are anchored to the interface by four Trps near the C terminus of each subunit in a membrane-spanning dimer. We masked the hydrogen bonding ability of individual or multiple Trps by 1-methylation of the indole ring and examined the structural and functional changes using circular dichroism spectroscopy, size exclusion chromatography, solid state (2)H NMR spectroscopy, and single channel analysis. N-Methylation causes distinct changes in the subunit conformational preference, channel-forming propensity, single channel conductance and lifetime, and average indole ring orientations within the membrane-spanning channels. The extent of the local ring dynamic wobble does not increase, and may decrease slightly, when the indole NH is replaced by the non-hydrogen-bonding and more bulky and hydrophobic N-CH(3) group. The changes in conformational preference, which are associated with a shift in the distribution of the aromatic residues across the bilayer, are similar to those observed previously with Trp --> Phe substitutions. We conclude that indole N-H hydrogen bonding is of major importance for the folding of gramicidin channels. The changes in ion permeability, however, are quite different for Trp --> Phe and Trp --> 1-methyl-tryptophan substitutions, indicating that the indole dipole moment and perhaps also ring size and are important for ion permeation through these channels.     

Stability of an ion channel in lipid bilayers: implicit solvent model calculations with gramicidin

Gramicidin is a helical peptide, 15 residues in length, which dimerizes to form ion-conducting channels in lipid bilayers. Here we report calculations of its free energy of transfer from the aqueous phase into bilayers of different widths. The electrostatic and nonpolar contributions to the desolvation free energy were calculated using implicit solvent models, in which gramicidin was described in atomic detail and the hydrocarbon region of the membrane was described as a slab of hydrophobic medium embedded in water. The free energy penalties from the lipid perturbation and membrane deformation effects, and the entropy loss associated with gramicidin immobilization in the bilayer, were estimated from a statistical thermodynamic model of the bilayer. The calculations were carried out using two classes of experimentally observed conformations: a head-to-head dimer of two single-stranded (SS) beta-helices and a double-stranded (DS) intertwined double helix. The calculations showed that gramicidin is likely to partition into the bilayer in all of these conformations. However, the SS conformation was found to be significantly more stable than the DS in the bilayer, in agreement with most of the experimental data. We tested numerous transmembrane and surface orientations of gramicidin in bilayers of various widths. Our calculations indicate that the most favorable orientation is transmembrane, which is indeed to be expected from a channel-forming peptide. The calculations demonstrate that gramicidin insertion into the membrane is likely to involve a significant deformation of the bilayer to match the hydrophobic width of the peptide (22 A), again in good agreement with experimental data. Interestingly, deformation of the bilayer was induced by all of the gramicidin conformations.     

Double-headed protease inhibitors from black-eyed peas. III. Subunit interactions of the native and half-site chemically modified proteins.

The chemical modification of two new double-headed-protease inhibitors from black-eyed peas, a trypsin-chymotrypsin inhibitor (BEPCI) and a trypsin inhibitor (BEPTI) with dansyl chloride was investigated under various conditions. The NH2-terminal serine of both BEPCI and BEPTI, the 4 lysyl residues of BEPCI, and 4 of the 5 lysyl residues of BEPTI, could not be dansylated in the absence of urea. The single tyrosine per subunit of BEPCI and BEPTI was unreactive even in the presence of urea but could be labeled with half-site reactivity by the Celite method. Lysine, NH2-terminal serine, and tyrosine were reactive in fully reduced, carbamidomethylated BEPCI and BEPTI. Gel filtration was used to study the subunit interactions of BEPCI and BEPTI. At pH 8 or pH 3.0 there is a complex set of multiple equilibria with widely differing rates of attainment. We have found evidence for a rapid dimer-tetramer equilibrium, a distinct moderate rate dimer-tetramer equilibrium, a very slow monomer-dimer equilibrium, and postulate slow isomerization of the two forms of dimer and the two forms of tetramer. The monomer-dimer equilibrium is quite unusual in that the dimer is stabilized by chaotropic ions and even slightly by guanidine HC1. In contrast to the complex pattern seen in native BEPCI, the half-site, dansylated BEPCI exists at similar concentration exclusively as a tetramer at neutral pH.     

Determination of the orientation of an integral membrane protein and sites of glycosylation by oligonucleotide-directed mutagenesis: influenza B virus NB glycoprotein lacks a cleavable signal sequence and has an extracellular NH2-terminal region.

The membrane orientation of the NB protein of influenza B virus, a small (Mr, approximately 18,000) glycoprotein with a single internal hydrophobic domain, was investigated by biochemical and genetic means. Cell fractionation and protein solubility studies indicate NB is an integral membrane protein, and NB has been shown to be a dimer under nonreducing conditions. Treatment of infected-cell surfaces with proteinase K and endoglycosidase F and immunoprecipitation with a site-specific antibody suggests that the 18-amino-acid NH2-terminal region of NB is exposed at the cell surface. Oligonucleotide-directed mutagenesis to eliminate each of the four potential sites of N-linked glycosylation and expression of the mutant NB proteins in eucaryotic cells suggest that the two sites adjacent to the NH2 terminus are glycosylated. This provides further evidence that NB, which lacks a cleavable NH2-terminal signal sequence, has an exposed NH2 terminus at the cell surface.