Structure and function of disk aggregates of the coat protein of tobacco mosaic virus

Experiments have been carried out on the coat protein of tobacco mosaic virus (TMVP) to test for the occurrence of the previously postulated RNA-induced direct switching, during in vitro assembly of tobacco mosaic virus (TMV), of the subunit packing from the cylindrical bilayer disk to the virus helical arrangement. No evidence was found for such RNA-induced switching and no evidence for the direct participation of the bilayer disk in either the nucleation or elongation phases of the in vitro virus assembly. Instead, virus assembly proceeds by an initiation step involving the binding of the RNA to the previously characterized two-plus turn helical aggregate that is formed from small oligomers of subunits. However, a bilayer disk, which has been characterized in high ionic strength crystals, has been observed in low ionic strength virus assembly solutions only as a transient species upon depolymerization of dimers of bilayer disks formed in solution at high ionic strength, and not as an equilibrium species of TMVP.     

Initial structural and dynamic characterization of the M2 protein transmembrane and amphipathic helices in lipid bilayers

Amphipathic helices in membrane proteins that interact with the hydrophobic/hydrophilic interface of the lipid bilayer have been difficult to structurally characterize. Here, the backbone structure and orientation of an amphipathic helix in the full-length M2 protein from influenza A virus has been characterized. The protein has been studied in hydrated DMPC/DMPG lipid bilayers above the gel to liquid-crystalline phase transition temperature by solid-state NMR spectroscopy. Characteristic PISA (Polar Index Slant Angle) wheels reflecting helical wheels have been observed in uniformly aligned bilayer preparations of both uniformly 15N labeled and amino acid specific labeled M2 samples. Hydrogen/deuterium exchange studies have shown the very slow exchange of some residues in the amphipathic helix and more rapid exchange for the transmembrane helix. These latter results clearly suggest the presence of an aqueous pore. A variation in exchange rate about the transmembrane helical axis provides additional support for this claim and suggests that motions occur about the helical axes in this tetramer to expose the entire backbone to the pore.     

Contribution of the γ-secretase subunits to the formation of catalytic pore of presenilin 1 protein

γ-Secretase is an intramembrane-cleaving protease related to the etiology of Alzheimer disease. γ-Secretase is a membrane protein complex composed of presenilin (PS) and three indispensable subunits: nicastrin, Aph-1, and Pen-2. PS functions as a protease subunit forming a hydrophilic catalytic pore structure within the lipid bilayer. However, it remains unclear how other subunits are involved in the pore formation. Here, we show that the hydrophilic pore adopted with an open conformation has already been formed by PS within the immature γ-secretase complex. The binding of the subunits induces the close proximity between transmembrane domains facing the catalytic pore. We propose a model in which the γ-secretase subunits restrict the arrangement of the transmembrane domains of PS during the formation of the functional structure of the catalytic pore.     

Charged multivesicular body protein 2B (CHMP2B) of the endosomal sorting complex required for transport-III (ESCRT-III) polymerizes into helical structures deforming the plasma membrane

The endosomal sorting complexes required for transport (ESCRT-0-III) allow membrane budding and fission away from the cytosol. This machinery is used during multivesicular endosome biogenesis, cytokinesis, and budding of some enveloped viruses. Membrane fission is catalyzed by ESCRT-III complexes made of polymers of charged multivesicular body proteins (CHMPs) and by the AAA-type ATPase VPS4. How and which of the ESCRT-III subunits sustain membrane fission from the cytoplasmic surface remain uncertain. In vitro, CHMP2 and CHMP3 recombinant proteins polymerize into tubular helical structures, which were hypothesized to drive vesicle fission. However, this model awaits the demonstration that such structures exist and can deform membranes in cellulo. Here, we show that depletion of VPS4 induces specific accumulation of endogenous CHMP2B at the plasma membrane. Unlike other CHMPs, overexpressed full-length CHMP2B polymerizes into long, rigid tubes that protrude out of the cell. CHMP4s relocalize at the base of the tubes, the formation of which depends on VPS4. Cryo-EM of the CHMP2B membrane tubes demonstrates that CHMP2B polymerizes into a tightly packed helical lattice, in close association with the inner leaflet of the membrane tube. This association is tight enough to deform the lipid bilayer in cases where the tubular CHMP2B helix varies in diameter or is closed by domes. Thus, our observation that CHMP2B polymerization scaffolds membranes in vivo represents a first step toward demonstrating its structural role during outward membrane deformation.     

Influenza hemagglutinin-mediated fusion pores connecting cells to planar membranes: flickering to final expansion

We have studied the fusion between voltage-clamped planar lipid bilayers and influenza virus infected MDCK cells, adhered to one side of the bilayer, using measurements of electrical admittance and fluorescence. The changes in currents in-phase and 90 degrees out-of- phase with respect to the applied sinusoidal voltage were used to monitor the addition of the cell membrane capacitance to that of the lipid bilayer through a fusion pore connecting the two membranes. When ethidium bromide was included in the solution of the cell-free side of the bilayer, increases in cell fluorescence accompanied tee admittance changes, independently confirming that these changes were due to formation of a fusion pore. Fusion required acidic pH on the cell- containing side and depended on temperature. For fusion to occur, the influenza hemagglutinin (HA) had to be cleaved into HA1 and HA2 subunits. The incorporation of gangliosides into the planar bilayers greatly augmented fusion. Fusion pores developed in four distinct stages after acidification: (a) a pre-pore, electrically quiescent stage; (b) a flickering stage, with 1-2 nS pores opening and closing repetitively; (c) an irreversibly opened stage, in which pore conductances varied between 2 and 100 nS and exhibited diverse kinetics; (d) a fully opened stage, initiated by an instantaneous, time- resolution limited, increase in conductance leveling at approximately 500 nS. The expansion of pores by stages has also been shown to occur during exocytosis in mast cells and fusion of HA-expressing cells and erythrocytes. We conclude that essential features of fusion pores are produced with proteins in just one of the two fusing membranes.     

SNARE Complex Zipping as a Driving Force in the Dilation of Proteinaceous Fusion Pores

The assembly of SNARE proteins into a tight complex has been hypothesized to drive membrane fusion. A model of the initial fusion pore as a proteinaceous channel formed by SNARE proteins places their membrane anchors in separate membranes. This leaves the possibility of a final assembly step that brings the membrane anchors together and drives fusion pore expansion. The present study develops a model for expansion in which the final SNARE complex zipping step drives a transition from a proteinaceous fusion pore to a lipidic fusion pore. An estimate of the energy released upon merger of the helical segments of the SNARE motifs with the helical segments of the membrane anchors indicates that completing the assembly of a few SNARE complexes can overcome the elastic energy that opposes lipid bilayer deformation into a narrow fusion pore. The angle between the helical axes of the membrane anchor and SNARE motif serves as a useful reaction coordinate for this transition. Energy was calculated as a function of this angle, incorporating contributions from membrane bending, SNARE complex assembly, membrane anchor flexing and hydrophobic interactions. The rate of this transition was evaluated as a process of diffusion over the barrier imposed by these combined energies, and the rates estimated were consistent with experimental measurements.     

Windex: a toolset for indexing helices

We describe here a set of procedures and algorithms that may be used as an aid in determining the indexing rule of a helical specimen. Crystallizing macromolecules into helical arrays has the potential to speed up and simplify the process of three-dimensional reconstruction of the macromolecular structure. The process of helical reconstruction has been largely automated except for the critical first step of indexing the helical diffraction pattern. This is quite often the rate-limiting step in the overall process, particularly in the case of large helical tubes, which have complicated helical diffraction patterns that may vary from tube to tube. We have developed a set of procedures, supported by a graphical user interface, that provide a straightforward and semi-automated approach to indexing a helical structure. The new procedures have been tested using a number of helical specimens, including TMV, acto-myosin, decorated microtubules, and a variety of helical tubes of a bacterial membrane protein.     

The exocytotic fusion pore interface: a model of the site of neurotransmitter release

infrastructurel techniques have shown that an early event in the exocytotic fusion of a secretory vesicle is the formation of a narrow, water-filled pore spanning both the vesicle and plasma membranes and connecting the lumen of the secretory vesicle to the extracellular environment. Smaller precursors of the exocytotic fusion pore have been detected using electrophysio-logical techniques, which reveal a dynamic fusion pore that quickly expands to the size of the pores seen with electron microscopy. While it is clear that in the latter stages of expansion, when the size of the fusion pore is several orders of magnitude bigger than any known macromolecule, the fusion pore must be mainly made of lipids, the structure of the smaller precursors is unknown. Patch-clamp measurements of the activity of individual fusion pores in mast cells have shown that the fusion pore has some unusual and unexpected properties, namely that there is a large flux of lipid through the pore and the rate of pore closure has a discontinuous temperature dependency, suggesting a purely lipidic fusion pore. Moreover, comparisons of experimental data with theoretical fusion pores and with breakdown pores support the view that the fusion pore is initially a pore through a single bilayer, as would be expected for membrane fusion proceeding through a hemifusion mechanism. Based on these observations we present a model where the fusion pore is initially a pore through a single bilayer. Fusion pore formation is regulated by a macromolecular scaffold of proteins that is responsible for bringing the plasma membrane into a highly curved dimple very close to a tense secretory granule membrane, creating the architecture where the strongly attractive hydrophobic force causes the membranes to form a ‘hemifusion’ intermediate. Membrane fusion is completed by the formation of an aqueous pore after rupture of the shared bilayer. We also propose that the microenvironment of the interface when the pore first opens, dominated by the charged groups on the secretory vesicle matrix and phospholipids, will greatly influence the release of secretory products.     

Distinct structural elements in the first membrane-spanning segment of the epithelial sodium channel

Epithelial Na+ channels (ENaCs) comprise three subunits that have been proposed to be arranged in either an alpha2betagamma or a higher ordered configuration. Each subunit has two putative membrane-spanning segments (M1 and M2), intracellular amino and carboxyl termini, and a large extracellular loop. We have used the TOXCAT assay (a reporter assay for transmembrane segment homodimerization) to identify residues within the transmembrane segments of ENaC that may participate in important structural interactions within ENaC, with which we identified a candidate site within alphaM1. We performed site-directed mutagenesis at this site and found that, although the mutants reduced channel activity, ENaC protein expression at the plasma membrane was unaffected. To deduce the role of alphaM1 in the pore structure of ENaC, we performed tryptophan-scanning mutagenesis throughout alphaM1 (residues 110-130). We found that mutations within the amino-terminal part of alphaM1 had effects on activity and selectivity with a periodicity consistent with a helical structure but no effect on channel surface expression. We also observed that mutations within the carboxyl-terminal part of alphaM1 had effects on activity and selectivity but with no apparent periodicity. Additionally, these mutants reduced channel surface expression. Our data support a model in which the amino-terminal half of alphaM1 is alpha-helical and packs against structural element(s) that contribute to the ENaC pore. Furthermore, these data suggest that the carboxyl-terminal half of alphaM1 may be helical or assume a different conformation and may be involved in tertiary interactions essential to proper channel folding or assembly. Together, our data suggest that alphaM1 is divided into two distinct regions.     

Molecular dynamics simulations of hydrophilic pores in lipid bilayers

Hydrophilic pores are formed in peptide free lipid bilayers under mechanical stress. It has been proposed that the transport of ionic species across such membranes is largely determined by the existence of such meta-stable hydrophilic pores. To study the properties of these structures and understand the mechanism by which pore expansion leads to membrane rupture, a series of molecular dynamics simulations of a dipalmitoylphosphatidylcholine (DPPC) bilayer have been conducted. The system was simulated in two different states; first, as a bilayer containing a meta-stable pore and second, as an equilibrated bilayer without a pore. Surface tension in both cases was applied to study the formation and stability of hydrophilic pores inside the bilayers. It is observed that below a critical threshold tension of approximately 38 mN/m the pores are stabilized. The minimum radius at which a pore can be stabilized is 0.7 nm. Based on the critical threshold tension the line tension of the bilayer was estimated to be approximately 3 x 10(-11) N, in good agreement with experimental measurements. The flux of water molecules through these stabilized pores was analyzed, and the structure and size of the pores characterized. When the lateral pressure exceeds the threshold tension, the pores become unstable and start to expand causing the rupture of the membrane. In the simulations the mechanical threshold tension necessary to cause rupture of the membrane on a nanosecond timescale is much higher in the case of the equilibrated bilayers, as compared with membranes containing preexisting pores.     

The structure of the M2 channel-lining segment from the nicotinic acetylcholine receptor

The structures of functional peptides corresponding to the predicted channel-lining M2 segment of the nicotinic acetylcholine (AChR) were determined using solution NMR experiments on micelle samples, and solid-state NMR experiments on bilayer samples. The AChR M2 peptide forms a straight transmembrane alpha-helix, with no kinks. M2 inserts in the lipid bilayer at an angle of 12 degrees relative to the bilayer normal, with a rotation about the helix long axis such that the polar residues face the N-terminus of the peptide, which is assigned to be intracellular. A molecular model of the AChR channel pore, constructed from the solid-state NMR 3-D structure of the AChR M2 helix in the membrane assuming a pentameric organization, results in a funnel-like architecture for the channel with the wide opening on the N-terminal intracellular side. A central narrow pore has a diameter ranging from about 3.0 A at its narrowest, to 8.6 A at its widest. Nonpolar residues are predominantly on the exterior of the bundle, while polar residues line the pore. This arrangement is in fair agreement with evidence collected from permeation, mutagenesis, affinity labeling and cysteine accessibility measurements. A pentameric M2 helical bundle may, therefore, represent the structural blueprint for the inner bundle that lines the channel of the nicotinic AChR.     

The internal cavity of the staphylococcal alpha-hemolysin pore accommodates approximately 175 exogenous amino acid residues

The cavity within the cap domain of the transmembrane staphylococcal alpha-hemolysin (alphaHL) pore is roughly a sphere of diameter approximately 45 A (molecular surface volume approximately 39,500 A(3)). We tested the ability of the cavity to accommodate exogenous polypeptide chains. Concatemerized Gly/Ser-containing sequences ("loops", L; number of repeats = n; number of residues = 10n + 5, n = 0-21) were inserted at a position located within the cavity of the fully assembled heptameric alphaHL pore. Homomeric pores containing 25 or less residues in each loop (n or= 7, only one L subunit was incorporated. As the inserted loop was lengthened, transient closures were observed in planar bilayer experiments with single pores. However, L(1)W(6) pores with very long loops (n = 14 and 21) had unitary conductance values close to those of W(7), suggesting that the loop is extruded through the opening in the cap of the pore into the external medium. Further analysis of bilayer recordings and electrophoretic migration patterns indicates that the upper capacity of the cavity is approximately 175 amino acids. The findings suggest that small functional peptides or proteins might be assembled within the alphaHL pore.     

Ruby-Helix: an implementation of helical image processing based on object-oriented scripting language

Helical image analysis in combination with electron microscopy has been used to study three-dimensional structures of various biological filaments or tubes, such as microtubules, actin filaments, and bacterial flagella. A number of packages have been developed to carry out helical image analysis. Some biological specimens, however, have a symmetry break (seam) in their three-dimensional structure, even though their subunits are mostly arranged in a helical manner. We refer to these objects as "asymmetric helices". All the existing packages are designed for helically symmetric specimens, and do not allow analysis of asymmetric helical objects, such as microtubules with seams. Here, we describe Ruby-Helix, a new set of programs for the analysis of "helical" objects with or without a seam. Ruby-Helix is built on top of the Ruby programming language and is the first implementation of asymmetric helical reconstruction for practical image analysis. It also allows easier and semi-automated analysis, performing iterative unbending and accurate determination of the repeat length. As a result, Ruby-Helix enables us to analyze motor-microtubule complexes with higher throughput to higher resolution.     

Melittin-induced bilayer leakage depends on lipid material properties: evidence for toroidal pores

The membrane-lytic peptide melittin has previously been shown to form pores in lipid bilayers that have been described in terms of two different structural models. In the "barrel stave" model the bilayer remains more or less flat, with the peptides penetrating across the bilayer hydrocarbon region and aggregating to form a pore, whereas in the "toroidal pore" melittin induces defects in the bilayer such that the bilayer bends sharply inward to form a pore lined by both peptides and lipid headgroups. Here we test these models by measuring both the free energy of melittin transfer (DeltaG degrees ) and melittin-induced leakage as a function of bilayer elastic (material) properties that determine the energetics of bilayer bending, including the area compressibility modulus (K(a)), bilayer bending modulus (k(c)), and monolayer spontaneous curvature (R(o)). The addition of cholesterol to phosphatidylcholine (PC) bilayers, which increases K(a) and k(c), decreases both DeltaG degrees and the melittin-induced vesicle leakage. In contrast, the addition to PC bilayers of molecules with either positive R(o), such as lysoPC, or negative R(o), such as dioleoylglycerol, has little effect on DeltaG degrees , but produces large changes in melittin-induced leakage, from 86% for 8:2 PC/lysoPC to 18% for 8:2 PC/dioleoylglycerol. We observe linear relationships between melittin-induced leakage and both K(a) and 1/R(o)(2). However, in contrast to what would be expected for a barrel stave model, there is no correlation between observed leakage and bilayer hydrocarbon thickness. All of these results demonstrate the importance of bilayer material properties on melittin-induced leakage and indicate that the melittin-induced pores are defects in the bilayer lined in part by lipid molecules.     

Models of beta-amyloid ion channels in the membrane suggest that channel formation in the bilayer is a dynamic process

Here we model the Alzheimer beta-peptide ion channel with the goal of obtaining insight into the mechanism of amyloid toxicity. The models are built based on NMR data of the oligomers, with the universal U-shaped (strand-turn-strand) motif. After 30-ns simulations in the bilayer, the channel dimensions, shapes and subunit organization are in good agreement with atomic force microscopy (AFM). The models use the Abeta(17-42) pentamer NMR-based coordinates. Extension and bending of the straight oligomers lead to two channel topologies, depending on the direction of the curvature: 1), the polar/charged N-terminal beta-strand of Abeta(17-42) faces the water-filled pore, and the hydrophobic C-terminal beta-strand faces the bilayer (CNpNC; p for pore); and 2), the C-terminal beta-strand faces the solvated pore (NCpCN). In the atomistic simulations in a fully solvated DOPC lipid bilayer, the first (CNpNC) channel preserves the pore and conducts solvent; by contrast, hydrophobic collapse blocks the NCpCN channel. AFM demonstrated open pores and collapsed complexes. The final averaged CNpNC pore dimensions (outer diameter 8 nm; inner diameter approximately 2.5 nm) are in the AFM range (8-12 nm; approximately 2 nm, respectively). Further, in agreement with high-resolution AFM images, during the simulations, the channels spontaneously break into ordered subunits in the bilayer; however, we also observe that the subunits are loosely connected by partially disordered inner beta-sheet, suggesting subunit mobility in the bilayer. The cationic channel has strong selective affinity for Ca(2+), supporting experimental calcium-selective beta-amyloid channels. Membrane permeability and consequent disruption of calcium homeostasis were implicated in cellular degeneration. Consequently, the CNpNC channel topology can sign cell death, offering insight into amyloid toxicity via an ion "trap-release" transport mechanism. The observed loosely connected subunit organization suggests that amyloid channel formation in the bilayer is a dynamic, fluid process involving subunit association, dissociation, and channel rearrangements.     

The structure of the M2 channel-lining segment from the nicotinic acetylcholine receptor

The structures of functional peptides corresponding to the predicted channel-lining M2 segment of the nicotinic acetylcholine (AChR) were determined using solution NMR experiments on micelle samples, and solid-state NMR experiments on bilayer samples. The AChR M2 peptide forms a straight transmembrane α-helix, with no kinks. M2 inserts in the lipid bilayer at an angle of 12° relative to the bilayer normal, with a rotation about the helix long axis such that the polar residues face the N-terminus of the peptide, which is assigned to be intracellular. A molecular model of the AChR channel pore, constructed from the solid-state NMR 3-D structure of the AChR M2 helix in the membrane assuming a pentameric organization, results in a funnel-like architecture for the channel with the wide opening on the N-terminal intracellular side. A central narrow pore has a diameter ranging from about 3.0 Å at its narrowest, to 8.6 Å at its widest. Nonpolar residues are predominantly on the exterior of the bundle, while polar residues line the pore. This arrangement is in fair agreement with evidence collected from permeation, mutagenesis, affinity labeling and cysteine accessibility measurements. A pentameric M2 helical bundle may, therefore, represent the structural blueprint for the inner bundle that lines the channel of the nicotinic AChR.     

An iterative Fourier-Bessel algorithm for reconstruction of helical structures with severe Bessel overlap

Classical Fourier-Bessel methodology fails when used to reconstruct helical structures with severe Bessel overlap on the layer lines. In the reconstruction of a peculiar type of double-layered helical tube of GDP-tubulin, we face the problem of Bessel overlap on all the layer lines due to the superposition of the Fourier components from the inner and outer layers of the tube. In order to decompose the Fourier terms of the inner and outer layers more than one image of the tubes must be combined and the orientations of their inner and outer layer helices must be determined. While there is no direct analytical method to determine these orientational parameters, we have devised an iterative Fourier-Bessel algorithm to calculate the correct orientations and thus allow us to obtain a reconstruction from multiple images of the double-layered tubes. The algorithm successfully works for the reconstruction of computer-modeled double-layered helical tubes as well as with real images obtained by cryo-electron microscopy. The algorithm has also been applied with very satisfactory results to the reconstruction of 13-protofilament microtubules, which is another helical structure that suffer Bessel overlap, suggesting its generality.     

Influence of the arrangement and secondary structure of melittin peptides on the formation and stability of toroidal pores

Melittin interactions with lipid bilayers and melittin formed pores are extensively studied to understand the mechanism of the toroidal pore formation. Early experimental studies suggested that melittin peptide molecules are anchored by their positively charged residues located next to the C-terminus to only one leaflet of the lipid bilayer (asymmetric arrangement). However, the recent non-linear spectroscopic experiment suggests a symmetric arrangement of the peptides with the C-terminus of the peptides anchored to both bilayers. Therefore, we present here a computational study that compares the effect of symmetric and asymmetric arrangements of melittin peptides in the toroidal pore formation. We also investigate the role of the peptide secondary structure during the pore formation. Two sets of the symmetric and asymmetric pores are prepared, one with a helical peptide from the crystal structure and the other set with a less helical peptide. We observe a stable toroidal pore being formed only in the system with a symmetric arrangement of the less helical peptides. Based on the simulation results we propose that the symmetric arrangement of the peptides might be more favorable than the asymmetric arrangement, and that the helical secondary structure is not a prerequisite for the formation of the toroidal pore.     

OmpA: a pore or not a pore? Simulation and modeling studies

The bacterial outer membrane protein OmpA is composed of an N-terminal 171-residue beta-barrel domain (OmpA(171)) that spans the bilayer and a periplasmic, C-terminal domain of unknown structure. OmpA has been suggested to primarily serve a structural role, as no continuous pore through the center of the barrel can be discerned in the crystal structure of OmpA(171). However, several groups have recorded ionic conductances for bilayer-reconstituted OmpA(171). To resolve this apparent paradox we have used molecular dynamics (MD) simulations on OmpA(171) to explore the conformational dynamics of the protein, in particular the possibility of transient formation of a central pore. A total of 19 ns of MD simulations of OmpA(171) have been run, and the results were analyzed in terms of 1) comparative behavior of OmpA(171) in different bilayer and bilayer-mimetic environments, 2) solvation states of OmpA(171), and 3) pore characteristics in different MD simulations. Significant mobility was observed for residues and water molecules within the beta-barrel. A simulation in which putative gate region side chains of the barrel interior were held in a non-native conformation led to an open pore, with a predicted conductance similar to experimental measurements. The OmpA(171) pore has been shown to be somewhat more dynamic than suggested by the crystal structure. A gating mechanism is proposed to explain its documented channel properties, involving a flickering isomerization of Arg138, forming alternate salt bridges with Glu52 (closed state) and Glu128 (open state).