OrganL: Dynamic triangulation of biomembranes using curved elements.

We describe a method for simulating biomembranes of arbitrary shape. In contrast to other dynamically triangulated surface (DTS) algorithms, our method provides a rich, quasi-tangent-continuous, yet local description of the surface. We use curved Nagata triangles, which we generalize to cubic order to achieve the requisite flexibility. The resulting interpolation can be constructed locally without iterations, at the cost of having only approximate tangent continuity away from the vertices. This allows us to provide a parallelized and fine-tuned Monte Carlo implementation. As a first example of the potential benefits of the enhanced description, our method supports inhomogeneous lipid distributions as well as lipid mixing. It also supports restraints and constraints of various types and is constructed to be as easily extensible as possible. We validate the approach by testing its numerical accuracy, followed by reproducing the known Helfrich solutions for shapes with rotational symmetry. Finally, we present some example applications, including curvature-driven demixing and stylized effects of proteins. Input files for these examples, as well as the implementation itself, are freely available for researchers under the name OrganL (https://zenodo.org/doi/10.5281/zenodo.11204709).

Cell Membrane State, Permeability, and Elasticity Assessment for Single Cells and Cell Ensembles.

The phase state and especially phase transitions of synthetic lipid membranes are known to drastically modulate mechanical membrane properties like permeability and bending modulus. Although the main transition of lipid membranes is typically detected employing differential scanning calorimetry (DSC), this technique is not suitable for many biological membranes. Moreover, often single cell data on the membrane state or order is of interest. We here first describe how to use a membrane polarity-sensitive dye, Laurdan, to optically determine the order of cell ensembles over a wide temperature range from T = -40 °C to +95 °C. This allows to quantify the position and width of biological membrane order-disorder transitions. Second, we show that the distribution of membrane order within a cell ensemble allows for correlation analysis of membrane order and permeability. Third, combining the technique with conventional atomic force spectroscopy allows for the quantitative correlation of an overall effective Young's modulus of living cells with the membrane order.

Dissecting the mechanisms of environment sensitivity of smart probes for quantitative assessment of membrane properties.

The plasma membrane, as a highly complex cell organelle, serves as a crucial platform for a multitude of cellular processes. Its collective biophysical properties are largely determined by the structural diversity of the different lipid species it accommodates. Therefore, a detailed investigation of biophysical properties of the plasma membrane is of utmost importance for a comprehensive understanding of biological processes occurring therein. During the past two decades, several environment-sensitive probes have been developed and become popular tools to investigate membrane properties. Although these probes are assumed to report on membrane order in similar ways, their individual mechanisms remain to be elucidated. In this study, using model membrane systems, we characterized the probes Pro12A, NR12S and NR12A in depth and examined their sensitivity to parameters with potential biological implications, such as the degree of lipid saturation, double bond position and configuration (cis versus trans), phospholipid headgroup and cholesterol content. Applying spectral imaging together with atomistic molecular dynamics simulations and time-dependent fluorescent shift analyses, we unravelled individual sensitivities of these probes to different biophysical properties, their distinct localizations and specific relaxation processes in membranes. Overall, Pro12A, NR12S and NR12A serve together as a toolbox with a wide range of applications allowing to select the most appropriate probe for each specific research question.

Membrane Remodeling Due to a Mixture of Multiple Types of Curvature Proteins.

We present an extension of the Monte Carlo based mesoscopic membrane model, where the membrane is represented as a dynamically triangulated surface and the proteins are modeled as anisotropic inclusions formulated as in-plane nematic field variables adhering to the deformable elastic sheet. In the extended model, we have augmented the Hamiltonian to study membrane deformation due to a mixture of multiple types of curvature generating proteins. This feature opens the door for understanding how multiple kinds of curvature-generating proteins may be working in a coordinated manner to induce desired membrane morphologies. For example, among other things, we study membrane deformations and tubulation due to a mixture of positive and negative curvature proteins as mimics of various proteins from BAR domain family. We also study the effect of membrane anisotropy that manifests as differential binding affinity and organization of curvature proteins, leading to insights into the tightly regulated cargo sorting and transport processes. Our simulation results show different morphology of deformed vesicles that depend on membrane tension, the curvatures and number of the participating proteins as well as on protein-protein and membrane-protein interactions.

Spatial Patterns of Light-Harvesting Antenna Complex Arrangements Tune the Transfer-to-Trap Efficiency of Excitons in Purple Bacteria.

In photosynthesis, the efficiency with which a photogenerated exciton reaches the reaction center is dictated by chromophore energies and the arrangement of chromophores in the supercomplex. Here, we explore the interplay between the arrangement of light-harvesting antennae and the efficiency of exciton transport in purple bacterial photosynthesis. Using a Miller-Abrahams-based exciton hopping model, we compare different arrangements of light-harvesting proteins on the intracytoplasmic membrane. We find that arrangements with aggregated LH1s have a higher efficiency than arrangements with randomly distributed LH1s in a wide range of physiological light fluences. This effect is robust to the introduction of defects on the intracytoplasmic membrane. Our result explains the absence of species with aggregated LH1 arrangements in low-light niches and the large increase seen in the expression of LH1 dimer complexes in high fluences. We suggest that the effect seen in our study is an adaptive strategy toward solar light fluence across different purple bacterial species.

A thermodynamic scaling law for electrically perturbed lipid membranes: Validation with steepest entropy ascent framework.

Experimental evidence has demonstrated the ability of transient pulses of electric fields to alter mammalian cell behavior. Strategies with these pulsed electric fields (PEFs) have been developed for clinical applications in cancer therapeutics, in-vivo decellularization, and tissue regeneration. Successful implementation of these strategies involve understanding how PEFs impact the cellular structures and, hence, cell behavior. The caveat, however, is that the PEF parameter space (i.e., comprising different pulse widths, amplitudes, number of pulses) is large, and design of experiments to explore all possible combinations of pulse parameters is prohibitive from a cost and time standpoint. In this study, a scaling law based on the Ising model is introduced to understand the impact of PEFs on the outer cell lipid membrane so that an understanding developed in one PEF pulse regime may be extended to another. Combining non-Markovian Monte Carlo techniques to determine density-of-states with a novel non-equilibrium thermodynamic framework based on the principle of steepest entropy ascent, the applicability of this scaling model to predict the behavior of both thermally quenched and electrically perturbed lipid membranes is demonstrated. A comparison of the predictions made by the steepest-entropy-ascent quantum thermodynamic (SEAQT) framework to experimental data is performed to validate the robustness of this computational methodology and the resulting scaling law.

Implications of an Improved Model of the TSH Receptor Transmembrane Domain (TSHR-TMD-TRIO).

The thyroid-stimulating hormone receptor (TSHR) is a G-protein-coupled receptor group A family member with 7 transmembrane helices. We generated 3 new models of its entire transmembrane region using a 600 ns molecular simulation. The simulation started from our previously published model, which we have now revised by also modeling the intracellular loops and the C-terminal tail, adding internal waters and embedding it into a lipid bilayer with a water layer and with ions added to complete the system. We have named this model TSHR-TMD-TRIO since 3 representative dominant structures were then extracted from the simulation trajectory and compared with the original model. These structures each showed small but significant changes in the relative positions of the helices. The 3 models were also used as targets to dock a set of small molecules that are known active compounds including a new TSHR antagonist (BT362), which confirmed the appropriateness of the model with some small molecules showing significant preference for one or other of the structures.

Computer simulations of a heterogeneous membrane with enhanced sampling techniques.

Computational determination of the equilibrium state of heterogeneous phospholipid membranes is a significant challenge. We wish to explore the rich phase diagram of these multi-component systems. However, the diffusion and mixing times in membranes are long compared to typical time scales of computer simulations. Here, we evaluate the combination of the enhanced sampling techniques molecular dynamics with alchemical steps and Monte Carlo with molecular dynamics with a coarse-grained model of membranes (Martini) to reduce the number of steps and force evaluations that are needed to reach equilibrium. We illustrate a significant gain compared to straightforward molecular dynamics of the Martini model by factors between 3 and 10. The combination is a useful tool to enhance the study of phase separation and the formation of domains in biological membranes.

Micropipette Aspiration-Based Assessment of Single Channel Water Permeability.

Measurements of the unitary hydraulic conductivity of membrane channels, pf , may be hampered by difficulties in producing sufficient quantities of purified and reconstituted proteins. Low yield expression, the purely empiric choice of detergents, as well as protein aggregation and misfolding during reconstitution may result in an average of less than one reconstituted channel per large unilamellar vesicle. This limits their applicability for pf measurements, independent of whether light scattering or fluorescence quenching of encapsulated dyes is monitored. Here the micropipette aspiration technique is adopted because its superb sensitivity allows resolving pf values for one order of magnitude smaller protein densities in sphingomyelin and cholesterol rich giant unilamellar vesicles (GUVs). Protein density is derived from intensity fluctuations that fluorescently labeled channels in the aspirated GUV induce by diffusing through the diffraction limited spot. A perfusion system minimizes unstirred layers in the immediate membrane vicinity as demonstrated by the distribution of both encapsulated and extravesicular aqueous dyes. pf amounted to 2.4 ± 0.1 × 10-13 cm³ s-1 for aquaporin-1 that served as a test case. The new assay paves the way for directly monitoring the effect that interaction of aquaporins with other proteins or inhibitors may have on pf on a single sample.

Effect of helical kink in antimicrobial peptides on membrane pore formation.

Every cell is protected by a semipermeable membrane. Peptides with the right properties, for example Antimicrobial peptides (AMPs), can disrupt this protective barrier by formation of leaky pores. Unfortunately, matching peptide properties with their ability to selectively form pores in bacterial membranes remains elusive. In particular, the proline/glycine kink in helical peptides was reported to both increase and decrease antimicrobial activity. We used computer simulations and fluorescence experiments to show that a kink in helices affects the formation of membrane pores by stabilizing toroidal pores but disrupting barrel-stave pores. The position of the proline/glycine kink in the sequence further controls the specific structure of toroidal pore. Moreover, we demonstrate that two helical peptides can form a kink-like connection with similar behavior as one long helical peptide with a kink. The provided molecular-level insight can be utilized for design and modification of pore-forming antibacterial peptides or toxins.

Need for speed: Super-resolving the dynamic nanoclustering of syntaxin-1 at exocytic fusion sites.

Communication between cells relies on regulated exocytosis, a multi-step process that involves the docking, priming and fusion of vesicles with the plasma membrane, culminating in the release of neurotransmitters and hormones. Key proteins and lipids involved in exocytosis are subjected to Brownian movement and constantly switch between distinct motion states which are governed by short-lived molecular interactions. Critical biochemical reactions between exocytic proteins that occur in the confinement of nanodomains underpin the precise sequence of priming steps which leads to the fusion of vesicles. The advent of super-resolution microscopy techniques has provided the means to visualize individual molecules on the plasma membrane with high spatiotemporal resolution in live cells. These techniques are revealing a highly dynamic nature of the nanoscale organization of the exocytic machinery. In this review, we focus on soluble N-ethylmaleimide-sensitive factor attachment receptor (SNARE) syntaxin-1, which mediates vesicular fusion. Syntaxin-1 is highly mobile at the plasma membrane, and its inherent speed allows fast assembly and disassembly of syntaxin-1 nanoclusters which are associated with exocytosis. We reflect on recent studies which have revealed the mechanisms regulating syntaxin-1 nanoclustering on the plasma membrane and draw inferences on the effect of synaptic activity, phosphoinositides, N-ethylmaleimide-sensitive factor (NSF), α-soluble NSF attachment protein (α-SNAP) and SNARE complex assembly on the dynamic nanoscale organization of syntaxin-1. This article is part of the special issue entitled 'Mobility and trafficking of neuronal membrane proteins'.

Nanoscale dynamics of cholesterol in the cell membrane.

Cholesterol constitutes ∼30-40% of the mammalian plasma membrane, a larger fraction than of any other single component. It is a major player in numerous signaling processes as well as in shaping molecular membrane architecture. However, our knowledge of the dynamics of cholesterol in the plasma membrane is limited, restricting our understanding of the mechanisms regulating its involvement in cell signaling. Here, we applied advanced fluorescence imaging and spectroscopy approaches on in vitro (model membranes) and in vivo (live cells and embryos) membranes as well as in silico analysis to systematically study the nanoscale dynamics of cholesterol in biological membranes. Our results indicate that cholesterol diffuses faster than phospholipids in live membranes, but not in model membranes. Interestingly, a detailed statistical diffusion analysis suggested two-component diffusion for cholesterol in the plasma membrane of live cells. One of these components was similar to a freely diffusing phospholipid analogue, whereas the other one was significantly faster. When a cholesterol analogue was localized to the outer leaflet only, the fast diffusion of cholesterol disappeared, and it diffused similarly to phospholipids. Overall, our results suggest that cholesterol diffusion in the cell membrane is heterogeneous and that this diffusional heterogeneity is due to cholesterol's nanoscale interactions and localization in the membrane.

Tubulation pattern of membrane vesicles coated with biofilaments.

Narrow membrane tubes are commonly pulled out from the surface of phospholipid vesicles using forces applied either through laser or magnetic tweezers or through the action of processive motor proteins. Recent examples have emerged in which an array of such tubes grows spontaneously from vesicles coated with bioactive cytoskeletal filaments (e.g., FtsZ, microtubule) in the presence GTP or ATP. We show how a soft vesicle deforms as a result of the interplay between its topology, local curvature, and the forces due to filament bundles. We present results from dynamically triangulated Monte Carlo simulations of a closed membrane vesicle coated with a nematic field (the filaments), and we show how the intrinsic curvature of the filaments and their bundling interactions drive membrane tubulation. We predict interesting patterns consisting of a large number of nematic defects that accompany tubulation. A common theme emerges: defect locations on vesicle surfaces are hot spots of membrane deformation activity, which could be useful for vesicle origami. Although our equilibrium model is not applicable to the nonequilibrium shape dynamics exhibited by active microtubule-coated vesicles, we show that some of the features, such as the size-dependent vesicle shape and the number of tubes, can still be understood from our equilibrium model.

Controlling Anomalous Diffusion in Lipid Membranes.

Diffusion in cell membranes is not just simple two-dimensional Brownian motion but typically depends on the timescale of the observation. The physical origins of this anomalous subdiffusion are unresolved, and model systems capable of quantitative and reproducible control of membrane diffusion have been recognized as a key experimental bottleneck. Here, we control anomalous diffusion using supported lipid bilayers containing lipids derivatized with polyethylene glycol (PEG) headgroups. Bilayers with specific excluded area fractions are formed by control of PEG lipid mole fraction. These bilayers exhibit a switch in diffusive behavior, becoming anomalous as bilayer continuity is disrupted. Using a combination of single-molecule fluorescence and interferometric imaging, we measure the anomalous behavior in this model over four orders of magnitude in time. Diffusion in these bilayers is well described by a power-law dependence of the mean-square displacement with observation time. Anomaleity in this system can be tailored by simply controlling the mole fraction of PEG lipid, producing bilayers with diffusion parameters similar to those observed for anomalous diffusion in biological membranes.

Secretory Expression of a Chimeric Peptide in Lactococcus lactis: Assessment of its Cytotoxic Activity and a Deep View on Its Interaction with Cell-Surface Glycosaminoglycans by Molecular Modeling.

Nowadays, cancer remains a major cause of death affecting millions of people. Currently, the antimicrobial peptides (AMPs) as potent anticancer therapeutic agents offer specificity and low levels of side effects in cancer therapy. In the present study, a cationic chimeric peptide (cLFchimera), derived from camel lactoferrin, was expressed as a secretory peptide using P170 expression system in L. lactis. Peptide purification was carried out using Ni-NTA agarose column from culture medium with 21 µ/mL concentration. The recombinant peptide was investigated for its activity against four tumor and one normal cell line. The cLFchimera was more active against two tumor cell lines (chondrosarcoma and colorectal cancer cells), but the activity against two other tumor cell lines (hepatoma and breast cancer cell line) and normal cells was low. Finally, to have better insight into the mode of action of the peptide on cytotoxic activity, we examined the interaction of cationic peptide with two glycosaminoglycans (GAGs), heparan sulfate (HS) and chondroitin sulfate (CS), as the two most anionic molecules on the cell surface by molecular dynamic simulation. The results of in silico analysis showed that the cLFchimera interacted with HS and CS with a totally different amino acid profile. Hydrogen bonding screening in GAGs-peptide complexes revealed K21, V23 and I3, R16 are the dominant amino acids involved in peptide-HS and CS interaction, respectively. Overall, the results of this investigation showed the P170 expression system successfully expressed a cationic peptide with potent anticancer activity. Moreover, molecular docking analysis revealed the pattern of peptide interaction with negatively charged membrane molecules.

Fluorimetric Techniques for the Assessment of Sperm Membranes.

Standard spermiograms describing sperm quality are mostly based on the physiological and visual parameters, such as ejaculate volume and concentration, motility and progressive motility, and sperm morphology and viability. However, none of these assessments is good enough to predict the semen quality. Given that maintenance of sperm viability and fertilization potential depends on membrane integrity and intracellular functionality, evaluation of these parameters might enable a better prediction of sperm fertilization competence. Here, we describe three feasible methods to evaluate sperm quality using specific fluorescent probes combined with fluorescence microscopy or flow cytometry analyses. Analyses assessed plasma membrane integrity using 4',6-diamidino-2-phenylindole (DAPI) and propidium iodide (PI), acrosomal membrane integrity using fluorescein isothiocyanate-conjugated Pisum sativum agglutinin (FITC-PSA) and mitochondrial membrane integrity using 5,5',6,6'-tetra-chloro-1,1',3,3'-tetraethylbenzimidazolyl carbocyanine iodide (JC-1). Combinations of these methods are also presented. For instance, use of annexin V combined with PI fluorochromes enables assessing apoptosis and calculating the proportion of apoptotic sperm (apoptotic index). We believe that these methodologies, which are based on examining spermatozoon membranes, are very useful for the evaluation of sperm quality.

LDL Cholesterol Uptake Assay Using Live Cell Imaging Analysis with Cell Health Monitoring.

The regulation of LDL cholesterol uptake through LDLR-mediated endocytosis is an important area of study in various major pathologies including metabolic disorder, cardiovascular disease, and kidney disease. Currently, there is no available method to assess LDL uptake while simultaneously monitoring for health of the cells. The current study presents a protocol, using a live cell imaging analysis system, to acquire serial measurements of LDL influx with concurrent monitoring for cell health. This novel technique is tested in three human cell lines (hepatic, renal tubular epithelial, and coronary artery endothelial cells) over a four-hour time course. Moreover, the sensitivity of this technique is validated with well-known LDL uptake inhibitors, Dynasore and recombinant PCSK9 protein, as well as by an LDL uptake promoter, Simvastatin. Taken together, this method provides a medium-to-high throughput platform for simultaneously screening pharmacological activity as well as monitoring of cell morphology, hence cytotoxicity of compounds regulating LDL influx. The analysis can be used with different imaging systems and analytical software.

Effect of the Organization of Rhodopsin on the Association between Transducin and a Photoactivated Receptor.

After photoactivation, rhodopsin (R), a G-protein-coupled receptor, rapidly activates multiple transducin G-proteins (G) in an initial amplification step of phototransduction. G-protein activation requires diffusion-mediated association with an active rhodopsin (R*) at the rod disk membrane. Different organizations of R within the membrane have been revealded by several microscopy studies, including static and freely diffusing situations. However, it is unclear how such different scenarios influence the activation rate of G proteins. Through Monte Carlo simulations, we study the association reaction between a photoactivated rhodopsin and transducin under different reported receptor organizations including (a) R monomers diffusing freely, (b) R forming static dispersed crystalline domains made of rows of dimers, and (c) R arranged in static tracks formed by two adjacent rows of dimers. A key parameter in our simulations is the probability of binding following a collision ( p). For high p, the association rate between R* and G is higher in the freely diffusive system than in the static organizations, but for low collision efficiencies, the static organizations can result in faster association rates than the mobile system. We also observe that with low p, increasing the concentration of R increases the association rate significantly in the dispersed crystals configuration and just slightly in the free diffusive system. In summary, the lateral organization of rhodopsin influences the association rate between R* and G in a manner strongly dependent on the collision efficiency.

Activity of Antimicrobial Peptide Aggregates Decreases with Increased Cell Membrane Embedding Free Energy Cost.

Antimicrobial peptides (AMPs) are a promising alternative to antibiotics for mitigating bacterial infections, in light of increasing bacterial resistance to antibiotics. However, predicting, understanding, and controlling the antibacterial activity of AMPs remain a significant challenge. While peptide intramolecular interactions are known to modulate AMP antimicrobial activity, peptide intermolecular interactions remain elusive in their impact on peptide bioactivity. Herein, we test the relationship between AMP intermolecular interactions and antibacterial efficacy by controlling AMP intermolecular hydrophobic and hydrogen bonding interactions. Molecular dynamics simulations and Gibbs free energy calculations in concert with experimental assays show that increasing intermolecular interactions via interpeptide aggregation increases the energy cost for the peptide to embed into the bacterial cell membrane, which in turn decreases the AMP antibacterial activity. Our findings provide a route for predicting and controlling the antibacterial activity of AMPs against Gram-negative bacteria via reductions of intermolecular AMP interactions.

Formation and stabilization of pores in bilayer membranes by peptide-like amphiphilic polymers.

We study pore formation in models of lipid bilayer membranes interacting with amphiphilic copolymers mimicking anti-microbial peptides using Monte Carlo simulations and we rationalize our results by a simple brush-model for the fluid membrane. In our study a weak tension on the membrane is required to observe pore-formation induced by the adsorption of flexible amphiphilic copolymers. The copolymers enhance the pore stability by decreasing the line tension due to weak adsorption along the rim of the pore. Pore formation is enhanced with increasing length of copolymers or stronger stretching of the membrane. Both solvent and copolymer permeability increase as the pore becomes stable. Pore-formation proceeds via a meta-stable pore-state according to a discontinuous phase transition scenario which lead to finite pore-sizes at once. Our generic model of copolymer-induced pore-formation does not require high polymer concentration at the pores nor any self-organization of the copolymers to open the pore.