Direct measurement of repulsive and attractive pair potentials using pairs of optical traps.

We present a technique for measuring the interactions between pairs of colloidal particles in two optical traps. This method is particularly suitable for measuring strongly attractive potentials, an otherwise challenging task. The interaction energy is calculated from the distribution of inter-particle separations by accounting for the contribution from the optical traps with arbitrary trap profiles. The method is simple to implement and applicable to different types of pair potentials and optical trapping geometries. We apply the method to measure dipolar pair interactions in experiments with paramagnetic colloids in external magnetic fields. We obtain consistent and accurate results in all regimes, from strongly attractive to repulsive potentials. By means of computer simulations, we demonstrate that the proposed method can be successfully applied to systems with complex pair interactions characterized by multiple attractive and repulsive regimes, which are ubiquitous in soft and biological matter.

Data-driven modelling makes quantitative predictions regarding bacteria surface motility.

In this work, we quantitatively compare computer simulations and existing cell tracking data of P. aeruginosa surface motility in order to analyse the underlying motility mechanism. We present a three dimensional twitching motility model, that simulates the extension, retraction and surface association of individual Type IV Pili (TFP), and is informed by recent experimental observations of TFP. Sensitivity analysis is implemented to minimise the number of model parameters, and quantitative estimates for the remaining parameters are inferred from tracking data by approximate Bayesian computation. We argue that the motility mechanism is highly sensitive to experimental conditions. We predict a TFP retraction speed for the tracking data we study that is in a good agreement with experimental results obtained under very similar conditions. Furthermore, we examine whether estimates for biologically important parameters, whose direct experimental determination is challenging, can be inferred directly from tracking data. One example is the width of the distribution of TFP on the bacteria body. We predict that the TFP are broadly distributed over the bacteria pole in both walking and crawling motility types. Moreover, we identified specific configurations of TFP that lead to transitions between walking and crawling states.

Erratum: Spontaneous Domain Formation in Spherically Confined Elastic Filaments [Phys. Rev. Lett. 123, 047801 (2019)].

This corrects the article DOI: 10.1103/PhysRevLett.123.047801.

Electrostatic interactions between charge regulated spherical macroions.

We study the interaction between two charge regulating spherical macroions with dielectric interior and dissociable surface groups immersed in a monovalent electrolyte solution. The charge dissociation is modelled via the Frumkin-Fowler-Guggenheim isotherm, which allows for multiple adsorption equilibrium states. The interactions are derived from the solutions of the mean-field Poisson-Boltzmann type theory with charge regulation boundary conditions. For a range of conditions we find symmetry breaking transitions from symmetric to asymmetric charge distribution exhibiting annealed charge patchiness, which results in like-charge attraction even in a univalent electrolyte-thus fundamentally modifying the nature of electrostatic interactions in charge-stabilized colloidal suspensions.

The role of surface topography in the self-assembly of polymeric surfactants.

We propose a classical density functional theory model to study the self-assembly of polymeric surfactants on curved surfaces. We use this model to investigate the thermodynamics of phase separation of a binary mixture of size asymmetric miscible surfactants on cylindrical and spherical surfaces, and observe that phase separation driven by size alone is thermodynamically unfavorable on both cylindrical and spherical surfaces. We use the theory, supplemented by dissipative particle dynamics (DPD) simulations, to predict pattern formation on a non-uniform surface with regions of positive and negative curvature. Our results suggest potential ways to couple surface topography and polymeric surfactants to design surfaces coated with non-uniform patterns.

Numerical Test of the Onsager Relations in a Driven System.

The Onsager reciprocity relations were formulated in the context of irreversible thermodynamics, but they are based on assumptions that have a wider applicability. Here, we present simulations testing the Onsager relations between surface-coupled diffusive and bulk fluxes in a system prepared in a nonequilibrium steady state. The system consists of a mixture of two identical species maintained at different temperatures inside a channel. In order to tune the friction of the two species with the walls independently, while keeping the particle-wall interaction potentials the same, we allow the kinematics of particle-wall collisions to be different: "bounce-back" (B) or "specular" (S). In the BB case, diffusio-capillary transport can only take place if the two species have different temperatures. We find that the Onsager reciprocity relations are obeyed in the linear regime, even in the BB case where all fluxes are the result of perturbing the system from a nonequilibrium steady state in a way that does not satisfy time-reversal symmetry. Our Letter provides a direct, numerical illustration of the validity of the Onsager relations outside their original range of application, and suggests their relevance for transport in driven or active systems.

Small Obstacle in a Large Polar Flock.

We show that arbitrarily large polar flocks are susceptible to the presence of a single small obstacle. In a wide region of parameter space, the obstacle triggers counterpropagating dense bands leading to reversals of the flow. In very large systems, these bands interact, yielding a never-ending chaotic dynamics that constitutes a new disordered phase of the system. While most of these results were obtained using simulations of aligning self-propelled particles, we find similar phenomena at the continuous level, not when considering the basic Toner-Tu hydrodynamic theory, but in simulations of truncations of the relevant Boltzmann equation.

Computational design of probes to detect bacterial genomes by multivalent binding.

Rapid methods for diagnosis of bacterial infections are urgently needed to reduce inappropriate use of antibiotics, which contributes to antimicrobial resistance. In many rapid diagnostic methods, DNA oligonucleotide probes, attached to a surface, bind to specific nucleotide sequences in the DNA of a target pathogen. Typically, each probe binds to a single target sequence; i.e., target-probe binding is monovalent. Here we show using computer simulations that the detection sensitivity and specificity can be improved by designing probes that bind multivalently to the entire length of the pathogen genomic DNA, such that a given probe binds to multiple sites along the target DNA. Our results suggest that multivalent targeting of long pieces of genomic DNA can allow highly sensitive and selective binding of the target DNA, even if competing DNA in the sample also contains binding sites for the same probe sequences. Our results are robust to mild fragmentation of the bacterial genome. Our conclusions may also be relevant for DNA detection in other fields, such as disease diagnostics more broadly, environmental management, and food safety.

Effect of the interaction strength and anisotropy on the diffusio-phoresis of spherical colloids.

Gradients in temperature, concentration or electrostatic potential cannot exert forces on a bulk fluid; they can, however, exert forces on a fluid in a microscopic boundary layer surrounding a (nano)colloidal solute, resulting in so-called phoretic flow. Here we present a simulation study of phoretic flow around a spherical colloid held fixed in a concentration gradient. We show that the resulting flow velocity depends non-monotonically on the strength of the colloid-fluid interaction. The reason for this non-monotonic dependence is that solute particles are effectively trapped in a shell around the colloid and cannot contribute to diffusio-phoresis. We also observe that the flow depends sensitively on the anisotropy of solute-colloid interaction.

The Lennard-Jones potential: when (not) to use it.

The Lennard-Jones 12-6 potential (LJ) is arguably the most widely used pair potential in molecular simulations. In fact, it is so popular that the question is rarely asked whether it is fit for purpose. In this paper, we argue that, whilst the LJ potential was designed for noble gases such as argon, it is often used for systems where it is not expected to be particularly realistic. Under those circumstances, the disadvantages of the LJ potential become relevant: most important among these is that in simulations the LJ potential is always modified such that it has a finite range. More seriously, there is by now a whole family of different potentials that are all called Lennard-Jones 12-6, and that are all different - and that may have very different macroscopic properties. In this paper, we consider alternatives to the LJ 12-6 potential that could be employed under conditions where the LJ potential is only used as a typical short-ranged potential with attraction. We construct a class of potentials that are, in many respects LJ-like but that are by construction finite ranged, vanishing quadratically at the cut-off distance, and that are designed to be computationally cheap. Below, we present this potential and report numerical data for its thermodynamic and transport properties, for the most important cases (cut-off distance rc = 2σ ("LJ-like") and rc = 1.2σ (a typical "colloidal" potential)).

Spontaneous Domain Formation in Spherically Confined Elastic Filaments.

Although the free energy of a genome packing into a virus is dominated by DNA-DNA interactions, ordering of the DNA inside the capsid is elasticity driven, suggesting general solutions with DNA organized into spool-like domains. Using analytical calculations and computer simulations of a long elastic filament confined to a spherical container, we show that the ground state is not a single spool as assumed hitherto, but an ordering mosaic of multiple homogeneously ordered domains. At low densities, we observe concentric spools, while at higher densities, other morphologies emerge, which resemble topological links. We discuss our results in the context of metallic wires, viral DNA, and flexible polymers.

Controlling the morphological evolution of a particle-stabilized binary-component system.

This work reveals the association between the morphological evolutions in particle-stabilized binary-component systems (i.e., binary liquids and polymer blends). It investigates the formation mechanism of bijels created via direct mixing, and proposes an empirical cost function to quantitatively evaluate the created structures.

Particle-stabilized Janus emulsions that exhibit pH-tunable stability.

By developing a deeper understanding of the formation mechanism and the origin of the stability, we report a simple and large-scale fabrication approach to create Janus emulsions that can be controlled in size, geometry and stability.

Active microrheology in two-dimensional magnetic networks.

We study active microrheology in two-dimensional (2D) magnetic networks. To this end, we use Langevin dynamics computer simulations where single non-magnetic or magnetic tracer particles are pulled through the network structures via a constant force f. Structural changes in the network around the pulled tracer particle are characterized in terms of pair correlation functions. These functions indicate that the non-magnetic tracer particles tend to strongly affect the network structure leading to the formation of channels at sufficiently high forces, while the magnetic tracer particles modify the network structure only slightly. At zero pulling force, f = 0, both non-magnetic and magnetic tracer particles are localized, i.e. they do not show diffusive behavior in the long-time limit. Nevertheless, the friction coefficient, as obtained from the steady-state velocity of the tracer particles, seems to indicate a linear-response regime at small values of f. Beyond the latter linear response regime, the diffusion dynamics of the tracer particles are anisotropic with superdiffusive behavior in force direction. This transport anomaly is investigated via van Hove correlation functions and residence time distributions.

Membrane potential drives direct translocation of cell-penetrating peptides.

Cell-penetrating peptides (CPPs) are frequently employed as drug delivery agents with rapid cellular uptake, however, the uptake mechanism and the detailed translocation pathway are at present not completely understood. Both endocytosis and direct translocation through membrane pores have been observed in experiments and simulations under different conditions. Here we report the molecular dynamics simulations providing evidence for the direct translocation of CPPs across the membrane driven by the membrane electrostatic potential. The local membrane potential can be produced by the ion concentration imbalance across the membrane, which is ubiquitous in biological environments. Moreover, if positively charged CPPs are adsorbed on the membrane, this further enhances the membrane potential, opening membrane pores through which CPPs can be instantly transported in a chain-like configuration. The classical nucleation theory is applied to estimate the translocation time by calculating the changes in the free energy upon transferring CPPs across the membrane at different potentials, showing good agreement with available experimental measurements. The revealed CPP translocation mechanism can be broadly relevant for cellular processes in biology.

Effect of Topographical Steps on the Surface Motility of the Bacterium Pseudomonas aeruginosa.

Bacteria traverse surfaces as part of colonizing solids, and it is of interest to hinder this motion to potentially thwart infections in humans. Here, we demonstrate that topographical steps hinder the ability of Pseudomonas aeruginosa PAO1 (P. aeruginosa) to traverse a solid-liquid interface. Using time-lapse fluorescence microscopy and image analysis, we analyzed the motion of P. aeruginosa that were challenged with steps ranging in height from 0.4 to 9.0 µm. Bacterial trajectories are sensitive to the height of the step, the curvature of the step face, and the direction of their motion relative to gravity. When the step height is ≥0.9 µm, which is similar to the cell diameter, there is a reduced probability of the cell crossing the step. For those bacteria that do cross a step, there is a time penalty for crossing steps of height 2-3 µm; this height is similar to the length of the bacterium. For higher steps, the bacteria reorient their cell body while traversing the step riser. Our findings elucidate how topography influences the motion of bacteria and informs strategies for hindering bacterial motion at surfaces.

Controlling Cargo Trafficking in Multicomponent Membranes.

Biological membranes typically contain a large number of different components dispersed in small concentrations in the main membrane phase, including proteins, sugars, and lipids of varying geometrical properties. Most of these components do not bind the cargo. Here, we show that such "inert" components can be crucial for the precise control of cross-membrane trafficking. Using a statistical mechanics model and molecular dynamics simulations, we demonstrate that the presence of inert membrane components of small isotropic curvatures dramatically influences cargo endocytosis, even if the total spontaneous curvature of such a membrane remains unchanged. Curved lipids, such as cholesterol, as well as asymmetrically included proteins and tethered sugars can, therefore, actively participate in the control of the membrane trafficking of nanoscopic cargo. We find that even a low-level expression of curved inert membrane components can determine the membrane selectivity toward the cargo size and can be used to selectively target membranes of certain compositions. Our results suggest a robust and general method of controlling cargo trafficking by adjusting the membrane composition without needing to alter the concentration of receptors or the average membrane curvature. This study indicates that cells can prepare for any trafficking event by incorporating curved inert components in either of the membrane leaflets.

What experiments on pinned nanobubbles can tell about the critical nucleus for bubble nucleation.

The process of homogeneous bubble nucleation is almost impossible to probe experimentally, except near the critical point or for liquids under large negative tension. Elsewhere in the phase diagram, the bubble nucleation barrier is so high as to be effectively insurmountable. Consequently, there is a severe lack of experimental studies of homogenous bubble nucleation under conditions of practical importance (e.g., cavitation). Here we use a simple geometric relation to show that we can obtain information about the homogeneous nucleation process from Molecular Dynamics studies of bubble formation in solvophobic nanopores on a solid surface. The free energy of pinned nanobubbles has two extrema as a function of volume: one state corresponds to a free-energy maximum ("the critical nucleus"), the other corresponds to a free-energy minimum (the metastable, pinned nanobubble). Provided that the surface tension does not depend on nanobubble curvature, the radius of the curvature of the metastable surface nanobubble is independent of the radius of the pore and is equal to the radius of the critical nucleus in homogenous bubble nucleation. This observation opens the way to probe the parameters that determine homogeneous bubble nucleation under experimentally accessible conditions, e.g. with AFM studies of metastable nanobubbles. Our theoretical analysis also indicates that a surface with pores of different sizes can be used to determine the curvature corrections to the surface tension. Our conclusions are not limited to bubble nucleation but suggest that a similar approach could be used to probe the structure of critical nuclei in crystal nucleation.

Crystallinity of Double-Stranded RNA-Antimicrobial Peptide Complexes Modulates Toll-Like Receptor 3-Mediated Inflammation.

Double-stranded RNA (dsRNA) induces production of pro-inflammatory cytokines in normal human epidermal keratinocytes (NHEK) by specific binding to endosomal Toll-like receptor-3 (TLR3). Recently, it has been shown that hyperactivation of TLR3 in psoriatic keratinocytes by dsRNA can occur in the presence of human antimicrobial peptide (AMP) LL37. Here, we combine synchrotron X-ray scattering, microscopy, computer simulations, and measurements of NHEK cytokine production to elucidate a previously unanticipated form of specific molecular pattern recognition. LL37 and similar α-helical AMPs can form pro-inflammatory nanocrystalline complexes with dsRNA that are recognized by TLR3 differently than dsRNA alone. dsRNA complexes that activate IL-6 production in NHEK and those that do not are both able to enter cells and co-localize with TLR3. However, the crystallinity of these AMP-dsRNA complexes, specifically the geometric spacing between parallel dsRNA and the repeat number of ordered dsRNA, strongly influences the level of TLR3 activation. Crystalline complexes that present dsRNA at a spacing that matches with the steric size of TLR3 can recruit and engage multiple TLR3 receptors, driving receptor clustering and immune amplification, whereas crystalline complexes that exhibit poor steric matching do not. Reverse-transcription quantitative PCR of IL-6 during siRNA knockdown of TLR3 confirms that cytokine production is due to TLR3: High levels of IL-6 transcription are observed for sterically matched complexes without TLR3 knockdown, whereas such activity is abrogated with TLR3 knockdown.

Optimal multivalent targeting of membranes with many distinct receptors.

Cells can often be recognized by the concentrations of receptors expressed on their surface. For better (targeted drug treatment) or worse (targeted infection by pathogens), it is clearly important to be able to target cells selectively. A good targeting strategy would result in strong binding to cells with the desired receptor profile and barely binding to other cells. Using a simple model, we formulate optimal design rules for multivalent particles that allow them to distinguish target cells based on their receptor profile. We find the following: (i) It is not a good idea to aim for very strong binding between the individual ligands on the guest (delivery vehicle) and the receptors on the host (cell). Rather, one should exploit multivalency: High sensitivity to the receptor density on the host can be achieved by coating the guest with many ligands that bind only weakly to the receptors on the cell surface. (ii) The concentration profile of the ligands on the guest should closely match the composition of the cognate membrane receptors on the target surface. And (iii) irrespective of all details, the effective strength of the ligand-receptor interaction should be of the order of the thermal energy [Formula: see text], where [Formula: see text] is the absolute temperature and [Formula: see text] is Boltzmann's constant. We present simulations that support the theoretical predictions. We speculate that, using the above design rules, it should be possible to achieve targeted drug delivery with a greatly reduced incidence of side effects.