The formation of biomolecular condensates has emerged as a new biophysical principle for subcellular compartmentalization within cells to facilitate the spatiotemporal regulation of a multitude of complex biomolecular reactions. In this Research Highlight, we summarize the findings that were published in Biophysical Journal during the past two years (2021 and 2022). These papers provided biophysical insights into the formation of biomolecular condensates via phase separation of proteins with or without nucleic acids.

Antibody-opsonized bacteria interact with Fc receptors in macrophages and trigger signaling cascades, which induce phagocytosis. The signaling pathways ultimately lead to actin polymerization that induces the protrusion of the membrane around the bacterium until it is completely engulfed. Although many proteins involved in the phagocytic cup formation have already been identified, it is still unclear how far the initial stimulus created by the bacterium-cell contact propagates in the cell. We hypothesize that this spreading distance is closely related to the spatial resolution limit of phagocytosis, the smallest distance in which two stimuli can be differentiated. Here, we probe this resolution limit by using holographic optical tweezers to attach pairs of immunoglobulin G-coated polystyrene microparticles (as models for opsonized bacteria) to murine macrophages in distances ranging from zero to several micrometers. By using 2-μm-sized particles, we found that the particles can be internalized jointly into one phagosome if they are attached to the cell very close together, but that they are taken up separately if they are attached far from each other. To explain this, we developed a model of the signaling process, which predicts the probabilities for separate uptake for different particle sizes and distances using cellular parameters including the average receptor distance. We tested the model by measuring the separate uptake probabilities for particles with a diameter of 1 to 3 μm and for cells with reduced numbers of Fcγ receptors and found very good agreement. Our model shows that the phagocytic uptake behavior can be explained by assuming an effective phagocytic signaling range of about 500 nm. Interestingly, this value corresponds to the lower size limit of phagocytosis. Our work provides quantitative access to spatial parameters of cellular signaling during phagocytosis and thereby contributes to a more quantitative understanding of cellular information processing.

DNA strands have to sample numerous states to find the alignment that maximizes Watson-Crick-Franklin base pairing. This process depends strongly on sequence, which affects the stability of the native duplex as well as the prevalence of non-native inter- and intramolecular helices. We present a theory that describes DNA hybridization as a three-stage process: diffusion, registry search, and zipping. We find that non-specific binding affects each of these stages in different ways. Mis-registered intermolecular binding in the registry search stage helps DNA strands sample different alignments and accelerates the hybridization rate. Non-native intramolecular structure affects all three stages by rendering portions of the molecule inert to intermolecular association, limiting mis-registered alignments to be sampled, and impeding the zipping process. Once in-register base pairs are formed, the stability of the native structure is important to hold the molecules together long enough for non-native contacts to break.

Small regulatory RNA molecules such as microRNA modulate gene expression through inhibiting the translation of messenger RNA (mRNA). Such posttranscriptional regulation has been recently hypothesized to reduce the stochastic variability of gene expression around average levels. Here, we quantify noise in stochastic gene expression models with and without such regulation. Our results suggest that silencing mRNA posttranscriptionally will always increase, rather than decrease, gene expression noise when the silencing of mRNA also increases its degradation, as is expected for microRNA interactions with mRNA. In that regime, we also find that silencing mRNA generally reduces the fidelity of signal transmission from deterministically varying upstream factors to protein levels. These findings suggest that microRNA binding to mRNA does not generically confer precision to protein expression.

Molecular dynamics (MD) simulation of biologically relevant processes at realistic time scale and atomistic precision is generally limited by prohibitively large computational cost, due to its restriction of using an ultrashort integration time step (1-2 fs). A popular numerical recipe to reduce the associated computational burden is adopting schemes that would allow relatively longer-time-step for MD propagation. Here, we explore the perceived potential of one of the most frequently used long-time-step protocols, namely the hydrogen mass repartitioning (HMR) approach, in alleviating the computational overhead associated with simulation of the kinetic process of protein-ligand recognition events. By repartitioning the mass of heavier atoms to their linked hydrogen atoms, HMR leverages around twofold longer time step than regular simulation, holding promise of significant performance boost. However, our probe into direct simulation of the protein-ligand recognition event, one of the computationally most challenging processes, shows that long-time-step HMR MD simulations do not necessarily translate to a computationally affordable solution. Our investigations spanning cumulative 176 μs in three independent proteins (T4 lysozyme, sensor domain of MopR, and galectin-3) show that long-time-step HMR-based MD simulations can catch the ligand in its act of recognizing the native cavity. But, as a major caveat, the ligand is found to require significantly longer time to identify buried native protein cavity in an HMR MD simulation than regular simulation, thereby defeating the purpose of its usage for performance upgrade. A molecular analysis shows that the longer time required by a ligand to recognize the protein in HMR is rooted in faster diffusion of the ligand, which reduces the survival probability of decisive on-pathway metastable intermediates, thereby slowing down the eventual recognition process at the native cavity. Together, the investigation stresses careful assessment of pitfalls of long-time-step algorithms before attempting to utilize them for higher performance for biomolecular recognition simulations.

Cell motility on flat substrates exhibits coexisting steady and oscillatory morphodynamics, the biphasic adhesion-velocity relation, and the universal correlation between speed and persistence (UCSP) as simultaneous observations common to many cell types. Their universality and concurrency suggest a unifying mechanism causing all three of them. Stick-slip models for cells on one-dimensional lanes suggest multistability to arise from the nonlinear friction of retrograde flow. This study suggests a mechanical mechanism controlled by integrin signaling on the basis of a biophysical model and analysis of trajectories of MDA-MB-231 cells on fibronectin lanes, which additionally explains the constitutive relations. The experiments exhibit cells with steady or oscillatory morphodynamics and either spread or moving with spontaneous transitions between the dynamic regimes, spread and moving, and spontaneous direction reversals. Our biophysical model is based on the force balance at the protrusion edge, the noisy clutch of retrograde flow, and a response function of friction and membrane drag to integrin signaling. The theory reproduces the experimentally observed cell states, characteristics of oscillations, and state probabilities. Analysis of experiments with the biophysical model establishes a stick-slip oscillation mechanism, and explains multistability of cell states and the statistics of state transitions. It suggests protrusion competition to cause direction reversal events, the statistics of which explain the UCSP. The effect of integrin signaling on drag and friction explains the adhesion-velocity relation and cell behavior at fibronectin density steps. The dynamics of our mechanism are nonlinear flow mechanics driven by F-actin polymerization and shaped by the noisy clutch of retrograde flow friction, protrusion competition via membrane tension, and drag forces. Integrin signaling controls the parameters of the mechanical system.

Members of the fatty acid binding protein (FABP) family function as intracellular transporters of long-chain fatty acids and other hydrophobic molecules to different cellular compartments. Brain FABP (FABP7) exhibits ligand-directed differences in cellular transport. For example, when FABP7 binds to docosahexaenoic acid (DHA), the complex relocates to the nucleus and influences transcriptional activity, whereas FABP7 bound with monosaturated fatty acids remains in the cytosol. Preferential binding of FABP7 to polyunsaturated fatty acids like DHA has been previously observed and is thought to play a role in differential localization. However, we find that at 37°C, FABP7 does not display strong selectivity, suggesting that the conformational ensemble of FABP7 and its perturbation upon binding may be important. We use molecular dynamics simulations, NMR, and a variety of biophysical techniques to better understand the conformational ensemble of FABP7, how it is perturbed by fatty acid binding, and how this may be related to ligand-directed transport. We find that FABP7 has high degree of conformational heterogeneity that is substantially reduced upon ligand binding. We also observe substantial heterogeneity in ligand binding poses, which is consistent with our finding that ligand binding is resistant to mutations in key polar residues in the binding pocket. Our NMR experiments show that DHA binding leads to chemical shift perturbations in residues near the nuclear localization signal, which may point toward a mechanism of differential transport.

Islets of Langerhans operate as multicellular networks in which several hundred β cells work in synchrony to produce secretory pulses of insulin, a hormone crucial for controlling metabolic homeostasis. Their collective rhythmic activity is facilitated by gap junctional coupling and affected by their functional heterogeneity, but the details of this robust and coordinated behavior are still not fully understood. Recent advances in multicellular imaging and optogenetic and photopharmacological strategies, as well as in network science, have led to the discovery of specialized β cell subpopulations that were suggested to critically determine the collective dynamics in the islets. In particular hubs, i.e., β cells with many functional connections, are believed to significantly enhance communication capacities of the intercellular network and facilitate an efficient spreading of intercellular Ca2+ waves, whereas wave-initiator cells trigger intercellular signals in their cohorts. Here, we determined Ca2+ signaling characteristics of these two β cell subpopulations and the relationship between them by means of functional multicellular Ca2+ imaging in mouse pancreatic tissue slices in combination with methods of complex network theory. We constructed network layers based on individual Ca2+ waves to identify wave initiators, and functional correlation-based networks to detect hubs. We found that both cell types exhibit a higher-than-average active time under both physiological and supraphysiological glucose concentrations, but also that they differ significantly in many other functional characteristics. Specifically, Ca2+ oscillations in hubs are more regular, and their role appears to be much more stable over time than for initiator cells. Moreover, in contrast to wave initiators, hubs transmit intercellular signals faster than other cells, which implies a stronger intercellular coupling. Our research indicates that hubs and wave-initiator cell subpopulations are both natural features of healthy pancreatic islets, but their functional roles in principle do not overlap and should thus not be considered equal.

In the mammalian cochlea, each longitudinal position of the basilar membrane (BM) has a nonlinear vibratory response in a limited frequency range around the location-dependent frequency of maximum response, known as the best frequency (BF). This nonlinear response arises from the electromechanical feedback from outer hair cells (OHCs). However, recent in vivo measurements have demonstrated that the mechanical response of other organ of Corti (OoC) structures, such as the reticular lamina (RL), and the electrical response of OHCs (measured in the local cochlear microphonic [LCM]) are nonlinear even at frequencies significantly below BF. In this work, a physiologically motivated model of the gerbil cochlea is used to demonstrate that the source of this discrepancy between the frequency range of the BM, RL, and LCM nonlinearities is greater compliance in the structures at the top of the OHCs. The predicted responses of the BM, RL, and LCM to pure tone and two-tone stimuli are shown to be in line with experimental evidence. Simulations then demonstrate that the sub-BF nonlinearity in the RL requires the structures at the top of the OHCs to be significantly more compliant than the BM. This same condition is also necessary for "optimal" gain near BF, i.e., high amplification that is in line with the experiment. This demonstrates that the conditions for OHCs to operate optimally at BF inevitably yield nonlinearity of the RL response over a broad frequency range.

The Young-Laplace, Kelvin, and Gibbs-Thomson equations form a cornerstone of colloidal and surface sciences and have found successful applications in many subfields of physics, chemistry, and biology. The Gibbs-Thomson effect, for example, predicts that small crystals are in equilibrium with their liquid melt at a lower temperature than large crystals and the positive interfacial energy increases the energy required to form small particles with a high curvature interface. In cases of liquids contained within porous media (confined geometry), the effect indicates decreasing the freezing/melting temperatures and the increment of the temperature is inversely proportional to the pore size. These phenomena can be reformulated for Gaussian maps of macromolecules and can be asked the following question: can one use the equations for predicting the melting temperature and shape of polymer chains in confined geometries? The answer is no, mainly because macromolecules form highly curved surfaces (Gaussian maps), and the equations hold only for simple geometries (sphere, plane, or cylinder). Here, we present general Young-Laplace, Kelvin, and Gibbs-Thomson equations for arbitrarily curved surfaces and apply them to predict temperature distribution on a few protein surfaces. Also, after increased interest toward liquid/liquid phase separation in biology, we derive generic Ostwald ripening and show that for shape-changing condensates, instead of a monotonic growing mechanism, a variety of processes are possible. Due to the generality of equations, we clarify that at appropriate internal/external pressure conditions systems, bounded by surfaces, may adopt any shape and thermal stability is strongly influenced by the geometries of confined spaces.

Necroptosis is a form of regulated cell death associated with degenerative disorders, autoimmune and inflammatory diseases, and cancer. To better understand the biochemical mechanisms regulating necroptosis, we constructed a detailed computational model of tumor necrosis factor-induced necroptosis based on known molecular interactions from the literature. Intracellular protein levels, used as model inputs, were quantified using label-free mass spectrometry, and the model was calibrated using Bayesian parameter inference to experimental protein time course data from a well-established necroptosis-executing cell line. The calibrated model reproduced the dynamics of phosphorylated mixed lineage kinase domain-like protein, an established necroptosis reporter. A subsequent dynamical systems analysis identified four distinct modes of necroptosis signal execution, distinguished by rate constant values and the roles of the RIP1 deubiquitinating enzymes A20 and CYLD. In one case, A20 and CYLD both contribute to RIP1 deubiquitination, in another RIP1 deubiquitination is driven exclusively by CYLD, and in two modes either A20 or CYLD acts as the driver with the other enzyme, counterintuitively, inhibiting necroptosis. We also performed sensitivity analyses of initial protein concentrations and rate constants to identify potential targets for modulating necroptosis sensitivity within each mode. We conclude by associating numerous contrasting and, in some cases, counterintuitive experimental results reported in the literature with one or more of the model-predicted modes of necroptosis execution. In all, we demonstrate that a consensus pathway model of tumor necrosis factor-induced necroptosis can provide insights into unresolved controversies regarding the molecular mechanisms driving necroptosis execution in numerous cell types under different experimental conditions.

In mammalian cells, all-trans farnesol, a 15-carbon isoprenol, is a product of the mevalonate pathway. It is the natural substrate of alcohol dehydrogenase and a substrate for CYP2E1, two enzymes implicated in ethanol metabolism. Studies have shown that farnesol is present in the human brain and inhibits voltage-gated Ca2+ channels at much lower concentrations than ethanol. Here we show that farnesol modulates the activity of γ-aminobutyric acid type A receptors (GABAARs), some of which also mediate the sedative activity of ethanol. Electrophysiology experiments performed in HEK cells expressing human α1β3γ2 or α6β3γ2 GABAARs revealed that farnesol increased chloride currents through positive allosteric modulation of these receptors and showed dependence on both the alcoholic functional group of farnesol and the length of the alkyl chain for activity. In silico studies using long-timescale unbiased all-atom molecular dynamics (MD) simulations of the human α1β3γ2 GABAA receptors revealed that farnesol modulates the channel by directly binding to the transmembrane neurosteroid-binding site, after partitioning into the surrounding membrane and reaching the receptor by lateral diffusion. Channel activation by farnesol was further characterized by several structural and dynamic variables, such as global twisting of the receptor's extracellular domain, tilting of the transmembrane M2 helices, radius, cross-sectional area, hydration status, and electrostatic potential of the channel pore. Our results expand the pharmacological activities of farnesol to yet another class of ion channels implicated in neurotransmission, thus providing a novel path for understanding and treatment of diseases involving GABAA receptor dysfunction.

The cytoplasm is a complex, crowded, actively driven environment whose biophysical characteristics modulate critical cellular processes such as cytoskeletal dynamics, phase separation, and stem cell fate. Little is known about the variance in these cytoplasmic properties. Here, we employed particle-tracking nanorheology on genetically encoded multimeric 40 nm nanoparticles (GEMs) to measure diffusion within the cytoplasm of individual fission yeast (Schizosaccharomyces pombe) cellscells. We found that the apparent diffusion coefficients of individual GEM particles varied over a 400-fold range, while the differences in average particle diffusivity among individual cells spanned a 10-fold range. To determine the origin of this heterogeneity, we developed a Doppelgänger simulation approach that uses stochastic simulations of GEM diffusion that replicate the experimental statistics on a particle-by-particle basis, such that each experimental track and cell had a one-to-one correspondence with their simulated counterpart. These simulations showed that the large intra- and inter-cellular variations in diffusivity could not be explained by experimental variability but could only be reproduced with stochastic models that assume a wide intra- and inter-cellular variation in cytoplasmic viscosity. The simulation combining intra- and inter-cellular variation in viscosity also predicted weak nonergodicity in GEM diffusion, consistent with the experimental data. To probe the origin of this variation, we found that the variance in GEM diffusivity was largely independent of factors such as temperature, the actin and microtubule cytoskeletons, cell-cyle stage, and spatial locations, but was magnified by hyperosmotic shocks. Taken together, our results provide a striking demonstration that the cytoplasm is not "well-mixed" but represents a highly heterogeneous environment in which subcellular components at the 40 nm size scale experience dramatically different effective viscosities within an individual cell, as well as in different cells in a genetically identical population. These findings carry significant implications for the origins and regulation of biological noise at cellular and subcellular levels.

For the past 50 years, evidence for the existence of functional lipid domains has been steadily accumulating. Although the notion of functional lipid domains, also known as "lipid rafts," is now widely accepted, this was not always the case. This ambiguity surrounding lipid domains could be partly attributed to the fact that they are highly dynamic, nanoscopic structures. Since most commonly used techniques are sensitive to microscale structural features, it is therefore, not surprising that it took some time to reach a consensus regarding their existence. In this review article, we will discuss studies that have used techniques that are inherently sensitive to nanoscopic structural features (i.e., neutron scatting, nuclear magnetic resonance, and Förster resonance energy transfer). We will also mention techniques that may be of use in the future (i.e., cryoelectron microscopy, droplet interface bilayers, inelastic x-ray scattering, and neutron reflectometry), which can further our understanding of the different and unique physicochemical properties of nanoscopic lipid domains.

The cationic lipid 1,2-dioleoyl-3-trimethylammonium propane (DOTAP) is one of the original synthetic cationic lipids used for the liposomal transfection of oligonucleotides in gene therapy. The key structural element of DOTAP is its quaternary ammonium headgroup that is responsible for interactions with both nucleic acids and target cell membranes. Because these interactions are fundamental to the design of a major class of transfection lipids, it is important to understand the structure of DOTAP and how it interacts with halide counterions. Here, we use x-ray and neutron diffraction techniques to examine the structure of DOTAP and how chloride (Cl-) and iodide (I-) counterions alter the hydration properties of the DOTAP headgroup. A problem of particular interest is the poor solubility of DOTAP/I- in water solutions. Our results show that the poor solubility results from very tight binding of the I- counterion to the headgroup and the consequent expulsion of water. The structural principles we report here are important for assessing the suitability of DOTAP and its quaternary ammonium derivatives for transfection.

Unfavorable lipid-lipid pairwise interactions between HiTm and LowTm lipids drive liquid-disordered (Ld) + liquid-ordered (Lo) phase separation. Large size of phase domains is opposed by lipid dipole repulsions, which are more significant compared with the pairwise interactions for naturally abundant LowTm lipids such as palmitoyl oleoyl phosphatidylcholine. During the nano-to-macro domain size transition, no lipid phase transition occurs, and measured properties of Ld + Lo nanodomains are found to be essentially the same as those of macrodomains. Use of macrodomains in mixtures to model cell plasma membranes (PM) is helpful, enabling study by optical microscopy. Use of asymmetric giant unilamellar vesicles to model a PM reveals that ordered phase domains in one leaflet induce ordered domains in an otherwise uniform phase in the apposing leaflet that models a cytoplasmic leaflet. Because macro and nano phase properties are so similar, we conclude that a cell PM that has nano-scale Ld + Lo phase domains in the exoplasmic leaflet is likely to induce nano-scale ordered domains in the cytoplasmic leaflet.

Long-chain polyunsaturated fatty acids (PUFAs) are prone to nonenzymatic oxidation in response to differing environmental stressors and endogenous cellular sources. There is increasing evidence that phospholipids containing oxidized PUFA acyl chains control the inflammatory response. However, the underlying mechanism(s) of action by which oxidized PUFAs exert their functional effects remain unclear. Herein, we tested the hypothesis that replacement of 1-palmitoyl-2-arachidonyl-phosphatidylcholine (PAPC) with oxidized 1-palmitoyl-2-arachidonyl-phosphatidylcholine (oxPAPC) regulates membrane architecture. Specifically, with solid-state 2H NMR of biomimetic membranes, we investigated how substituting oxPAPC for PAPC modulates the molecular organization of liquid-ordered (Lo) domains. 2H NMR spectra for bilayer mixtures of 1,2-dipalmitoylphosphatidylcholine-d62 (an analog of DPPC deuterated throughout sn-1 and -2 chains) and cholesterol to which PAPC or oxPAPC was added revealed that replacing PAPC with oxPAPC disrupted molecular organization, indicating that oxPAPC does not mix favorably in a tightly packed Lo phase. Furthermore, unlike PAPC, adding oxPAPC stabilized 1,2-dipalmitoylphosphatidylcholine-d6-rich/cholesterol-rich Lo domains formed in mixtures with 1,2-dioleoylphosphatidylcholine while decreasing the molecular order within 1,2-dioleoylphosphatidylcholine-rich liquid-disordered regions of the membrane. Collectively, these results suggest a mechanism in which oxPAPC stabilizes Lo domains-by disordering the surrounding liquid-disordered region. Changes in the structure, and thereby functionality, of Lo domains may underly regulation of plasma membrane-based inflammatory signaling by oxPAPC.

Bilayer membranes composed of cholesterol and phospholipids exhibit diverse forms of nonideal mixing. In particular, many previous studies document macroscopic liquid-liquid phase separation as well as nanometer-scale heterogeneity in membranes of phosphatidylcholine (PC) lipids and cholesterol. Here, we present experimental measurements of cholesterol chemical potential (μc) in binary membranes containing dioleoyl PC (DOPC), 1-palmitoyl-2-oleoyl PC (POPC), or dipalmitoyl PC (DPPC), and in ternary membranes of DOPC and DPPC, referenced to crystalline cholesterol. μc is the thermodynamic quantity that dictates the availability of cholesterol to bind other factors, and notably must be equal between coexisting phases of a phase separated mixture. It is simply related to concentration under conditions of ideal mixing, but is far from ideal for the majority of lipid mixtures investigated here. Measurements of μc can vary with phospholipid composition by 1.5 kBT at constant cholesterol mole fraction implying a more than fivefold change in its availability for binding receptors and other reactions. Experimental measurements are fit to thermodynamic models including cholesterol-DPPC complexes or pairwise interactions between lipid species to provide intuition about the magnitude of interactions. These findings reinforce that μc depends on membrane composition overall, suggesting avenues for cells to alter the availability of cholesterol without varying cholesterol concentration.

Dopamine (DA) is a neurotransmitter that also acts as a neuromodulator, with both functions being essential to brain function. Here, we present the first experimental measurement of DA location in lipid bilayers using x-ray diffuse scattering, solid-state deuterium NMR, and electron paramagnetic resonance. We find that the association of DA with lipid headgroups as seen in electron density profiles leads to an increase of intermembrane repulsion most likely due to electrostatic charging. DA location in the lipid headgroup region also leads to an increase of the cross-sectional area per lipid without affecting the bending rigidity significantly. The order parameters measured by solid-state deuterium NMR decrease in the presence of DA for the acyl chains of PC and PS lipids, consistent with an increase in the area per lipid due to DA. Most importantly, these results support the hypothesis that three-dimensional diffusion of DA to target membranes could be followed by relatively more efficient two-dimensional diffusion to receptors within those membranes.

Lipid membrane viscosity is critical to biological function. Bacterial cells grown in different environments alter their lipid composition in order to maintain a specific viscosity, and membrane viscosity has been linked to the rate of cellular respiration. To understand the factors that determine the viscosity of a membrane, we ran equilibrium all-atom simulations of single component lipid bilayers and calculated their viscosities. The viscosity was calculated via a Green-Kubo relation, with the stress-tensor autocorrelation function modeled by a stretched exponential function. By simulating a series of lipids at different temperatures, we establish the dependence of viscosity on several aspects of lipid chemistry, including hydrocarbon chain length, unsaturation, and backbone structure. Sphingomyelin is found to have a remarkably high viscosity, roughly 20 times that of DPPC. Furthermore, we find that inclusion of the entire range of the dispersion interaction increases viscosity by up to 140%. The simulated viscosities are similar to experimental values obtained from the rotational dynamics of small chromophores and from the diffusion of integral membrane proteins but significantly lower than recent measurements based on the deformation of giant vesicles.

"Membrane order" is a term commonly used to describe the elastic and mechanical properties of the lipid bilayer, though its exact meaning is somewhat context- and method dependent. These mechanical properties of the membrane control many cellular functions and are measured using various biophysical techniques. Here, we ask if the results obtained from various techniques are mutually consistent. Such consistency cannot be assumed a priori because these techniques probe different spatial locations and different spatial and temporal scales. We evaluate the change of membrane order induced by serotonin using nine different techniques in lipid bilayers of three different compositions. Serotonin is an important neurotransmitter present at 100s of mM concentrations in neurotransmitter vesicles, and therefore its interaction with the lipid bilayer is biologically relevant. Our measurement tools include fluorescence of lipophilic dyes (Nile Red, Laurdan, TMA-DPH, DPH), whose properties are a function of membrane order; atomic force spectroscopy, which provides a measure of the force required to indent the lipid bilayer; 2H solid-state NMR spectroscopy, which measures the molecular order of the lipid acyl chain segments; fluorescence correlation spectroscopy, which provides a measure of the diffusivity of the probe in the membrane; and Raman spectroscopy, where spectral intensity ratios are affected by acyl chain order. We find that different measures often do not correlate with each other and sometimes even yield conflicting results. We conclude that no probe provides a general measure of membrane order and that any inference based on the change of membrane order measured by a particular probe may be unreliable.

We monitored the effect on function of the G-protein-coupled receptor (GPCR) rhodopsin from small, stepwise changes in bilayer thickness induced by cholesterol. Over a range of phosphatidylcholine bilayers with hydrophobic thickness from ≈21 Å to 38 Å, the metarhodopsin-I (MI)/metarhodopsin-II (MII) equilibrium was monitored with UV-visible spectroscopy while ordering of hydrocarbon chains was probed by 2H-NMR. Addition of cholesterol shifted equilibrium toward MII for bilayers thinner than the average length of hydrophobic transmembrane helices (27 Å) and to MI for thicker bilayers, while small bilayer thickness changes within the range of the protein hydrophobic thickness drastically up- or downregulated MII formation. The cholesterol-induced shifts toward MII for thinner membranes correlated with the cholesterol-induced increase of bilayer hydrophobic thickness measured by NMR, consistent with continuum elastic modeling. The energetic penalty of adding cholesterol to thick bilayers caused rhodopsin oligomerization and a shift toward MI. In membranes of physiological thickness, changes in bilayer mechanical properties induced by cholesterol potentiated the interplay between bilayer and protein thickness resulting in large swings of the MI-MII equilibrium. In membrane containing cholesterol, elastic deformations near the protein are a dominant energetic contribution to the functional equilibrium of the model GPCR rhodopsin.

Cyclic lipopeptides (CLiPs) have many biological functions, including the selective permeabilization of target membranes, and technical and medical applications. We studied the anionic CLiP viscosin from Pseudomonas along with a neutral analog, pseudodesmin A, and the cationic viscosin-E2K to better understand electrostatic effects on target selectivity. Calcein leakage from liposomes of anionic phosphatidylglycerol (PG) and phosphatidylethanolamine (PE) is measured in comparison with net-neutral phosphatidylcholine by time-resolved fluorescence. By contrast to the typical selectivity of cationic peptides against anionic membranes, we find viscosin more active against PG/PE at 30 μM lipid than viscosin-E2K. At very low lipid concentration, the selectivity is reversed. An equi-activity analysis reveals the reciprocal partition coefficients, 1/K, and the CLiP-to-lipid mole ratio within the membrane as leakage after 1 h reaches 50%, Re50. As expected, 1/K to PG/PE is much lower (higher affinity) for viscosin-E2K (3 μM) than viscosin (15 μM). However, the local damage to the PG/PE membrane caused by a viscosin molecule is much stronger than that of viscosin-E2K. This can be explained by the strong membrane expansion due to PG/viscosin repulsion inducing asymmetry stress between the two leaflets and, ultimately, transient limited leakage at Re50 = 0.08. PG/viscosin-E2K attraction opposes expansion and leakage starts only as the PG charges in the outer leaflet are essentially compensated by the cationic peptide (Re50 = 0.32). In the high-lipid regime (at lipid concentrations cL ≫ 1/K), virtually all CLiP is membrane bound anyway and Re50 governs selectivity, favoring viscosin. In the low-lipid regime at cL ≪ 1/K, virtually all CLiP is in solution, 1/K becomes important and the "cation attacks anionic membrane" selectivity gets restored. Overall, activity and selectivity data can only properly be interpreted if the lipid regime is known and predictions for other lipid concentrations or cell counts require knowledge of 1/K and Re50.

High-resolution x-ray data are reported for the ordered phases of long-chain di-monounsaturated C22:1 phosphocholine lipid bilayers. Similar to PC lipids that have saturated chains, diC22:1PC has a subgel phase and a gel phase, but dissimilarly, we find no ripple phase. Our quantitative focus is on the structure of the gel phase. We have recorded 17 lamellar orders, indicating a very well-ordered structure. Fitting to a model provides the phases of the orders. The Fourier construction of the electron density profile has two well-defined headgroup peaks and a very sharp and deep methyl trough. The wide-angle scattering exhibits two Bragg rods that provide the area per molecule. They have an intensity pattern quite different than that of lipids with saturated chains. Models of chain packing indicate that ground state chain configurations are tilted primarily toward next nearest neighbors with an angle that is also consistent with the modeling of the electron density profile. Wide-angle modeling also indicates broken mirror symmetry between the monolayers. Our wide-angle results and our electron density profile together leads to the hypothesis that the sn-1 and sn-2 chains have equivalent penetration depths in contrast to the gel phase structure of lipids with saturated hydrocarbon chains.

Krokinobacter eikastus rhodopsin 2 (KR2) is a light-driven pentameric sodium pump. Its ability to translocate cations other than protons and to create an electrochemical potential makes it an attractive optogenetic tool. Tailoring its ion-pumping characteristics by mutations is therefore of great interest. In addition, understanding the functional and structural consequences of certain mutations helps to derive a functional mechanism of ion selectivity and transfer of KR2. Based on solid-state NMR spectroscopy, we report an extensive chemical shift resonance assignment of KR2 within lipid bilayers. This data set was then used to probe site-resolved allosteric effects of sodium binding, which revealed multiple responsive sites including the Schiff base nitrogen and the NDQ motif. Based on this data set, the consequences of the H180A mutation are probed. The mutant is silenced in the presence of sodium while in its absence proton pumping is observed. Our data reveal specific long-range effects along the sodium transfer pathway. These experiments are complemented by time-resolved optical spectroscopy. Our data suggest a model in which sodium uptake by the mutant can still take place, while sodium release and backflow control are disturbed.

Lipid bilayers form the main matrix of functional cell membranes, and their dynamics underlie a host of physical and biological processes. Here we show that elastic membrane properties and collective molecular dynamics (MD) are related by the mean-square amplitudes (order parameters) and relaxation rates (correlation times) of lipid acyl chain motions. We performed all-atom MD simulations of liquid-crystalline bilayers that allow direct comparison with carbon-hydrogen (CH) bond relaxations measured with NMR spectroscopy. Previous computational and theoretical approaches have assumed isotropic relaxation, which yields inaccurate description of lipid chain dynamics and incorrect data interpretation. Instead, the new framework includes a fixed bilayer normal (director axis) and restricted anisotropic motion of the CH bonds in accord with their segmental order parameters, enabling robust validation of lipid force fields. Simulated spectral densities of thermally excited CH bond fluctuations exhibited well-defined spin-lattice (Zeeman) relaxations analogous to those in NMR measurements. Their frequency signature could be fit to a simple power-law function, indicative of nematic-like collective dynamics. Moreover, calculated relaxation rates scaled as the squared order parameters yielding an apparent κC modulus for bilayer bending. Our results show a strong correlation with κC values obtained from solid-state NMR studies of bilayers without and with cholesterol as validated by neutron spin-echo measurements of membrane elasticity. The simulations uncover a critical role of interleaflet coupling in membrane mechanics and thus provide important insights into molecular sites of emerging elastic properties within lipid bilayers.

Antimicrobial peptides are an important class of membrane-active peptides that can provide alternatives or complements to classic antibiotics. Among the many classes of AMPs, the histidine-rich family is of particular interest since they may induce pH-sensitive interactions with cell membranes. The AMP caerin 1.1 (Cae-1), from Australian tree frogs, has three histidine residues, and thus we studied the pH dependence of its interactions with model cell membranes. Using NMR spectroscopy and molecular dynamics simulations, we showed that Cae-1 induced greater perturbation of the lipid dynamics and water penetrations within the membrane interior in an acidic environment compared with physiological conditions. Using 31P solid-state NMR, the packing, chemical environment, and dynamics of the lipid headgroup were monitored. 2H solid-state NMR showed that Cae-1 ordered the acyl chains of the hydrophobic core of the bilayer. These results supported the molecular dynamics data, which showed that Cae-1 was mainly inserted within the lipid bilayer for both neutral and negatively charged membranes, with the charged residues pulling the water and phosphate groups inward. This could be an early step in the mechanism of membrane disruption by histidine-rich antimicrobial peptides and indicated that Cae-1 acts via a transmembrane mechanism in bilayers of neutral and anionic phospholipid membranes, especially in acidic conditions.

The triple glutamine (Q) mutant (QQQ) structure of a Cl-/H+ antiporter from Escherichia coli (ClC-ec1) displaying a novel backbone arrangement has been used to challenge the long-held notion that Cl-/H+ antiporters do not operate through large conformational motions. The QQQ mutant substitutes the glutamine residue for an external glutamate E148, an internal glutamate E203, and a third glutamate E113 that hydrogen-bonds with E203. However, it is unknown if QQQ represents a physiologically relevant state, as well as how the protonation of the wild-type glutamates relates to the global dynamics. We herein apply continuous constant-pH molecular dynamics to investigate the H+-coupled dynamics of ClC-ec1. Although any large-scale conformational rearrangement upon acidification would be due to the accumulation of excess charge within the protein, protonation of the glutamates significantly impacts mainly the local structure and dynamics. Despite the fact that the extracellular pore enlarges at acidic pHs, an occluded ClC-ec1 within the active pH range of 3.5-7.5 requires a protonated E148 to facilitate extracellular Cl- release. E203 is also involved in the intracellular H+ transfer as an H+ acceptor. The water wire connection of E148 with the intracellular solution is regulated by the charge states of the E113/E203 dyad with coupled proton titration. However, the dynamics extracted from our simulations are not QQQ-like, indicating that the QQQ mutant does not represent the behavior of the wild-type ClC-ec1. These findings reinforce the necessity of having a protonatable residue at the E203 position in ClC-ec1 and suggest that a higher level of complexity exists for the intracellular H+ transfer in Cl-/H+ antiporters.

The fusion peptide (FP) domain is necessary for the fusogenic activity of spike proteins in a variety of enveloped viruses, allowing the virus to infect the host cell, and is the only part of the protein that interacts directly with the target membrane lipid tails during fusion. There are consistent findings of poration by this domain in experimental model membrane systems, and, in certain conditions, the isolated FPs can generate pores. Here, we use molecular dynamics simulations to investigate the specifics of how these FP-induced pores form in membranes with different compositions of lysolipid and POPC. The simulations show that pores form spontaneously at high lysolipid concentrations via hybrid intermediates, where FP aggregates in the cis leaflet tilt to form a funnel-like structure that spans the leaflet and locally reduces the hydrophobic thickness that must be traversed by water to form a pore. By restraining a single FP within an FP aggregate to this tilted conformation, pores can be formed in lower-lysolipid-content membranes, including pure POPC, on the 100-ns timescale, much more rapidly than in unbiased simulations in bilayers with the same composition. The pore formation pathway is similar to the spontaneous formation in high lysolipid concentrations. Depending on the membrane composition, the pores can be metastable (as seen in POPC) or lead to membrane rupture.

Upon nutrient limitation, budding yeast of Saccharomyces cerevisiae shift from fast growth (the log stage) to quiescence (the stationary stage). This shift is accompanied by liquid-liquid phase separation in the membrane of the vacuole, an endosomal organelle. Recent work indicates that the resulting micrometer-scale domains in vacuole membranes enable yeast to survive periods of stress. An outstanding question is which molecular changes might cause this membrane phase separation. Here, we conduct lipidomics of vacuole membranes in both the log and stationary stages. Isolation of pure vacuole membranes is challenging in the stationary stage, when lipid droplets are in close contact with vacuoles. Immuno-isolation has previously been shown to successfully purify log-stage vacuole membranes with high organelle specificity, but it was not previously possible to immuno-isolate stationary-stage vacuole membranes. Here, we develop Mam3 as a bait protein for vacuole immuno-isolation, and demonstrate low contamination by non-vacuolar membranes. We find that stationary-stage vacuole membranes contain surprisingly high fractions of phosphatidylcholine lipids (∼40%), roughly twice as much as log-stage membranes. Moreover, in the stationary stage, these lipids have higher melting temperatures, due to longer and more saturated acyl chains. Another surprise is that no significant change in sterol content is observed. These lipidomic changes, which are largely reflected on the whole-cell level, fit within the predominant view that phase separation in membranes requires at least three types of molecules to be present: lipids with high melting temperatures, lipids with low melting temperatures, and sterols.