A continuum membrane model can predict curvature sensing by helix insertion

Charge distribution and helicity tune the binding of septin's amphipathic helix domain to membranes

Amphipathic helices (AHs) are secondary structures that can facilitate binding of proteins to the membrane by folding into a helix with hydrophobic and hydrophilic faces that interact with the same surfaces in the lipid membrane. Septins are cytoskeletal proteins that preferentially bind to domains of micron-scale curvature on the cell membrane. Studies have shown that AH domains in septin are essential for curvature sensing. We present the first computational study of septin AH interactions with lipid bilayers. Using all-atom simulations and metadynamics-enhanced sampling, we study the effect of charge distribution at the flanking ends of septin AH on the energy for helical folding and its consequences on the binding configuration and affinity to the membrane. This is relevant to septins, since the net positive charge on the flanking C-terminal amino acids is a conserved property across several organisms. Simulations revealed that the energy barrier for folding in the neutral-capped AH is much larger than the charge-capped AH, leading to a small fraction of AH folding and integration to the membrane compared to a significantly folded configuration in the bound charge-capped AH. These observations are consistent with the binding measurements of synthetic AH constructs with variable helicity to lipid vesicles. Additionally, we examined an extended AH sequence including eight amino acids upstream and downstream of the AH to mimic the native protein. Again, simulations and experiments show that the extended peptide, with a net positive charge at C-terminus, adopts a strong helical configuration in solution, giving rise to a higher membrane affinity. Altogether, these results identify the energy cost for folding of AHs as a regulator of AH binding configuration and affinity and provide a basic template for parameterizing AH-membrane interactions as a starting point for the future multiscale simulations for septin-membrane interactions.

Membrane Remodeling Driven by Shallow Helix Insertions via a Cooperative Mechanism

Helix-membrane interactions are key to membrane deformation and play significant biological roles. However, systematic studies on the mechanisms behind these interactions are limited. This study uses a continuum membrane model to investigate how shallowly inserted helices interact with biological membranes, focusing on membrane deformation and the cooperative effects of multiple helices. Our findings show that even short helices (2 nm in length) can induce anisotropic membrane deformation. Longer helices and deeper insertions result in more significant deformations, and the spatial arrangement of helices affects the nature of these deformations. The perturbation area (PA) and perturbation extent (PE) are quantified to describe membrane deformation, revealing stronger cooperative effects in parallel insertions and more complex deformations in other arrangements. Additionally, membrane properties, such as lipid composition, influence the extent of deformation. In multi-helix systems, we observe local clustering behavior when perturbations are strong enough, with cooperativity varying based on helix length, insertion depth, and membrane composition. This study provides criteria for helix cooperativity, advancing our understanding of helix-membrane interactions and their biological significance in processes like membrane remodeling.

Predicting protein curvature sensing across membrane compositions with a bilayer continuum model

Cytoplasmic proteins must recruit to membranes to function in processes such as endocytosis and cell division. Many of these proteins recognize not only the chemical structure of the membrane lipids, but the curvature of the surface, binding more strongly to more highly curved surfaces, or 'curvature sensing'. Curvature sensing by amphipathic helices is known to vary with membrane bending rigidity, but changes to lipid composition can simultaneously alter membrane thickness, spontaneous curvature, and leaflet symmetry, thus far preventing a systematic characterization of lipid composition on such curvature sensing through either experiment or simulation. Here we develop and apply a bilayer continuum membrane model that can tractably address this gap, quantifying how controlled changes to each material property can favor or disfavor protein curvature sensing. We evaluate both energetic and structural changes to vesicles upon helix insertion, with strong agreement to new in vitro experiments and all-atom MD simulations, respectively. Our membrane model builds on previous work to include both monolayers of the bilayer via representation by continuous triangular meshes. We introduce a coupling energy that captures the incompressibility of the membrane and the established energetics of lipid tilt. In agreement with experiment, our model predicts stronger curvature sensing in membranes with distinct tail groups (POPC vs DOPC vs DLPC), despite having identical head-group chemistry; the model shows that the primary driving force for weaker curvature sensing in DLPC is that it is thinner, and more wedge shaped. Somewhat surprisingly, asymmetry in lipid shape composition between the two leaflets has a negligible contribution to membrane mechanics following insertion. Our multi-scale approach can be used to quantitatively and efficiently predict how changes to membrane composition in flat to highly curved surfaces alter membrane energetics driven by proteins, a mechanism that helps proteins target membranes at the correct time and place.

Membrane curvature sensing and symmetry breaking of the M2 proton channel from Influenza A

The M2 proton channel aids in the exit of mature influenza viral particles from the host plasma membrane through its ability to stabilize regions of high negative Gaussian curvature (NGC) that occur at the neck of budding virions. The channels are homo-tetramers that contain a cytoplasm-facing amphipathic helix (AH) that is necessary and sufficient for NGC generation; however, constructs containing the transmembrane spanning helix, which facilitates tetramerization, exhibit enhanced curvature generation. Here, we used all-atom molecular dynamics (MD) simulations to explore the conformational dynamics of M2 channels in lipid bilayers revealing that the AH is dynamic, quickly breaking the fourfold symmetry observed in most structures. Next, we carried out MD simulations with the protein restrained in four- and twofold symmetric conformations to determine the impact on the membrane shape. While each pattern was distinct, all configurations induced pronounced curvature in the outer leaflet, while conversely, the inner leaflets showed minimal curvature and significant lipid tilt around the AHs. The MD-generated profiles at the protein-membrane interface were then extracted and used as boundary conditions in a continuum elastic membrane model to calculate the membrane-bending energy of each conformation embedded in different membrane surfaces characteristic of a budding virus. The calculations show that all three M2 conformations are stabilized in inward-budding, concave spherical caps and destabilized in outward-budding, convex spherical caps, the latter reminiscent of a budding virus. One of the C2-broken symmetry conformations is stabilized by 4 kT in NGC surfaces with the minimum energy conformation occurring at a curvature corresponding to 33 nm radii. In total, our work provides atomistic insight into the curvature sensing capabilities of M2 channels and how enrichment in the nascent viral particle depends on protein shape and membrane geometry.

Physical limits to membrane curvature sensing by a single protein

Membrane curvature sensing is essential for a diverse range of biological processes. Recent experiments have revealed that a single nanometer-sized septin protein has different binding rates to membrane-coated glass beads of 1-µm and 3-µm diameters, even though the septin is orders of magnitude smaller than the beads. This sensing ability is especially surprising since curvature-sensing proteins must deal with persistent thermal fluctuations of the membrane, leading to discrepancies between the bead's curvature and the local membrane curvature sensed instantaneously by a protein. Using continuum models of fluctuating membranes, we investigate whether it is feasible for a protein acting as a perfect observer of the membrane to sense micron-scale curvature either by measuring local membrane curvature or by using bilayer lipid densities as a proxy. To do this, we develop algorithms to simulate lipid density and membrane shape fluctuations. We derive physical limits to the sensing efficacy of a protein in terms of protein size, membrane thickness, membrane bending modulus, membrane-substrate adhesion strength, and bead size. To explain the experimental protein-bead association rates, we develop two classes of predictive models: (i) for proteins that maximally associate to a preferred curvature and (ii) for proteins with enhanced association rates above a threshold curvature. We find that the experimentally observed sensing efficacy is close to the theoretical sensing limits imposed on a septin-sized protein. Protein-membrane association rates may depend on the curvature of the bead, but the strength of this dependence is limited by the fluctuations in membrane height and density.

Structure of the HIV immature lattice allows for essential lattice remodeling within budded virions

For HIV virions to become infectious, the immature lattice of Gag polyproteins attached to the virion membrane must be cleaved. Cleavage cannot initiate without the protease formed by the homo-dimerization of domains linked to Gag. However, only 5% of the Gag polyproteins, termed Gag-Pol, carry this protease domain, and they are embedded within the structured lattice. The mechanism of Gag-Pol dimerization is unknown. Here, we use spatial stochastic computer simulations of the immature Gag lattice as derived from experimental structures, showing that dynamics of the lattice on the membrane is unavoidable due to the missing 1/3 of the spherical protein coat. These dynamics allow for Gag-Pol molecules carrying the protease domains to detach and reattach at new places within the lattice. Surprisingly, dimerization timescales of minutes or less are achievable for realistic binding energies and rates despite retaining most of the large-scale lattice structure. We derive a formula allowing extrapolation of timescales as a function of interaction free energy and binding rate, thus predicting how additional stabilization of the lattice would impact dimerization times. We further show that during assembly, dimerization of Gag-Pol is highly likely and therefore must be actively suppressed to prevent early activation. By direct comparison to recent biochemical measurements within budded virions, we find that only moderately stable hexamer contacts (-12kBT<∆G<-8kBT) retain both the dynamics and lattice structures that are consistent with experiment. These dynamics are likely essential for proper maturation, and our models quantify and predict lattice dynamics and protease dimerization timescales that define a key step in understanding formation of infectious viruses.

Modeling membrane reshaping driven by dynamic protein assemblies

Remodeling of membranes in living systems is almost universally coupled to self-assembly of soluble proteins. Proteins assemble into semi-rigid shells that reshape attached membranes, and into filaments that protrude membranes. These assemblies are temporary, building from reversible protein and membrane interactions that must nucleate in the proper location. The interactions are strongly influenced by the nonequilibrium environment of the cell, such as gradients of components or active modifications by kinases. From a modeling perspective, understanding mechanisms and control thus requires 1. time-dependent approaches that ideally incorporate 2. macromolecular structure, 3. out-of-equilibrium processes, and 4. deformable membranes over microns and seconds. Realistically, tradeoffs must be made with these last three features. However, we see recent developments from the highly coarsened molecule-based scale, the continuum reaction-diffusion scale, and hybrid approaches as stimulating efforts in diverse applications. We discuss here methodological advances and progress towards simulating these processes as they occur physiologically.

Membrane curvature as a signal to ensure robustness of diverse cellular processes

An increasing corpus of research has demonstrated that membrane shape, generated either by the external environment of the cell or by intrinsic mechanisms such as cytokinesis and vesicle or organelle formation, is an important parameter in the control of diverse cellular processes. In this review we discuss recent findings that demonstrate how membrane curvature (from nanometer to micron length-scales) alters protein function. We describe an expanding toolkit for experimentally modulating membrane curvature to reveal effects on protein function, and discuss how membrane curvature - far from being a passive consequence of the physical environment and the internal protein activity of a cell - is an important signal that controls protein affinity and enzymatic activity to ensure robust forward progression of key processes within the cell.