An in-silico study on the mechanical behavior of colorectal cancer cell lines in the micropipette aspiration process.

Cancer alters the structural integrity and morphology of cells. Consequently, the cell function is overshadowed. In this study, the micropipette aspiration process is computationally modeled to predict the mechanical behavior of the colorectal cancer cells. The intended cancer cells are modeled as an incompressible Neo-Hookean visco-hyperelastic material. Also, the micropipette is assumed to be rigid with no deformation. The proposed model is validated with an in-vitro study. To capture the equilibrium and time-dependent behaviors of cells, ramp, and creep tests are respectively performed using the finite element method. Through the simulations, the effects of the micropipette geometry and the aspiration pressure on the colorectal cancer cell lines are investigated. Our findings indicate that, as the inner radius of the micropipette increases, despite the increase in deformation rate and aspirated length, the time to reach the equilibrium state increases. Nevertheless, it is obvious that increasing the tip curvature radius has a small effect on the change of the aspirated length. But, due to the decrease in the stress concentration, it drastically reduces the equilibrium time and increases the deformation rate significantly. Interestingly, our results demonstrate that increasing the aspiration pressure somehow causes the cell stiffening, thereby reducing the upward trend of deformation rate, equilibrium time, and aspirated length. Our findings provide valuable insights for researchers in cell therapy and cancer treatment and can aid in developing more precise microfluidic.

Disease-related metabolites affect protein-nanoparticle interactions.

Once in biological fluids, the surface of nanoparticles (NPs) is rapidly covered with a layer of biomolecules (i.e., the "protein corona") whose composition strongly determines their biological identity, regulates interactions with biological entities including cells and the immune system, and consequently directs the biological fate and pharmacokinetics of nanoparticles. We recently introduced the concept of a "personalized protein corona" which refers to the formation of different biological identities of the exact same type of NP after being exposed to extract plasmas from individuals who have various types of diseases. As different diseases have distinct metabolomic profiles and metabolites can interact with proteins, it is legitimate to hypothesize that metabolomic profiles in plasma may have the capacity to, at least partially, drive the formation of a personalized protein corona. To test this hypothesis, we employed a multi-scale approach composed of coarse-grained (CG) and all atom (AA) molecular dynamics (MD) simulations to probe the role of glucose and cholesterol (model metabolites in diabetes and hypercholesterolemia patients) in the interaction of fibrinogen protein and polystyrene NPs. Our results revealed that glucose and cholesterol had the capacity to induce substantial changes in the binding site of fibrinogen to the surface of NPs. More specifically, the simulation results demonstrated that increasing the metabolite amount could change the profiles of fibrinogen adsorption and replacement, what is known as the Vroman effect, on the NP surface. In addition, we also found out that metabolites can substantially determine the immune triggering potency of the fibrinogen-NP complex. Our proof-of-concept outcomes further emphasize the need for the development of patient-specific NPs in a disease type-specific manner for high yielding and safe clinical applications.

The Combined Effect of Hydrophobic Mismatch and Bilayer Local Bending on the Regulation of Mechanosensitive Ion Channels.

The hydrophobic mismatch between the lipid bilayer and integral membrane proteins has well-defined effect on mechanosensitive (MS) ion channels. Also, membrane local bending is suggested to modulate MS channel activity. Although a number of studies have already shown the significance of each individual factor, the combined effect of these physical factors on MS channel activity have not been investigated. Here using finite element simulation, we study the combined effect of hydrophobic mismatch and local bending on the archetypal mechanosensitive channel MscL. First we show how the local curvature direction impacts on MS channel modulation. In the case of MscL, we show inward (cytoplasmic) bending can more effectively gate the channel compared to outward bending. Then we indicate that in response to a specific local curvature, MscL inserted in a bilayer with the same hydrophobic length is more expanded in the constriction pore region compared to when there is a protein-lipid hydrophobic mismatch. Interestingly in the presence of a negative mismatch (thicker lipids), MscL constriction pore is more expanded than in the presence of positive mismatch (thinner lipids) in response to an identical membrane curvature. These results were confirmed by a parametric energetic calculation provided for MscL gating. These findings have several biophysical consequences for understanding the function of MS channels in response to two major physical stimuli in mechanobiology, namely hydrophobic mismatch and local membrane curvature.

Influence of Global and Local Membrane Curvature on Mechanosensitive Ion Channels: A Finite Element Approach.

Mechanosensitive (MS) channels are ubiquitous molecular force sensors that respond to a number of different mechanical stimuli including tensile, compressive and shear stress. MS channels are also proposed to be molecular curvature sensors gating in response to bending in their local environment. One of the main mechanisms to functionally study these channels is the patch clamp technique. However, the patch of membrane surveyed using this methodology is far from physiological. Here we use continuum mechanics to probe the question of how curvature, in a standard patch clamp experiment, at different length scales (global and local) affects a model MS channel. Firstly, to increase the accuracy of the Laplace's equation in tension estimation in a patch membrane and to be able to more precisely describe the transient phenomena happening during patch clamping, we propose a modified Laplace's equation. Most importantly, we unambiguously show that the global curvature of a patch, which is visible under the microscope during patch clamp experiments, is of negligible energetic consequence for activation of an MS channel in a model membrane. However, the local curvature (RL < 50) and the direction of bending are able to cause considerable changes in the stress distribution through the thickness of the membrane. Not only does local bending, in the order of physiologically relevant curvatures, cause a substantial change in the pressure profile but it also significantly modifies the stress distribution in response to force application. Understanding these stress variations in regions of high local bending is essential for a complete understanding of the effects of curvature on MS channels.

The effect of local bending on gating of MscL using a representative volume element and finite element simulation.

Many physiological processes such as cell division, endocytosis and exocytosis cause severe local curvature of the cell membrane. Local curvature has been shown experimentally to modulate numerous mechanosensitive (MS) ion channels. In order to quantify the effects of local curvature we introduced a coarse grain representative volume element for the bacterial mechanosensitive ion channel of large conductance (MscL) using continuum elasticity. Our model is designed to be consistent with the channel conformation in the closed and open states to capture its major continuum rheological behavior in response to the local membrane curvature. Herein we show that change in the local curvature of the lipid bilayer can modulate MscL activity considerably by changing both bilayer thickness and lateral pressure profile. Intriguingly, although bending in any direction results in almost the same free-energy cost, inward (cytoplasmic) bending favors channel opening, whereas outward (periplasmic) bending facilitates closing of the narrowest part of the MscL pore. This quantitative study using MscL as a model channel may have wide reaching consequences for the effect of local curvature on the physiological function of other types of prokaryotic and eukaryotic membrane proteins.

Modeling and simulation of the water gradient within a Nafion membrane.

The perfluorinated sulfonic acid membrane (Nafion) shows among ionomers high water uptake and cationic conductivity. These properties allow Nafion to be used in nanocomposite actuators, sensors and fuel cells. In situ experiments have shown that there is a water gradient within the Nafion membrane. The water gradient causes the alteration of other physical properties within the thickness of the membrane and has a drastic impact on the performance of the devices made of the Nafion membrane. Deriving closed-form equations and using molecular dynamics (MD) simulation results, we bridge Nafion properties at the atomic scale and the macroscopic behavior of the membrane within a hierarchical multi-scale model. Multiple discrete simulation cells are selected across the thickness of the membrane with a wide range of water contents as representative volume elements (RVEs). The present framework is able to quantitatively predict the macroscopic properties of Nafion with the nanometric resolution regarding the water gradient across the membrane.

Nonlinear membrane model for large amplitude vibration of single layer graphene sheets.

The nonlinear vibrational properties of single layer graphene sheets (SLGSs) are investigated using a membrane model. The nonlinear equation of motion is considered for the SLGSs by including the effects of stretching due to large amplitudes. The equation of motion is numerically solved utilizing the finite difference method for SLGSs with different initial and boundary conditions, sizes and pretensions. It is concluded that the nonlinear fundamental frequency of SLGSs increases by increasing the pretension and initial velocity. In addition, it is observed that an increase in the pretension weakens the effects of the initial velocity on the fundamental frequency, such that the fundamental frequency approximately becomes independent of the initial velocity. This is an important feature of the vibrating systems consisting of SLGSs which are used in the nano-electromechanical systems (NEMS), where resonators with a specific fundamental frequency and independent of the initial velocity are of interest.