Transient pores in hemifusion diaphragms.

Exchange of material across two membranes, as in the case of synaptic neurotransmitter release from a vesicle, involves the formation and poration of a hemifusion diaphragm (HD). The nontrivial geometry of the HD leads to environment-dependent control, regarding the stability and dynamics of the pores required for this kind of exocytosis. This work combines particle simulations, field-based calculations, and phenomenological modeling to explore the factors influencing the stability, dynamics, and possible control mechanisms of pores in HDs. We find that pores preferentially form at the HD rim, and that their stability is sensitive to a number of factors, including the three line tensions, membrane tension, HD size, and the ability of lipids to "flip-flop" across leaflets. Along with a detailed analysis of these factors, we discuss ways that vesicles or cells may use them to open and close pores and thereby quickly and efficiently transport material.

Thermodynamically reversible paths of the first fusion intermediate reveal an important role for membrane anchors of fusion proteins.

Biological membrane fusion proceeds via an essential topological transition of the two membranes involved. Known players such as certain lipid species and fusion proteins are generally believed to alter the free energy and thus the rate of the fusion reaction. Quantifying these effects by theory poses a major challenge since the essential reaction intermediates are collective, diffusive and of a molecular length scale. We conducted molecular dynamics simulations in conjunction with a state-of-the-art string method to resolve the minimum free-energy path of the first fusion intermediate state, the so-called stalk. We demonstrate that the isolated transmembrane domains (TMDs) of fusion proteins such as SNARE molecules drastically lower the free energy of both the stalk barrier and metastable stalk, which is not trivially explained by molecular shape arguments. We relate this effect to the local thinning of the membrane (negative hydrophobic mismatch) imposed by the TMDs which favors the nearby presence of the highly bent stalk structure or prestalk dimple. The distance between the membranes is the most crucial determinant of the free energy of the stalk, whereas the free-energy barrier changes only slightly. Surprisingly, fusion enhancing lipids, i.e., lipids with a negative spontaneous curvature, such as PE lipids have little effect on the free energy of the stalk barrier, likely because of its single molecular nature. In contrast, the lipid shape plays a crucial role in overcoming the hydration repulsion between two membranes and thus rather lowers the total work required to form a stalk.

Mechanics of membrane fusion/pore formation.

Lipid bilayers play a fundamental role in many biological processes, and a considerable effort has been invested in understanding their behavior and the mechanism of topological changes like fusion and pore formation. Due to the time- and length-scale on which these processes occur, computational methods have proven to be an especially useful tool in their study. With their help, a number of interesting findings about the shape of fusion intermediates could be obtained, and novel hypotheses about the mechanism of topological changes and the involvement of peptides therein were suggested. In this work, we try to present a summary of these developments together with some hitherto unpublished results, featuring, among others, the shape of stalks and fusion pores, possible modes of action of the influenza HA fusion peptide and the SNARE protein complex, the mechanism of supported lipid bilayer formation by vesicle spreading, and the free energy and transition pathway of the fusion process.

Transition path from two apposed membranes to a stalk obtained by a combination of particle simulations and string method.

The formation of an hourglass-shaped passage (stalk) connecting two apposed membranes is an essential initial step in membrane fusion. The most probable transition path from two separate membranes to a stalk, i.e., the minimum free-energy path (MFEP), is constructed using a combination of particle simulations and string method. For the reversible transition path in the coarse-grained membrane model, a collective order parameter, m, can be identified as the local difference of hydrophilic and hydrophobic densities. In particle simulations, the free energy F[m] as a functional of m is not readily available. This difficulty is overcome by an equation-free approach, where the morphology and the excess free energy along the MFEP are obtained by an on-the-fly string method. The transition state is confirmed by diagonalization of order-parameter fluctuations and by the probability of reaching either stalk or bilayer morphology from different positions along the MFEP.

Line-tension controlled mechanism for influenza fusion.

Our molecular simulations reveal that wild-type influenza fusion peptides are able to stabilize a highly fusogenic pre-fusion structure, i.e. a peptide bundle formed by four or more trans-membrane arranged fusion peptides. We rationalize that the lipid rim around such bundle has a non-vanishing rim energy (line-tension), which is essential to (i) stabilize the initial contact point between the fusing bilayers, i.e. the stalk, and (ii) drive its subsequent evolution. Such line-tension controlled fusion event does not proceed along the hypothesized standard stalk-hemifusion pathway. In modeled influenza fusion, single point mutations in the influenza fusion peptide either completely inhibit fusion (mutants G1V and W14A) or, intriguingly, specifically arrest fusion at a hemifusion state (mutant G1S). Our simulations demonstrate that, within a line-tension controlled fusion mechanism, these known point mutations either completely inhibit fusion by impairing the peptide's ability to stabilize the required peptide bundle (G1V and W14A) or stabilize a persistent bundle that leads to a kinetically trapped hemifusion state (G1S). In addition, our results further suggest that the recently discovered leaky fusion mutant G13A, which is known to facilitate a pronounced leakage of the target membrane prior to lipid mixing, reduces the membrane integrity by forming a 'super' bundle. Our simulations offer a new interpretation for a number of experimentally observed features of the fusion reaction mediated by the prototypical fusion protein, influenza hemagglutinin, and might bring new insights into mechanisms of other viral fusion reactions.

Solvent-exposed tails as prestalk transition states for membrane fusion at low hydration.

Membrane fusion is a key step in intracellular trafficking and viral infection. The underlying molecular mechanism is poorly understood. We have used molecular dynamics simulations in conjunction with a coarse grained model to study early metastable and transition states during the fusion of two planar palmitoyl-oleoyl-phosphatidylcholine (POPC) bilayers separated by five waters per lipid in the cis leaflets at zero tension. This system mimics the contact area between two vesicles with large diameters compared to the membrane thickness at conditions where fusion may start in the core of the contact area. At elevated temperatures, the two proximal leaflets become connected via multiple lipid molecules and form a stalklike structure. At room temperature, this structure has a free energy of 3k(B)T and is separated from the unconnected state by a significant free energy barrier of 20k(B)T. Stalk formation is initiated by the establishment of a localized hydrophobic contact between the bilayers. This contact is either formed by two partially splayed lipids or a single fully splayed one leading to the formation of a (metastable) splayed lipid bond intermediate. These findings indicate that, for low hydration, early membrane fusion kinetics is not determined by the stalk energy but by the energy of prestalk transition states involving solvent-exposed lipid tails.

Nonmonotonic incommensurability effects in lamellar-in-lamellar self-assembled multiblock copolymers.

Using the self-consistent-field theory numerical procedure we find that the period D of the lamellar-in-lamellar morphology formed in symmetric multiblock copolymer melts A(mN/2)(B(N/2)A(N/2))(n)B(mN/2) at intermediate segregations changes nonmonotonically with an increase in the relative tail length m. Therewith D reveals, as a function of the Flory chi-parameter, a drastic change in the vicinity of the internal structure formation, which can be both a drop and a rise, depending on the value of m. It is argued that the unusual behavior found is a particular case of a rather general effect of the incommensurability between the two length scales that characterize the system under consideration.

Nonconventional morphologies in two-length scale block copolymer systems beyond the weak segregation theory.

The order-disorder and order-order transitions (ODT and OOT) in the linear multiblock copolymers with two-length scale architecture A(fmN)(B(N2)A(N2))(n)B((1-f)mN) are studied under intermediate cooling below the ODT critical point where a nonconventional sequence of the OOTs was predicted previously [Smirnova et al., J. Chem. Phys. 124, 054907 (2006)] within the weak segregation theory (WST). To describe the ordered morphologies appearing in block copolymers (BCs) under cooling, we use the pseudospectral version of the self-consistent field theory (SCFT) with some modifications providing a good convergence speed and a high precision of the solution due to using the Ng iterations [J. Chem. Phys. 61, 2680 (1974)] and a reasonable choice of the predefined symmetries of the computation cell as well as initial guess for the iterations. The WST predicted sequence of the phase transitions is found to hold if the tails of the BCs under consideration are symmetric enough (mid R:0.5-fmid R:0.05, a large region of the face-centered cubic phase stability is found (up to our knowledge, first within the SCFT framework) inside of the body-centered cubic phase stability region. Occurrence of the two-dimensional and three-dimensional phases with the micelles formed, unlike the conventional diblock copolymers, by the longer (rather than shorter) tails, and its relationship to the BC architecture is first described in detail. The calculated spectra of the ordered phases show that nonmonotonous temperature dependence of the secondary peak scattering intensities accompanied by their vanishing and reappearance is rather a rule than exception.

Microphase separation in multiblock copolymer melts: Nonconventional morphologies and two-length-scale switching.

The phase behavior of AfmN(BN2AN2)B(1-fmN) multiblock copolymer melts is studied within the weak segregation theory. The interplay between ordering on different length scales is shown to cause dramatic changes both in the ordered phase symmetry and periodicity upon small variation of the architectural parameters of the macromolecules. Phase diagrams are presented in the (f,chiN) plane (chi is the Flory-Huggins parameter) for various values of the architecture parameters n and m. Near the critical surface, i.e., for (f-0.5)2<<1, such nonconventional cubic phases as the face-centered cubic (FCC), simple cubic (SC), (double) gyroid, and the so-called BCC(2) (single gyroid) are found to be stable. The lamellar morphology is shown to be replaced by BCC2, FCC, or SC (depending on the structural parameters) as the most stable low-temperature phase.