PubMed HealthSearch

Biomedical subjects

M S Sansom

Publications and source records attributed to M S Sansom.

At least 19 recordsLinked to original sources

Simulation of voltage-dependent interactions of alpha-helical peptides with lipid bilayers.

Pore formation in lipid bilayers by channel-forming peptides and toxins is thought to follow voltage-dependent insertion of amphipathic alpha-helices into lipid bilayers. We have developed an approximate potential for use within the CHARMm molecular mechanics program which enables one to simulate voltage-dependent interaction of such helices with a lipid bilayer. Two classes of helical peptides which interact with lipid bilayers have been studied: (a) delta-toxin, a 26 residue channel-forming peptide from Staphylococcus aureus; and (b) synthetic peptides corresponding to the alpha 5 and alpha 7 helices of the pore-forming domain of Bacillus thuringiensis CryIIIA delta-endotoxin. Analysis of delta-toxin molecular dynamics (MD) simulations suggested that the presence of a transbilayer voltage stabilized the inserted location of delta-toxin helices, but did not cause insertion per se. A series of simulations for the alpha 5 and alpha 7 peptides revealed dynamic switching of the alpha 5 helix between a membrane-associated and a membrane-inserted state in response to a transbilayer voltage. In contrast the alpha 7 helix did not exhibit such switching but instead retained a membrane associated state. These results are in agreement with recent experimental studies of the interactions of synthetic alpha 5 and alpha 7 peptides with lipid bilayers.

Amino Acid Sequence

Structure and orientation of the mammalian antibacterial peptide cecropin P1 within phospholipid membranes.

Cecropins are positively charged antibacterial peptides that act by permeating the membrane of susceptible bacteria. To gain insight into the mechanism of membrane permeation, the secondary structure and the orientation within phospholipid membranes of the mammalian cecropin P1 (CecP) was studied using attenuated total reflectance Fourier-transform infrared (ATR-FTIR) spectroscopy and molecular dynamics simulations. The shape and frequency of the amide I and II absorption peaks of CecP within acidic PE/PG multibilayers (phosphatidylethanolamine/phosphatidylglycerol) in a 7:3 (w/w) ratio (a phospholipid composition similar to that of many bacterial membranes), indicated that the peptide is predominantly alpha-helical. Polarized ATR-FTIR spectroscopy was used to determine the orientation of the peptide relative to the bilayer normal of phospholipid multibilayers. The ATR dichroic ratio of the amide I band of CecP peptide reconstituted into oriented PE/PG phospholipid membranes indicated that the peptide is preferentially oriented nearly parallel to the surface of the lipid membranes. A similar secondary structure and orientation were found when zwitterionic phosphatidylcholine phospholipids were used. The incorporation of CecP did not significantly change the order parameters of the acyl chains of the multibilayer, further suggesting that CecP does not penetrate the hydrocarbon core of the membranes. Molecular dynamics simulations were used to gain insight into possible effects of transmembrane potential on the orientation of CecP relative to the membrane. The simulations appear to confirm that CecP adopts an orientation parallel to the membrane surface and does not insert into the bilayer in response to a cis positive transmembrane voltage difference. Taken together, the results further support a "carpet-like" mechanism, rather than the formation of transmembrane pores, as the mode of action of CecP. According to this model, formation of a layer of peptide monomers on the membrane surface destablizes the phospholipid packing of the membrane leading to its eventual disintegration.

Amino Acid Sequence

Engineering stabilized ion channels: covalent dimers of alamethicin.

The peptide alamethicin forms channels with a variety of conductance states. Selective stabilization of a particular state should simplify the task of understanding conductance in terms of channel structure. We synthesized two different covalent dimers of alamethicin in which peptides were linked at their C-terminal ends by flexible tethers. Both dimeric peptides formed channels with conductances that matched those of alamethicin channels. Particular conductance states were selectively stabilized, however, with lifetimes up to 170-fold longer than the same states observed with monomers. In addition, tethering appeared to limit the size of the structures formed so that, even at higher peptide concentrations, a single predominant conductance state was obtained. We suggest this state corresponds to a channel made from six alamethicin molecules (three dimers).

Alamethicin

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

The ion channel of the nicotinic acetylcholine receptor is a water-filled pore formed by five M2 helix segments, one from each subunit. Molecular dynamics simulations on bundles of five M2 alpha 7 helices surrounding a central column of water and with caps of water molecules at either end of the pore have been used to explore the effects of intrapore water on helix packing. Interactions of water molecules with the N-terminal polar sidechains lead to a conformational transition from right- to left-handed supercoils during these stimulations. These studies reveal that the pore formed by the bundle of M2 helices is flexible. A structural role is proposed for water molecules in determining the geometry of bundles of isolated pore-forming helices.

Amino Acid Sequence

Modelling packing interactions in parallel helix bundles: pentameric bundles of nicotinic receptor M2 helices.

The transbilayer pore of the nicotinic acetylcholine receptor (nAChR) is formed by a pentameric bundle of M2 helices. Models of pentameric bundles of M2 helices have been generated using simulated annealing via restrained molecular dynamics. The influence of: (a) the initial C alpha template; and (b) screening of sidechain electrostatic interactions on the geometry of the resultant M2 helix bundles is explored. Parallel M2 helices, in the absence of sidechain electrostatic interactions, pack in accordance with simple ridges-in-grooves considerations. This results in a helix crossing angle of ca. +12 degrees, corresponding to a left-handed coiled coil structure for the bundle as a whole. Tilting of M2 helices away from the central pore axis at their C-termini and/or inclusion of sidechain electrostatic interactions may perturb such ridges-in-grooves packing. In the most extreme cases right-handed coiled coils are formed. An interplay between inter-helix H-bonding and helix bundle geometry is revealed. The effects of changes in electrostatic screening on the dimensions of the pore mouth are described and the significance of these changes in the context of models for the nAChR pore domain is discussed.

Amino Acid Sequence

Ion channel formation by synthetic analogues of staphylococcal delta-toxin.

Ion channel formation by three analogues of staphylococcal delta-toxin, an amphipathic and alpha-helical channel-forming peptide, has been evaluated by measurement of ionic currents across planar lipid bilayers. Replacement of beta-branched, hydrophobic residues by leucine and movement of a tryptophan residue from the hydrophilic to the hydrophobic face of the helix does not significantly alter ion channel activity. Removal of the N-terminal blocking group combined with the substitution of glycine-10 by leucine changes the single channel properties of delta-toxin, without altering macroscopic conductance/voltage behaviour. Truncation of the N-terminus by three residues results in complete loss of channel-forming activity. These changes in channel-forming properties upon altering the peptide sequence do not mirror changes in haemolytic activity. The results lend support to the proposal that channel formation and haemolysis are distinct events. Channel properties are discussed in the context of a model in which the pore is formed by a bundle of approximately parallel transbilayer helices.

Amino Acid Sequence

Ion-channel gating. Twist to open.

Three-dimensional images of the open state of the nicotinic acetylcholine receptor have been obtained at 9 A resolution. Comparison with the closed state reveals the structural basis of channel gating.

Acetylcholine

Packing interactions of Aib-containing helices: molecular modeling of parallel dimers of simple hydrophobic helices and of alamethicin.

alpha-Aminoisobutyric acid (Aib) is a helicogenic alpha, alpha-dimethyl amino acid found in channel-forming peptaibols such as alamethicin. Possible effects of Aib on helix-helix packing are analyzed. Simulated annealing via restrained molecular dynamics is used to generate ensembles of approximately parallel helix dimers. Analysis of variations in geometrical and energetic parameters within ensembles defines how tightly a pair of helices interact. Simple hydrophobic helix dimers are compared: Ala20, Leu20, Aib20, and P20, the latter a simple channel-forming peptide [G. Menestrina, K.P. Voges, G. Jung, and G. Boheim (1986) Journal of Membrane Biology, Vol. 93, pp. 111-132]. Ala20 and Leu20 dimers exhibit well-defined ridges-in-grooves packing with helix crossing angles (omega) of the order of +20 degrees. Aib20 alpha-helix dimers are much more loosely packed, as evidenced by a wide range of omega values and small helix-helix interaction energies. However, when in a 3(10) conformation Aib20 helices pack in three well-defined parallel modes, with omega ca. -15 degrees, +5 degrees, and 10 degrees. Comparison of helix-helix interaction energies suggests that dimerization may favor the 3(10) conformation. P20, with 8 Aib residues, also shows looser packing of alpha-helices. The results of these studies of hydrophobic helix dimers are analyzed in the context of the ridges-in-grooves packing model. Simulations are extended to dimers of alamethicin, and of an alamethicin derivative in which all Aib residues are replaced by Leu. This substitution has little effect on helix-helix packing. Rather, such interactions appear to be sensitive to interactions between polar side chains. Overall, the results suggest that Aib may modulate the packing of simple hydrophobic helices, in favor of looser interactions. For more complex amphipathic helices, interactions between polar side chains may be more critical.

Alamethicin

Structural features of isolated M2 helices of nicotinic receptors. Simulated annealing via molecular dynamics studies.

The nicotinic acetylcholine receptor is an integral membrane protein and a ligand-gated cation channel. It has stoichiometry alpha 2 beta gamma delta, the subunits arranged symmetrically around an approximate five-fold axis. Five M2 helices, one from each subunit, form a parallel helix bundle surrounding a central pore. Simulated annealing via restrained molecular dynamics (SA/MD) has been employed to generate ensembles of isolated M2 transmembrane helices. Four ensembles of two different M2 helix sequences, M2 delta and M2 gamma, have been generated by SA/MD. The ensembles differed in their treatment of electrostatic interactions. Analysis of the simulated structures showed that intra-helical H-bonds were more strongly conserved in the C-terminal (and more hydrophobic) segment of M2 helices. Conformations of polar sidechains have been analyzed, placing particular emphasis on EK (and QK) pairs at the N-termini of M2 delta (and M2 gamma) helices. Conformations of EK sidechain pairs were obtained for the high resolution structures in the protein database in order to guide our analysis of simulated structures. Serine and threonine sidechain conformations in the M2 models also have been determined. Implications of studies of isolated M2 helices for models of the intact pore region of the nicotinic receptor are discussed.

Amino Acid Sequence

Transbilayer pores formed by beta-barrels: molecular modeling of pore structures and properties.

Transmembrane beta-barrels, first observed in bacterial porins, are possible models for a number of membrane channels. Restrained molecular dynamics simulations based on idealized C alpha beta templates have been used to generate models of such beta-barrels. Model beta-barrels have been analyzed in terms of their conformational, energetic, and pore properties. Model beta-barrels formed by N = 4, 8, 12 and 16 anti-parallel Ala10 strands have been developed. For each N, beta-barrels with shear numbers S = N to 2N have been modeled. In all beta-barrel models the constituent beta-strands adopt a pronounced right-handed twist. Interstrand interactions are of approximately equal stability for all models with N > or = 8, whereas such interactions are weaker for the N = 4 beta-barrels. In N = 4 beta-barrels the pore is too narrow (minimum radius approximately 0.6 A) to allow ion permeation. For N > or = 8, the pore radius depends on both N and S; for a given value of N an increase in S from N to 2N is predicted to result in an approximately threefold increase in pore conductance. Calculated maximal conductances for the beta-barrel models are compared with experimental values for porins and for K+ channels.

Alanine

Seven-helix bundles: molecular modeling via restrained molecular dynamics.

Simulated annealing via restrained molecular dynamics (SA/MD) has been used to model compact bundles of seven approximately (anti)parallel alpha-helices. Seven such helix bundles occur, e.g., in bacteriorhodopsin, in rhodopsin, and in the channel-forming N-terminal domain of Bacillus thuringiensis delta-endotoxin. Two classes of model are considered: (a) those consisting of seven Ala20 peptide chains; and (b) those containing a single polypeptide chain, made up of seven Ala20 helices linked by GlyN interhelix loops (where N = 5 or 10). Three different starting C alpha templates for SA/MD are used, in which the seven helices are arranged (a) on a left-handed circular template, (b) on a bacteriorhodopsin-like template, or (c) on a zig-zag template. The ensembles of models generated by SA/MD are analyzed in terms of their geometry and energetics, and the most stable structures from each ensemble are examined in greater detail. Structures resembling bacteriorhodopsin and structures resembling delta-endotoxin are both represented among the most stable structures. delta-Endotoxin-like structures arise from both circular and bacteriorhodopsin-like C alpha templates. A third helix-packing mode occurs several times among the stable structures, regardless of the C alpha template and of the presence or absence of interhelix loops. It is characterized by a "4 + 1" core, in which four helices form a distorted left-handed supercoil around a central, buried helix. The remaining two helices pack onto the outside of the core. This packing mode is comparable with that proposed for rhodopsin on the basis of two-dimensional electron crystallographic and sequence analysis studies.

Amino Acid Sequence

Modelling membrane proteins using structural restraints.

Here we present a procedure for modelling membrane proteins which employs molecular dynamics simulations incorporating target restraints derived from low-resolution structures alongside distance restraints derived from mutagenesis data. The application of the modelling procedure to the closed conformation of the pore domain of the nicotinic acetylcholine receptor is described. This domain is formed by a parallel bundle of five M2 helices. Each M2 helix is kinked due to cumulative distortions of backbone (phi, psi) values. The central region of M2 may adopt a more distorted conformation. This would enable a ring of conserved leucine residues (one from each M2 helix) to pack together, occluding the central pore and thus preventing ion permeation. Molecular dynamics simulations on isolated helices that kink formation is not an inherent property of M2.

Amino Acid Sequence