Lipid asymmetry in plasma membrane of eukaryotes is ubiquitous. The first measurements reported compositional asymmetry: phosphatidylethanolamine and phosphatidylserine are mostly on the cytoplasmic leafet, while phosphatidylcholine and sphingomyelin are mostly on the exoplasmic leaflet. More recent experiments using lipidomics have evidenced the presence of saturation asymmetry between the two leaflets. A question that naturally arises is why such an asymmetry? To complicate matters, it is still largely unknown in which leaflet cholesterol lies. Here, we use chemical potentials to mimic flippase proteins responsible for maintenance of compositional asymmetry in silico. We show that saturation asymmetry naturally arises as a byproduct of phospholipid number asymmetry and sphingomyelin contents, thereby showing that some reported asymmetries may naturally result from others and do not necessarily require being externally driven. We also show that plasmalogen lipids' tendency to be highly unsaturated is also natural. Additionally, we tackle the problem of cholesterol and show that, while it is influenced by all asymmetries, the resulting cholesterol asymmetry tends to be fairly mild.

Computer simulation techniques form a versatile tool, a computational microscope, for exploring biological processes. This tool has been particularly effective in exploring different features of biological membranes. In recent years, thanks to elegant multiscale simulation schemes, some fundamental limitations of investigations by distinct simulation techniques have been resolved. As a result, we are now capable of exploring processes spanning multiple scales beyond the capacity of any single technique. In this perspective, we argue that mesoscale simulations require more attention and must be further developed to fill evident gaps in a quest toward simulating and modeling living cell membranes.

The mammalian cell membrane consists of thousands of different lipid species, and this variety is critical for biological function. Alterations to this balance can be dangerous as they can lead to permanent disruption of lipid metabolism, a hallmark in several viral diseases. The Flaviviridae family is made up of positive single-stranded RNA viruses that assemble at or near the location of lipid droplet formation in the endoplasmic reticulum. These viruses are known to interfere with lipid metabolism during the onset of liver disease, albeit to different extents. Pathogenesis of these infections involves specific protein-lipid interactions that alter lipid sorting and metabolism to sustain propagation of the viral infection. Recent experimental studies identify a correlation between viral proteins and lipid content or location in the cell, but these do not assess membrane-embedded interactions. Molecular modeling, specifically molecular dynamics simulations, can provide molecular-level spatial and temporal resolution for characterization of biomolecular interactions. This review focuses on recent advancements and current knowledge gaps in the molecular mechanisms of lipid-mediated liver disease preceded by viral infection. We discuss three viruses from the Flaviviridae family: dengue, zika, and hepatitis C, with a particular focus on lipid interactions with their respective ion channels, known as viroporins.

Gram-negative bacteria are equipped with a cell wall that contains a complex matrix of lipids, proteins, and glycans, which form a rigid layer protecting bacteria from the environment. Major components of this outer membrane are the high-molecular weight and amphiphilic lipopolysaccharides (LPSs). They form the extracellular part of a heterobilayer with phospholipids. Understanding LPS properties within the outer membrane is therefore important to develop new antimicrobial strategies. Model systems, such as giant unilamellar vesicles (GUVs), provide a suitable platform for exploring membrane properties and interactions. However, LPS molecules contain large polysaccharide parts that confer high water solubility, which makes LPS incorporation in artificial membranes difficult; this hindrance is exacerbated for LPS with long polysaccharide chains, i.e., the smooth LPS. Here, a novel emulsification step of the inverted emulsion method is introduced to incorporate LPS in the outer or the inner leaflet of GUVs, exclusively. We developed an approach to determine the LPS content on individual GUVs and quantify membrane asymmetry. The asymmetric membranes with outer leaflet LPS show incorporations of 1-16 mol % smooth LPS (corresponding to 16-79 wt %), while vesicles with inner leaflet LPS reach coverages of 2-7 mol % smooth LPS (28-60 wt %). Diffusion coefficient measurements in the obtained GUVs showed that increasing LPS concentrations in the membranes resulted in decreased diffusivity.

Lateral phase heterogeneity in biomembranes can govern cellular functions and may serve as a platform for enrichment or depletion of membrane-anchored molecules. In this work, we address the question of how the process of membrane fusion is affected by the membrane phase state (fluid or gel) and by phase coexistence, as well as the effects of fusion-mediated incorporation of exogeneous lipids on phase separation. Our system is based on the fusion of cationic fluid large unilamellar vesicles (LUVs) composed of dioleoyl trimethylammonium propane (DOTAP) and dioleoyl phosphoethanolamine (DOPE) with neutral and anionic giant unilamellar vesicles (GUVs) composed of phosphatidylcholine and phosphatidylglycerol. By changing the lipid composition of the GUVs, we modulated the phase state and charge of the different phases (charged or neutral, fluid or gel) and identified systems in which we can target fusion to specific domains on phase-separated membranes. Fusion efficiency was quantified using fluorescence microscopy-based lipid and content mixing assays, and flow chamber devices were used to assess the real-time sequence of events of the fusion process. To investigate the bilayer thermal behavior, differential scanning calorimetry (DSC) experiments were performed on LUVs. The results show that fusion is extensive in single-component GUVs only for fluid and negatively charged acceptor membranes. On the other hand, in phase-separated GUVs, high fusion efficiency was observed even when the gel phase was anionic and phase separation somewhat increased the fusion efficiency. Extensive fusion led to dissolution of the gel domains as a result of extensive incorporation of lipids in the fluid state from the fusogenic liposomes. Altogether, these findings have the potential to unravel the important role of membrane phase state, phase separation, charge, and the effects of extensive fusion on membrane organization and may give insights in the regulation of the interactions between cells and liposomes that are used in drug delivery systems.

Functionally meaningful reversible protein-membrane interactions mediate many biological events. Fluorescence correlation spectroscopy (FCS) is increasingly used to quantitatively study the non-reversible binding of proteins to membranes using lipid vesicles in solution. However, the lack of a complete description of the phase and statistical equilibria in the case of reversible protein-membrane partitioning has hampered the application of FCS to quantify the partition coefficient (Kx). In this work, we further extend the theory that describes membrane-protein partitioning to account for spontaneous protein-membrane dissociation and reassociation to the same or a different lipid vesicle. We derive the probability distribution of proteins on lipid vesicles for reversible binding and demonstrate that FCS is a suitable technique for accurate Kx quantification of membrane-protein reversible association. We also establish the limits to Kx determination by FCS studying the Cramer-Rao bound on the variance of the retrieved parameters. We validate the mathematical formulation against reaction-diffusion simulations to study phase and statistical equilibria and compare the Kx obtained from a computational FCS titration experiment with the experimental ground truth. Finally, we demonstrate the application of our methodology studying the association of anti-HIV broadly neutralizing antibody (10E8-3R) to the membrane.

We investigate the effects on the distribution of lipids in the plasma membrane that are caused by the insertion of a protein, Piezo1, that significantly distorts the membrane toward the cytosol. From coarse-grained molecular dynamics simulations, we find that the major effects occur in the outer, extracellular, leaflet. The mol fraction of cholesterol increases significantly in the curved region of the membrane close to Piezo1, while those of phosphatidylcholine and of sphingomyelin decrease. In the inner leaflet, mol fractions of cholesterol and of phosphatidylethanolamine decrease slightly as the protein is approached, while that of phosphatidylserine increases slightly. The mol fraction of phosphatidylcholine decreases markedly as the protein is approached. Most of these results are understood in the context of a theoretical model that utilizes two elements: 1) a coupling between the leaflets' actual curvatures and their compositionally dependent spontaneous curvatures and 2) the dependence of the spontaneous curvatures not only on the mol fractions of the phospholipids, but also on the effect that cholesterol has on the spontaneous curvatures of the phospholipids.

Tubular networks of the endoplasmic reticulum (ER) are dynamic structures whose steady-state conformations are maintained by a balance between the persistent generation and vanishing of the network elements. While factors producing the ER tubules and intertubular junctions have been investigated, the mechanisms behind their elimination remained unknown. Here, we addressed the ER ring closure, the process resulting in the tubule and junction removal through constriction of the network unit cells into junctional knots followed by the knot remodeling into regular junctions. We considered the ring closure to be driven by the tension existing in ER membranes. We based our consideration on the notion of Gibbs' thermodynamic tension and reviewed its relationship to other tension definitions used in the literature. We modeled, computationally, the structures of the junctional knots containing internal nanopores and analyzed their tension dependence. We analyzed the process of the pore sealing through membrane fission resulting in the formation of regular junctions. Considering the hemi-fission as the rate-limiting stage of the fission reaction, we evaluated the membrane tensions guaranteeing the spontaneous character of the pore sealing. We concluded that feasible membrane tensions explain all stages of the ER ring closure.

The twin arginine translocase (Tat) exports folded proteins across bacterial membranes. The putative pore-forming or membrane-weakening component (TatAd in B. subtilis) is anchored to the lipid bilayer via an unusually short transmembrane α-helix (TMH), with less than 16 residues. Its tilt angle in different membranes was analyzed under hydrophobic mismatch conditions, using synchrotron radiation circular dichroism and solid-state NMR. Positive mismatch (introduced either by reconstitution in short-chain lipids or by extending the hydrophobic TMH length) increased the helix tilt of the TMH as expected. Negative mismatch (introduced either by reconstitution in long-chain lipids or by shortening the TMH), on the other hand, led to protein aggregation. These data suggest that the TMH of TatA is just about long enough for stable membrane insertion. At the same time, its short length is a crucial factor for successful translocation, as demonstrated here in native membrane vesicles using an in vitro translocation assay. Furthermore, when reconstituted in model membranes with negative spontaneous curvature, the TMH was found to be aligned parallel to the membrane surface. This intrinsic ability of TatA to flip out of the membrane core thus seems to play a key role in its membrane-destabilizing effect during Tat-dependent translocation.

An approach for the efficient simulation of phase-separated lipid bilayers, for use in the calculation of equilibrium free energies of partitioning between lipid domains, is proposed. The methodology exploits restraint potentials and rectangular aspect ratios that enforce lipid phase separation, allowing for the simulation of smaller systems that approximately reproduce bulk behavior. The utility of this approach is demonstrated through the calculation of potentials of mean force for the translation of a transmembrane protein between lipid domains. The impact of the imposed restraints on lipid tail ordering and lipid packing are explored, providing insight into how restraints can best be employed to compute accurate free-energy surfaces. This approach should be useful in the accurate calculation of equilibrium partition coefficients for transmembrane protein partitioning in heterogeneous membranes, providing insight into the thermodynamic driving forces that control this fundamental biophysical phenomenon.

In cell membranes, proteins and lipids are organized into submicrometric nanodomains of varying sizes, shapes, and compositions, performing specific functions. Despite their biological importance, the detailed morphology of these nanodomains remains unknown. Not only can they hardly be observed by conventional microscopy due to their small size, but there is no full consensus on the theoretical models to describe their structuring and their shapes. Here, we use a combination of analytical calculations and Monte Carlo simulations based upon a model coupling membrane composition and shape to show that increasing protein concentration leads to an elongation of membrane nanodomains. The results are corroborated by single-particle tracking measurements on HIV receptors, whose level of expression in the membrane of specifically designed living cells can be tuned. These findings highlight that protein abundance can modulate nanodomain shape and potentially their biological function. Beyond biomembranes, this mesopatterning mechanism is of relevance in several soft-matter systems because it relies on generic physical arguments.

"Rafts" are nanometer-size inhomogeneities in the plasma membrane that, in the outer leaflet, are enriched in sphingomyelin and cholesterol. They are thought to provide a platform for proteins to carry out biological processes. Here, we employ a model asymmetric plasma membrane to address the question of the range of sphingomyelin and cholesterol compositions in which one would expect the formation of rafts. We define a weight for the likelihood of raft formation and evaluate it as a function of the sphingomyelin mole fraction in the outer leaflet for three bilayers with total cholesterol mole fractions of 0.30, 0.40, and 0.50. Not surprisingly, the weight decreases when there is little sphingomyelin. Less expected, we find that the weight also decreases when there is a large mole fraction of sphingomyelin. The weight is largest in the bilayer with a total cholesterol mole fraction of 0.30 and decreases rapidly with increasing total cholesterol. We explicate the reasons for these behaviors. In the 0.30 cholesterol bilayer, the largest weight occurs at a sphingomyelin mole fraction in the outer leaflet of approximately 0.23. The weight falls to one half its maximum value at sphingomyelin mole fractions of 0.15 and 0.33. In terms of the sphingomyelin mole fraction of the asymmetric bilayer, the maximum weight occurs at 0.12 and falls to half maximum at 0.08 and 0.17.

Plasma membranes as well as their simplified model systems show an inherent nanoscale heterogeneity. As a result of strong interleaflet interactions, these nanoheterogeneities (called here lipid nanodomains) can be found in perfect registration (i.e., nanodomains in the inner leaflet are registered with the nanodomains in the outer leaflet). Alternatively, they might be interleaflet independent, antiregistered, or located asymmetrically in one bilayer leaflet only. To distinguish these scenarios from each other appears to be an experimental challenge. In this work, we analyzed the potential of Förster resonance energy transfer to characterize interleaflet organization of nanodomains. We generated in silico time-resolved fluorescence decays for a large set of virtual as well as real donor/acceptor pairs distributed over the bilayer containing registered, independent, antiregistered, or asymmetrically distributed nanodomains. In this way, we were able to identify conditions that gave satisfactory or unsatisfactory resolution. Overall, Förster resonance energy transfer appears as a robust method that, when using donor/acceptor pairs with good characteristics, yields otherwise difficult-to-reach characteristics of membrane lipid nanodomains.

Next-generation sequencing of human genomes reveals millions of missense variants, some of which may lead to loss of protein function and ultimately disease. Here, we investigate missense variants in membrane proteins-key drivers in cell signaling and recognition. We find enrichment of pathogenic variants in the transmembrane region across 19,000 functionally classified variants in human membrane proteins. To accurately predict variant consequences, one fundamentally needs to understand the underlying molecular processes. A key mechanism underlying pathogenicity in missense variants of soluble proteins has been shown to be loss of stability. Membrane proteins, however, are widely understudied. Here, we interpret variant effects on a larger scale by performing structure-based estimations of changes in thermodynamic stability using a membrane-specific energy function and analyses of sequence conservation during evolution of 15 transmembrane proteins. We find evidence for loss of stability being the cause of pathogenicity in more than half of the pathogenic variants, indicating that this is a driving factor also in membrane-protein-associated diseases. Our findings show how computational tools aid in gaining mechanistic insights into variant consequences for membrane proteins. To enable broader analyses of disease-related and population variants, we include variant mappings for the entire human proteome.

Neutralizing antibodies against influenza have generally been classified according to their recognition sites, with antibodies against the head domain of hemagglutinin thought to inhibit attachment and antibodies against the stalk region thought to inhibit fusion. Here, we report the development of a microfluidic assay to measure neutralization of viral entry that can clearly differentiate between effects on attachment and fusion. Testing multiple broadly neutralizing antibodies against the hemagglutinin stalk domain, we obtain a surprising result: some broadly neutralizing antibodies inhibit fusion only, while others inhibit both fusion and viral attachment. Antibodies binding the globular head domain primarily inhibit attachment but can also reduce the fusogenic capability of viral particles that nonetheless bind the receptor. These findings shed light on the unexpectedly heterogeneous mechanisms of antibody neutralization even within similar recognition sites. The assay we have developed also provides a tool to optimize vaccine design by permitting assessment of the elicited antibody response with greater mechanistic resolution.

Many phospholipid membranes in the cell have a high curvature; for instance, in caveolae, mitochondrial crystae, nanotubes, membrane pearls, small liposomes, or exosomes. Molecular dynamics (MD) simulations are a computational tool to gain insight in the transport behavior at the atomic scale. Membrane permeability is a key kinetic property that might be affected in these highly curved membranes. Unfortunately, the geometry of highly curved membranes creates ambiguity in the permeability value, even with an arbitrarily large factor purely based on geometry, caused by the radial flux not being a constant value in steady state. In this contribution, the ambiguity in permeability for liposomes is countered by providing a new permeability definition. First, the inhomogeneous solubility diffusion model based on the Smoluchowski equation is solved analytically under radial symmetry, from which the entrance and escape permeabilities are defined. Next, the liposome permeability is defined guided by the criterion that a flat and curved membrane should have equal permeability, in case these were to be carved out from an imaginary homogeneous medium. With this criterion, our new definition allows for a fair comparison of flat and curved membranes. The definition is then transferred to the counting method, which is a practical computational approach to derive permeability by counting complete membrane crossings. Finally, the usability of the approach is illustrated with MD simulations of diphosphatidylcholine (DPPC) bilayers, without or with some cholesterol content. Our new liposome permeability definition allows us to compare a spherically shaped membrane with its flat counterpart, thus showcasing how the curvature effect on membrane transport may be assessed.

Actin networks polymerize and depolymerize to construct highly organized structures, thereby endowing the mechanical phenotypes found in a cell. It is generally believed that the amount of filamentous actin and actin network architecture determine cytoplasmic viscoelasticity of the whole cell. However, the intrinsic complexity of a cell and the presence of endogenous cellular components make it difficult to study the differential roles of distinct actin networks in regulating cell mechanics. Here, we model a cell by using giant unilamellar vesicles (GUVs) encapsulating actin filaments and networks assembled by various actin cross-linker proteins. Perturbation of these cytoskeletal vesicles using alternating current electric fields revealed that deformability depends on actin network architecture. While actin-free vesicles exhibited large electromechanical deformations, deformations of GUVs encapsulating actin filaments were significantly dampened. The suppression of electrodeformation of actin-GUVs can be similarly recapitulated by using aqueous poly(ethylene glycol) 8000 solutions at different concentrations to modulate solution viscoelasticity. Furthermore, networks cross-linked by alpha actinin resulted in decreased GUV deformability compared with actin-filament-encapsulating GUVs, and membrane-associated actin networks, through the formation of the dendritic actin cortex, greatly dampened electrodeformation of GUVs. These results highlight that the organization of actin networks regulates the mechanics of GUVs and shed insights into the origin of differential deformability of cells.

The actin cortex is a complex cytoskeletal machinery that drives and responds to changes in cell shape. It must generate or adapt to plasma membrane curvature to facilitate diverse functions such as cell division, migration, and phagocytosis. Due to the complex molecular makeup of the actin cortex, it remains unclear whether actin networks are inherently able to sense and generate membrane curvature, or whether they rely on their diverse binding partners to accomplish this. Here, we show that curvature sensing is an inherent capability of branched actin networks nucleated by Arp2/3 and VCA. We develop a robust method to encapsulate actin inside giant unilamellar vesicles (GUVs) and assemble an actin cortex at the inner surface of the GUV membrane. We show that actin forms a uniform and thin cortical layer when present at high concentration and distinct patches associated with negative membrane curvature at low concentration. Serendipitously, we find that the GUV production method also produces dumbbell-shaped GUVs, which we explain using mathematical modeling in terms of membrane hemifusion of nested GUVs. We find that branched actin networks preferentially assemble at the neck of the dumbbells, which possess a micrometer-range convex curvature comparable with the curvature of the actin patches found in spherical GUVs. Minimal branched actin networks can thus sense membrane curvature, which may help mammalian cells to robustly recruit actin to curved membranes to facilitate diverse cellular functions such as cytokinesis and migration.

Remorins are a family of multigenic plasma membrane phosphoproteins involved in biotic and abiotic plant interaction mechanisms, partnering in molecular signaling cascades. Signaling activity of remorins depends on their phosphorylation states and subsequent clustering into nanosized membrane domains. The presence of a coiled-coil domain and a C-terminal domain is crucial to anchor remorins to negatively charged membrane domains; however, the exact role of the N-terminal intrinsically disordered domain (IDD) on protein clustering and lipid interactions is largely unknown. Here, we combine chemical biology and imaging approaches to study the partitioning of group 1 remorin into anionic model membranes mimicking the inner leaflet of the plant plasma membrane. Using reconstituted membranes containing a mix of saturated and unsaturated phosphatidylcholine, phosphatidylinositol phosphates, and sterol, we investigate the clustering of remorins to the membrane and monitor the formation of nanosized membrane domains. REM1.3 promoted membrane nanodomain organization on the exposed external leaflet of both spherical lipid vesicles and flat supported lipid bilayers. Our results reveal that REM1.3 drives a mechanism allowing lipid reorganization, leading to the formation of remorin-enriched nanodomains. Phosphorylation of the N-terminal IDD by the calcium protein kinase CPK3 influences this clustering and can lead to the formation of smaller and more disperse domains. Our work reveals the phosphate-dependent involvement of the N-terminal IDD in the remorin-membrane interaction process by driving structural rearrangements at lipid-water interfaces.

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.

Mycobacteria, such as Mycobacterium tuberculosis, are characterized by a uniquely thick and waxy cell envelope that consists of two membranes, with a variety of mycolates comprising their outer membrane (OM). The protein Mycobacterial membrane protein Large 3 (MmpL3) is responsible for the transport of a primary OM component, trehalose monomycolate (TMM), from the inner (cytoplasmic) membrane (IM) to the periplasmic space, a process driven by the proton gradient. Although multiple structures of MmpL3 with bound substrates have been solved, the exact pathway(s) for TMM or proton transport remains elusive. Here, employing molecular dynamics simulations we investigate putative pathways for either transport species. We hypothesized that MmpL3 will cycle through similar conformational states as the related transporter AcrB, which we used as targets for modeling the conformation of MmpL3. A continuous water pathway through the transmembrane region was found in one of these states, illustrating a putative pathway for protons. Additional equilibrium simulations revealed that TMM can diffuse from the membrane into a binding pocket in MmpL3 spontaneously. We also found that acetylation of TMM, which is required for transport, makes it more stable within MmpL3's periplasmic cavity compared with the unacetylated form.

AlphaFold2 (AF2) has revolutionized the field of protein structural prediction. Here, we test its ability to predict the tertiary and quaternary structure of a previously undescribed scaffold with new folds and unusual architecture, the monotopic membrane protein caveolin-1 (CAV1). CAV1 assembles into a disc-shaped oligomer composed of 11 symmetrically arranged protomers, each assuming an identical new fold, and contains the largest parallel β-barrel known to exist in nature. Remarkably, AF2 predicts both the fold of the protomers and the interfaces between them. It also assembles between seven and 15 copies of CAV1 into disc-shaped complexes. However, the predicted multimers are energetically strained, especially the parallel β-barrel. These findings highlight the ability of AF2 to correctly predict new protein folds and oligomeric assemblies at a granular level while missing some elements of higher-order complexes, thus positing a new direction for the continued development of deep-learning protein structure prediction approaches.

The function of the actin cytoskeleton in cellular motility and trafficking has been widely studied. However, reorganization of the actin cytoskeleton upon modulation of membrane cholesterol and its consequences on membrane dynamics are addressed only rarely. In a recent work, we reported that chronic cholesterol depletion using statins leads to significant polymerization of the actin cytoskeleton. In this work, we explore the effect of reorganization of the actin cytoskeleton on membrane dynamics under cholesterol-depleted condition. Specifically, we explore the role of actin cytoskeleton in regulating the dynamics of the serotonin1A receptor, a crucial neurotransmitter G protein-coupled receptor (GPCR) that plays a major role in the generation and modulation of cognitive and behavioral functions. For this, we analyzed the lateral dynamics of the serotonin1A receptor in cholesterol-depleted cells (using statins) by fluorescence recovery after photobleaching (FRAP) and fluorescence loss in photobleaching (FLIP) measurements. Our results indicate that lateral diffusion parameters of serotonin1A receptors in normal cells are consistent with models describing the diffusion of molecules in a homogeneous membrane. Interestingly, these parameters are altered in cholesterol-depleted cells and the receptor exhibits dynamic confinement. Notably, our results show that statin-induced dynamic confinement could be reversed by destabilization of the actin cytoskeleton. On a broader perspective, these results assume significance in understanding the modulatory role of the membrane environment on the organization and dynamics of GPCRs in diseases caused by altered cholesterol biosynthesis.

The mechanical properties of red blood cells (RBCs) play key roles in their biological functions in microcirculation. In particular, RBCs must deform significantly to travel through microcapillaries with sizes comparable with or even smaller than their own. Although the dynamics of RBCs in microcapillaries have received considerable attention, the effect of membrane viscoelasticity has been largely overlooked. In this work, we present a computational study based on the boundary integral method and thin-shell mechanics to examine how membrane viscoelasticity influences the dynamics of RBCs flowing through straight and constricted microcapillaries. Our results reveal that the cell with a viscoelastic membrane undergoes substantially different motion and deformation compared with results based on a purely elastic membrane model. Comparisons with experimental data also suggest the importance of accounting for membrane viscoelasticity to properly capture the transient dynamics of an RBC flowing through a microcapillary. Taken together, these findings demonstrate the significant effects of membrane viscoelasticity on RBC dynamics in different microcapillary environments. The computational framework also lays the groundwork for more accurate quantitative modeling of the mechanical response of RBCs in their mechanotransduction process in subsequent investigations.

The precise spatiotemporal control of nanoscale membrane shape and composition is the result of a complex interplay of individual and collective molecular behaviors. Here, we employed single-molecule localization microscopy and computational simulations to observe single-lipid diffusion and sorting in model membranes with varying compositions, phases, temperatures, and curvatures. Supported lipid bilayers were created over 50-nm-radius nanoparticles to mimic the size of naturally occurring membrane buds, such as endocytic pits and the formation of viral envelopes. The curved membranes recruited liquid-disordered lipid phases while altering the diffusion and sorting of tracer lipids. Disorder-preferring fluorescent lipids sorted to and experienced faster diffusion on the nanoscale curvature only when embedded in a membrane capable of sustaining lipid phase separation at low temperatures. The curvature-induced sorting and faster diffusion even occurred when the sample temperature was above the miscibility temperature of the planar membrane, implying that the nanoscale curvature could induce phase separation in otherwise homogeneous membranes. Further confirmation and understanding of these results are provided by continuum and coarse-grained molecular dynamics simulations with explicit and spontaneous curvature-phase coupling, respectively. The curvature-induced membrane compositional heterogeneity and altered dynamics were achieved only with a coupling of the curvature with a lipid phase separation. These cross-validating results demonstrate the complex interplay of lipid phases, molecular diffusion, and nanoscale membrane curvature that are critical for membrane functionality.

Biological and model membranes are frequently subjected to fluid shear stress. However, membrane mechanical responses to flow remain incompletely described. This is particularly true of membranes supported on a solid substrate, and the influences of membrane composition and substrate roughness on membrane flow responses remain poorly understood. Here, we combine microfluidics, fluorescence microscopy, and neutron reflectivity to explore how supported lipid bilayer patches respond to controlled shear stress. We demonstrate that lipid membranes undergo a significant, passive, and partially reversible increase in membrane area due to flow. We show that these fluctuations in membrane area can be constrained, but not prevented, by increasing substrate roughness. Similar flow-induced changes to membrane structure may contribute to the ability of living cells to sense and respond to flow.

Discoidal lipid-protein nanoparticles known as nanodiscs are widely used tools in structural and membrane biology. Amphipathic, synthetic copolymers have recently become an attractive alternative to membrane scaffold proteins for the formation of nanodiscs. Such copolymers can directly intercalate into, and form nanodiscs from, intact membranes without detergents. Although these copolymer nanodiscs can extract native membrane lipids, it remains unclear whether native membrane properties are also retained. To determine the extent to which bilayer lipid packing is retained in nanodiscs, we measured the behavior of packing-sensitive fluorescent dyes in various nanodisc preparations compared with intact lipid bilayers. We analyzed styrene-maleic acid (SMA), diisobutylene-maleic acid (DIBMA), and polymethacrylate (PMA) as nanodisc scaffolds at various copolymer-to-lipid ratios and temperatures. Measurements of Laurdan spectral shifts revealed that dimyristoyl-phosphatidylcholine (DMPC) nanodiscs had increased lipid headgroup packing compared with large unilamellar vesicles (LUVs) above the lipid melting temperature for all three copolymers. Similar effects were observed for DMPC nanodiscs stabilized by membrane scaffolding protein MSP1E1. Increased lipid headgroup packing was also observed when comparing nanodiscs with intact membranes composed of binary mixtures of 1-palmitoyl-2-oleoyl-phosphocholine (POPC) and di-palmitoyl-phosphocholine (DPPC), which show fluid-gel-phase coexistence. Similarly, Laurdan reported increased headgroup packing in nanodiscs for biomimetic mixtures containing cholesterol, most notable for relatively disordered membranes. The magnitudes of these ordering effects were not identical for the various copolymers, with SMA being the most and DIBMA being the least perturbing. Finally, nanodiscs derived from mammalian cell membranes showed similarly increased lipid headgroup packing. We conclude that nanodiscs generally do not completely retain the physical properties of intact membranes.

Many transient processes in cells arise from the binding of cytosolic proteins to membranes. Quantifying this membrane binding and its associated diffusion in the living cell is therefore of primary importance. Dynamic photonic microscopies, e.g., single/multiple particle tracking, fluorescence recovery after photobleaching, and fluorescence correlation spectroscopy (FCS), enable non-invasive measurement of molecular mobility in living cells and their plasma membranes. However, FCS with a single beam waist is of limited applicability with complex, non-Brownian, motions. Recently, the development of FCS diffusion laws methods has given access to the characterization of these complex motions, although none of them is applicable to the membrane binding case at the moment. In this study, we combined computer simulations and FCS experiments to propose an FCS diffusion law for membrane binding. First, we generated computer simulations of spot-variation FCS (svFCS) measurements for a membrane binding process combined to 2D and 3D diffusion at the membrane and in the bulk/cytosol, respectively. Then, using these simulations as a learning set, we derived an empirical diffusion law with three free parameters: the apparent binding constant KD, the diffusion coefficient on the membrane D2D, and the diffusion coefficient in the cytosol, D3D. Finally, we monitored, using svFCS, the dynamics of retroviral Gag proteins and associated mutants during their binding to supported lipid bilayers of different lipid composition or at plasma membranes of living cells, and we quantified KD and D2D in these conditions using our empirical diffusion law. Based on these experiments and numerical simulations, we conclude that this new approach enables correct estimation of membrane partitioning and membrane diffusion properties (KD and D2D) for peripheral membrane molecules.

Pancratistatin (PST) and narciclasine (NRC) are natural therapeutic agents that exhibit specificity toward the mitochondria of cancerous cells and initiate apoptosis. Unlike traditional cancer therapeutic agents, PST and NRC are effective, targeted, and have limited adverse effects on neighboring healthy, noncancerous cells. Currently, the mechanistic pathway of action for PST and NRC remains elusive, which in part inhibits PST and NRC from becoming efficacious therapeutic alternatives. Herein, we use neutron and x-ray scattering in combination with calcein leakage assays to characterize the effects of PST, NRC, and tamoxifen (TAM) on a biomimetic model membrane. We report an increase in lipid flip-flop half-times (t1/2) (≈12.0%, ≈35.1%, and a decrease of ≈45.7%) with 2 mol percent PST, NRC, and TAM respectively. An increase in bilayer thickness (≈6.3%, ≈7.8%, and ≈7.8%) with 2 mol percent PST, NRC, and TAM, respectively, was also observed. Lastly, increases in membrane leakage (≈31.7%, ≈37.0%, and ≈34.4%) with 2 mol percent PST, NRC, and TAM, respectively, were seen. Considering the maintenance of an asymmetric lipid composition across the outer mitochondrial membrane (OMM) is crucial to eukaryotic cellular homeostasis and survival, our results suggest PST and NRC may play a role in disrupting the native distribution of lipids within the OMM. A possible mechanism of action for PST- and NRC-induced mitochondrial apoptosis is proposed via the redistribution of the native OMM lipid organization and through OMM permeabilization.

Sphingolipids are a structurally diverse class of lipids predominantly found in the plasma membrane of eukaryotic cells. These lipids can laterally segregate with other rigid lipids and cholesterol into liquid-ordered domains that act as organizing centers within biomembranes. Owing the vital role of sphingolipids for lipid segregation, controlling their lateral organization is of utmost significance. Hence, we made use of the light-induced trans-cis isomerization of azobenzene-modified acyl chains to develop a set of photoswitchable sphingolipids with different headgroups (hydroxyl, galactosyl, phosphocholine) and backbones (sphingosine, phytosphingosine, tetrahydropyran-blocked sphingosine) that are able to shuttle between liquid-ordered and liquid-disordered regions of model membranes upon irradiation with UV-A (λ = 365 nm) and blue (λ = 470 nm) light, respectively. Using combined high-speed atomic force microscopy, fluorescence microscopy, and force spectroscopy, we investigated how these active sphingolipids laterally remodel supported bilayers upon photoisomerization, notably in terms of domain area changes, height mismatch, line tension, and membrane piercing. Hereby, we show that the sphingosine-based (Azo-β-Gal-Cer, Azo-SM, Azo-Cer) and phytosphingosine-based (Azo-α-Gal-PhCer, Azo-PhCer) photoswitchable lipids promote a reduction in liquid-ordered microdomain area when in the UV-adapted cis-isoform. In contrast, azo-sphingolipids having tetrahydropyran groups that block H-bonding at the sphingosine backbone (lipids named Azo-THP-SM, Azo-THP-Cer) induce an increase in the liquid-ordered domain area when in cis, accompanied by a major rise in height mismatch and line tension. These changes were fully reversible upon blue light-triggered isomerization of the various lipids back to trans, pinpointing the role of interfacial interactions for the formation of stable liquid-ordered domains.

Synaptophysin (syp) is a major protein of secretory vesicles with four transmembrane domains (TMDs) and a large cytoplasmic C-terminus. Syp has been shown to regulate exocytosis, vesicle cycling, and synaptic plasticity through its C-terminus. However, the roles of its TMDs remain unclear. The TMDs of soluble N-ethylmaleimide-sensitive factor attachment protein receptor (SNARE) proteins are thought to line initial fusion pores, and structural work together with sequence analysis suggest that TMD III of syp may play a similar role. To test this hypothesis, we performed tryptophan scanning experiments of TMD III in chromaffin cells and used amperometry to evaluate fusion pores. In contrast to SNARE TMDs, tryptophan substitutions in syp TMD III had no effect on the flux through initial fusion pores. However, a number of these mutants increased the fraction of kiss-and-run events and decreased the initial fusion pore lifetime. These results indicate that TMD III stabilizes the initial fusion pore and controls the initial choice between kiss and run and full fusion. Late-stage fusion pores were not impacted by TMD III mutations. These results indicate that syp TMD III does not line the initial fusion pore. However, its impact on pore dynamics suggests that it interacts with a SNARE protein implicated as a part of the fusion pore that forms at the onset of exocytosis.

Galectin-3 (Gal-3) is a β-galactosidase-binding protein involved in various biological processes, including neuronal growth and adhesion. The pairing of Gal-3 with ganglioside GM1's pentasaccharide chain at the outer leaflet of the plasma membrane, which triggers downstream cell-signaling cascades, seems to be involved in these processes. A crucial feature of Gal-3 is its ability to form oligomers and supramolecular assemblies that connect various carbohydrate-decorated molecules. Although we know the atomistic structure of Gal-3 bound to small carbohydrate ligands, it remains unclear how Gal-3 binds GM1 in a membrane. Furthermore, the influence of this interaction on Gal-3's structure and oligomeric assembly has to be elucidated. In this study, we used X-ray reflectivity (XR) from a model membrane to determine the structure and surface coverage of Gal-3 bound to a membrane containing GM1. We observed that the carbohydrate recognition domain interacts with GM1's pentasaccharide, while the N-terminal domain is pointed away from the membrane, likely to facilitate protein-protein interactions. In a membrane containing 20 mol % GM1, Gal-3 covered ∼50% of the membrane surface with one Gal-3 molecule bound per 2130 Å2. We used molecular dynamics simulations and Voronoi tessellation algorithms to build an atomistic model of membrane-bound Gal-3, which is supported by the XR results. Overall, this work provides structural information describing how Gal-3 can bind GM1's pentasaccharide chain, a prerequisite for triggering regulatory processes in neuronal growth and adhesion.

The interplay and communication between cells build the foundation of life. Many signaling processes at the cell surface and inside the cell, as well as the cellular function itself, depend on protein-protein interactions and the oligomerization of proteins. In the past, we presented an approach to single out interactions of fluorescently labeled membrane proteins by combining photobleaching and single-molecule microscopy. With this approach, termed "thinning out clusters while conserving the stoichiometry of labeling" (TOCCSL), oligomerization can be detected even at physiologically high surface densities of fluorescently labeled proteins. In TOCCSL, an aperture-restricted region of the plasma membrane is irreversibly photobleached by applying a high-intensity laser pulse. During a recovery time, in which illumination is turned off, nonphotobleached molecules from the nonilluminated area of the plasma membrane re-populate the aperture-restricted region. At the onset of this recovery process, these molecules can be detected as well-separated, diffraction-limited signals and their oligomerization state can be quantified. Here, we used extensive Monte Carlo simulations to provide a theoretical framework for quantitative interpretation of TOCCSL measurements. We determined the influence of experimental parameters and intrinsic characteristics of the investigated system on the outcome of a TOCCSL experiment. We identified the diffraction-affected laser intensity profile and the diffusion of molecules at the aperture edges during photobleaching as major sources of generating partially photobleached oligomers. They are falsely detected as lower-order oligomers and, hence, higher-order oligomers might be prevented from detection. The amount of partially photobleached oligomers that are analyzed depends on the photobleaching and the recovery time, on the mobility of molecules and-for mixed populations of oligomers-on mobility differences between different kinds of oligomers. Moreover, we quantified random colocalizations of molecules after recovery, which are falsely detected as higher-order oligomers.

Membrane fusion is a tool to increase the complexity of model membrane systems. Here, we use silica nanoparticles to fuse liquid-disordered DOPC giant unilamellar vesicles (GUVs) and liquid-ordered DPPC:cholesterol (7:3) GUVs. After fusion, GUVs display large membrane domains as confirmed by fluorescence confocal microscopy. Laurdan spectral imaging of the membrane phases in the fused GUVs shows differences compared with the initial vesicles indicating some lipid redistribution between phase domains as dictated by the tie lines of the phase diagram. Remarkably, using real-time confocal microscopy we were able to record the dynamics of formation of asymmetric membrane domains in hemifused GUVs and detected interleaflet coupling phenomena by which the DOPC-rich liquid-disordered domains in outer monolayer modulates the phase state of the DPPC:cholesterol inner membrane leaflet which transitions from liquid-ordered to liquid-disordered phase. We find that internal membrane stresses generated by membrane asymmetry enhance the efficiency of full fusion compared with our previous studies on symmetric vesicle fusion. Furthermore, under these conditions, the liquid-disordered monolayer dictates the bilayer phase state of asymmetric membrane domains in >90% of observed cases. By comparison to the findings of previous literature, we suggest that the monolayer phase that dominates the bilayer properties could be a mechanoresponsive signaling mechanism sensitive to the local membrane environment.

Previous studies have documented the formation of a heterodimer between the two protein kinases PDK1 and PKCα on a lipid bilayer containing their target lipids. This work investigates the association-dissociation kinetics of this PDK1:PKCα heterodimer. The approach monitors the two-dimensional diffusion of single, membrane-associated PDK1 molecules for diffusivity changes as PKCα molecules bind and unbind. In the absence of PKCα, a membrane-associated PDK1 molecule exhibits high diffusivity (or large diffusion constant, D) because its membrane-contacting PH domain binds the target PIP3 lipid headgroup with little bilayer penetration, yielding minimal frictional drag against the bilayer. In contrast, membrane-associated PKCα contacts the bilayer via its C1A, C1B, and C2 domains, which each bind at least one target lipid with significant bilayer insertion, yielding a large frictional drag and low diffusivity. The present findings reveal that individual fluor-PDK1 molecules freely diffusing on the membrane surface undergo reversible switching between distinct high and low diffusivity states, corresponding to the PDK1 monomer and the PDK1:PKCα heterodimer, respectively. The observed single-molecule diffusion trajectories are converted to step length time courses, then subjected to two-state, hidden Markov modeling and dwell time analysis. The findings reveal that both the PDK1 monomer state and the PDK1:PKCα heterodimer state decay via simple exponential kinetics, yielding estimates of rate constants for state switching in both directions. Notably, the PDK1:PKCα heterodimer has been shown to competitively inhibit PDK1 phosphoactivation of AKT1, and is believed to play a tumor suppressor role by limiting excess activation of the highly oncogenic PDK1/AKT1/mTOR pathway. Thus, the present elucidation of the PDK1:PKCα association-dissociation kinetics has important biological and medical implications. More broadly, the findings illustrate the power of single-molecule diffusion measurements to reveal the kinetics of association-dissociation events in membrane signaling reactions that yield a large change in diffusive mobility.

Molecular dynamics (MD) simulations have become the predominant computational analysis method in membrane biophysics, as this technique is uniquely suited for investigations of complex molecular systems through the relevant physical principles. Owing to continued improvements in scope and performance, the trajectories generated through this approach contain ever-increasing amounts of information, which must be synthesized and simplified in post-analysis using tools that are not only mechanistically insightful but also computationally efficient and highly scalable. Here, we introduce MOSAICS, a self-contained high-performance suite of C++ software tools designed for advanced analyses of lipid bilayer structure and dynamics from MD trajectories. MOSAICS is to our knowledge the most comprehensive software suite of this kind, enabling analysis of a wide array of morphological and kinetic properties, for both simple and complex membranes, irrespective of system size or resolution. Importantly, MOSAICS is designed to provide spatial distributions of all computed quantities, with built-in masking tools, noise filtering, and statistical significance metrics to facilitate quantitative interpretations of the trajectory data; it is also fully parallelized and can therefore leverage the capabilities of supercomputing facilities. Despite its technical sophistication, MOSAICS is user-friendly and requires minimal computational expertise, making it accessible to researchers of all skill levels. This sofware suite can be freely downloaded at https://github.com/MOSAICS-NIH/.

Voltage-sensing phosphatase (VSP) consists of the voltage sensor domain (VSD) similar to that of voltage-gated ion channels and the cytoplasmic phosphatase region with remarkable similarity to the phosphatase and tensin homolog deleted on chromosome 10 (PTEN). Membrane depolarization activates VSD, leading to dephosphorylation of three species of phosphoinositides (phosphatidylinositol phosphates (PIPs)), PI(3,4,5)P3, PI(4,5)P2, and PI(3,4)P2. VSP dephosphorylates 3- and 5-phosphate of PIPs, unlike PTEN, which shows rigid 3-phosphate specificity. In this study, a bioinformatics search showed that some mammals have VSP orthologs with amino acid diversity in the active center motif, Cx5R, which is highly conserved among protein tyrosine phosphatases and PTEN-related phosphatases; lysine next to the active site cysteine in the Cx5R motif was substituted for methionine in VSP orthologs of Tasmanian devil, koala, and prairie deer mouse, and leucine in opossum. Since lysine at the corresponding site in PTEN is known to be critical for enzyme activities, we attempted to address the significance of amino acid diversity among VSP orthologs at this site. K364 was changed to different amino acids in sea squirt VSP (Ci-VSP), and voltage-dependent phosphatase activity in Xenopus oocyte was studied using fluorescent probes for PI(4,5)P2 and PI(3,4)P2. All mutants retained both 5-phosphatase and 3-phosphatase activity, indicating that lysine at this site is dispensable for 3-phosphatase activity, unlike PTEN. Notably, K364M mutant showed increased activity both of 5-phosphatase and 3-phosphatase compared with the wild type (WT). It also showed slower kinetics of voltage sensor motion. Malachite green assay of K364M mutant did not show significant difference of phosphatase activity from WT, suggesting tighter interaction between substrate binding and voltage sensing. Mutation corresponding to K364M in the zebrafish VSP led to enhanced voltage-dependent dephosphorylation of PI(4,5)P2. Further studies will provide clues to understanding of substrate preference in PIPs phosphatases as well as to customization of a molecular tool.

In this article, we review contemporary evidence that GluD receptors are functional ion channels whose depolarizing currents contribute to their biological functions, akin to all other members of the ionotropic glutamate receptor (iGluR) family.

Gene expression states are often stably sustained in cis despite massively disruptive events like DNA replication. This is achieved by on-going enzymatic activity that maintains parts of the DNA in either heterochromatic (packed) or euchromatic (free) states, each of which is stabilized by both positive feedback and bridging interactions between individual nucleosomes. In contrast to condensed matter, however, the dynamics is not only governed by equilibrium binding interactions but is also mediated by enzymes that recognize and act on specific amino acid tails of the nucleosomes. The mechanical result is that some nucleosomes can bind to one another and form tightly packed polymer configurations, whereas others remain unbound and form free, noncompact polymer configurations. Here, we study the consequences of such an asymmetric interaction pattern on the dynamics of epigenetic switching. We develop a 3D polymer model and show that traits associated with epigenetic switching, such as bistability and epigenetic memory, are permitted by such a model. We find, however, that the experimentally observed burst-like nature of some epigenetic switches is difficult to reproduce by this biologically motivated interaction. Instead, the behavior seen in experiments can be explained by introducing partial confinement, which particularly affects the euchromatic regions of the chromosome.

Addition of amyloid seeds to aggregation-prone monomers allows for amyloid fiber growth (elongation) omitting slow nucleation. We here combine Thioflavin T fluorescence (probing formation of amyloids) and solution-state NMR spectroscopy (probing disappearance of monomers) to assess elongation kinetics of the amyloidogenic protein, α-synuclein, for which aggregation is linked to Parkinson's disease. We found that both spectroscopic detection methods give similar kinetic results, which can be fitted by applying double exponential decay functions. When the origin of the two-phase behavior was analyzed by mathematical modeling, parallel paths as well as stop-and-go behavior were excluded as possible explanations. Instead, supported by previous theory, the experimental elongation data reveal distinct kinetic regimes that depend on instantaneous monomer concentration. At low monomer concentrations (toward end of experiments), amyloid growth is limited by conformational changes resulting in β-strand alignments. At the higher monomer concentrations (initial time points of experiments), growth occurs rapidly by incorporating monomers that have not successfully completed the conformational search. The presence of a fast disordered elongation regime at high monomer concentrations agrees with coarse-grained simulations and theory but has not been detected experimentally before. Our results may be related to the wide range of amyloid folds observed.

Lyophilization is promising for tackling degradation during the drying and storage of protein-based drugs. Tardigrade cytosolically abundant heat soluble (CAHS) proteins are necessary and sufficient for desiccation-tolerance in vivo and protein protection in vitro. Hydrated CAHS proteins form coiled-coil-based fine-stranded, cold-setting hydrogels, but the dried protein remains largely uncharacterized. Here, we show that dried CAHS D gels (i.e., aerogels) retain the structural units of their hydrogels, but the details depend on prelyophilization CAHS concentrations. Low concentration samples (<10 g/L) form thin (<0.2 μm) tangled fibrils lacking regular structure on the micron scale. Upon increasing the concentration, the fibers thicken and form slabs comprising the walls of the aerogel pores. These changes in morphology are associated with a loss in disorder and an increase in large β sheets and a decrease in α helices and random coils. This disorder-to-order transition is also seen in hydrated gels as a function of concentration. These results suggest a mechanism for pore formation and indicate that using CAHS proteins as excipients will require attention to initial conditions because the starting concentration impacts the lyophilized product.

The spike protein of the SARS-CoV-2 virus (the causative agent of COVID-19) recognizes the host cell by binding to the peptidase domain (PD) of the extracellular receptor angiotensin-converting enzyme 2 (ACE2). A variety of carbohydrates could be attached to the six asparagines in the PD, resulting in a heterogeneous population of ACE2 glycoforms. Experiments have shown that the binding affinity of glycosylated and deglycosylated ACE2 to the virus is virtually identical. In most cases, the reduction in glycan size correlates with stronger binding, which suggests that volume exclusion, and hence entropic forces, determine the binding affinity. Here, we quantitatively test the entropy-based hypothesis by developing a lattice model for the complex between ACE2 and the SARS-CoV-2 spike protein receptor-binding domain (RBD). Glycans are treated as branched polymers with only volume exclusion, which we justify using all-atom molecular dynamics simulations in explicit water. We show that the experimentally measured changes in the ACE2-RBD dissociation constants for a variety of engineered ACE2 glycoforms are in reasonable agreement with our theory, thus supporting our hypothesis. However, a quantitative recovery of all the experimental data could require weak attractive interactions.

We studied the mechanical leaflet coupling of prototypic mammalian plasma membranes using neutron spin-echo spectroscopy. In particular, we examined a series of asymmetric phospholipid vesicles with phosphatidylcholine and sphingomyelin enriched in the outer leaflet and inner leaflets composed of phosphatidylethanolamine/phosphatidylserine mixtures. The bending rigidities of most asymmetric membranes were anomalously high, exceeding even those of symmetric membranes formed from their cognate leaflets. Only asymmetric vesicles with outer leaflets enriched in sphingolipid displayed bending rigidities in conformity with these symmetric controls. We performed complementary small-angle neutron and x-ray experiments on the same vesicles to examine possible links to structural coupling mechanisms, which would show up in corresponding changes in membrane thickness. In addition, we estimated differential stress between leaflets originating either from a mismatch of their lateral areas or spontaneous curvatures. However, no correlation with asymmetry-induced membrane stiffening was observed. To reconcile our findings, we speculate that an asymmetric distribution of charged or H-bond forming lipids may induce an intraleaflet coupling, which increases the weight of hard undulatory modes of membrane fluctuations and hence the overall membrane stiffness.

Transmembrane assemblies of the peptaibol alamethicin (ALM) are among the most extensively studied ion channels not only because of their antimicrobial activity but also as models for channel structure and aggregation. In this study, several oligomeric states of ALM are investigated with molecular dynamics simulations to establish properties of the channel and obtain free energy profiles for ion transport and the corresponding values of conductance. The hexamer, heptamer, and octamer of ALM in phospholipid membrane are found to be stable but highly dynamic in barrel-stave structures, with calculated conductance equal to 18, 195, and 1270 pS, respectively, in 1 M KCl ion solution. The corresponding free energy profiles, reported for the first time, are reconstructed from simulations at applied voltage of 200 mV with the aid of the electrodiffusion model both with and without the knowledge of diffusivity. The calculated free energy barriers are equal to 2.5, 1.5, and 0.5 kcal/mol for K+ and 4.0, 2.2, and 1.5 kcal/mol for Cl-, for hexamer, heptamer, and octamer, respectively. The calculated conductance and the ratio between conductance in consecutive states are in good agreement with those measured experimentally. This suggests that the hexamer is the lowest conducting state, with measured conductance equal to 19 pS. The selectivity of K+ over Cl- is calculated as 1.5 and 2.3 for the octameric and heptameric channels, close to the selectivity measured for high-conductance states. Selectivity increases to 13 in the hexameric channel in which the narrowest Gln7 site has a pore radius of only ∼1.6 Å, again in accord with experiment. A good agreement found between calculated and measured conductance through a hexamer templated on cyclodextrin lands additional support for the results of our simulations, and the comparison with ALM reveals the dependence of conductance on the nature of phospholipid membrane.

The cell elastic modulus (Ec) is widely used as the mechanics-based marker to analyze the biological effects of substrates on cells. However, the employment of the Hertz model to extract the apparent Ec can cause errors due to the disobedience of the small deformation assumption and the infinite half-space assumption, as well as an inability to deduct the deformation of the substrate. So far, no model can effectively solve the errors caused by the above-mentioned aspects simultaneously. In response to this, herein, we propose an active learning model to extract Ec. The numerical calculation with finite element suggests the good prediction accuracy of the model. The indentation experiments on both hydrogel and cell indicate that the established model can efficiently reduce the error caused by the method of extracting Ec. The application of this model may facilitate our understanding about the role of Ec in correlating the stiffness of substrate and the biological behavior of cell.

Single-molecule data are of great significance in biology, chemistry, and medicine. However, new experimental tools to characterize, in a multiplexed manner, protein bond rupture under force are still needed. Acoustic force spectroscopy is an emerging manipulation technique which generates acoustic waves to apply force in parallel on multiple microbeads tethered to a surface. We here exploit this configuration in combination with the recently developed modular junctured-DNA scaffold that has been designed to study protein-protein interactions at the single-molecule level. By applying repetitive constant force steps on the FKBP12-rapamycin-FRB complex, we measure its unbinding kinetics under force at the single-bond level. Special efforts are made in analyzing the data to identify potential pitfalls. We propose a calibration method allowing in situ force determination during the course of the unbinding measurement. We compare our results with well-established techniques, such as magnetic tweezers, to ensure their accuracy. We also apply our strategy to study the force-dependent rupture of a single-domain antibody with its antigen. Overall, we get a good agreement with the published parameters that have been obtained at zero force and population level. Thus, our technique offers single-molecule precision for multiplexed measurements of interactions of biotechnological and medical interest.

The physical stability of peptide-based drugs is of great interest to the pharmaceutical industry. Glucagon-like peptide 1 (GLP-1) is a 31-amino acid peptide hormone, the analogs of which are frequently used in the treatment of type 2 diabetes. We investigated the physical stability of GLP-1 and its C-terminal amide derivative, GLP-1-Am, both of which aggregate into amyloid fibrils. While off-pathway oligomers have been proposed to explain the unusual aggregation kinetics observed previously for GLP-1 under specific conditions, these oligomers have not been studied in any detail. Such states are important as they may represent potential sources of cytotoxicity and immunogenicity. Here, we identified and isolated stable, low-molecular-weight oligomers of GLP-1 and GLP-1-Am, using size-exclusion chromatography. Under the conditions studied, isolated oligomers were shown to be resistant to fibrillation or dissociation. These oligomers contain between two and five polypeptide chains and they have a highly disordered structure as indicated by a variety of spectroscopic techniques. They are highly stable with respect to time, temperature, or agitation despite their noncovalent character, which was established using liquid chromatography-mass spectrometry and sodium dodecyl sulfate-polyacrylamide gel electrophoresis. These results provide evidence of stable, low-molecular-weight oligomers that are formed by an off-pathway mechanism which competes with amyloid fibril formation.

The cytoplasmic proteins of some halophilic organisms remain stable and functional at multimolar concentrations of KCl, i.e., under conditions that most mesophilic proteins cannot withstand. Their stability arises from their unusual amino acid composition. The most dramatic difference between halophilic and mesophilic proteins is that the former are rich in acidic amino acids. It has been proposed that one of the evolutionary driving forces for this difference is the occurrence of synergistic interactions between multiple acidic amino acids at the surface of the protein, the potassium cations in solution, and water. We investigate this possibility with molecular dynamics simulations, using high-quality force fields for the protein-water, protein-ion, and ion-ion interactions. We create a rigorous thermodynamic definition of interactions between acidic amino acids on proteins that can be used to distinguish between synergistic, noninteracting and interfering interactions. Our results demonstrate that synergistic interactions between neighboring acidic amino acids in halophilic proteins are frequent at multimolar KCl concentration. Synergistic interactions have an electrostatic origin, and are associated with stronger water-to-carboxylate hydrogen bonds than for acidic amino acids without synergistic interactions. Synergistic interactions are not observed in minimal systems of carboxylates, indicating that the protein environment is critical for their emergence. Our results demonstrate that synergistic interactions are neither associated with rigid amino acid orientations nor with highly structured and slow moving water networks, as had been originally proposed. Moreover, synergistic interactions can also be found in unfolded protein conformations. However, because these conformations are only a small subset of the unfolded state ensemble, synergistic interactions should contribute to the net stabilization of the folded state.

The cell-cell adhesion cadherin-catenin complexes recruit vinculin to the adherens junction (AJ) to modulate the mechanical couplings between neighboring cells. However, it is unclear how vinculin influences the AJ structure and function. Here, we identified two patches of salt bridges that lock vinculin in the head-tail autoinhibited conformation and reconstituted the full-length vinculin activation mimetics bound to the cadherin-catenin complex. The cadherin-catenin-vinculin complex contains multiple disordered linkers and is highly dynamic, which poses a challenge for structural studies. We determined the ensemble conformation of this complex using small-angle x-ray and selective deuteration/contrast variation small-angle neutron scattering. In the complex, both α-catenin and vinculin adopt an ensemble of flexible conformations, but vinculin has fully open conformations with the vinculin head and actin-binding tail domains well separated from each other. F-actin binding experiments show that the cadherin-catenin-vinculin complex binds and bundles F-actin. However, when the vinculin actin-binding domain is removed from the complex, only a minor fraction of the complex binds to F-actin. The results show that the dynamic cadherin-catenin-vinculin complex employs vinculin as the primary F-actin binding mode to strengthen AJ-cytoskeleton interactions.

Erythrophagocytosis occurring in the spleen is a critical process for removing senescent and diseased red blood cells (RBCs) from the microcirculation. Although some progress has been made in understanding how the biological signaling pathways mediate the phagocytic processes, the role of the biophysical interaction between RBCs and macrophages, particularly under pathological conditions such as sickle cell disease, has not been adequately studied. Here, we combine computational simulations with microfluidic experiments to quantify RBC-macrophage adhesion dynamics under flow conditions comparable to those in the red pulp of the spleen. We also investigate the RBC-macrophage interaction under normoxic and hypoxic conditions. First, we calibrate key model parameters in the adhesion model using microfluidic experiments for normal and sickle RBCs under normoxia and hypoxia. We then study the adhesion dynamics between the RBC and the macrophage. Our simulation illustrates three typical adhesion states, each characterized by a distinct dynamic motion of the RBCs, namely firm adhesion, flipping adhesion, and no adhesion (either due to no contact with macrophages or detachment from the macrophages). We also track the number of bonds formed when RBCs and macrophages are in contact, as well as the contact area between the two interacting cells, providing mechanistic explanations for the three adhesion states observed in the simulations and microfluidic experiments. Furthermore, we quantify, for the first time to our knowledge, the adhesive forces between RBCs (normal and sickle) and macrophages under different oxygenated conditions. Our results show that the adhesive forces between normal cells and macrophages under normoxia are in the range of 33-58 pN and 53-92 pN for sickle cells under normoxia and 155-170 pN for sickle cells under hypoxia. Taken together, our microfluidic and simulation results improve our understanding of the biophysical interaction between RBCs and macrophages in sickle cell disease and provide a solid foundation for investigating the filtration function of the splenic macrophages under physiological and pathological conditions.

Could the phenomenon of catch bonding-force-strengthened cellular adhesion-play a role in sickle cell disease, where abnormal red blood cell (RBC) adhesion obstructs blood flow? Here, we investigate the dynamics of sickle RBCs adhering to a surface functionalized with the protein laminin (a component of the extracellular matrix around blood vessels) under physiologically relevant microscale flow. First, using total internal reflectance microscopy we characterize the spatial fluctuations of the RBC membrane above the laminin surface before detachment. The complex dynamics we observe suggest the possibility of catch bonding, where the mean detachment time of the cell from the surface initially increases to a maximum and then decreases as a function of shear force. We next conduct a series of shear-induced detachment experiments on blood samples from 25 sickle cell disease patients, quantifying the number and duration of adhered cells under both sudden force jumps and linear force ramps. The experiments reveal that a subset of patients does indeed exhibit catch bonding. By fitting the data to a theoretical model of the bond dynamics, we can extract the mean bond lifetime versus force for each patient. The results show a striking heterogeneity among patients, both in terms of the qualitative behavior (whether or not there is catch bonding) and in the magnitudes of the lifetimes. Patients with large bond lifetimes at physiological forces are more likely to have certain adverse clinical features, like a diagnosis of pulmonary arterial hypertension and intracardiac shunts. By introducing an in vitro platform for fully characterizing RBC-laminin adhesion dynamics, our approach could contribute to the development of patient-specific antiadhesive therapies for sickle cell disease. The experimental setup is also easily generalizable to studying adhesion dynamics in other cell types, for example, leukocytes or cancer cells, and can incorporate disease-relevant environmental conditions like oxygen deprivation.

Mechanisms behind the fluctuations in the ionic current through single acetylcholine receptor (AChR) channels have remained elusive. In a recent study of muscle AChR we showed that mutation of a conserved intramembrane salt bridge in the β- and δ-subunits markedly increased fluctuations in the open channel current that extended from low to high frequency. Here, we show that extracellular divalent cations reduce the high-frequency fluctuations and increase the low-frequency fluctuations. The low-frequency fluctuations are shown to arise from steps between two current levels, with the ratio of the time at each level changing e-fold for a 70 mV increase in membrane potential, indicating modulation by a charged element within the membrane field. Increasing the charge on the ion selectivity filter biases the ratio of current levels equivalent to a 50 mV increase in membrane potential but does not alter the voltage dependence of the ratio. The magnitudes of the voltage dependence and voltage bias allow estimates of the distance between the ion selectivity filter and the voltage-sensing element. Studies with either calcium or magnesium show that the two divalent cations synergize to increase the low-frequency fluctuations, whereas they act independently to decrease the high-frequency fluctuations, indicating multiple divalent cation binding sites. Molecular dynamics simulations applied to the structure of the Torpedo AChR reveal that mutation of the salt bridge alters the equilibrium positions and dynamics of residues local to the site of the mutation and within the adjacent ion selectivity filter in a calcium-dependent manner. Thus, disruption of a conserved intramembrane salt bridge in the muscle AChR induces fluctuations in open channel current that are sensitive to divalent cation binding at multiple sites and modulated by a charged element within the membrane field.