Diffusive transport of nanoscale objects through cell membranes: a computational perspective

Diffusion coefficients and MSD measurements on curved membranes and porous media

We study some geometric aspects that influence the transport properties of particles that diffuse on curved surfaces. We compare different approaches to surface diffusion based on the Laplace-Beltrami operator adapted to predict concentration along entire membranes, confined subdomains along surfaces, or within porous media. Our goal is to summarize, firstly, how diffusion in these systems results in different types of diffusion coefficients and mean square displacement measurements, and secondly, how these two factors are affected by the concavity of the surface, the shape of the possible barriers or obstacles that form the available domains, the sinuosity, tortuosity, and constrictions of the trajectories and even how the observation plane affects the measurements of the diffusion. In addition to presenting a critical and organized comparison between different notions of MSD, in this review, we test the correspondence between theoretical predictions and numerical simulations by performing finite element simulations and illustrate some situations where diffusion theory can be applied. We briefly reviewed computational schemes for understanding surface diffusion and finally, discussed how this work contributes to understanding the role of surface diffusion transport properties in porous media and their relationship to other transport processes.

Transition between Different Diffusion Modes of Individual Lipids during the Membrane-Specific Action of As-CATH4 Peptides

The cell membrane permeabilization ability of immune defense antimicrobial peptides (AMPs) is widely applied in biomedicine. Although the mechanisms of peptide-membrane interactions have been widely investigated, analyses at the molecular level are still lacking. Herein, the membrane-specific action of a native AMP, As-CATH4, is investigated using a single-lipid tracking method in combination with live cell and model membrane assays conducted at different scales. The peptide-membrane interaction process is characterized by analyzing single-lipid diffusion behaviors. As-CATH4 exhibits potent antimicrobial activity through bacterial membrane permeabilization, with moderate cytotoxicity against mammalian cells. In-plane diffusion analyses of individual lipids show that the lipid molecules exhibit non-Gaussian and heterogeneous diffusion behaviors in both pristine and peptide-treated membranes, which can be decomposed into two Gaussian subgroups corresponding to normal- and slow-diffusive lipids. Assessment of the temporal evolution of these two diffusion modes of lipids reveal that the peptide action states of As-CATH4 include surface binding, transmembrane defect formation, and dynamic equilibrium. The action mechanisms of As-CATH4 at varying concentrations and against different membranes are distinguished. This work resolves the simultaneous mixed diffusion mechanisms of single lipids in biomimetic cell membranes, especially during dynamic membrane permeabilization by AMPs.

Effects of Shape on Interaction Dynamics of Tetrahedral Nanoplastics and the Cell Membrane

Cellular uptake of nanoplastics is instrumental in their environmental accumulation and transfer to humans through the food chain. Despite extensive studies using spherical plastic nanoparticles, the influence of the morphological characteristics of environmentally released nanoplastics is understudied. Using dissipative particle dynamics simulations, we modeled the interactions between a cell membrane and hydrophobic nanotetrahedra, which feature high shape anisotropy and large surface curvature seen for environmental nanoplastics. We observe robust uptake of nanotetrahedra with sharp vertices and edges by the lipid membrane. Two local energy minimum configurations of nanotetrahedra embedded in the membrane bilayer were identified for particles of large sizes. Further analysis of particle dynamics within the membrane shows that the two interaction states exhibit distinct translational and rotational dynamics in the directions normal and parallel to the plane of the membrane. The membrane confinement significantly arrests the out-of-plane motion, resulting in caged translation and subdiffusive rotation. While the in-plane diffusion remains Brownian, we find that the translational and rotational modes decouple from each other as the particle size increases. The rotational diffusion decreases by a greater extent compared to the translational diffusion, deviating from the continuum theory predictions. These results provide fundamental insights into the shape effect on the nanoparticle dynamics in crowded lipid membranes.

Large stress asymmetries of lipid bilayers and nanovesicles generate lipid flip-flops and bilayer instabilities

Much effort has been devoted to lipid bilayers and nanovesicles with a compositional asymmetry between the two leaflets of the bilayer membranes. Here, we address another fundamental asymmetry related to lipid densities and membrane tensions. To avoid membrane rupture, the osmotic conditions must be adjusted in such a way that the bilayer membranes are subject to a relatively low bilayer tension. However, even for vanishing bilayer tension, the individual leaflets can still experience significant leaflet tensions if one leaflet is stretched whereas the other leaflet is compressed. Such a stress asymmetry between the two leaflets can be directly controlled in molecular dynamics simulations by the initial assembly of the lipid bilayers. This stress asymmetry is varied here over a wide range to determine the stability and instability regimes of the asymmetric bilayers. The stability regime shrinks with decreasing size and increasing membrane curvature of the nanovesicle. In the instability regimes, the lipids undergo stress-induced flip-flops with a flip-flop rate that increases with increasing stress asymmetry. The onset of flip-flops can be characterized by a cumulative distribution function that is well-fitted by an exponential function for planar bilayers but has a sigmoidal shape for nanovesicles. In addition, the bilayer membranes form transient non-bilayer structures that relax back towards ordered bilayers with a reduced stress asymmetry. Our study reveals intrinsic limits for the possible magnitude of the transbilayer stress asymmetry and shows that the leaflet tensions represent key parameters for the flip-flop rates.

Copolyelectrolyte-Based Nanocapsules for Oral Monoclonal Antibody Therapy: A Mesoscale Modeling Survey

Antibody therapy generally requires parenteral injection to attain the required bioavailability and pharmacokinetics, but improved formulations may slow enzymatic degradation of the antibody in the gastrointestinal tract, permitting the use of noninvasive oral delivery. Rationally designed carrier materials can potentially improve therapeutic activity both by shielding fragile biopharmaceuticals from proteolytic degradation and targeting specific receptors in vivo. One potentially useful class of protein carriers is block copolyelectrolytes, one polyelectrolyte plus one neutral hydrophilic polymer block, that self-assemble into stable micelles, providing a simple and biocompatible nanocapsule separating the protein from the outer medium. Here, we develop and implement an integrated mesoscale model to design molecular structures for block copolyelectrolyte nanocapsules predicted to protect Trastuzumab, an antibody used to treat breast cancer, in the low pH gastrointestinal tract and to selectively release this antibody in the more neutral intestinal environment. The simulations show a tightly packed self-assembled core-shell structure at pH = 3 that is ruptured and dynamically reassembled into a weaker structure at pH = 7. Our model identifies that the designed block copolyelectrolyte characteristics, such as block length ratio, can control the level of drug protection and release in vivo, providing simple design rules for engineering polyelectrolyte-based formulations that may allow oral administration of targeted antibody chemotherapies.

Comparative study of SARS-CoV-2 infection in different cell types: Biophysical-computational approach to the role of potential receptors

Cellular susceptibility to SARS-CoV-2 infection in the respiratory tract has been associated with the ability of the virus to interact with potential receptors on the host membrane. We have modeled viral dynamics by simulating various cellular systems and artificial conditions, including macromolecular crowding, based on experimental and transcriptomic data to infer parameters associated with viral growth and predict cell susceptibility. We have accomplished this based on the type, number and level of expression of the angiotensin-converting enzyme 2 (ACE2), transmembrane serine 2 (TMPRSS2), basigin2 (CD147), FURIN protease, neuropilin 1 (NRP1) or other less studied candidate receptors such as heat shock protein A5 (HSPA5) and angiotensin II receptor type 2 (AGTR2). In parallel, we studied the effect of simulated artificial environments on the accessibility to said proposed receptors. In addition, viral kinetic behavior dependent on the degree of cellular susceptibility was predicted. The latter was observed to be more influenced by the type of proteins and expression level, than by the number of potential proteins associated with the SARS CoV-2 infection. We predict a greater theoretical propensity to susceptibility in cell lines such as NTERA-2, SCLC-21H, HepG2 and Vero6, and a lower theoretical propensity in lines such as CaLu3, RT4, HEK293, A549 and U-251MG. An important relationship was observed between expression levels, protein diffusivity, and thermodynamically favorable interactions between host proteins and the viral spike, suggesting potential sites of early infection other than the lungs. This research is expected to stimulate future quantitative experiments and promote systematic investigation of the effect of crowding presented here.

RNA Secondary Structures Regulate Adsorption of Fragments onto Flat Substrates

RNA is a functionally rich molecule with multilevel, hierarchical structures whose role in the adsorption to molecular substrates is only beginning to be elucidated. Here, we introduce a multiscale simulation approach that combines a tractable coarse-grained RNA structural model with an interaction potential of a structureless flat adsorbing substrate. Within this approach, we study the specific role of stem-hairpin and multibranch RNA secondary structure motifs on its adsorption phenomenology. Our findings identify a dual regime of adsorption for short RNA fragments with and without the secondary structure and underline the adsorption efficiency in both cases as a function of the surface interaction strength. The observed behavior results from an interplay between the number of contacts formed at the surface and the conformational entropy of the RNA molecule. The adsorption phenomenology of RNA seems to persist also for much longer RNAs as qualitatively observed by comparing the trends of our simulations with a theoretical approach based on an ideal semiflexible polymer chain.

Durotaxis of Passive Nanoparticles on Elastic Membranes

The transport of macromolecules and nanoscopic particles to a target cellular site is a crucial aspect in many physiological processes. This directional motion is generally controlled via active mechanical and chemical processes. Here we show, by means of molecular dynamics simulations and an analytical theory, that completely passive nanoparticles can exhibit directional motion when embedded in nonuniform mechanical environments. Specifically, we study the motion of a passive nanoparticle adhering to a mechanically nonuniform elastic membrane. We observe a nonmonotonic affinity of the particle to the membrane as a function of the membrane's rigidity, which results in the particle transport. This transport can be both up or down the rigidity gradient, depending on the absolute values of the rigidities that the gradient spans across. We conclude that rigidity gradients can be used to direct average motion of passive macromolecules and nanoparticles on deformable membranes, resulting in the preferential accumulation of the macromolecules in regions of certain mechanical properties.

Durotaxis of Passive Nanoparticles on Elastic Membranes

The transport of macromolecules and nanoscopic particles to a target cellular site is a crucial aspect in many physiological processes. This directional motion is generally controlled via active mechanical and chemical processes. Here we show, by means of molecular dynamics simulations and an analytical theory, that completely passive nanoparticles can exhibit directional motion when embedded in nonuniform mechanical environments. Specifically, we study the motion of a passive nanoparticle adhering to a mechanically nonuniform elastic membrane. We observe a nonmonotonic affinity of the particle to the membrane as a function of the membrane's rigidity, which results in the particle transport. This transport can be both up or down the rigidity gradient, depending on the absolute values of the rigidities that the gradient spans across. We conclude that rigidity gradients can be used to direct average motion of passive macromolecules and nanoparticles on deformable membranes, resulting in the preferential accumulation of the macromolecules in regions of certain mechanical properties.

Tracking of Nanoparticle Diffusion at a Liquid-Liquid Interface Adsorbed by Nonionic Surfactants

Emulsions stabilized by both nanoparticles and surfactants often display longer shelf life than those stabilized by nanoparticles or surfactants alone. Although numerous works have been conducted to understand the effect of nanoparticles and surfactants on the variation of interfacial tension, little is known about interfacial diffusion when both nanoparticles and surfactants are present at interfaces. In this work, we used single-particle fluorescence tracking to study the lateral diffusion of individual hydrophobic nanoparticles at hexane-glycerol interfaces adsorbed by different amounts of nonionic surfactants. When the surfactant concentration is over a threshold, we found that the nanoparticle diffusion exhibits a two-regime behavior involving short-time Brownian and the emergence of subdiffusive, non-Gaussian, and dynamically anticorrelated diffusion in the long lag time regime. A stepwise analysis rationalized diffusion in different lag time regimes, leading to a mechanistic interpretation regarding the two-regime behavior. These results could provide insight into the understanding of the synergistic effect for the surfactant-assistant Pickering emulsion.

Experimental and Theoretical Analysis of Metal Complex Diffusion through Cell Monolayer

The study of drugs diffusion through different biological membranes constitutes an essential step in the development of new pharmaceuticals. In this study, the method based on the monolayer cell culture of CHO-K1 cells has been developed in order to emulate the epithelial cells barrier in permeability studies by laser interferometry. Laser interferometry was employed for the experimental analysis of nickel(II) and cobalt(II) complexes with 1-allylimidazole or their chlorides’ diffusion through eukaryotic cell monolayers. The amount (mol) of nickel(II) and cobalt(II) chlorides transported through the monolayer was greater than that of metals complexed with 1-allylimidazole by 4.34-fold and 1.45-fold, respectively, after 60 min. Thus, laser interferometry can be used for the quantitative analysis of the transport of compounds through eukaryotic cell monolayers, and the resulting parameters can be used to formulate a mathematical description of this process.

Sliding dynamics of multi-rings on a semiflexible polymer in poly[n]catenanes

The sliding dynamics of one- or multi-ring structures along a semiflexible cyclic polymer in radial poly[n]catenanes is investigated using molecular dynamics simulations. The fixed and fluctuating (non-fixed) semiflexible central cyclic polymers are considered, respectively. With increasing bending energy of the central cyclic polymer, for the fixed case, the diffusion coefficient increases monotonically due to the reduction of the tortuous sliding path, while for the fluctuating case, the diffusion coefficient decreases. This indicates that the contribution of the polymer fluctuation is suppressed by a further increase in the stiffness of the central cyclic chain. Compared with the one ring case, the mean-square displacement of the multiple rings exhibits a unique sub-diffusive behavior at intermediate time scales due to the repulsion between two neighboring rings. In addition, for the multi-ring system, the whole set of rings exhibit relatively slower diffusion, but faster local dynamics of threading rings and rotational diffusion of the central cyclic polymer arise. These results may help us to understand the diffusion motion of rings in radial poly[n]catenanes from a fundamental point of view and control the sliding dynamics in molecular designs.

Regulation Mechanism of Bubbling Deformation and Fracture Toughness of the Membrane by Asymmetric Phospholipids: A Model System Study

Dynamic regulation of the deformation modulus and fracture toughness of a membrane is critical to organelles and cells for matching their conflicting needs of resilient and fractured behaviors. These properties implement the protection of the function in the normal condition and the fission function in the endocytosis condition of a membrane. Naturally, a membrane contains phospholipids that have different hydrophilic and hydrophobic group length. The diffusion and aggregation of the phospholipids with asymmetry of the hydrophilic-hydrophobic ratio on the membrane play a key role in regulating the mechanical behaviors passively to the external force. In present work, the effects of the asymmetry of phospholipids on the bubbling deformation and fracture toughness of the membrane to external stretching are investigated in a model system. A disk-shaped micelle formed from the blend of symmetric and asymmetric diblock copolymers in a selective solvent is considered as the membrane sheet. Its mechanically responsive behaviors are investigated by self-consistent field theory. By analyzing the evolution of different components during the stretching process, the mechanism of formation of the bubbling structure is revealed. Moreover, the fracture toughness depending on the asymmetry of the phospholipids is determined quantitatively.

Entropy at Bio-Nano Interfaces

Entropy, one of the central concepts of thermodynamics, can be a predominant contribution to structural formation and transition. Although it is well-known that diverse forces and energies can significantly contribute to the structures and activities at bio-nano interfaces, the potential entropic contribution remains less well understood. Therefore, this review article seeks to provide a conceptual framework demonstrating that entropy can be exploited to shape the physicochemical properties of bio-nano interfaces and thereby regulate the structures, responses, and functions of biological systems. We introduce the typical types of entropy that matter at bio-nano interfaces. Moreover, some key characteristics featuring entropy at bio-nano interfaces, such as the difference between entropic force and energetic interaction and the associated implications for biomimetic research, are discussed. We expect that this review could stimulate further effort in the fundamental research of entropy in biology and in the biological applications of entropic effects in designer biomaterials.

Correlation between Single-Molecule Dynamics and Biological Functions of Antimicrobial Peptide Melittin

Many fundamental biological processes occur on cell membranes, and a typical example is the membrane permeabilization by peptides for an antimicrobial purpose. Previous studies of the underlying mechanism mostly focus on structural changes of membranes and peptides during their interactions. Herein, from a new perspective of single-molecule dynamics, the real-time three-dimensional motions of individual phospholipid and peptide molecules were monitored, and specifically, their correlation with the membrane poration function of melittin, a most representative natural antimicrobial peptide, was studied. We found that the adsorption and accumulation of melittin on the membrane surface significantly sped up the lateral diffusion of lipids surrounding the peptides, which in turn facilitated the peptide insertion at such heterogeneous regions. A unique "U"-bending pathway of melittin during membrane insertion and the ultimate formation of toroidal pores with dynamical translocations of peptides and lipids with several metastable states between the two leaflets of bilayer were observed.