Effects of Stretching Speed on Mechanical Rupture of Phospholipid/Cholesterol Bilayers: Molecular Dynamics Simulation

Microfluidic Electroporation Arrays for Investigating Electroporation-Induced Cellular Rupture Dynamics

Electroporation is pivotal in bioelectrochemistry for cellular manipulation, with prominent applications in drug delivery and cell membrane studies. A comprehensive understanding of pore generation requires an in-depth analysis of the critical pore size and the corresponding energy barrier at the onset of cell rupture. However, many studies have been limited to basic models such as artificial membranes or theoretical simulations. Challenging this paradigm, our study pioneers using a microfluidic electroporation chip array. This tool subjects live breast cancer cell species to a diverse spectrum of alternating current electric field conditions, driving electroporation-induced cell rupture. We conclusively determined the rupture voltages across varying applied voltage loading rates, enabling an unprecedented characterization of electric cell rupture dynamics encompassing critical pore radius and energy barrier. Further bolstering our investigation, we probed cells subjected to cholesterol depletion via methyl-β-cyclodextrin and revealed a strong correlation with electroporation. This work not only elucidates the dynamics of electric rupture in live cell membranes but also sets a robust foundation for future explorations into the mechanisms and energetics of live cell electroporation.

How to Build Plasmon-Driven Molecular Jackhammers that Disassemble Cell Membranes and Cytoskeletons in Cancer

Plasmon-driven molecular machines with ultrafast motion at the femtosecond scale are effective for the treatment of cancer and other diseases. It is recently shown that cyanine dyes act as molecular jackhammers (MJH) through vibronic (vibrational and electronic mode coupling) driven activation that causes the molecule to stretch longitudinally and axially through concerted whole molecule vibrations. However, the theoretical and experimental underpinnings of these plasmon-driven motions in molecules are difficult to assess. Here the use of near-infrared (NIR) light-activated plasmons in a broad array of MJH that mechanically disassemble membranes and cytoskeletons in human melanoma A375 cells is described. The characteristics of plasmon-driven molecular mechanical disassembly of supramolecular biological structures are observed and recorded using real-time fluorescence confocal microscopy. Molecular plasmon resonances in MJH are quantified through a new experimental plasmonicity index method. This is done through the measurement of the UV-vis-NIR spectra in various solvents, and quantification of the optical response as a function of the solvent polarity. Structure-activity relationships are used to optimize the synthesis of plasmon-driven MJH, applying them to eradicate human melanoma A375 cells at low lethal concentrations of 75 nm and 80 mW cm-2 of 730 nm NIR-light for 10 min.

Effect of Force Field Resolution on Membrane Mechanical Response and Mechanoporation Damage under Deformation Simulations

Damage induced by transient disruption and mechanoporation in an intact cell membrane is a vital nanoscale biomechanical mechanism that critically affects cell viability. To complement experimental studies of mechanical membrane damage and disruption, molecular dynamics (MD) simulations have been performed at different force field resolutions, each of which follows different parameterization strategies and thus may influence the properties and dynamics of membrane systems. Therefore, the current study performed tensile deformation MD simulations of bilayer membranes using all-atom (AA), united-atom (UA), and coarse-grained Martini (CG-M) models to investigate how the damage biomechanics differs across atomistic and coarse-grained (CG) simulations. The mechanical response and mechanoporation damage were qualitatively similar but quantitatively different in the three models, including some progressive changes based on the coarse-graining level. The membranes yielded and reached ultimate strength at similar strains; however, the coarser systems exhibited lower average yield stresses and failure strains. The average failure strain in the UA model was approximately 7% lower than the AA, and the CG-M was 20% lower than UA and 27% lower than AA. The CG systems also nucleated a higher number of pores and larger pores, which resulted in higher damage during the deformation process. Overall, the study provides insight on the impact of force field-a critical factor in modeling biomolecular systems and their interactions-in inspecting membrane mechanosensitive responses and serves as a reference for justifying the appropriate force field for future studies of more complex membranes and more diverse biomolecular assemblies.

Effects of Nitro-Oxidative Stress on Biomolecules: Part 1-Non-Reactive Molecular Dynamics Simulations

Plasma medicine, or the biomedical application of cold atmospheric plasma (CAP), is an expanding field within plasma research. CAP has demonstrated remarkable versatility in diverse biological applications, including cancer treatment, wound healing, microorganism inactivation, and skin disease therapy. However, the precise mechanisms underlying the effects of CAP remain incompletely understood. The therapeutic effects of CAP are largely attributed to the generation of reactive oxygen and nitrogen species (RONS), which play a crucial role in the biological responses induced by CAP. Specifically, RONS produced during CAP treatment have the ability to chemically modify cell membranes and membrane proteins, causing nitro-oxidative stress, thereby leading to changes in membrane permeability and disruption of cellular processes. To gain atomic-level insights into these interactions, non-reactive molecular dynamics (MD) simulations have emerged as a valuable tool. These simulations facilitate the examination of larger-scale system dynamics, including protein-protein and protein-membrane interactions. In this comprehensive review, we focus on the applications of non-reactive MD simulations in studying the effects of CAP on cellular components and interactions at the atomic level, providing a detailed overview of the potential of CAP in medicine. We also review the results of other MD studies that are not related to plasma medicine but explore the effects of nitro-oxidative stress on cellular components and are therefore important for a broader understanding of the underlying processes.

Changes in free energy barrier for water permeation by stretch-induced phase transitions in phospholipid/cholesterol bilayers

Water permeation through phospholipid/cholesterol bilayers is the key to understanding tension-induced rupture of biological cell membranes. We performed molecular dynamics simulations of stretched phospholipid/cholesterol bilayers to investigate changes in the free energy profile of water molecules across the bilayer and the lipid structure responsible for water permeation. We modeled stretching of the bilayer by applying areal strain. In stretched phospholipid/cholesterol bilayers, the hydrophobic tail of the phospholipids became disordered and the free energy barrier to water permeation decreased. Upon exceeding the critical areal strain, a phase transition to an interdigitated gel phase occurred before rupture, and the hydrophobic tail ordering as well as the free energy barrier were restored. In pure phospholipid bilayers, we did not observe such recoveries. These transient recoveries in the phospholipid/cholesterol bilayer suppressed water permeation and membrane rupture, followed by an increase in the critical areal strain at which the bilayer ruptured. This result agrees with experimental results and provides a reasonable molecular mechanism for the toughness of phospholipid/cholesterol bilayers under tension.Communicated by Ramaswamy H. Sarma.

The ATP-dependent Pathways and Human Diseases

Adenosine triphosphate (ATP) is one of the most important molecules of life, present both inside the cells and extracellularly. It is an essential building block for nucleic acids biosynthesis and crucial intracellular energy storage. However, one of the most interesting functions of ATP is the role of a signaling molecule. Numerous studies indicate the involvement of ATP-dependent pathways in maintaining the proper functioning of individual tissues and organs. Herein, the latest data indicating the ATP function in the network of intra- and extracellular signaling pathways including purinergic signaling, MAP kinase pathway, mTOR and calcium signaling are collected. The main ATP-dependent processes maintaining the proper functioning of the nervous, cardiovascular and immune systems, as well as skin and bones, are summarized. The disturbances in the ATP amount, its cellular localization, or interaction with target elements may induce pathological changes in signaling pathways leading to the development of serious diseases. The impact of an ATP imbalance on the development of dangerous health dysfunctions such as neurodegeneration diseases, cardiovascular diseases (CVDs), diabetes mellitus, obesity, cancers and immune pathogenesis are discussed here.

Extreme Membrane Tensile Loads Induce Half-Activation of the Thermosensitive TRPV1 Channel

Transient receptor potential (TRP) channels are sensors for a wide range of cellular and environmental signals, but elucidating how these channels convert a wide range of physical and chemical stimuli into channel opening is essential to understanding their normal function. Here, half-activation of thermosensitive TRPV1 channel under extreme membrane stretching from blast loads was provided by molecular dynamics simulations. The results show that such extreme membrane stretch loading will only lead to half-activation of the TRPV1 channel: that is, the upper gate is open for high-speed stretching (>15m/s), but the lower gate is still closed. The corresponding activation threshold also depends on both the tensile speed and the area strain. This means that the direct mechanical gating of TRP channels in one step is unlikely to occur.

Effect of the cholesterol on electroporation of planar lipid bilayer

Electroporation threshold depends on the membrane composition, with cholesterol being one of its key components already studied in the past, but the results were inconclusive. The aim of our study was to determine behaviour of planar lipid bilayers with varying cholesterol concentrations under electric field. This would give us a better insight into cholesterol effect on membrane properties during electroporation process, since cholesterol is one of the major components of biological membranes and plays a crucial role in membrane organisation, dynamics, and function. Planar lipid bilayers were prepared from phosphatidylcholine lipids with 0, 20, 30, 50 and 80 mol% cholesterol. Capacitance was measured using the discharge method. Results show no statistical difference of cBLM between the cholesterol concentrations. Breakdown voltage Ubr of planar lipid bilayers was measured by means of linear rising voltage with seven different slopes. Obtained results were fitted to a strength-duration curve, where parameter Ubrmin represents minimal breakdown voltage, and parameter τRC represents the inclination of the strength-duration curve. Adding cholesterol to planar lipid bilayer gradually increased its Ubrmin until 50 mol% cholesterol concentration. Afterwards at 80 mol% Ubrmin does not further increase, in fact it reduces by 20% of the Ubrmin at 50 mol% cholesterol concentration.

How the spleen reshapes and retains young and old red blood cells: A computational investigation

The spleen, the largest secondary lymphoid organ in humans, not only fulfils a broad range of immune functions, but also plays an important role in red blood cell’s (RBC) life cycle. Although much progress has been made to elucidate the critical biological processes involved in the maturation of young RBCs (reticulocytes) as well as removal of senescent RBCs in the spleen, the underlying mechanisms driving these processes are still obscure. Herein, we perform a computational study to simulate the passage of RBCs through interendothelial slits (IES) in the spleen at different stages of their lifespan and investigate the role of the spleen in facilitating the maturation of reticulocytes and in clearing the senescent RBCs. Our simulations reveal that at the beginning of the RBC life cycle, intracellular non-deformable particles in reticulocytes can be biomechanically expelled from the cell upon passage through IES, an insightful explanation of why this peculiar “pitting” process is spleen-specific. Our results also show that immature RBCs shed surface area by releasing vesicles after crossing IES and progressively acquire the biconcave shape of mature RBCs. These findings likely explain why RBCs from splenectomized patients are significantly larger than those from nonsplenectomized subjects. Finally, we show that at the end of their life span, senescent RBCs are not only retained by IES due to reduced deformability but also become susceptible to mechanical lysis under shear stress. This finding supports the recent hypothesis that transformation into a hemolyzed ghost is a prerequisite for phagocytosis of senescent RBCs. Altogether, our computational investigation illustrates critical biological processes in the spleen that cannot be observed in vivo or in vitro and offer insights into the role of the spleen in the RBC physiology.

The spleen, the largest secondary lymphoid organ in humans, not only fulfils a broad range of immune functions, but also plays an important role in red blood cell (RBC) life cycle. In this study, we perform a computational study to simulate the passage of RBCs through interendothelial slits (IES) in the spleen at different stages of their lifespan, a critical biological process that cannot be observed in humans. Our simulation results illustrate a specific role of spleen in shaping young RBCs, which points to a probable missing step in current in vitro RBC culture protocols that fail to generate a majority of typical biconcave RBCs. Our results also reveal that intra-splenic mechanical constraints likely contribute to the final clearance and elimination of aged RBCs. Altogether, we demonstrate that our computational model can provide mechanistic rationales for experimental studies, offer insights into the role of the spleen in the RBC physiology and help the optimization of in vitro RBC culture techniques.

Development of a mesoscopic framework spanning nanoscale protofibril dynamics to macro-scale fibrin clot formation

Thrombi form a micro-scale fibrin network consisting of an interlinked structure of nanoscale protofibrils, resulting in haemostasis. It is theorized that the mechanical effect of the fibrin clot is caused by the polymeric protofibrils between crosslinks, or to their dynamics on a nanoscale order. Despite a number of studies, however, it is still unknown, how the nanoscale protofibril dynamics affect the formation of the macro-scale fibrin clot and thus its mechanical properties. A mesoscopic framework would be useful to tackle this multi-scale problem, but it has not yet been established. We thus propose a minimal mesoscopic model for protofibrils based on Brownian dynamics, and performed numerical simulations of protofibril aggregation. We also performed stretch tests of polymeric protofibrils to quantify the elasticity of fibrin clots. Our model results successfully captured the conformational properties of aggregated protofibrils, e.g., strain-hardening response. Furthermore, the results suggest that the bending stiffness of individual protofibrils increases to resist extension.

Impact of Lipid Peroxidation on the Response of Cell Membranes to High-Speed Equibiaxial Stretching: A Computational Study

The difference between diseased and healthy cellular membranes in response to mechanical stresses is crucial for biology, as well as in the development of medical devices. However, the biomolecular mechanisms by which mechanical stresses interact with diseased cellular components remain largely unknown. In this work, we focus on the response of diseased cellular membranes with lipid peroxidation to high-speed tensile loadings. We find that the critical areal strain (ξ_c, when the pore forms) is highly sensitive to lipid peroxidation. For example, ξ_c of a fully oxidized bilayer is only 64 and 69% of the nonoxidized one at the stretching speed of 0.1 and 0.6 m/s, respectively. ξ_c decreases with the increase in the oxidized lipid ratio, regardless of the speeds. Also, the critical rupture tension of membranes exhibits a similar change. It is obvious that the oxidized membranes are more easily damaged than normal ones by high-speed stretching, which coincides with experimental findings. The reason is that peroxidation introduces a polar group to the tail of lipids, increases the hydrophilicity of tails, and warps the tails to the membrane-water interface, which causes loose accumulation and disorder of lipid tails. This can be deduced from the variation in the area per lipid and order parameter. In addition, the lowering stretching modulus and line tension of membranes (i.e., softening) after lipid peroxidation is also a significant factor. We reveal the difference between the peroxidized (diseased) and normal membrane in response to high-speed stretching, give the ξ_c value in the pore formation of membranes and analyze the influence of the stretching speed, peroxidation ratio, and molecular structure of phospholipids. We hope that the molecular-level information will be useful for the development of biological and medical devices in the future.

Contact probing of prestressed adhesive membranes of living cells

Atomic force microscopy (AFM) studies of living biological cells is one of main experimental tools that enable quantitative measurements of deformation of the cells and extraction of information about their structural and mechanical properties. However, proper modelling of AFM probing and related adhesive contact problems are of crucial importance for interpretation of experimental data. The Johnson-Kendall-Roberts (JKR) theory of adhesive contact has often been used as a basis for modelling of various phenomena including cell-cell interactions. However, strictly speaking the original JKR theory is valid only for contact of isotropic linearly elastic spheres, while the cell membranes are often prestressed. For the first time, effects caused by molecular adhesion for living cells are analytically studied taking into account the mechanical properties of cell membranes whose stiffness depends on the level of the tensile prestress. Another important question is how one can extract the work of adhesion between the probe and the cell. An extended version of the Borodich-Galanov method for non-direct extraction of elastic and adhesive properties of contacted materials is proposed to apply to experiments of cell probing. Evidently, the proposed models of adhesive contact for cells with prestressed membranes do not cover all types of biological cells because the structure and properties of the cells may vary considerably. However, the obtained results can be applied to many types of smooth cells and can be used to describe initial stages of contact and various other processes when effects of adhesion are of crucial importance. This article is part of a discussion meeting issue 'A cracking approach to inventing new tough materials: fracture stranger than friction'.

A new purification technique to obtain specific size distribution of giant lipid vesicles using dual filtration

A new purification technique is developed for obtaining distribution of giant unilamellar vesicles (GUVs) within a specific range of sizes using dual filtration. The GUVs were prepared using well known natural swelling method. For filtration, different combinations of polycarbonate membranes were implemented in filter holders. In our experiment, the combinations of membranes were selected with corresponding pore sizes-(i) 12 and 10 μm, (ii) 12 and 8 μm, and (iii) 10 and 8 μm. By these filtration arrangements, obtained GUVs size distribution were in the ranges of 6-26 μm, 5-38 μm and 5-30 μm, respectively. In comparison, the size distribution range was much higher for single filtration technique, for example, 6-59 μm GUVs found for a membrane with 12 μm pores. Using this technique, the water-soluble fluorescent probe, calcein, can be removed from the suspension of GUVs successfully. The size distributions were analyzed with lognormal distribution. The skewness became smaller (narrow size distribution) when a dual filtration was used instead of single filtration. The mode of the size distribution obtained in dual filtration was also smaller to that of single filtration. By continuing this process of purification for a second time, the GUVs size distribution became even narrower. After using an extra filtration with dual filtration, two different size distributions of GUVs were obtained at a time. This experimental observation suggests that different size specific distributions of GUVs can be obtained easily, even if GUVs are prepared by different other methods.

Mechanics of Bacterial Interaction and Death on Nanopatterned Surfaces

Nanopatterned surfaces are believed to kill bacteria through physical deformation, a mechanism that has immense potential against biochemical resistance. Because of its elusive nature, this mechanism is mostly understood through biophysical modeling. Problematically, accurate descriptions of the contact mechanics and various boundary conditions involved in the bacteria-nanopattern interaction remain to be seen. This may underpin conflicting predictions, found throughout the literature, regarding two important aspects of the mechanism-that is, its critical action site and relationship with geometry. Herein, a robust computational analysis of bacteria-nanopattern interaction is performed using a three-dimensional finite element modeling that incorporates relevant continuum mechanical properties, multilayered envelope structure, and adhesion interaction conditions. The model is applied to more accurately study the elusory mechanism and its enhancement via nanopattern geometry. Additionally, micrographs of bacteria adhered on a nanopatterned cicada wing are examined to further inform and verify the major modeling predictions. Together, the results indicate that nanopatterned surfaces do not kill bacteria predominantly by rupture in between protruding pillars as previously thought. Instead, nondevelopable deformation about pillar tips is more likely to create a critical site at the pillar apex, which delivers significant in-plane strains and may locally rupture and penetrate the cell. The computational analysis also demonstrates that envelope deformation is increased by adhesion to nanopatterns with smaller pillar radii and spacing. These results further progress understanding of the mechanism of nanopatterned surfaces and help guide their design for enhanced bactericidal efficiency.

Topological analysis of single-stranded DNA with an alpha-hederin nanopore

Nanopores have been emerged as a powerful tool for analyzing the structural information and interactional properties of a range of biomolecules. The spatial resolution of nanopore is determined by the diameter and effective thickness of its constriction region, but the presence of vestibule or stem structure in protein-based nanopore could negatively affect the sensitivity of the nanopore when applied for genome sequencing and topological analysis of DNA. Recently, alpha-hederin (Ah) has been reported to form a sub-nanometer scale pore structure in lipid membrane. With the simple structure and extremely small effective thickness, the Ah nanopore was shown to discriminate four different types of nucleotides. However, identification of a certain nucleotide in a strand of DNA, which is essential for genome sequencing, remains challenging. Here, we investigated the resolving capability of Ah nanopore to discriminate few nucleotides in a strand of single-stranded DNA, and the factors determining the sensitivity of Ah nanopore. The Ah nanopore was shown to be able to identify as few as three adenosine nucleotides in a strand of poly cytidine, in which the dwell time of the additional current blockade that represents the adenosine residue was in good agreement with their physical length. We also found that the lateral tension and chain pressure generated around the nanopore were influenced by pore's diameter and played as a dependent variables to determine the geometry of nanopore's constriction as well as the spatial resolution of the Ah nanopore.

Coarse-Grained Modeling of Pore Dynamics on the Red Blood Cell Membrane under Large Deformations

Transient pore formation on the membrane of red blood cells (RBCs) under high mechanical tensions is of great importance in many biomedical applications, such as RBC damage (hemolysis) and mechanoporation-based drug delivery. The dynamic process of pore formation, growth, and resealing is hard to visualize in experiments. We developed a mesoscale coarse-grained model to study the characteristics of transient pores on a patch of the lipid bilayer that is strengthened by an elastic meshwork representing the cytoskeleton. Unsteady molecular dynamics was used to study the pore formation and reseal at high strain rates close to the physiological ranges. The critical strain for pore formation, pore characteristics, and cytoskeleton effects were studied. Results show that the presence of the cytoskeleton increases the critical strain of pore formation and confines the pore growth. Moreover, the pore recovery process under negative strain rates (compression) is analyzed. Simulations show that pores can remain open for a long time during the high-speed tank-treading induced stretching and compression process that a patch of the RBC membrane usually experiences under high shear flow. Furthermore, complex loading conditions can affect the pore characteristics and result in denser pores. Finally, the effects of strain rate on pore formation are analyzed. Higher rate stretching of membrane patch can result in a significant increase in the critical areal strain and density of pores. Such a model reveals the dynamic molecular process of RBC damage in biomedical devices and mechanoporation that, to our knowledge, has not been reported before.

Mechanobiology predicts raft formations triggered by ligand-receptor activity across the cell membrane

Clustering of ligand-binding receptors of different types on thickened isles of the cell membrane, namely lipid rafts, is an experimentally observed phenomenon. Although its influence on cell's response is deeply investigated, the role of the coupling between mechanical processes and multiphysics involving the active receptors and the surrounding lipid membrane during ligand-binding has not yet been understood. Specifically, the focus of this work is on G-protein-coupled receptors (GPCRs), the widest group of transmembrane proteins in animals, which regulate specific cell processes through chemical signalling pathways involving a synergistic balance between the cyclic Adenosine Monophosphate (cAMP) produced by active GPCRs in the intracellular environment and its efflux, mediated by the Multidrug Resistance Proteins (MRPs) transporters. This paper develops a multiphysics approach based on the interplay among energetics, multiscale geometrical changes and mass balance of species, i.e. active GPCRs and MRPs, including diffusion and kinetics of binding and unbinding. Because the obtained energy depends upon both the kinematics and the changes of species densities, balance of mass and of linear momentum are coupled and govern the space-time evolution of the cell membrane. The mechanobiology involving remodelling and change of lipid ordering of the cell membrane allows to predict dynamics of transporters and active receptors -in full agreement with experimentally observed cAMP levels- and how the latter trigger rafts formation and cluster on such sites. Within the current scientific debate on Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2) and on the basis of the ascertained fact that lipid rafts often serve as an entry port for viruses, it is felt that approaches accounting for strong coupling among mechanobiological aspects could even turn helpful in better understanding membrane-mediated phenomena such as COVID-19 virus-cell interaction.

Cell-free protein synthesis: The transition from batch reactions to minimal cells and microfluidic devices

Thanks to the synthetic biology, the laborious and restrictive procedure for producing a target protein in living microorganisms by biotechnological approaches can now experience a robust, pliant yet efficient alternative. The new system combined with lab-on-chip microfluidic devices and nanotechnology offers a tremendous potential envisioning novel cell-free formats such as DNA brushes, hydrogels, vesicular particles, droplets, as well as solid surfaces. Acting as robust microreactors/microcompartments/minimal cells, the new platforms can be tuned to perform various tasks in a parallel and integrated manner encompassing gene expression, protein synthesis, purification, detection, and finally enabling cell-cell signaling to bring a collective cell behavior, such as directing differentiation process, characteristics of higher order entities, and beyond. In this review, we issue an update on recent cell-free protein synthesis (CFPS) formats. Furthermore, the latest advances and applications of CFPS for synthetic biology and biotechnology are highlighted. In the end, contemporary challenges and future opportunities of CFPS systems are discussed.

Localized Axolemma Deformations Suggest Mechanoporation as Axonal Injury Trigger

Traumatic brain injuries are a leading cause of morbidity and mortality worldwide. With almost 50% of traumatic brain injuries being related to axonal damage, understanding the nature of cellular level impairment is crucial. Experimental observations have so far led to the formulation of conflicting theories regarding the cellular primary injury mechanism. Disruption of the axolemma, or alternatively cytoskeletal damage has been suggested mainly as injury trigger. However, mechanoporation thresholds of generic membranes seem not to overlap with the axonal injury deformation range and microtubules appear too stiff and too weakly connected to undergo mechanical breaking. Here, we aim to shed a light on the mechanism of primary axonal injury, bridging finite element and molecular dynamics simulations. Despite the necessary level of approximation, our models can accurately describe the mechanical behavior of the unmyelinated axon and its membrane. More importantly, they give access to quantities that would be inaccessible with an experimental approach. We show that in a typical injury scenario, the axonal cortex sustains deformations large enough to entail pore formation in the adjoining lipid bilayer. The observed axonal deformation of 10–12% agree well with the thresholds proposed in the literature for axonal injury and, above all, allow us to provide quantitative evidences that do not exclude pore formation in the membrane as a result of trauma. Our findings bring to an increased knowledge of axonal injury mechanism that will have positive implications for the prevention and treatment of brain injuries.

Molecular dynamics simulations of mechanical stress on oxidized membranes

Biomembranes are under constant attack of free radicals that may lead to lipid oxidation in conditions of oxidative stress. The products generated during lipid oxidation are responsible for structural and dynamical changes which may jeopardize the membrane function. For instance, the local rearrangements of oxidized lipid molecules may induce membrane rupture. In this study, we investigated the effects of mechanical stress on oxidized phospholipid bilayers (PLBs). Model bilayers were stretched until pore formation (or poration) using non-equilibrium molecular dynamics simulations. We studied single-component homogeneous membranes composed of lipid oxidation products, as well as two-component heterogeneous membranes with coexisting native and oxidized domains. In homogeneous membranes, the oxidation products with -OH and -OOH groups reduced the areal strain required for pore formation, whereas the oxidation product with O group behaved similarly to the native membrane. In heterogeneous membranes composed of oxidized and non-oxidized domains, we tested the hypothesis according to which poration may be facilitated at the domain interface region. However, results were inconclusive due to their large statistical variance and sensitivity to simulation setup parameters. We pointed out important technical issues that need to be considered in future simulations of mechanically-induced poration of heterogeneous membranes. This research is of interest for photodynamic therapy and plasma medicine, because ruptured and intact plasma membranes are experimentally considered hallmarks of necrotic and apoptotic cell death.

Self-Assembly and Biogenesis of the Cellular Membrane are Dictated by Membrane Stretch and Composition

The cell plasma membrane is a highly dynamic organelle governing a wide range of cellular activities including ion transport, secretion, cell division, growth, and development. The fundamental process involved in the addition of new membranes to pre-existing plasma membranes, however, is unclear. Here, we report, using biophysical, morphological, biochemical, and molecular dynamic simulations, the selective incorporation of proteins and lipids from the cytosol into the cell plasma membrane dictated by membrane stretch and composition. Stretching of the cell membrane as a consequence of volume increase following incubation in a hypotonic solution and results in the incorporation of cytosolic proteins and lipids into the existing plasma membrane. Molecular dynamic simulations further confirm that increased membrane stretch results in the rapid insertion of lipids into the existing plasma membrane. Similarly, depletion of cholesterol from the cell plasma membrane selectively alters the incorporation of lipids into the membrane.

The role of membrane tension in the action of antimicrobial peptides and cell-penetrating peptides in biomembranes

For antimicrobial peptides (AMPs) with antimicrobial and bactericidal activities and cell-penetrating peptides (CPPs) with activity to permeate through plasma membrane, their interactions with lipid bilayer region in plasma membrane play important roles in these functions. However, the elementary processes and mechanisms of their functions have not been clear. The single giant unilamellar vesicle (GUV) method has revealed the details of elementary processes of interaction of some AMPs and CPPs with lipid vesicles. In this review, we summarize the mode of action of AMPs such as magainin 2 (Mag) and CPPs such as transportan 10 (TP10), revealed by the single GUV methods, and especially we focus on the role of membrane tension in actions of Mag and TP10 and the mechanisms of their actions. First, we explain the characteristics of the single GUV method briefly. Next, we summarize the recent view on the effect of tension on physical properties of lipid bilayers and describe the role of tension in actions of Mag and TP10. Some experimental results indicate that Mag-induced pore is a stretch-activated pore. The effect of packing of transbilayer asymmetric lipid on Mag-induced pore formation is described. On the other hand, entry of fluorescent dye, carboxyfluorescein (CF)-labeled TP10 (i.e., CF-TP10), into single GUVs without pore formation is affected by tension and high concentration of cholesterol. Pre-pore model for translocation of CF-TP10 across lipid bilayer is described. The experimental methods and their analysis described here are useful for investigation of functions of the other types of AMPs, CPPs, and proteins.

Response of Single Cells to Shock Waves and Numerically Optimized Waveforms for Cancer Therapy

Shock waves are used clinically for breaking kidney stones and treating musculoskeletal indications. The mechanisms by which shock waves interact with tissue are still not well understood. Here, ultra-high-speed imaging was used to visualize the deformation of individual cells embedded in a tissue-mimicking phantom when subject to shock-wave exposure from a clinical source. Three kidney epithelial cell lines were considered to represent normal healthy (human renal epithelial), cancer (CAKI-2), and virus-transformed (HK-2) cells. The experimental results showed that during the compressive phase of the shock waves, there was a small (<2%) decrease in the projected cell area, but during the tensile phase, there was a relatively large (∼10%) increase in the projected cell area. The experimental observations were captured by a numerical model with a constitutive material framework consisting of an equation of state for the volumetric response and hyper-viscoelasticity for the deviatoric response. To model the volumetric cell response, it was necessary to change from a higher bulk modulus during the compression to a lower bulk modulus during the tensile shock loading. It was discovered that cancer cells showed a smaller deformation but faster response to the shock-wave tensile phase compared to their noncancerous counterparts. Cell viability experiments, however, showed that cancer cells suffered more damage than other cell types. These data suggest that the cell response to shock waves is specific to the type of cell and waveforms that could be tailored to an application. For example, the model predicts that a shock wave with a tensile stress of 4.59 MPa would increase cell membrane permeability for cancer cells with minimal impact on normal cells.

Mechanosensitive ATP release in the lungs: New insights from real-time luminescence imaging studies

Extracellular ATP and other nucleotides are important autocrine/paracrine mediators that stimulate purinergic receptors and regulate diverse processes in the normal lungs. They are also associated with pathogenesis of a number of respiratory diseases and clinical complications including acute respiratory distress syndrome and ventilator induced lung injury. Mechanical forces are major stimuli for cellular ATP release but precise mechanisms responsible for this release are still debated. The present review intends to provide the current state of knowledge of the mechanisms of ATP release in the lung. Putative pathways of the release, including the contribution of cell membrane injury and cell lysis are discussed addressing their strength, weaknesses and missing evidence that requires future study. We also provide an overview of the recent technical advances in studying cellular ATP release in vitro and ex vivo. Special attention is given to new insights into lung ATP release obtained with the real-time luminescence ATP imaging. This includes recent data on stretch-induced mechanosensitive ATP release in a model and primary cells of lung alveoli in vitro as well as inflation-induced ATP release in airspaces and pulmonary blood vessels of lungs, ex vivo.

3D multicellular model of shock wave-cell interaction

Understanding the interaction between shock waves and tissue is critical for advancing the use of shock waves for medical applications, such as cancer therapy. This work aims to study shock wave-cell interaction in a more realistic environment, relevant to in vitro and in vivo studies, by using 3D computational models of healthy and cancerous cells. The results indicate that for a single cell embedded in an extracellular environment, the cellular geometry does not influence significantly the membrane strain but does influence the von Mises stress. On the contrary, the presence of neighbouring cells has a strong effect on the cell response, by increasing fourfold both quantities. The membrane strain response of a cell converges with more than three neighbouring cell layers, indicating that a cluster of four layers of cells is sufficient to model the membrane strain in a large domain of tissue. However, a full 3D tissue model is needed if the stress evaluation is of main interest. A tumour mimicking multicellular spheroid model is also proposed to study mutual interaction between healthy and cancer cells and shows that cancer cells can be specifically targeted in an early stage tumour-mimicking environment.

STATEMENT OF SIGNIFICANCE

This work presents 3D computational models of shock-wave/cell interaction in a biophysically realistic environment using real cell morphology in tissue-mimicking phantoms and multicellular spheroids. Results show that cell morphology does not strongly influence the membrane strain but influences the von Mises stress. While the presence of neighbouring cells significantly increases the cell response, four cell layers are enough to capture the membrane strain change in tissue. However, a full tissue model is necessary if accurate stress analysis is needed. The work also shows that cancer cells can be specifically targeted in early stage tumour mimicking environment. This work is a step towards realistic modelling of shock-wave/cell interactions in tissues and provides insight on the use of 3D models for different scenarios.

Utilizing a Structural Mechanics Approach to Assess the Primary Effects of Injury Loads Onto the Axon and Its Components

Diffuse axonal injury (DAI) occurs as a result of the transmission of rapid dynamic loads from the head to the brain and specifically to its neurons. Despite being one of the most common and most deleterious types of traumatic brain injury (TBI), the inherent cell injury mechanism is yet to be understood. Experimental observations have led to the formulation of different hypotheses, such as mechanoporation of the axolemma and microtubules (MTs) breakage. With the goal of singling out the mechanical aspect of DAI and to resolve the ambiguity behind its injury mechanism, a composite finite element (FE) model of a representative volume of an axon was developed. Once calibrated and validated against published experimental data, the axonal model was used to simulate injury scenarios. The resulting strain distributions along its components were then studied in dependence of strain rate and of typical modeling choices such as the applied MT constraints and tau proteins failure. Results show that oversimplifying the MT bundle kinematic entails a systematic attenuation (cf = 2.33) of the computed maximum MT strain. Nevertheless, altogether, results support the hypothesis of axolemma mechanoporation as a cell-injury trigger. Moreover, for the first time the interconnection between the axolemma and the MT bundle is shown to play a role in damage localization. The proposed FE axonal model is a valuable tool to understand the axonal injury mechanism from a mechanical perspective and could be used in turn for the definition of well-informed injury criteria at the tissue and organ level.

Modeling thermal inkjet and cell printing process using modified pseudopotential and thermal lattice Boltzmann methods

Pseudopotential lattice Boltzmann methods (LBMs) can simulate a phase transition in high-density ratio multiphase flow systems. If coupled with thermal LBMs through equation of state, they can be used to study instantaneous phase transition phenomena with a high-temperature gradient where only one set of formulations in an LBM system can handle liquid, vapor, phase transition, and heat transport. However, at lower temperatures an unrealistic spurious current at the interface introduces instability and limits its application in real flow system. In this study, we proposed new modifications to the LBM system to minimize a spurious current which enables us to study nucleation dynamic at room temperature. To demonstrate the capabilities of this approach, the thermal ejection process is modeled as one example of a complex flow system. In an inkjet printer, a thermal pulse instantly heats up the liquid in a microfluidic chamber and nucleates bubble vapor providing the pressure pulse necessary to eject droplets at high speed. Our modified method can present a more realistic model of the explosive vaporization process since it can also capture a high-temperature/density gradient at nucleation region. Thermal inkjet technology has been successfully applied for printing cells, but cells are susceptible to mechanical damage or death as they squeeze out of the nozzle head. To study cell deformation, a spring network model, representing cells, is connected to the LBM through the immersed boundary method. Looking into strain and stress distribution of a cell membrane at its most deformed state, it is found that a high stretching rate effectively increases the rupture tension. In other words, membrane deformation energy is released through creation of multiple smaller nanopores rather than big pores. Overall, concurrently simulating multiphase flow, phase transition, heat transfer, and cell deformation in one unified LB platform, we are able to provide a better insight into the bubble dynamic and cell mechanical damage during the printing process.

Molecular machines open cell membranes

Beyond the more common chemical delivery strategies, several physical techniques are used to open the lipid bilayers of cellular membranes. These include using electric and magnetic fields, temperature, ultrasound or light to introduce compounds into cells, to release molecular species from cells or to selectively induce programmed cell death (apoptosis) or uncontrolled cell death (necrosis). More recently, molecular motors and switches that can change their conformation in a controlled manner in response to external stimuli have been used to produce mechanical actions on tissue for biomedical applications. Here we show that molecular machines can drill through cellular bilayers using their molecular-scale actuation, specifically nanomechanical action. Upon physical adsorption of the molecular motors onto lipid bilayers and subsequent activation of the motors using ultraviolet light, holes are drilled in the cell membranes. We designed molecular motors and complementary experimental protocols that use nanomechanical action to induce the diffusion of chemical species out of synthetic vesicles, to enhance the diffusion of traceable molecular machines into and within live cells, to induce necrosis and to introduce chemical species into live cells. We also show that, by using molecular machines that bear short peptide addends, nanomechanical action can selectively target specific cell-surface recognition sites. Beyond the in vitro applications demonstrated here, we expect that molecular machines could also be used in vivo, especially as their design progresses to allow two-photon, near-infrared and radio-frequency activation.

Sequential bottom-up assembly of mechanically stabilized synthetic cells by microfluidics

Compartments for the spatially and temporally controlled assembly of biological processes are essential towards cellular life. Synthetic mimics of cellular compartments based on lipid-based protocells lack the mechanical and chemical stability to allow their manipulation into a complex and fully functional synthetic cell. Here, we present a high-throughput microfluidic method to generate stable, defined sized liposomes termed 'droplet-stabilized giant unilamellar vesicles (dsGUVs)'. The enhanced stability of dsGUVs enables the sequential loading of these compartments with biomolecules, namely purified transmembrane and cytoskeleton proteins by microfluidic pico-injection technology. This constitutes an experimental demonstration of a successful bottom-up assembly of a compartment with contents that would not self-assemble to full functionality when simply mixed together. Following assembly, the stabilizing oil phase and droplet shells are removed to release functional self-supporting protocells to an aqueous phase, enabling them to interact with physiologically relevant matrices.

Nanoparticle-Mediated Mechanical Destruction of Cell Membranes: A Coarse-Grained Molecular Dynamics Study

The effects of binding mode, shape, binding strength, and rotational speed of actively rotating nanoparticles on the integrity of cell membranes have been systematically studied using dissipative particle dynamics simulations. With theoretical analyses of lipid density, surface tension, stress distribution, and water permeation, we demonstrate that the rotation of nanoparticles can provide a strong driving force for membrane rupture. The results show that nanoparticles embedded inside a cell membrane via endocytosis are more capable of producing large membrane deformations under active rotation than nanoparticles attached on the cell membrane surface. Nanoparticles with anisotropic shapes produce larger deformation and have a higher rupture efficiency than those with symmetric shapes. Our findings provide useful design guidelines for a general strategy based on utilizing mechanical forces to rupture cell membranes and therefore destroy the integrity of cells.

Molecular dynamics simulations of heterogeneous cell membranes in response to uniaxial membrane stretches at high loading rates

The chemobiomechanical signatures of diseased cells are often distinctively different from that of healthy cells. This mainly arises from cellular structural/compositional alterations induced by disease development or therapeutic molecules. Therapeutic shock waves have the potential to mechanically destroy diseased cells and/or increase cell membrane permeability for drug delivery. However, the biomolecular mechanisms by which shock waves interact with diseased and healthy cellular components remain largely unknown. By integrating atomistic simulations with a novel multiscale numerical framework, this work provides new biomolecular mechanistic perspectives through which many mechanosensitive cellular processes could be quantitatively characterised. Here we examine the biomechanical responses of the chosen representative membrane complexes under rapid mechanical loadings pertinent to therapeutic shock wave conditions. We find that their rupture characteristics do not exhibit significant sensitivity to the applied strain rates. Furthermore, we show that the embedded rigid inclusions markedly facilitate stretch-induced membrane disruptions while mechanically stiffening the associated complexes under the applied membrane stretches. Our results suggest that the presence of rigid molecules in cellular membranes could serve as “mechanical catalysts” to promote the mechanical destructions of the associated complexes, which, in concert with other biochemical/medical considerations, should provide beneficial information for future biomechanical-mediated therapeutics.

The Optimization of a Microfluidic CTC Filtering Chip by Simulation

The detection and separation of circulating tumor cells (CTCs) are crucial in early cancer diagnosis and cancer prognosis. Filtration through a thin film is one of the size and deformability based separation methods, which can isolate rare CTCs from the peripheral blood of cancer patients regardless of their heterogeneity. In this paper, volume of fluid (VOF) multiphase flow models are employed to clarify the cells’ filtering processes. The cells may deform significantly when they enter a channel constriction, which will induce cell membrane stress and damage if the area strain is larger than the critical value. Therefore, the cellular damage criterion characterized by membrane area strain is presented in our model, i.e., the lysis limit of the lipid bilayer is taken as the critical area strain. Under this criterion, we discover that the microfilters with slit-shaped pores do less damage to cells than those with circular pores. The influence of contact angle between the microfilters and blood cells on cellular injury is also discussed. Moreover, the optimal film thickness and flux in our simulations are obtained as 0.5 μm and 0.375 mm/s, respectively. These findings will provide constructive guidance for the improvement of next generation microfilters with higher throughput and less cellular damage.

Effect of Boundary Edge in DOPC/DPPC/Cholesterol Liposomes on Acceleration of l-Histidine Preferential Adsorption

In order to investigate the interaction of hydrophilic molecules with liposomal membranes, we employed 1-(4-(trimethylamino)phenyl)-6-phenyl-1,3,5-hexatriene and 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine-N-(5-dimethylamino-1-naphthalenesulfonyl) as fluorescent probes to monitor the surface regions of the membrane, and the results for various liposomes were plotted in correlation diagrams. According to the formation of a variety of phase states, different tendencies of decreasing surface hydrophobicity were observed in the liposomes that were modified with high concentrations of cholesterol or in the liposomes that were composed of ternary components. These liposomes, with hydrophobic surfaces, also showed preferential adsorption of l-histidine (l-His), and the hydrophobicity of the liposomal membrane at the surface changed during l-His adsorption regardless of the initial liposomal properties. Furthermore, we revealed that accelerated adsorption of l-His and preferential binding was induced in ternary liposomes forming boundaries between two separate phases.