Source author record

Oded Farago

Oded Farago appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

30works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

30 published item(s)

preprint2020arXiv

'Rocket propulsion' of Janus micro-swimmers

We report simulations of a spherical Janus particle undergoing exothermic surface reactions around one pole only. Our model excludes self-phoretic transport by design. Nevertheless, net motion occurs from direct momentum transfer between solvent and colloid, with speed scaling as the square root of the energy released during the reaction. We find that such propulsion is dominated by the system's short-time response, when neither the time dependence of the flow around the colloid nor the solvent compressibility can be ignored. Our simulations agree reasonably well with previous experiments.

preprint2020arXiv

A Langevin dynamics approach for multi-layer mass transfer problems

We use Langevin dynamics simulations to study the mass diffusion problem across two adjacent porous layers of different transport property. At the interface between the layers, we impose the Kedem-Katchalsky (KK) interfacial boundary condition that is well suited in a general situation. A detailed algorithm for the implementation of the KK interfacial condition in the Langevin dynamics framework is presented. As a case study, we consider a two-layer diffusion model of a drug-eluting stent. The simulation results are compared with those obtained from the solution of the corresponding continuum diffusion equation, and an excellent agreement is shown.

preprint2020arXiv

Defining velocities for accurate kinetic statistics in the GJF thermostat

We expand on two previous developments in the modeling of discrete-time Langevin systems. One is the well-documented Grønbech-Jensen Farago (GJF) thermostat, which has been demonstrated to give robust and accurate configurational sampling of the phase space. Another is the recent discovery that also kinetics can be accurately sampled for the GJF method. Through a complete investigation of all possible finite difference approximations to the velocity, we arrive at two main conclusions:~1) It is not possible to define a so-called on-site velocity such that kinetic temperature will be correct and independent of the time step, and~2) there exists a set of infinitely many possibilities for defining a two-point (leap-frog) velocity that measures kinetic energy correctly for linear systems in addition to the correct configurational statistics obtained from the GJF algorithm. We give explicit expressions for the possible definitions, and we incorporate these into convenient and practical algorithmic forms of the normal Verlet-type algorithms along with a set of suggested criteria for selecting a useful definition of velocity.

preprint2016arXiv

Formation of dilute adhesion domains driven by weak elasticity-mediated interactions

Cell-cell adhesion is established by specific binding of receptor and ligand proteins. The adhesion bonds attract each other and often aggregate into large clusters that are central to many biological processes. One possible origin of attractive interactions between adhesion bonds is the elastic response of the membranes to their deformation by the bonds. Here, we analyze these elasticity-mediated interactions using a novel mean-field approach. Analysis of systems at different densities of bonds, $ϕ$, reveals that the phase diagram exhibits a nearly-universal behavior when the temperature $T$ is plotted vs. the scaled density $x=ϕξ^2$, where $ξ$ is the linear size of the membrane's region affected by a single bond. The critical point $(ϕ_c,T_c)$ is located at very low densities and slightly below $T_c$ we identify phase coexistence between two low-density phases. Dense domains are observed only when the height by which the bonds deform the membranes, $h_0$, is much larger than their thermal roughness, $Δ$, which occurs at very low temperatures $T\ll T_c$. We conclude that the elasticity-mediated interactions are weak and cannot be regarded as responsible for the formation of dense adhesion domains. The weakness of the elasticity-mediated effect and its relevance to dilute systems only can be attributed to the fact that the membrane's elastic energy saturates in the semi-dilute regime, when the typical spacing between the bonds $r\gtrsim ξ$, i.e., for $x\lesssim 1$. Therefore, at higher densities, only the mixing entropy of the bonds (which favors uniform distributions) is thermodynamically relevant. We discuss the implications of our results to the question of immunological synapse formation, and demonstrate the elasticity-mediated interactions may be involved in the aggregation of these semi-dilute membrane domains.

preprint2016arXiv

Isothermal Langevin dynamics in systems with power-law spatially-dependent friction

We study the dynamics of Brownian particles in a heterogeneous one-dimensional medium with a spatially-dependent diffusion coefficient of the form $D(x)\sim |x|^c$, at constant temperature. The particle's probability distribution function (PDF) is calculated both analytically, by solving Fick's diffusion equation, and from numerical simulations of the underdamped Langevin equation. At large times, the PDFs calculated by both approaches yield identical results, corresponding to subdiffusion for $c<0$, and superdiffusion for $0<c<1$. For $c>1$, the diffusion equation predicts that the particles accelerate. Here, we show that this phenomenon, previously considered in several works as an illustration for the possible dramatic effects of spatially-dependent thermal noise, is unphysical. We argue that in an isothermal medium, the motion cannot exceed the ballistic limit ($\left\langle x^2\right\rangle \sim t^2$). The ballistic limit is reached when the friction coefficient drops sufficiently fast at large distances from the origin, and is correctly captured by Langevin's equation.

preprint2016arXiv

On the connection between dissipative particle dynamics and the Itô-Stratonovich dilemma

Dissipative Particle Dynamics (DPD) is a popular simulation model for investigating hydrodynamic behavior of systems with non-negligible equilibrium thermal fluctuations. DPD employs soft core repulsive interactions between the system particles, thus allowing them to overlap. This supposedly permits relatively large integration time steps, which is an important feature for simulations on large temporal scales. In practice, however, an increase in the integration time step leads to increasingly larger systematic errors in the sampling statistics. Here, we demonstrate that the prime origin of these systematic errors is the multiplicative nature of the thermal noise term in Langevin's equation, i.e., the fact that it depends on the instantaneous coordinates of the particles. This lead to an ambiguity in the interpretation of the stochastic differential Langevin equation, known as the Itô-Stratonovich dilemma. Based on insights from previous studies of the dilemma, we propose a novel algorithm for DPD simulations exhibiting almost an order of magnitude improvement in accuracy, and nearly twice the efficiency of commonly used DPD Langevin thermostats.

preprint2016arXiv

The G-JF Thermostat for Accurate Configurational Sampling in Soft-Matter Simulations

We implement the statistically sound G-JF thermostat for Langevin Dynamics simulations into the ESPREesSo molecular package for large-scale simulations of soft matter systems. The implemented integration method is tested against the integrator currently used by the molecular package in simulations of a fluid bilayer membrane. While the latter exhibits deviations in the sampling statistics that increase with the integration time step dt, the former reproduces near-correct configurational statistics for all dt within the stability range of the simulations. We conclude that, with very modest revisions to existing codes, one can significantly improve the performance of statistical sampling using Langevin thermostats.

preprint2016arXiv

The Thermodynamics of Endosomal Escape and DNA Release from Lipoplexes

Complexes of cationic and neutral lipids and DNA (lipoplexes) are emerging as promising vectors for gene therapy applications. Their appeal stems from their non pathogenic nature and the fact that they self-assemble under conditions of thermal equilibrium. Lipoplex adhesion to the cell plasma membrane initiates a three-stage process termed transfection, consisting of (i) endocytosis, (ii) lipoplex breakdown, and (iii) DNA release followed by gene expression. As successful transfection requires lipoplex degradation, it tends to be hindered by the lipoplex thermodynamic stability; nevertheless, it is known that the transfection process may proceed spontaneously. Here, we use a simple model to study the thermodynamic driving forces governing transfection. We demonstrate that after endocytosis [stage (i)], the lipoplex becomes inherently unstable. This instability, which is triggered by interactions between the cationic lipids of the lipoplex and the anionic lipids of the enveloping plasma membrane, is entropically controlled involving both remixing of the lipids and counterions release. Our detailed calculation shows that the free energy gain during stage (ii) is approximately linear in $Φ_+$, the mole fraction of cationic lipids in the lipoplex. This free energy gain, $ΔF$, reduces the barrier for fusion between the enveloping and the lipoplex bilayers, which produces a hole allowing for DNA release [stage (iii)]. The linear relationship between $ΔF$ and the fraction of cationic lipids explains the experimentally observed exponential increase of transfection efficiency with $Φ_+$ in lamellar lipoplexes.

preprint2015arXiv

Formation of adhesion domains in stressed and confined membranes

The adhesion bonds connecting a lipid bilayer to an underlying surface may undergo a condensation transition resulting from an interplay between a short range attractive potential between them, and a long range fluctuation-induced potential of mean force. Here, we use computer simulations of a coarse-grained molecular model of supported lipid bilayers to study this transition in confined membranes, and in membranes subjected to a non-vanishing surface tension. Our results show that confinement may alter significantly the condensation transition of the adhesion bonds, whereas the application of surface tension has a very minor effect on it. We also investigate domain formation in membranes under negative tension which, in free membranes, causes enhancement of the amplitude of the membrane thermal undulations. Our results indicate that in supported membranes, this effect of a negative surface tension on the fluctuation spectrum is largely eliminated by the pressure resulting from the mixing entropy of the adhesion bonds.

preprint2015arXiv

Small membranes under negative surface tension

We use computer simulations and a simple free energy model to study the response of a bilayer membrane to the application of a negative (compressive) mechanical tension. Such a tension destabilizes the long wavelength undulation modes of giant vesicles, but it can be sustained when small membranes and vesicles are considered. Our negative tension simulation results reveal two regimes - (i) a weak negative tension regime characterized by stretching-dominated elasticity, and (ii) a strong negative tension regime featuring bending-dominated elastic behavior. This resembles the findings of the classic Evans and Rawicz micropipette aspiration experiment in giant unilamellar vesicles (GUVs) [Phys, Rev. Lett. {\bf 64}, 2094 (1990)]. However, while in GUVs the crossover between the two elasticity regimes occurs at a small positive surface tension, in smaller membranes it takes place at a moderate negative tension. Another interesting observation concerning the response of a small membrane to negative surface tension is related to the relationship between the mechanical and fluctuation tensions, which are equal to each other for non-negative values. When the tension decreases to negative values, the fluctuation tension $γ$ drops somewhat faster than $τ$ in the small negative tension regime, before it saturates (and becomes larger than $τ$) for large negative tensions. The bending modulus exhibit an "opposite" trend. It remains almost unchanged in the stretching-dominated elastic regime, and decreases in the bending-dominated regime. Both the amplitudes of the thermal height undulations and the projected area variations diverge at the onset of mechanical instability.

preprint2014arXiv

Coarse grained molecular simulations of membrane adhesion domains

We use a coarse grained molecular model of supported lipid bilayers to study the formation of adhesion domains. We find that this process is a first order phase transition, triggered by a combination of pairwise short range attractive interactions between the adhesion bonds and many-body Casimir-like interactions, mediated by the membrane thermal undulations. The simulation results display an excellent agreement with the recently proposed Weil-Farago 2D lattice model, in which the occupied and empty sites represent, respectively, the adhesion bonds and unbound segments of the membrane. A second phase transition, into a hexatic phase, is observed when the attraction between the adhesion bonds is further strengthened.

preprint2014arXiv

Constant pressure and temperature discrete-time Langevin molecular dynamics

We present a new and improved method for simultaneous control of temperature and pressure in molecular dynamics simulations with periodic boundary conditions. The thermostat-barostat equations are build on our previously developed stochastic thermostat, which has been shown to provide correct statistical configurational sampling for any time step that yields stable trajectories. Here, we extend the method and develop a set of discrete-time equations of motion for both particle dynamics and system volume in order to seek pressure control that is insensitive to the choice of the numerical time step. The resulting method is simple, practical, and efficient. The method is demonstrated through direct numerical simulations of two characteristic model systems - a one dimensional particle chain for which exact statistical results can be obtained and used as benchmarks, and a three dimensional system of Lennard-Jones interacting particles simulated in both solid and liquid phases. The results, which are compared against the method of Kolb & Dunweg, show that the new method behaves according to the objective, namely that acquired statistical averages and fluctuations of configurational measures are accurate and robust against the chosen time step applied to the simulation.

preprint2014arXiv

Elasticity and mechanical instability of charged lipid bilayers in ionic solutions

We use coarse-grained Monte Carlo simulations to study the elastic properties of charged membranes in solutions of monovalent and pentavalent counterions. The simulation results of the two cases reveal trends opposite to each other. The bending rigidity and projected area increase with the membrane charge density for monovalent counterions, while they decrease for the pentavalent ions. These observations can be related to the counterion screening of the lipid charges. While the monovalent counterions only weakly screen the Coulomb interactions, which implies a repulsive Coulomb system, the multivalent counterions condense on the membrane and, through spatial charge correlations, make the effective interactions due the charged lipids attractive. The differences in the elastic properties of the charged membranes in monovalent and multivalent counterion solutions are reflected in the mechanisms leading to their mechanical instability at high charge densities. In the former case, the membranes develop pores to relieve the electrostatic tensile stresses, while in the latter case, the membrane exhibit large wavelength bending instability.

preprint2014arXiv

Fluctuation-dissipation relation for systems with spatially varying friction

When a particle diffuses in a medium with spatially dependent friction coefficient $α(r)$ at constant temperature $T$, it drifts toward the low friction end of the system even in the absence of any real physical force $f$. This phenomenon, which has been previously studied in the context of non-inertial Brownian dynamics, is termed "spurious drift", although the drift is real and stems from an inertial effect taking place at the short temporal scales. Here, we study the diffusion of particles in inhomogeneous media within the framework of the inertial Langevin equation. We demonstrate that the quantity which characterizes the dynamics with non-uniform $α(r)$ is not the displacement of the particle $Δr=r-r^0$ (where $r^0$ is the initial position), but rather $ΔA(r)=A(r)-A(r^0)$, where $A(r)$ is the primitive function of $α(r)$. We derive expressions relating the mean and variance of $ΔA$ to $f$, $T$, and the duration of the dynamics $Δt$. For a constant friction coefficient $α(r)=α$, these expressions reduce to the well known forms of the force-drift and fluctuation-dissipation relations. We introduce a very accurate method for Langevin dynamics simulations in systems with spatially varying $α(r)$, and use the method to validate the newly derived expressions.

preprint2013arXiv

A simple and effective Verlet-type algorithm for simulating Langevin dynamics

We present a revision to the well known Stormer-Verlet algorithm for simulating second order differential equations. The revision addresses the inclusion of linear friction with associated stochastic noise, and we analytically demonstrate that the new algorithm correctly reproduces diffusive behavior of a particle in a flat potential. For a harmonic oscillator, our algorithm provides the exact Boltzmann distribution for any value of damping, frequency, and time step for both underdamped and over damped behavior within the usual the stability limit of the Verlet algorithm. Given the structure and simplicity of the method we conclude this approach can trivially be adapted for contemporary applications, including molecular dynamics with extensions such as molecular constraints.

preprint2013arXiv

Application of the G-JF Discrete-Time Thermostat for Fast and Accurate Molecular Simulations

A new Langevin-Verlet thermostat that preserves the fluctuation-dissipation relationship for discrete time steps, is applied to molecular modeling and tested against several popular suites (AMBER, GROMACS, LAMMPS) using a small molecule as an example that can be easily simulated by all three packages. Contrary to existing methods, the new thermostat exhibits no detectable changes in the sampling statistics as the time step is varied in the entire numerical stability range. The simple form of the method, which we express in the three common forms (Velocity-Explicit, Stormer-Verlet, and Leap-Frog), allows for easy implementation within existing molecular simulation packages to achieve faster and more accurate results with no cost in either computing time or programming complexity.

preprint2013arXiv

Langevin dynamics in inhomogeneous media: Re-examining the Itô-Stratonovich Dilemma

The diffusive dynamics of a particle in a medium with space-dependent friction coefficient is studied within the framework of the inertial Langevin equation. In this description, the ambiguous interpretation of the stochastic integral, known as the Itô-Stratonovich dilemma, is avoided since all interpretations converge to the same solution in the limit of small time steps. We use a newly developed method for Langevin simulations to measure the probability distribution of a particle diffusing in a flat potential. Our results reveal that both the Itô and Stratonovich interpretations converge very slowly to the uniform equilibrium distribution for vanishing time step sizes. Three other conventions exhibit significantly improved accuracy: (i) the "isothermal" (Hänggi) convention, (ii) the Stratonovich convention corrected by a drift term, and (iii) a new convention employing two different effective friction coefficients representing two different averages of the friction function during the time step. We argue that the most physically accurate dynamical description is provided by the third convention, in which the particle experiences a drift originating from the dissipation instead of the fluctuation term. This feature is directly related to the fact that the drift is a result of an inertial effect that cannot be well understood in the Brownian, overdamped limit of the Langevin equation.

preprint2013arXiv

Muscle contraction and the elasticity-mediated crosstalk effect

Cooperative action of molecular motors is essential for many cellular processes. One possible regulator of motor coordination is the elasticity-mediated crosstalk (EMC) coupling between myosin II motors whose origin is the tensile stress that they collectively generate in actin filaments. Here, we use a statistical mechanical analysis to investigate the influence of the EMC effect on the sarcomere - the basic contractile unit of skeletal muscles. We demonstrate that the EMC effect leads to an increase in the attachment probability of motors located near the end of the sarcomere while, simultaneously, decreasing the attachment probability of the motors in the central part. Such a polarized attachment probability would impair the motors ability to cooperate efficiently. Interestingly, this undesired phenomenon becomes significant only when the system size exceeds that of the sarcomere in skeletal muscles, which provides an explanation for the remarkable lack of sarcomere variability in vertebrates. Another phenomenon that we investigate is the recently observed increase in the duty ratio of the motors with the tension in muscle. We reveal that the celebrated Hill's equation for muscle contraction is very closely related to this observation.

preprint2012arXiv

Duty-ratio of cooperative molecular motors

Molecular motors are found throughout the cells of the human body, and have many different and important roles. These micro-machines move along filament tracks, and have the ability to convert chemical energy into mechanical work that powers cellular motility. Different types of motors are characterized by different duty-ratios, which is the fraction of time that a motor is attached to its filament. In the case of myosin II - a non-processive molecular machine with a low duty ratio - cooperativity between several motors is essential to induce motion along its actin filament track. In this work we use statistical mechanical tools to calculate the duty ratio of cooperative molecular motors. The model suggests that the effective duty ratio of non-processive motors that work in cooperation is lower than the duty ratio of the individual motors. The origin of this effect is the elastic tension that develops in the filament which is relieved when motors detach from the track.

preprint2012arXiv

Hierarchical self-organization of cytoskeletal active networks

The structural reorganization of the actin cytoskeleton is facilitated through the action of motor proteins that crosslink the actin filaments and transport them relative to each other. Here, we present a combined experimental-computational study that probes the dynamic evolution of mixtures of actin filaments and clusters of myosin motors. While on small spatial and temporal scales the system behaves in a very noisy manner, on larger scales it evolves into several well distinct patterns such as bundles, asters, and networks. These patterns are characterized by junctions with high connectivity, whose formation is possible due to the organization of the motors in "oligoclusters" (intermediate-size aggregates). The simulations reveal that the self-organization process proceeds through a series of hierarchical steps, starting from local microscopic moves and ranging up to the macroscopic large scales where the steady-state structures are formed. Our results shed light into the mechanisms involved in processes like cytokinesis and cellular contractility, where myosin motors organized in clusters operate cooperatively to induce structural organization of cytoskeletal networks.

preprint2011arXiv

Crosstalk between non-processive myosin motors mediated by the actin filament elasticity

Many biological processes involve the action of molecular motors that interact with the cell cytoskeleton. Some processes, such as the transport of cargoes is achieved mainly by the action of individual motors. Other, such as cell motility and division, require the cooperative work of many motors. Collective motor dynamics can be quite complex and unexpected. One beautiful example is the bidirectional ("back and forth") motion of filaments which is induced when the motors within a group exert forces in opposite directions. This review tackles the puzzle emerging from a recent experimental work in which it has been shown that the characteristic reversal times of the bidirectional motion are practically independent of the number of motors. This result is in a striking contradiction with existing theoretical models that predict an exponential growth of the reversal times with the size of the system. We argue that the solution to this puzzle may be the crosstalk between the motors which is mediated by the elastic tensile stress that develops in the cytoskeleton track. The crosstalk does not directly correlate the attachment and detachment of the motors, which work independently of each other. However, it highly accelerates their detachments by making the detachment rates system size dependent.

preprint2011arXiv

Entropic attraction of adhesion bonds toward cell boundaries

Adhesion bonds between membranes and surfaces are attracted to each other via effective interactions whose origin the entropy loss due to the reduction in the amplitude of the membrane thermal fluctuations in the vicinity of the adhesion bonds. These fluctuation-induced interactions are also expected to drive the adhesion bonds toward the rim of the cell, as well as toward the surfaces of membrane inclusions. In this paper, we analyze the attraction of adhesion bonds to the cell inner and outer boundaries. Our analysis shows that the probability distribution function of a single (diffusing) adhesion bond decay algebraically with the distance from the boundaries. Upon increasing the concentration of the adhesion bonds, the attraction to the boundaries becomes strongly self-screened.

preprint2011arXiv

Mechanical surface tension governs membrane thermal fluctuations

Motivated by the still ongoing debate about the various possible meanings of the term surface tension of bilayer membranes, we present here a detailed discussion that explains the differences between the "intrinsic", "renormalized", and "mechanical" tensions. We use analytical considerations and computer simulations to show that the membrane spectrum of thermal fluctuations is governed by the mechanical and not the intrinsic tension. Our study highlights the fact that the commonly used quadratic approximation of Helfrich effective Hamiltonian is not rotationally invariant. We demonstrate that this non-physical feature leads to a calculated mechanical tension that differs dramatically from the correct mechanical tension. Specifically, our results suggest that the mechanical and intrinsic tensions vanish simultaneously, which contradicts recent theoretical predictions derived for the approximated Hamiltonian.

preprint2011arXiv

Molecular simulation analysis of structural variations in lipoplexes

We use a coarse-grained molecular model to study the self-assembly process of complexes of cationic and neutral lipids with DNA molecules ("lipoplexes") - a promising nonviral carrier of DNA for gene therapy. We identify the resulting structures through direct visualization of the molecular arrangements and through calculations of the corresponding scattering plots. The latter approach provides a means for comparison with published data from X-ray scattering experiments. Consistent with experimental results, we find that upon increasing the stiffness of the lipid material, the system tends to form lamellar structures. Two characteristic distances can be extracted from the scattering plots of lamellar complexes - the lamellar (interlayer) spacing and the DNA-spacing within each layer. We find a remarkable agreement between the computed values of these two quantities and the experimental data [J. O. Rädler, I. Koltover, T. Salditt and C. R. Safinya, Science Vol. 275, 810 (1997)] over the entire range of mole fractions of charged lipids (CLs) studied experimentally. A visual inspection of the simulated systems reveals that, for very high fractions of CLs, disordered structures consisting of DNA molecules bound to small membrane fragments are spontaneously formed. The diffraction plots of these non-lamellar disordered complexes appear very similar to that of the lamellar structure, which makes the interpretation of the X-ray data ambiguous. The loss of lamellar order may be the origin of the observed increase in the efficiency of lipoplexes as gene delivery vectors at high charge densities.

preprint2011arXiv

Statistical thermodynamics of adhesion points in supported membranes

Supported lipid membranes are useful and important model systems for studying cell membrane properties and membrane mediated processes. One attractive application of supported membranes is the design of phantom cells exhibiting well defined adhesive properties and receptor densities. Adhesion of membranes may be achieved by specific and non-specific interactions, and typically requires the clustering of many adhesion bonds into "adhesion domains". One potential mediator of the early stages of the aggregation process is the Casimir-type forces between adhesion sites induced by the membrane thermal fluctuations. In this review, I will present a theoretical analysis of fluctuation induced aggregation of adhesion sites in supported membranes. I will first discuss the influence of a single attachment point on the spectrum of membrane thermal fluctuations, from which the free energy cost of the attachment point will be deduced. I will then analyze the problem of a supported membrane with two adhesion points. Using scaling arguments and Monte Carlo simulations, I will demonstrate that two adhesion points attract each other via an infinitely long range effective potential that grows logarithmically with the pair distance. Finally, I will discuss the many-body nature of the fluctuation induced interactions. I will show that while these interactions alone are not sufficient to allow the formation of aggregation clusters, they greatly reduce the strength of the residual interactions required to facilitate cluster formation. Specifically, for adhesion molecules interacting via a short range attractive potential, the strength of the direct interactions required for aggregation is reduced by about a factor of two to below the thermal energy $k_BT$.

preprint2010arXiv

Biased transport of elastic cytoskeletal filaments with alternating polarities by molecular motors

We present a simple model for the bidirectional dynamics of actin bundles with alternating polarities in gliding assays with non-processive myosin motors. In the model, the bundle is represented as an elastic chain consisting of monomers with positive and negative polarities. The motion of the bundle is induced by the pulling forces of the underlying motors which stochastically attach to the monomers and, depending on the polarity of the monomers, pull them in the right or left direction. We demonstrate that perfectly a-polar chains consisting of equal numbers of monomers with positive and negative polarities may exhibit biased bidirectional motion with non-zero drift. This effect is attributed to the elastic tension developed in the chain due to the action of the myosin motors. We also show that as a result of this tension, the attachment probability of the motors is greatly reduced and becomes strongly dependent on the length of the chain. These surprising effects point to the necessity of considering the elasticity of the cytoskeleton in theoretical studies of cooperative dynamics of molecular motors.

preprint2010arXiv

Entropy driven aggregation of adhesion sites of supported membranes

We study, by means of mean field calculations and Monte Carlo simulations of a lattice-gas model, the distribution of adhesion sites of a bilayer membrane and a supporting flat surface. Our model accounts for the many-body character of the attractive interactions between adhesion points induced by the membrane thermal fluctuations. We show that while the fluctuation-mediated interactions alone are not sufficient to allow the formation of aggregation domains, they greatly reduce the strength of the direct interactions required to facilitate cluster formation. Specifically, for adhesion molecules interacting via a short range attractive potential, the strength of the direct interactions required for aggregation is reduced by about a factor of two to below the thermal energy $k_BT$.

preprint2010arXiv

Fluctuation induced attraction between adhesion sites of supported membranes

We use scaling arguments and coarse grained Monte Carlo simulations to study the fluctuation mediated interactions between a pair of adhesion sites of a bilayer membrane and a supporting surface. We find that the potential of mean force is an infinitely long range attractive potential that grows logarithmically with the pair distance $r$: $ϕ(\vec{r})/k_BT=c\ln r$, where the constant $c=2$ and $c=1$ for non-stressed and stressed membranes, respectively. When, in addition to excluded volume repulsion, the membrane also interacts with the underlying surface through a height-dependent attractive potential, the potential $ϕ(\vec{r})$ is screened at large pair distances.

preprint2009arXiv

Cooperative molecular motors moving back and forth

We use a two-state ratchet model to study the cooperative bidirectional motion of molecular motors on cytoskeletal tracks with randomly alternating polarities. Our model is based on a previously proposed model [Badoual et al., {\em Proc. Natl. Acad. Sci. USA} {\bf 99}, 6696 (2002)] for collective motor dynamics and, in addition, takes into account the cooperativity effect arising from the elastic tension that develops in the cytoskeletal track due to the joint action of the walking motors. We show, both computationally and analytically, that this additional cooperativity effect leads to a dramatic reduction in the characteristic reversal time of the bidirectional motion, especially in systems with a large number of motors. We also find that bidirectional motion takes place only on (almost) a-polar tracks, while on even slightly polar tracks the motion is unidirectional. We argue that the origin of these observations is the sensitive dependence of the cooperative dynamics on the difference between the number of motors typically working in and against the instantaneous direction of motion.

preprint2008arXiv

Transitions between distinct compaction regimes in complexes of multivalent cationic lipids and DNA

We use X-ray scattering and molecular simulations to investigate the structural properties of complexes of multivalent cationic lipids and DNA molecules. At low mole fraction of neutral lipids (NLs), $Φ_{\rm NL}$, the complexes show dramatic DNA compaction down to essentially close packed DNA arrays with a DNA interaxial spacing $d_{\rm DNA}=25Å$. A gradual increase in $Φ_{\rm NL}$ does not lead to a continuous increase in $d_{\rm DNA}$ as observed for DNA complexes of monovalent cationic lipids (CLs). Instead, distinct spacing regimes exist, with sharp transitions between them. Three packing states have been identified: (i) close packed, (ii) condensed, but not close packed, with $d_{\rm DNA}=27-28Å$, and (iii) an expanded state, where $d_{\rm DNA}$ increases gradually with $Φ_{\rm NL}$. Based on our experimental and computational results, we conclude that the DNA condensation is mediated by the multivalent cationic lipids, which assemble between the negatively charged DNA rods. Quite remarkably, the computational results show that the less tightly packed structure in regime (ii) is thermodynamically more stable than the close packed structure in regime (i). Accordingly, the constant DNA spacing observed in regime (ii) is attributed to lateral phase coexistence between this stable CL-DNA complex and neutral membranes. This finding may explain the reduced transfection efficiency measured for such complexes: Transfection involves endosomal escape and disassembly of the complex, and these processes are inhibited by the high thermodynamic stability. Our results, which demonstrate the existence of an inverse correlation between the stability and transfection activity of lamellar CL-DNA complexes are, therefore, consistent with a recently proposed model of cellular entry.