Membrane Curvature Sensing by Amphipathic Helices Is Modulated by the Surrounding Protein Backbone

Membrane curvature is involved in numerous biological pathways like vesicle trafficking, endocytosis or nuclear pore complex assembly. In addition to its topological role, membrane curvature is sensed by specific proteins, enabling the coordination of biological processes in space and time. Amongst membrane curvature sensors are the ALPS (Amphipathic Lipid Packing Sensors). ALPS motifs are short peptides with peculiar amphipathic properties. They are found in proteins targeted to distinct curved membranes, mostly in the early secretory pathway. For instance, the ALPS motif of the golgin GMAP210 binds trafficking vesicles, while the ALPS motif of Nup133 targets nuclear pores. It is not clear if, besides curvature sensitivity, ALPS motifs also provide target specificity, or if other domains in the surrounding protein backbone are involved. To elucidate this aspect, we studied the subcellular localization of ALPS motifs outside their natural protein context. The ALPS motifs of GMAP210 or Nup133 were grafted on artificial fluorescent probes. Importantly, ALPS motifs are held in different positions and these contrasting architectures were mimicked by the fluorescent probes. The resulting chimeras recapitulated the original proteins localization, indicating that ALPS motifs are sufficient to specifically localize proteins. Modulating the electrostatic or hydrophobic content of Nup133 ALPS motif modified its avidity for cellular membranes but did not change its organelle targeting properties. In contrast, the structure of the backbone surrounding the helix strongly influenced targeting. In particular, introducing an artificial coiled-coil between ALPS and the fluorescent protein increased membrane curvature sensitivity. This coiled-coil domain also provided membrane curvature sensitivity to the amphipathic helix of Sar1. The degree of curvature sensitivity within the coiled-coil context remains correlated to the natural curvature sensitivity of the helices. This suggests that the chemistry of ALPS motifs is a key parameter for membrane curvature sensitivity, which can be further modulated by the surrounding protein backbone.     

Protein Partitioning into Ordered Membrane Domains: Insights from Simulations

Cellular membranes are laterally organized into domains of distinct structures and compositions by the differential interaction affinities between various membrane lipids and proteins. A prominent example of such structures are lipid rafts, which are ordered, tightly packed domains that have been widely implicated in cellular processes. The functionality of raft domains is driven by their selective recruitment of specific membrane proteins to regulate their interactions and functions; however, there have been few general insights into the factors that determine the partitioning of membrane proteins between coexisting liquid domains. In this work, we used extensive coarse-grained and atomistic molecular dynamics simulations, potential of mean force calculations, and conceptual models to describe the partitioning dynamics and energetics of a model transmembrane domain from the linker of activation of T cells. We find that partitioning between domains is determined by an interplay between protein-lipid interactions and differential lipid packing between raft and nonraft domains. Specifically, we show that partitioning into ordered domains is promoted by preferential interactions between peptides and ordered lipids, mediated in large part by modification of the peptides by saturated fatty acids (i.e., palmitoylation). Ordered phase affinity is also promoted by elastic effects, specifically hydrophobic matching between the membrane and the peptide. Conversely, ordered domain partitioning is disfavored by the tight molecular packing of the lipids therein. The balance of these dominant drivers determines partitioning. In the case of the wild-type linker of activation of T cells transmembrane domain, these factors combine to yield enrichment of the peptide at L_o/L_d interfaces. These results define some of the general principles governing protein partitioning between coexisting membrane domains and potentially explain previous disparities among experiments and simulations across model systems.     

Interplay between Membrane Curvature and Cholesterol: Role of Palmitoylated Caveolin-1

Caveolin-1 (cav-1) is an important player in cell signaling and endocytosis that has been shown to colocalize with cholesterol-rich membrane domains. Experimental studies with varying cav-1 constructs have suggested that it can induce both cholesterol clustering and membrane curvature. Here, we probe the molecular origin of membrane curvature and cholesterol clustering by cav-1 by using coarse-grain molecular dynamics simulations. We have performed a series of simulations of a functionally important cav-1 construct, comprising the membrane-interacting domains and a C-terminal palmitoyl tail. Our results suggest that cav-1 is able to induce cholesterol clustering in the membrane leaflet to which it is bound as well as the opposing leaflet. A positive membrane curvature is observed upon cav-1 binding in cholesterol-containing bilayers. Interestingly, we observe an interplay between cholesterol clustering and membrane curvature such that cav-1 is able to induce higher membrane curvature in cholesterol-rich membranes. The role of the cav-1 palmitoyl tail is less clear and appears to increase the membrane contacts. Further, we address the importance of the secondary structure of cav-1 domains and show that it could play an important role in membrane curvature and cholesterol clustering. Our work is an important step toward a molecular picture of caveolae and vesicular endocytosis.     

Synergism of Antimicrobial Frog Peptides Couples to Membrane Intrinsic Curvature Strain

Mixtures of the frog peptides magainin 2 and PGLa are well-known for their pronounced synergistic killing of Gram-negative bacteria. We aimed to gain insight into the underlying biophysical mechanism by interrogating the permeabilizing efficacies of the peptides as a function of stored membrane curvature strain. For Gram-negative bacterial-inner-membrane mimics, synergism was only observed when the anionic bilayers exhibited significant negative intrinsic curvatures imposed by monounsaturated phosphatidylethanolamine. In contrast, the peptides and their mixtures did not exhibit significant activities in charge-neutral mammalian mimics, including those with negative curvature, which is consistent with the requirement of charge-mediated peptide binding to the membrane. Our experimental findings are supported by computer simulations showing a significant decrease of the peptide-insertion free energy in membranes upon shifting intrinsic curvatures toward more positive values. The physiological relevance of our model studies is corroborated by a remarkable agreement with the peptide’s synergistic activity in Escherichia coli. We propose that synergism is related to a lowering of a membrane-curvature-strain-mediated free-energy barrier by PGLa that assists membrane insertion of magainin 2, and not by strict pairwise interactions of the two peptides as suggested previously.     

Modeling the Effect of Curvature on the Collective Behavior of Cells Growing New Tissue

The growth of several biological tissues is known to be controlled in part by local geometrical features, such as the curvature of the tissue interface. This control leads to changes in tissue shape that in turn can affect the tissue’s evolution. Understanding the cellular basis of this control is highly significant for bioscaffold tissue engineering, the evolution of bone microarchitecture, wound healing, and tumor growth. Although previous models have proposed geometrical relationships between tissue growth and curvature, the role of cell density and cell vigor remains poorly understood. We propose a cell-based mathematical model of tissue growth to investigate the systematic influence of curvature on the collective crowding or spreading of tissue-synthesizing cells induced by changes in local tissue surface area during the motion of the interface. Depending on the strength of diffusive damping, the model exhibits complex growth patterns such as undulating motion, efficient smoothing of irregularities, and the generation of cusps. We compare this model with in vitro experiments of tissue deposition in bioscaffolds of different geometries. By including the depletion of active cells, the model is able to capture both smoothing of initial substrate geometry and tissue deposition slowdown as observed experimentally.     

The Structural Basis of IKs Ion-Channel Activation: Mechanistic Insights from Molecular Simulations

Relating ion channel (iCh) structural dynamics to physiological function remains a challenge. Current experimental and computational techniques have limited ability to explore this relationship in atomistic detail over physiological timescales. A framework associating iCh structure to function is necessary for elucidating normal and disease mechanisms. We formulated a modeling schema that overcomes the limitations of current methods through applications of artificial intelligence machine learning. Using this approach, we studied molecular processes that underlie human IKs voltage-mediated gating. IKs malfunction underlies many debilitating and life-threatening diseases. Molecular components of IKs that underlie its electrophysiological function include KCNQ1 (a pore-forming tetramer) and KCNE1 (an auxiliary subunit). Simulations, using the IKs structure-function model, reproduced experimentally recorded saturation of gating-charge displacement at positive membrane voltages, two-step voltage sensor (VS) movement shown by fluorescence, iCh gating statistics, and current-voltage relationship. Mechanistic insights include the following: 1) pore energy profile determines iCh subconductance; 2) the entire protein structure, not limited to the pore, contributes to pore energy and channel subconductance; 3) interactions with KCNE1 result in two distinct VS movements, causing gating-charge saturation at positive membrane voltages and current activation delay; and 4) flexible coupling between VS and pore permits pore opening at lower VS positions, resulting in sequential gating. The new modeling approach is applicable to atomistic scale studies of other proteins on timescales of physiological function.     

Understanding the Role of Amphipathic Helices in N-BAR Domain Driven Membrane Remodeling

Endophilin N-BAR (N-terminal helix and Bin/amphiphysin/Rvs) domain tubulates and vesiculates lipid membranes in vitro via its crescent-shaped dimer and four amphipathic helices that penetrate into membranes as wedges. Like F-BAR domains, endophilin N-BAR also forms a scaffold on membrane tubes. Unlike F-BARs, endophilin N-BARs have N-terminal H0 amphipathic helices that are proposed to interact with other N-BARs in oligomer lattices. Recent cryo-electron microscopy reconstructions shed light on the organization of the N-BAR lattice coats on a nanometer scale. However, because of the resolution of the reconstructions, the precise positioning of the amphipathic helices is still ambiguous. In this work, we applied a coarse-grained model to study various membrane remodeling scenarios induced by endophilin N-BARs. We found that H0 helices of N-BARs prefer to align in an antiparallel manner at two ends of the protein to form a stable lattice. The deletion of H0 helices causes disruption of the lattice. In addition, we analyzed the persistence lengths of the protein-coated tubes and found that the stiffness of endophilin N-BAR-coated tubules qualitatively agrees with previous experimental work studying N-BAR-coated tubules. Large-scale simulations on membrane liposomes revealed a systematic relation between H0 helix density and local membrane curvature fluctuations. The data also suggest that the H0 helix is required for BARs to form organized structures on the liposome, further illustrating its important function.     

Molecular Insights on the Release of Avibactam from the Acyl-Enzyme Complex

Avibactam is a non-β-lactam β-lactamase inhibitor for treating complicated urinary tract and respiratory infections caused by multidrug-resistant bacterial pathogens, a serious public health threat. Despite its importance, the release mechanism of avibactam from the enzyme-inhibitor complex has been scarcely studied from first principles, considering the total protein environment. This information at the molecular level is essential for the rational design of new antibiotics and inhibitors. In this article, we addressed the release of avibactam from the complex CTX-M-15 by means of molecular dynamics simulations and quantum mechanics/molecular mechanics calculations. This study provides molecular information not available earlier, including exploration of the potential energy surfaces, characterization of the observed intermediate, and their critical points, as well. Our results show that unlike that observed in the acylation reaction, the residues Glu166 and Lys73 would be in their neutral forms. Release of avibactam follows a stepwise mechanism in which the first stage corresponds to the formation of a tetrahedral intermediate, whereas the second stage corresponds to the cleavage of the Ser70-C7 bond, mediated by Lys73, either directly or through Ser130.     

Competitive Binding of Mg2+ and Na+ Ions to Nucleic Acids: From Helices to Tertiary Structures

Nucleic acids generally reside in cellular aqueous solutions with mixed divalent/monovalent ions, and the competitive binding of divalent and monovalent ions is critical to the structures of nucleic acids because of their polyanionic nature. In this work, we first proposed a general and effective method for simulating a nucleic acid in mixed divalent/monovalent ion solutions with desired bulk ion concentrations via molecular dynamics (MD) simulations and investigated the competitive binding of Mg²⁺/Na⁺ ions to various nucleic acids by all-atom MD simulations. The extensive MD-based examinations show that single MD simulations conducted using the proposed method can yield desired bulk divalent/monovalent ion concentrations for various nucleic acids, including RNA tertiary structures. Our comprehensive analyses show that the global binding of Mg²⁺/Na⁺ to a nucleic acid is mainly dependent on its structure compactness, as well as Mg²⁺/Na⁺ concentrations, rather than the specific structure of the nucleic acid. Specifically, the relative global binding of Mg²⁺ over Na⁺ is stronger for a nucleic acid with higher effective surface charge density and higher relative Mg²⁺/Na⁺ concentrations. Furthermore, the local binding of Mg²⁺/Na⁺ to a phosphate of a nucleic acid mainly depends on the local phosphate density in addition to Mg²⁺/Na⁺ concentrations.     

Roles of Amphipathic Helices and the Bin/Amphiphysin/Rvs (BAR) Domain of Endophilin in Membrane Curvature Generation

Control of membrane curvature is required in many important cellular processes, including endocytosis and vesicular trafficking. Endophilin is a bin/amphiphysin/rvs (BAR) domain protein that induces vesicle formation by promotion of membrane curvature through membrane binding as a dimer. Using site-directed spin labeling and EPR spectroscopy, we show that the overall BAR domain structure of the rat endophilin A1 dimer determined crystallographically is maintained under predominantly vesiculating conditions. Spin-labeled side chains on the concave surface of the BAR domain do not penetrate into the acyl chain interior, indicating that the BAR domain interacts only peripherally with the surface of a curved bilayer. Using a combination of EPR data and computational refinement, we determined the structure of residues 63–86, a region that is disordered in the crystal structure of rat endophilin A1. Upon membrane binding, residues 63–75 in each subunit of the endophilin dimer form a slightly tilted, amphipathic α-helix that directly interacts with the membrane. In their predominant conformation, these helices are located orthogonal to the long axis of the BAR domain. In this conformation, the amphipathic helices are positioned to act as molecular wedges that induce membrane curvature along the concave surface of the BAR domain.     

Membrane Diffusion Occurs by Continuous-Time Random Walk Sustained by Vesicular Trafficking

Diffusion in cellular membranes is regulated by processes that occur over a range of spatial and temporal scales. These processes include membrane fluidity, interprotein and interlipid interactions, interactions with membrane microdomains, interactions with the underlying cytoskeleton, and cellular processes that result in net membrane movement. The complex, non-Brownian diffusion that results from these processes has been difficult to characterize, and moreover, the impact of factors such as membrane recycling on membrane diffusion remains largely unexplored. We have used a careful statistical analysis of single-particle tracking data of the single-pass plasma membrane protein CD93 to show that the diffusion of this protein is well described by a continuous-time random walk in parallel with an aging process mediated by membrane corrals. The overall result is an evolution in the diffusion of CD93: proteins initially diffuse freely on the cell surface but over time become increasingly trapped within diffusion-limiting membrane corrals. Stable populations of freely diffusing and corralled CD93 are maintained by an endocytic/exocytic process in which corralled CD93 is selectively endocytosed, whereas freely diffusing CD93 is replenished by exocytosis of newly synthesized and recycled CD93. This trafficking not only maintained CD93 diffusivity but also maintained the heterogeneous distribution of CD93 in the plasma membrane. These results provide insight into the nature of the biological and biophysical processes that can lead to significantly non-Brownian diffusion of membrane proteins and demonstrate that ongoing membrane recycling is critical to maintaining steady-state diffusion and distribution of proteins in the plasma membrane.     

Real-Time Insights into Biological Events: In-Cell Processes and Protein-Ligand Interactions

FlowNMR has the aim of continuously monitoring processes that occur in conditions that are not compatible with being carried out within a closed tube. However, it is sample intensive and not suitable for samples, such as proteins or living cells, that are often available in limited volumes and possibly low concentrations. We here propose a dialysis-based modification of a commercial flowNMR setup that allows for recycling the medium while confining the sample (proteins and cells) within the active volume of the tube. This approach is demonstrated in the specific cases of in-cell NMR and protein-based ligand studies.     

Membrane Morphologies Induced by Arc-Shaped Scaffolds Are Determined by Arc Angle and Coverage

The intricate shapes of biological membranes such as tubules and membrane stacks are induced by proteins. In this article, we systematically investigate the membrane shapes induced by arc-shaped scaffolds such as proteins and protein complexes with coarse-grained modeling and simulations. We find that arc-shaped scaffolds induce membrane tubules at membrane coverages larger than a threshold of ～40%, irrespective of their arc angle. The membrane morphologies at intermediate coverages below this tubulation threshold, in contrast, strongly depend on the arc angle. Scaffolds with arc angles of about 60°, akin to N-BAR domains, do not change the membrane shape at coverages below the tubulation threshold, whereas scaffolds with arc angles larger than about 120° induce double-membrane stacks at intermediate coverages. The scaffolds stabilize the curved membrane edges that connect the membrane stacks, as suggested for complexes of reticulon proteins. Our results provide general insights on the determinants of membrane shaping by arc-shaped scaffolds.     

Gaussian Curvature Directs Stress Fiber Orientation and Cell Migration

We show that substrates with nonzero Gaussian curvature influence the organization of stress fibers and direct the migration of cells. To study the role of Gaussian curvature, we developed a sphere-with-skirt surface in which a positive Gaussian curvature spherical cap is seamlessly surrounded by a negative Gaussian curvature draping skirt, both with principal radii similar to cell-length scales. We find significant reconfiguration of two subpopulations of stress fibers when fibroblasts are exposed to these curvatures. Apical stress fibers in cells on skirts align in the radial direction and avoid bending by forming chords across the concave gap, whereas basal stress fibers bend along the convex direction. Cell migration is also strongly influenced by the Gaussian curvature. Real-time imaging shows that cells migrating on skirts repolarize to establish a leading edge in the azimuthal direction. Thereafter, they migrate in that direction. This behavior is notably different from migration on planar surfaces, in which cells typically migrate in the same direction as the apical stress fiber orientation. Thus, this platform reveals that nonzero Gaussian curvature not only affects the positioning of cells and alignment of stress fiber subpopulations but also directs migration in a manner fundamentally distinct from that of migration on planar surfaces.     

Active Transport of Membrane Components by Self-Organization of the Min Proteins

Heterogeneous distribution of components in the biological membrane is critical in the process of cell polarization. However, little is known about the mechanisms that can generate and maintain the heterogeneous distribution of the membrane components. Here, we report that the propagating wave patterns of the bacterial Min proteins can impose steric pressure on the membrane, resulting in transport and directional accumulation of the component in the membrane. Therefore, the membrane component waves represent transport of the component in the membrane that is caused by the steric pressure gradient induced by the differential levels of binding and dissociation of the Min proteins in the propagating waves on the membrane surface. The diffusivity, majorly influenced by the membrane anchor of the component, and the repulsed ability, majorly influenced by the steric property of the membrane component, determine the differential spatial distribution of the membrane component. Thus, transportation of the membrane component by the Min proteins follows a simple physical principle, which resembles a linear peristaltic pumping process, to selectively segregate and maintain heterogeneous distribution of materials in the membrane.

Modulation of DEG/ENaCs by Amphiphiles Suggests Sensitivity to Membrane Alterations

The bile acid-sensitive ion channel is activated by amphiphilic substances such as bile acids or artificial detergents via membrane alterations; however, the mechanism of membrane sensitivity of the bile acid-sensitive ion channel is not known. It has also not been systematically investigated whether other members of the degenerin/epithelial Na⁺ channel (DEG/ENaC) gene family are affected by amphiphilic compounds. Here, we show that DEG/ENaCs ASIC1a, ASIC3, ENaC, and the purinergic receptor P2X2 are modulated by a large number of different, structurally unrelated amphiphilic substances, namely the detergents N-lauroylsarcosine, Triton X-100, and β-octylglucoside; the fenamate flufenamic acid; the antipsychotic drug chlorpromazine; the natural phenol resveratrol; the chili pepper compound capsaicin; the loop diuretic furosemide; and the antiarrythmic agent verapamil. We determined the modification of membrane properties using large-angle x-ray diffraction experiments on model lipid bilayers, revealing that the amphiphilic compounds are positioned in a characteristic fashion either in the lipid tail group region or in the lipid head group region, demonstrating that they perturbed the membrane structure. Collectively, our results show that DEG/ENaCs and structurally related P2X receptors are modulated by diverse amphiphilic molecules. Furthermore, they suggest alterations of membrane properties by amphiphilic compounds as a mechanism contributing to modulation.