Simulated dynamic cholesterol redistribution favors membrane fusion pore constriction.

Endo- and exocytosis proceed through a highly strained membrane fusion pore topology regardless of the aiding protein machinery. The membrane's lipid components bias fusion pores toward expansion or closure, modifying the necessary work done by proteins. Cholesterol, a key component of plasma membranes, promotes both inverted lipid phases with concave leaflets (i.e., negative total curvature, which thins the leaflet) and flat bilayer phases with thick, ordered hydrophobic interiors. We demonstrate by theory and simulation that both leaflets of nascent catenoidal fusion pores have negative total curvature. Furthermore, the hydrophobic core of bilayers with strong negative Gaussian curvature is thinned. Therefore, it is an open question whether cholesterol will be enriched in these regions because of the negative total curvature or depleted because of the membrane thinning. Here, we compare all-atom molecular dynamics simulations (built using a procedure to create specific fusion pore geometries) and theory to understand the underlying reasons for lipid redistribution on fusion pores. Our all-atom molecular dynamics simulations resolve this question by showing that cholesterol is strongly excluded from the thinned neck of fusion and fission pores, revealing that thickness (and/or lipid order) influences cholesterol distributions more than curvature. The results imply that cholesterol exclusion can drive fusion pore closure by creating a small, cholesterol-depleted zone in the neck. This model agrees with literature evidence that membrane reshaping is connected to cholesterol-dependent lateral phase separation.

Molecular mechanisms of spontaneous curvature and softening in complex lipid bilayer mixtures.

Membrane reshaping is an essential biological process. The chemical composition of lipid membranes determines their mechanical properties and thus the energetics of their shape. Hundreds of distinct lipid species make up native bilayers, and this diversity complicates efforts to uncover what compositional factors drive membrane stability in cells. Simplifying assumptions, therefore, are used to generate quantitative predictions of bilayer dynamics based on lipid composition. One assumption commonly used is that "per lipid" mechanical properties are both additive and constant-that they are an intrinsic property of lipids independent of the surrounding composition. Related to this is the assumption that lipid bulkiness, or "shape," determines its curvature preference, independently of context. In this study, all-atom molecular dynamics simulations on three separate multilipid systems were used to explicitly test these assumptions, applying methodology recently developed to isolate properties of single lipids or nanometer-scale patches of lipids. The curvature preference experienced by populations of lipid conformations were inferred from their redistribution on a dynamically fluctuating bilayer. Representative populations were extracted by both structural similarity and semi-automated hidden Markov model analysis. The curvature preferences of lipid dimers were then determined and compared with an additive model that combines the monomer curvature preference of both the individual lipids. In all three systems, we identified conformational subpopulations of lipid dimers that showed non-additive curvature preference, in each case mediated by a special chemical interaction (e.g., hydrogen bonding). Our study highlights the importance of specific chemical interactions between lipids in multicomponent bilayers and the impact of interactions on bilayer stiffness. We identify two mechanisms of bilayer softening: diffusional softening, driven by the dynamic coupling between lipid distributions and membrane undulations, and conformational softening, driven by the inter-conversion between distinct dimeric conformations.

Biophysical and functional characterization of Norrin signaling through Frizzled4.

Wnt signaling is initiated by Wnt ligand binding to the extracellular ligand binding domain, called the cysteine-rich domain (CRD), of a Frizzled (Fzd) receptor. Norrin, an atypical Fzd ligand, specifically interacts with Fzd4 to activate ß-catenin-dependent canonical Wnt signaling. Much of the molecular basis that confers Norrin selectivity in binding to Fzd4 was revealed through the structural study of the Fzd4CRD-Norrin complex. However, how the ligand interaction, seemingly localized at the CRD, is transmitted across full-length Fzd4 to the cytoplasm remains largely unknown. Here, we show that a flexible linker domain, which connects the CRD to the transmembrane domain, plays an important role in Norrin signaling. The linker domain directly contributes to the high-affinity interaction between Fzd4 and Norrin as shown by ∼10-fold higher binding affinity of Fzd4CRD to Norrin in the presence of the linker. Swapping the Fzd4 linker with the Fzd5 linker resulted in the loss of Norrin signaling, suggesting the importance of the linker in ligand-specific cellular response. In addition, structural dynamics of Fzd4 associated with Norrin binding investigated by hydrogen/deuterium exchange MS revealed Norrin-induced conformational changes on the linker domain and the intracellular loop 3 (ICL3) region of Fzd4. Cell-based functional assays showed that linker deletion, L430A and L433A mutations at ICL3, and C-terminal tail truncation displayed reduced ß-catenin-dependent signaling activity, indicating the functional significance of these sites. Together, our results provide functional and biochemical dissection of Fzd4 in Norrin signaling.

Gramicidin A Channel Formation Induces Local Lipid Redistribution II: A 3D Continuum Elastic Model.

To change conformation, a protein must deform the surrounding bilayer. In this work, a three-dimensional continuum elastic model for gramicidin A in a lipid bilayer is shown to describe the sensitivity to thickness, curvature stress, and the mechanical properties of the lipid bilayer. A method is demonstrated to extract the gramicidin-lipid boundary condition from all-atom simulations that can be used in the three-dimensional continuum model. The boundary condition affects the deformation dramatically, potentially much more than typical variations in the material stiffness do as lipid composition is changed. Moreover, it directly controls the sensitivity to curvature stress. The curvature stress and hydrophobic surfaces of the all-atom and continuum models are found to be in excellent agreement. The continuum model is applied to estimate the enrichment of hydrophobically matched lipids near the channel in a mixture, and the results agree with single-channel experiments and extended molecular dynamics simulations from the companion article by Beaven et al. in this issue of Biophysical Journal.

Gramicidin A Channel Formation Induces Local Lipid Redistribution I: Experiment and Simulation.

Integral membrane protein function can be modulated by the host bilayer. Because biological membranes are diverse and nonuniform, we explore the consequences of lipid diversity using gramicidin A channels embedded in phosphatidylcholine (PC) bilayers composed of equimolar mixtures of di-oleoyl-PC and di-erucoyl-PC (dC18:1+dC22:1, respectively), di-palmitoleoyl-PC and di-nervonoyl-PC (dC16:1+dC24:1, respectively), and di-eicosenoyl-PC (pure dC20:1), all of which have the same average bilayer chain length. Single-channel lifetime experiments, molecular dynamics simulations, and a simple lipid compression model are used in tandem to gain insight into lipid redistribution around the channel, which partially alleviates the bilayer deformation energy associated with channel formation. The average single-channel lifetimes in the two-component bilayers (95 ± 10 ms for dC18:1+dC22:1 and 195 ± 20 ms for dC16:1+dC24:1) were increased relative to the single-component dC20:1 control bilayer (65 ± 10 ms), implying lipid redistribution. Using a theoretical treatment of thickness-dependent changes in channel lifetimes, the effective local enrichment of lipids around the channel was estimated to be 58 ± 4% dC18:1 and 66 ± 2% dC16:1 in the dC18:1+dC22:1 and dC16:1+dC24:1 bilayers, respectively. 3.5-µs molecular dynamics simulations show 66 ± 2% dC16:1 in the first lipid shell around the channel in the dC16:1+dC24:1 bilayer, but no significant redistribution (50 ± 4% dC18:1) in the dC18:1+dC22:1 bilayer; these simulated values are within the 95% confidence intervals of the experimental averages. The strong preference for the better matching lipid (dC16:1) near the channel in the dC16:1+dC24:1 mixture and lesser redistribution in the dC18:1+dC22:1 mixture can be explained by the energetic cost associated with compressing the lipids to match the channel's hydrophobic length.

An extensive simulation study of lipid bilayer properties with different head groups, acyl chain lengths, and chain saturations.

Previous MD simulations of six phosphocholine (PC) lipid bilayers demonstrated the accuracy of the CHARMM36 force field (C36FF) for PC bilayer simulation at varied temperatures (BBA-Biomembranes, 1838 (2014): 2520-2529). In this work, we further examine the accuracy of C36FF over a wide temperature range for a broader range of lipid types such as various head groups (phosphatidic acid (PA), PC, phosphoethanolamine (PE), phosphoglycerol (PG), and phosphoserine (PS)), and tails (saturated, mono-, mixed- and poly-unsaturated acyl chains with varied chain lengths). The structural properties (surface area per lipid (SA/lip), overall bilayer thickness, hydrophobic thickness, headgroup-to-headgroup thickness, deuterium order parameter (SCD), and spin-lattice relaxation time (T1)) obtained from simulations agree well with nearly all available experimental data. Our analyses indicate that PS lipids have the most inter-lipid hydrogen bonds, while PG lipids have the most intra-lipid hydrogen bonds, which play the main role in their low SA/lip in PS lipids and low thicknesses in PG lipids, respectively. PS, PE, and PA lipids have the largest contact clusters with on average 5-8 lipids per cluster, while PC and PG have clusters of 4 lipids based on a cutoff distance of 6.5Å. PS lipids have much slower lipid wobble (i.e., higher correlation time) than other head groups at a given temperature as the hydrogen bonded network significantly reduces a lipid's mobility, and the rate of lipid wobble increases dramatically as temperature increases. These in-depth analyses facilitate further understanding of lipid bilayers at the atomic level.

CHARMM-GUI 10 years for biomolecular modeling and simulation.

CHARMM-GUI, http://www.charmm-gui.org, is a web-based graphical user interface that prepares complex biomolecular systems for molecular simulations. CHARMM-GUI creates input files for a number of programs including CHARMM, NAMD, GROMACS, AMBER, GENESIS, LAMMPS, Desmond, OpenMM, and CHARMM/OpenMM. Since its original development in 2006, CHARMM-GUI has been widely adopted for various purposes and now contains a number of different modules designed to set up a broad range of simulations: (1) PDB Reader & Manipulator, Glycan Reader, and Ligand Reader & Modeler for reading and modifying molecules; (2) Quick MD Simulator, Membrane Builder, Nanodisc Builder, HMMM Builder, Monolayer Builder, Micelle Builder, and Hex Phase Builder for building all-atom simulation systems in various environments; (3) PACE CG Builder and Martini Maker for building coarse-grained simulation systems; (4) DEER Facilitator and MDFF/xMDFF Utilizer for experimentally guided simulations; (5) Implicit Solvent Modeler, PBEQ-Solver, and GCMC/BD Ion Simulator for implicit solvent related calculations; (6) Ligand Binder for ligand solvation and binding free energy simulations; and (7) Drude Prepper for preparation of simulations with the CHARMM Drude polarizable force field. Recently, new modules have been integrated into CHARMM-GUI, such as Glycolipid Modeler for generation of various glycolipid structures, and LPS Modeler for generation of lipopolysaccharide structures from various Gram-negative bacteria. These new features together with existing modules are expected to facilitate advanced molecular modeling and simulation thereby leading to an improved understanding of the structure and dynamics of complex biomolecular systems. Here, we briefly review these capabilities and discuss potential future directions in the CHARMM-GUI development project. © 2016 Wiley Periodicals, Inc.

Softening in Two-Component Lipid Mixtures by Spontaneous Curvature Variance.

The bending modulus of a lipid bilayer quantifies its mechanical resistance to curvature. It is typically understood in terms of thickness; e.g., thicker bilayers are usually stiffer. Here, we describe an additional and powerful molecular determinant of stiffness─the variance in the distribution of curvature sensitivity of lipids and lipid conformations. Zwitterionic choline and ethanolamine headgroups of glycerophospholipids dynamically explore inter- and intraspecies interactions, leading to transient clustering. We demonstrate that these clusters couple strongly to negative curvature, exciting undulatory membrane modes and reducing the apparent bending modulus. Three force fields (Martini 2, Martini 3, and all-atom CHARMM C36) each show the effect to a different extent, with the coarse-grained Martini models showing the most clustering and thus the most softening. The theory is a guide to understanding the stiffness of biological membranes with their complex composition, as well as how choices of force field parameterization are translated into mechanical stiffness.

Biomimetic Vesicles with Designer Phospholipids Can Sense Environmental Redox Cues.

Cell-like materials that sense environmental cues can serve as next-generation biosensors and help advance the understanding of intercellular communication. Currently, bottom-up engineering of protocell models from molecular building blocks remains a grand challenge chemists face. Herein, we describe giant unilamellar vesicles (GUVs) with biomimetic lipid membranes capable of sensing environmental redox cues. The GUVs employ activity-based sensing through designer phospholipids that are fluorescently activated in response to specific reductive (hydrogen sulfide) or oxidative (hydrogen peroxide) conditions. These synthetic phospholipids are derived from 1,2-dipalmitoyl-rac-glycero-3-phosphocholine and they possess a headgroup with heterocyclic aromatic motifs. Despite their structural deviation from the phosphocholine headgroup, the designer phospholipids (0.5-1.0 mol %) mixed with natural lipids can vesiculate, and the resulting GUVs (7-20 µm in diameter) remain intact over the course of redox sensing. All-atom molecular dynamics simulations gave insight into how these lipids are positioned within the hydrophobic core of the membrane bilayer and at the membrane-water interface. This work provides a purely chemical method to investigate potential redox signaling and opens up new design opportunities for soft materials that mimic protocells.

Synaptotagmin 7 C2 domains induce membrane curvature stress via electrostatic interactions and the wedge mechanism.

Synaptotagmin 7 (Syt-7) is part of the synaptotagmin protein family that regulates exocytotic lipid membrane fusion. Among the family, Syt-7 stands out by its membrane binding strength and stabilization of long-lived membrane fusion pores. Given that Syt-7 vesicles form long-lived fusion pores, we hypothesize that its interactions with the membrane stabilize the specific curvatures, thicknesses, and lipid compositions that support a metastable fusion pore. Using all-atom molecular dynamics simulations and FRET-based assays of Syt-7's membrane-binding C2 domains (C2A and C2B), we found that Syt-7 C2 domains sequester anionic lipids, are sensitive to cholesterol, thin membranes, and generate lipid membrane curvature by two competing, but related mechanisms. First, Syt-7 forms strong electrostatic contacts with the membrane, generating negative curvature stress. Second, Syt-7's calcium binding loops embed in the membrane surface, acting as a wedge to thin the membrane and induce positive curvature stress. These curvature mechanisms are linked by the protein insertion depth as well as the resulting protein tilt. Simplified quantitative models of the curvature-generating mechanisms link simulation observables to their membrane-reshaping effectiveness.

Softening in two-component lipid mixtures by spontaneous curvature variance.

The bending modulus of a lipid bilayer quantifies its mechanical resistance to curvature. It is typically understood in terms of thickness, e.g., thicker bilayers are stiffer. Here, we describe an additional and powerful molecular determinant of stiffness - the variance in the distribution of curvature sensitivity of lipids and lipid conformations. Zwitterionic choline and ethanolamine head-groups of glycero-phospholipids dynamically explore inter- and intra-species interactions, leading to transient clustering. We demonstrate that these clusters couple strongly to negative curvature, exciting undulatory membrane modes and reducing the apparent bending modulus. Three forcefields (Martini 2, Martini 3, and all-atom CHARMM C36) each show the effect to a different extent, with the coarse-grained Martini models showing the most clustering and thus the most softening. The theory is a guide to understanding the stiffness of biological membranes with their complex composition, as well as how choices of forcefield parameterization are translated into mechanical stiffness.

Cholesterol Binds the Amphipathic Helix of IFITM3 and Regulates Antiviral Activity.

The interferon-induced transmembrane (IFITM) proteins broadly inhibit the entry of diverse pathogenic viruses, including Influenza A virus (IAV), Zika virus, HIV-1, and SARS coronaviruses by inhibiting virus-cell membrane fusion. IFITM3 was previously shown to disrupt cholesterol trafficking, but the functional relationship between IFITM3 and cholesterol remains unclear. We previously showed that inhibition of IAV entry by IFITM3 is associated with its ability to promote cellular membrane rigidity, and these activities are functionally linked by a shared requirement for the amphipathic helix (AH) found in the intramembrane domain (IMD) of IFITM3. Furthermore, it has been shown that the AH of IFITM3 alters lipid membranes in vitro in a cholesterol-dependent manner. Therefore, we aimed to elucidate the relationship between IFITM3 and cholesterol in more detail. Using a fluorescence-based in vitro binding assay, we found that a peptide derived from the AH of IFITM3 directly interacted with the cholesterol analog, NBD-cholesterol, while other regions of the IFITM3 IMD did not, and native cholesterol competed with this interaction. In addition, recombinant full-length IFITM3 protein also exhibited NBD-cholesterol binding activity. Importantly, previously characterized mutations within the AH of IFITM3 that strongly inhibit antiviral function (F63Q and F67Q) disrupted AH structure in solution, inhibited cholesterol binding in vitro, and restricted bilayer insertion in silico. Our data suggest that direct interactions with cholesterol may contribute to the inhibition of membrane fusion pore formation by IFITM3. These findings may facilitate the design of therapeutic peptides for use in broad-spectrum antiviral therapy.

Spatial extent of a single lipid's influence on bilayer mechanics.

To what spatial extent does a single lipid affect the mechanical properties of the membrane that surrounds it? The lipid composition of a membrane determines its mechanical properties. The shapes available to the membrane depend on its compositional material properties, and therefore, the lipid environment. Because each individual lipid species' chemistry is different, it is important to know its range of influence on membrane mechanical properties. This is defined herein as the lipid's mechanical extent. Here, a lipid's mechanical extent is determined by quantifying lipid redistribution and the average curvature that lipid species experience on fluctuating membrane surfaces. A surprising finding is that, unlike unsaturated lipids, saturated lipids have a complicated, nonlocal effect on the surrounding surface, with the interaction strength maximal at a finite length-scale. The methodology provides the means to substantially enrich curvature-energy models of membrane structures, quantifying what was previously only conjecture.

Curvature Energetics Determined by Alchemical Simulation on Four Topologically Distinct Lipid Phases.

The relative curvature energetics of two lipids are tested using thermodynamic integration (TI) on four topologically distinct lipid phases. Simulations use TI to switch between choline headgroup lipids (POPC; that prefers to be flat) and ethanolamine headgroup lipids (POPE; that prefer, for example, the inner monolayer of vesicles). The thermodynamical moving of the lipids between planar, inverse hexagonal (HII), cubic (QII; Pn3m space group), and vesicle topologies reveals differences in material parameters that were previously challenging to access. The methodology allows for predictions of two important lipid material properties: the difference in POPC/POPE monolayer intrinsic curvature (ΔJ0) and the difference in POPC/POPE monolayer Gaussian curvature modulus (Δκ̅m), both of which are connected to the energetics of topological variation. Analysis of the TI data indicates that, consistent with previous experiment and simulation, the J0 of POPE is more negative than POPC (ΔJ0 = -0.018 ± 0.001 Å-1). The theoretical framework extracts significant differences in κ̅m of which POPE is less negative than POPC by 2.0 to 4.0 kcal/mol. The range of these values is determined by considering subsets of the simulations, and disagreement between these subsets suggests separate mechanical parameters at very high curvature. Finally, the fit of the TI data to the model indicates that the position of the pivotal plane of curvature is not constant across topologies at high curvature. Overall, the results offer insights into lipid material properties, the limits of a single HC model, and how to test them using simulation.

A combined molecular/continuum-modeling approach to predict the small-angle neutron scattering of curved membranes.

This paper develops a framework to compute the small-angle neutron scattering (SANS) from highly curved, dynamically fluctuating, and potentially inhomogeneous membranes. This method is needed to compute the scattering from nanometer-scale membrane domains that couple to curvature, as predicted by molecular modeling. The detailed neutron scattering length density of a small planar bilayer patch is readily available via molecular dynamics simulation. A mathematical, mechanical transformation of the planar scattering length density is developed to predict the scattering from curved bilayers. By simulating a fluctuating, curved, surface-continuum model, long time- and length-scales can be reached while, with the aid of the planar-to-curved transformation, the molecular features of the scattering length density can be retained. A test case for the method is developed by constructing a coarse-grained lipid vesicle following a protocol designed to relieve both the osmotic stress inside the vesicle and the lipid-number stress between the leaflets. A question was whether the hybrid model would be able to replicate the scattering from the highly deformed inner and outer leaflets of the small vesicle. Matching the scattering of the full (molecular vesicle) and hybrid (continuum vesicle) models indicated that the inner and outer leaflets of the full vesicle were expanded laterally, consistent with previous simulations of the Martini forcefield that showed thinning in small vesicles. The vesicle structure is inconsistent with a zero-tension leaflet deformed by a single set of elastic parameters, and the results show that this is evident in the scattering. The method can be applied to translate observations of any molecular model's neutron scattering length densities from small patches to large length and timescales.

Characterizing Residue-Bilayer Interactions Using Gramicidin A as a Scaffold and Tryptophan Substitutions as Probes.

Previous experiments have shown that the lifetime of a gramicidin A dimer channel (which forms from two nonconducting monomers) in a lipid bilayer is modulated by mutations of the tryptophan (Trp) residues at the bilayer-water interface. We explore this further using extensive molecular dynamics simulations of various gA dimer and monomer mutants at the Trp positions in phosphatidylcholine bilayers with different tail lengths. gA interactions with the surrounding bilayer are strongly modulated by mutating these Trp residues. There are three principal effects: eliminating residue hydrogen bonding ability (i.e., reducing the channel-monolayer coupling strength) reduces the extent of the bilayer deformation caused by the assembled dimeric channel; a residue's size and geometry affects its orientation, leading to different hydrogen bonding partners; and increasing a residue's hydrophobicity increases the depth of gA monomer insertion relative to the bilayer center, thereby increasing the lipid bending frustration.

How Tolerant are Membrane Simulations with Mismatch in Area per Lipid between Leaflets?

Difficulties in estimating the correct number of lipids in each leaflet of complex bilayer membrane simulation systems make it inevitable to introduce a mismatch in lipid packing (i.e., area per lipid) and thus alter the lateral pressure of each leaflet. To investigate potential impacts of such mismatch on simulation results, we performed molecular dynamics simulations of saturated and monounsaturated lipid bilayers with and without gramicidin A or WALP23 at various mismatches by adjusting the number of lipids in the lower leaflet from no mismatch to a 25% reduction compared to that in the upper leaflet. All simulations were stable under the constant pressure barostat, but the mismatch induces asymmetric lipid packing between the leaflets, so that the upper leaflet becomes more ordered, and the lower leaflet becomes less ordered. The mismatch impacts on various bilayer properties are mild up to 5-10% mismatch, and bilayers with fully saturated chains appear to be more prone to these impacts than those with unsaturated tails. The nonvanishing leaflet surface tensions and the free energy derivatives with respect to the bilayer curvature indicate that the bilayer would be energetically unstable in the presence of mismatch. We propose a quantitative criterion for allowable mismatch based on the energetics derived from a continuum elastic model, which grows as a square root of the number of the lipids in the system. On the basis of this criterion, we infer that the area per lipid mismatch up to 5% would be tolerable in various membrane simulations of reasonable all-atom system sizes (40-160 lipids per leaflet).