Search

The oblique insertion of type 1 viral fusion peptides into the cell membrane of the host cell has been shown previously to be an essential element of viral fusion. The actual physical explanation of the cause of the oblique insertion has been the subject of speculation. In this study the physical properties of the fusion peptide surface have been determined computationally and compared to the tilt angles determined both experimentally and by the use of molecular dynamics. It has been shown that the relationship between the distribution of lipophilic potential over the peptide surface and the peptide geometry control the tilt angle of the peptide in a biomimetic DMPC bilayer whereas the depth of penetration into the bilayer appears to be determined by the electrostatic potential and hydrogen bonding at the C-terminus.

The energetic profile of an ion translated along the axis of an ion channel should reveal whether the structure corresponds to a functionally open or closed state of the channel. In this study, we explore the combined use of Poisson-Boltzmann electrostatic calculations and evaluation of van der Waals interactions between ion and pore to provide an initial appraisal of the gating state of a channel. This approach is exemplified by its application to the bacterial inward rectifier potassium channel KirBac3.1, where it reveals the closed gate to be formed by a ring of leucine (L124) side chains. We have extended this analysis to a comparative survey of gating profiles, including model hydrophobic nanopores, the nicotinic acetylcholine receptor, and a number of potassium channel structures and models. This enables us to identify three gating regimes, and to show the limitation of this computationally inexpensive method. For a (closed) gate radius of 0.4 nm < R < 0.8 nm, a hydrophobic gate may be present. For a gate radius of 0.2 nm < R < 0.4 nm, both electrostatic and van der Waals interactions will contribute to the barrier height. Below R = 0.2 nm, repulsive van der Waals interactions are likely to dominate, resulting in a sterically occluded gate. In general, the method is more useful when the channel is wider; for narrower channels, the flexibility of the protein may allow otherwise-unsurmountable energetic barriers to be overcome.

Despite a growing repertoire of membrane protein structures (currently approximately 120 unique structures), considerations of low resolution and crystallization in the absence of a lipid bilayer require the development of techniques to assess the global quality of membrane protein folds. This is also the case for assessment of, e.g. homology models of human membrane proteins based on structures of (distant) bacterial homologues. Molecular dynamics (MD) simulations may be used to help evaluate the quality of a membrane protein structure or model. We have used a structure of the bacterial ABC transporter MsbA which has the correct transmembrane helices but an incorrect handedness and topology of their packing to test simulation methods of quality assessment. An MD simulation of the MsbA model in a lipid bilayer is compared to a simulation of another bacterial ABC transporter, BtuCD. The latter structure has demonstrated good conformational stability in the same bilayer environment and over the same timescale (20 ns) as for the MsbA model simulation. A number of comparative analyses of the two simulations were performed to assess changes in the structural integrity of each protein. The results show a significant difference between the two simulations, chiefly due to the dramatic structural deformations of MsbA. We therefore propose that MD could become a useful quality control tool for membrane protein structural biology. In particular, it provides a way in which to explore the global conformational stability of a model membrane protein fold.

The structures of three bacterial outer membrane proteins (OmpA, OmpX and PagP) have been determined by both X-ray diffraction and NMR. We have used multiple (7 x 15 ns) MD simulations to compare the conformational dynamics resulting from the X-ray versus the NMR structures, each protein being simulated in a lipid (DMPC) bilayer. Conformational drift was assessed via calculation of the root mean square deviation as a function of time. On this basis the 'quality' of the starting structure seems mainly to influence the simulation stability of the transmembrane beta-barrel domain. Root mean square fluctuations were used to compare simulation mobility as a function of residue number. The resultant residue mobility profiles were qualitatively similar for the corresponding X-ray and NMR structure-based simulations. However, all three proteins were generally more mobile in the NMR-based than in the X-ray simulations. Principal components analysis was used to identify the dominant motions within each simulation. The first two eigenvectors (which account for >50% of the protein motion) reveal that such motions are concentrated in the extracellular loops and, in the case of PagP, in the N-terminal alpha-helix. Residue profiles of the magnitude of motions corresponding to the first two eigenvectors are similar for the corresponding X-ray and NMR simulations, but the directions of these motions correlate poorly reflecting incomplete sampling on a approximately 10 ns timescale.

The structures of the mammalian water transport protein Aqp1 and of its bacterial homologue GlpF enables us to test whether homology models can be used to explore relationships between structure, dynamics and function in mammalian transport proteins. Molecular dynamics simulations (totalling almost 40 ns) were performed starting from: the X-ray structure of Aqp1; a homology model of Aqp1 based on the GlpF structure; and intermediate resolution structures of Aqp1 derived from electron microscopy. Comparisons of protein RMSDs vs. time suggest that the homology models are of comparable conformational stability to the X-ray structure, whereas the intermediate resolution structures exhibit significant conformation drift. For simulations based on the X-ray structure and on homology models, the flexibility profile vs. residue number correlates well with the crystallographic B-values for each residue. In the simulations based on intermediate resolution structures, mobility of the highly conserved NPA loops is substantially higher than in the simulations based on the X-ray structure or the homology models. Pore radius profiles remained relatively constant in the X-ray and homology model simulations but showed substantial fluctuations (reflecting the higher NPA loop mobility) in the intermediate resolution simulations. The orientation of the dipoles of water molecules within the pore is of key importance in maintaining low proton permeability through Aqp1. This property seems to be quite robust to the starting model used in the simulation. These simulations suggest that homology models based on bacterial homologues may be used to derive functionally relevant information on the structural dynamics of mammalian transport proteins.

Bacterial homologues of mammalian potassium channels provide structures of two states of a gated K channel. Thus, the crystal structure of KcsA represents a closed state whilst that of MthK represents an open state. Using homology modelling and molecular dynamics simulations we have built a model of the transmembrane domain of KcsA in an open state and have compared its conformational stability with that of the same domain of KcsA in a closed state. Approximate Born energy calculations of monovalent cations within the two KcsA channel states suggest that the intracellular hydrophobic gate in the closed state provides a barrier of height approximately 5 kT to ion permeation, whilst in the open state the barrier is absent. Simulations (10 ns duration) in an octane slab (a simple membrane mimetic) suggest that closed- and open-state models are of comparable conformational stability, both exhibiting conformational drifts of approximately 3.3 A Calpha RMSD relative to the respective starting models. Substantial conformational fluctuations are observed in the intracellular gate region during both simulations (closed state and open state). In the simulation of open-state KcsA, rapid ( < 5 ns) exit of all three K+ ions occurs through the intracellular mouth of the channel. Helix kink and swivel motion is observed at the molecular hinge formed by residue G99 of the M2 helix. This motion is more substantial for the open- than for the closed-state model of the channel.

We use molecular dynamics simulations to investigate the position-dependent free energy of a potassium ion in a model of an ion channel formed by the synthetic amphipathic leucine-serine peptide, LS3. The channel model is a parallel bundle of six LS3 helices around which are packed 146 methane-like spheres in order to mimic a membrane. At either end of and within the channel are 1051 water molecules, plus four ions (two potassium and two chloride). The free energy of a potassium ion in the channel was estimated using the weighted histogram analysis (WHAM) method. This is the first time to our knowledge that such a calculation has been carried out as a function of the position of an ion in three dimensions within a channel. The results indicate that for this channel, which is lined by hydrophilic serine sidechains, there is a relatively weak dependence of the free energy on the axial/off-axial position of the ion. There are some off-axis local minima, especially in the C-terminal half of the channel. Using the free energy results, a single channel current-voltage curve was estimated using a one-dimensional Nernst-Planck equation. Although reasonable agreement with experiment is achieved for K(+) ions flowing from the N-terminal to the C-terminal mouth, in the opposite direction the current is underestimated. This underestimation may be a consequence of under-sampling of the conformational dynamics of the channel. We suggest that our simulations may have captured, for example, a sub-conductance level (i.e. an incompletely open state) of the LS3 channel.

The gating mechanism of a bacterial potassium channel, KcsA, has been investigated via multi-nanosecond molecular dynamic simulations of the channel molecules embedded in a fully solvated palmitoyloleoylphosphatidylcholine bilayer. Four events are seen in which a cation (K(+) or, in one case, Na(+)) initially present in the central cavity exits through the intracellular mouth (the presumed gate) of the channel. Whilst in the cavity a cation interacts with the sidechain T107 O gamma atom of one of the subunits prior to its exit from the channel. Secondary structure analysis as a function of time reveals a break in the helicity of one of the M2 helices. This break is expected to lend flexibility to the helices, enabling them to "open" (minimum pore radius >0.13 nm) and "close" (minimum pore radius <0.13 nm) the channel. Fluctuations in the pore radius at the intracellular gate region are of the order of 0.05 nm, with an average radius in the region of the gate of ca. 0.1 nm. However, around the time of exit of a cation, the pore widens to about 0.15 nm. The distances between the C alpha atoms of the inner helices M2 reveal a coupled increase and decrease between the opposite pair of helices at about the time of exit of the ion. This suggests a breathing motion of the M2 helices that may form the basis for a gating mechanism.

In this paper we describe a method for setting up an atomistic simulation of a membrane protein in a hydrated lipid bilayer and report the effect of differing electrostatic parameters on the drift in the protein structure during the subsequent simulation. The method aims to generate a suitable cavity in the interior of a lipid bilayer, using the solvent-accessible surface of the protein as a template, during the course of a short steered molecular dynamics simulation of a solvated lipid membrane. This is achieved by a two-stage process: firstly, lipid molecules whose headgroups are inside a cylindrical volume equivalent to that defined by the protein surface are removed; then the protein-lipid interface is optimized by applying repulsive forces perpendicular to the protein surface, and of gradually increased magnitude, to the remaining lipid atoms inside the volume occupied by the protein surface until it is emptied. The protein itself may then be inserted. Using the bacterial membrane proteins KcsA and FhuA as test cases, we show how the method achieves the formation of a suitable cavity in the interior of a dimyristoylphosphatidylcholine lipid bilayer without perturbing the configuration of the non-interfacial regions of the previously equilibrated lipid bilayer, even in cases of membrane proteins with irregular geometrical shapes. In addition, we compare subsequent simulations in which the long-range electrostatic interactions are treated via either a cut-off or a particle-mesh Ewald method. The results show that the drift from the initial structure is less in the latter case, especially for KcsA and for the non-core secondary structural elements (i.e. surface loops) of both proteins.

NB is short auxiliary protein with ca. 100 amino acids, encoded in the viral genome of influenza B. It is believed to be similar to M2 from influenza A and Vpu from HIV-1 in that it demonstrates ion channel activity. Channels formed by the protein can be blocked by amantadine. We have synthesized the putative transmembrane segment of NB (IRG S20 IIITICVSL I30 VILIVFGCI A40 KIFI (NB, Lee)). Reconstituted in a lipid bilayer, the peptide shows channel activity. The addition of amantadine leads to dose-dependent loss of channel activity. Channel blocking is reversible. Channel behaviour of the peptide in the presence of amantadine is in accordance with findings for the intact channel. Thus, the synthetic transmembrane peptide captures the ion channel activity of the intact NB protein.

A homology model of the pore domain of the Shaker K+ channel has been constructed using a bacterial K+ channel, KcsA, as a template structure. The model is in agreement with mutagenesis and sequence variability data. A number of structural features are conserved between the two channels, including a ring of tryptophan sidechains on the outer surface of the pore domain at the extracellular end of the helix bundle, and rings of acidic sidechains close to the extracellular mouth of the channel. One of these rings, that formed by four Asp447 sidechains at the mouth of the Shaker pore, is shown by pK(A) calculations to be incompletely ionized at neutral pH. The potential energy profile for a K+ ion moved along the central axis of the Shaker pore domain model selectivity filter reveals a shallow well, the depth of which is modulated by the ionization state of the Asp447 ring. This is more consistent with the high cation flux exhibited by the channel in its conductance value of 19 pS.

Bacteriorhodopsin (BR) is a membrane protein which pumps protons through the plasma membrane. Transmembrane BR helical segments are subjected to simulation studies in order to investigate the effect of bilayer environment in various simulation conditions. A bilayer potential is introduced to the system to mimic the lipid membrane. The structures from the simulations are compared with the experimentally determined structures in terms of geometrical properties. Electrostatic contribution to the helix packing is also investigated. The simulation results show that the packing geometry of the transmembrane helices is highly affected by the bilayer potential. The results obtained from the simulations may be used for further simulation studies and analysis in investigating transmembrane helix packing.

Bacteriorhodopsin (BR) is a membrane protein which pumps protons through the plasma membrane. Seven transmembrane BR helical segments are subjected to simulation studies in order to investigate the packing process of transmembrane helices. A Monte Carlo simulated annealing protocol is employed to optimize the helix bundle system. Helix packing is optimized according to a semi-empirical potential mainly composed of six components: a bilayer potential, a crossing angle potential, a helix dipole potential, a helix-helix distance potential, a helix orientation potential and a helix-helix distance restraint potential (a loop potential). Necessary parameters are derived from theoretical studies and statistical analysis of experimentally determined protein structures. The structures from the simulations are compared with the experimentally determined structures in terms of geometry. The structures generated show similar shapes to the experimentally suggested structure even without the helix-helix distance restraint potential. However, the relative locations of individual helices were reproduced only when the helix-helix distance restraint potential was used with restraint conditions. Our results suggest that transmembrane helix bundles resembling those observed experimentally may be generated by simulations using simple potentials.

The nicotinic acetylcholine receptor (nAChR) is the archetypal ligand-gated ion channel. A model of the alpha7 homopentameric nAChR is described in which the pore-lining M2 helix bundle is treated atomistically and the remainder of the molecule is treated as a "low resolution" cylinder. The surface charge on the cylinder is derived from the distribution of charged amino acids in the amino acid sequence (excluding the M2 segments). This model is explored in terms of its predicted single-channel properties. Based on electrostatic potential profiles derived from the model, the one-dimensional Poisson-Nernst-Planck equation is used to calculate single-channel current/voltage curves. The predicted single-channel conductance is three times higher (ca. 150 pS) than that measured experimentally, and the predicted ion selectivity agrees with the observed cation selectivity of nAChR. Molecular dynamics (MD) simulations are used to estimate the self-diffusion coefficients (D) of water molecules within the channel. In the narrowest region of the pore, D is reduced ca. threefold relative to that of bulk water. Assuming that the diffusion of ions scales with that of water, this yields a revised prediction of the single-channel conductance (ca. 50 pS) in good agreement with the experimental value. We conclude that combining atomistic (MD) and continuum electrostatics calculations is a promising approach to bridging the gap between structure and physiology of ion channels.

The aim of this study is to investigate if the packing motifs of native transmembrane helices can be produced by simulations with simple potentials and to develop a method for the rapid generation of initial candidate models for integral membrane proteins composed of bundles of transmembrane helices. Constituent residues are mapped along the helix axis in order to maintain the amino acid sequence-dependent properties of the helix. Helix packing is optimized according to a semi-empirical potential mainly composed of four components: a bilayer potential, a crossing angle potential, a helix dipole potential and a helix-helix distance potential. A Monte Carlo simulated annealing protocol is employed to optimize the helix bundle system. Necessary parameters are derived from theoretical studies and statistical analysis of experimentally determined protein structures. Preliminary testing of the method has been conducted with idealized seven Ala(20) helix bundles. The structures generated show a high degree of compactness. It was observed that both bacteriorhodopsin-like and delta-endotoxin-like structures are generated in seven-helix bundle simulations, within which the composition varies dependent upon the cooling rate. The simulation method has also been employed to explore the packing of N = 4 and N = 12 transmembrane helix bundles. The results suggest that seven and 12 transmembrane helix bundles resembling those observed experimentally (e.g., bacteriorhodopsin, rhodopsin and cytochrome c oxidase subunit I) may be generated by simulations using simple potentials.

Ion channels may be formed by bundles of amphipathic alpha-helices aligned parallel to one another and spanning a lipid bilayer membrane, with the hydrophilic faces of the helices lining a central pore. In order to provide insight into the packing of such helices in bundles, a method has been developed to evaluate hydrophilic surface maps of amphipathic alpha-helices and to display these surfaces in a readily interpretable form. The procedure is based upon empirical energy calculations of interactions of a water molecule with an amphipathic alpha-helix. The method has been applied to three channel-forming peptides: Staphylococcal delta-toxin; alamethicin; and a synthetic leucine- and serine-containing peptide. Particular emphasis is placed upon the effects of sidechain conformational flexibility on hydrophilic surface maps. A family of models of the delta-toxin helix is generated by a simulated annealing procedure. The results of hydrophilic surface map analyses provide more exact definition of the centre of the hydrophilic face of amphipathic helices, and of the variation of the position of the centre in response to changes in sidechain conformation. This information is used to define families of preliminary models for a given ion channel, as is illustrated for delta-toxin.

Peptaibols are considered as models of those ion channels which consist of a bundle of transbilayer helices surrounding a central pore. X-Ray diffraction and NMR studies have yielded high resolution structures for several peptaibols. In conjunction with other spectroscopic investigations and molecular dynamics simulations, these studies suggest that peptaibols form proline-kinked alpha-helices, and that there may be "hinge-bending" movement of the helix in the region of the central proline residue. The amphipathicity of peptaibol helices is analyzed in relation to their channel-forming properties. Studies of the interactions of peptaibols with lipid bilayers suggest that they are helical when in a membrane-like environment, and that the helix orientation relative to the bilayer is sensitive to the peptaibol:lipid ratio, and to the degree of hydration of the bilayer. Electrical studies reveal that many peptaibols form multiple-conductance level channels in a voltage-dependent fashion. Analysis of conductance levels provides support for the "barrel stave" model of channel formation, whereby different conductance levels correspond to different numbers of monomers in a helix bundle. Alternative models for voltage-activation are discussed, and the roles of molecular dipoles and of hinge-bending in this process are considered. Two molecular models for an N = 6 bundle of alamethicin helices are presented and their electrostatic properties analyzed. The relevance of studies of peptaibols to channel and transport proteins in general is considered.

Zervamicin-IIB (Zrv-IIB) is a 16 residue peptaibol which forms voltage-activated, multiple conductance level channels in planar lipid bilayers. A molecular model of Zrv-IIB channels is presented. The structure of monomeric Zrv-IIB is based upon the crystal structure of Zervamicin-Leu. The helical backbone is kinked by a hydroxyproline residue at position 10. Zrv-IIB channels are modelled as helix bundles of from 4 to 8 parallel helices surrounding a central pore. The monomers are packed with their C-terminal helical segments in close contact, and the bundles are stabilized by hydrogen bonds between glutamine 11 and hydroxyproline 10 of adjacent helices. Interaction energy profiles for movement of three different probes species (K+, Cl- and water) through the central pore are analyzed. The conformations of: (a) the sidechain of glutamine 3; (b) the hydroxyl group of hydroxyproline 10; and (c) the C-terminal hydroxyl group are "optimized" in order to maximize favourable interactions between the channel and the probes, resulting in favourable interaction energy profiles for all three. This suggests that conformational flexibility of polar sidechains enables the channel lining to mimic an aqueous environment.

The zervamicins (Zrv) are a family of 16 residue peptaibol channel formers, related to the 20 residue peptaibol alamethicin (Alm), but containing a higher proportion of polar sidechains. Zrv-IIB forms multi-level channels in planar lipid (diphytanoyl phosphatidylcholine) bilayers in response to cis positive voltages. Analysis of the voltage and concentration dependence of macroscopic conductances induced by Zrv-IIB suggests that, on average, channels contain ca. 13 peptide monomers. Analysis of single channel conductance levels suggests a similar value. The pattern of successive conductance levels is consistent with a modified helix bundle model in which the higher order bundle are distorted within the plane of the bilayer towards a "torpedo" shaped cross-section. The kinetics of intra-burst switching between adjacent conductance levels are shown to be approximately an order of magnitude faster for Zrv-IIB than for Alm. The channel forming properties of the related naturally occurring peptaibols, Zrv-Leu and Zrv-IC, have also been demonstrated, as have those of the synthetic apolar analogue Zrv-Al-16. The experimental studies on channel formation are combined with the known crystallographic structures of Zrv-Al-16 and Zrv-Leu to develop a molecular model of Zrv-IIB channels.

The ion channel of the nicotinic acetylcholine receptor (nAChR) is believed to be lined by transmembrane M2 helices. A "4-8-12" sequence motif, comprising serine (S) or threonine (T) residues at positions 4, 8 and 12 of M2, is conserved between different members, anion and cation selective, of the nAChR superfamily. Parallel bundles of 4-8-12 motif-containing helices are considered as simplified models of ion channels. The relationship between S and T sidechain conformations and channel-ion interactions is explored via evaluation of interaction energies of K+ and of Cl- ions with channel models. Energy calculations are used to determine optimal chi 2 (C alpha-C beta-O gamma-H gamma) values in the presence of K+ or Cl- ions. 4-8-12 motif-containing bundles may form favourable interactions with either cations or anions, dependent upon the chi 2 values adopted. Parallel-helix and tilted-helix bundles are considered, as are heteromeric models designed to mimic the Torpedo nAChR. The main conclusion of the study is that conformational flexibility at chi 2 enables both S and T residues to form favourable interactions with anions or cations. Consequently, there is apparently no difference between S and T residues in their interactions with permeant ions, which suggests that the presence of T vs. S residues within the 4-8-12 motif is not a major mechanism whereby anion/cation selectivity may be generated. The implications of these studies with respect to more elaborate models of nAChR and related receptors are considered.

Channel forming peptides (CFPs) are amphipathic peptides, of length ca. 20 residues, which adopt an alpha-helical conformation in the presence of lipid bilayers and form ion channels with electrophysiological properties comparable to those of ion channel proteins. We have modelled CFP channels as bundles of parallel trans-bilayer helices surrounding a central ion-permeable pore. Ion-channel interactions have been explored via accessible surface area calculations, and via evaluation of changes in van der Waals and electrostatic energies as a K+ ion is translated along the length of the pore. Two CFPs have been modelled: (a) zervamicin-A1-16, a synthetic apolar peptaibol related to alamethicin, and (b) delta-toxin from Staphylococcus aureus. Both of these CFPs have previously been shown to form ion channels in planar lipid bilayers, and have been shown to have predominantly helical conformations. Zervamicin-A1-16 channels were modelled as bundles of 4 to 8 parallel helices. Two related helix bundle geometries were explored. K(+)-channel interactions have been shown to involve exposed backbone carbonyl oxygen atoms. delta-Toxin channels were modelled as bundles of 6 parallel helices. Residues Q3, D11 and D18 generate favourable K(+)-channel interactions. Rotation of W15 about its C beta-C gamma bond has been shown to be capable of occluding the central pore, and is discussed as a possible model for sidechain conformational changes in relation to ion channel gating.