Using Liquid Crystals to Reveal How Mechanical Anisotropy Changes Interfacial Behaviors of Motile Bacteria

Bacteria often inhabit and exhibit distinct dynamical behaviors at interfaces, but the physical mechanisms by which interfaces cue bacteria are still poorly understood. In this work, we use interfaces formed between coexisting isotropic and liquid crystal (LC) phases to provide insight into how mechanical anisotropy and defects in LC ordering influence fundamental bacterial behaviors. Specifically, we measure the anisotropic elasticity of the LC to change fundamental behaviors of motile, rod-shaped Proteus mirabilis cells (3 μm in length) adsorbed to the LC interface, including the orientation, speed, and direction of motion of the cells (the cells follow the director of the LC at the interface), transient multicellular self-association, and dynamical escape from the interface. In this latter context, we measure motile bacteria to escape from the interfaces preferentially into the isotropic phase, consistent with the predicted effects of an elastic penalty associated with strain of the LC about the bacteria when escape occurs into the nematic phase. We also observe boojums (surface topological defects) present at the interfaces of droplets of nematic LC (tactoids) to play a central role in mediating the escape of motile bacteria from the LC interface. Whereas the bacteria escape the interface of nematic droplets via a mechanism that involved nematic director-guided motion through one of the two boojums, for isotropic droplets in a continuous nematic phase, the elasticity of the LC generally prevented single bacteria from escaping. Instead, assemblies of bacteria piled up at boojums and escape occurred through a cooperative, multicellular phenomenon. Overall, our studies show that the dynamical behaviors of motile bacteria at anisotropic LC interfaces can be understood within a conceptual framework that reflects the interplay of LC elasticity, surface-induced order, and topological defects.     

GPU-Accelerated Molecular Dynamics Simulation to Study Liquid Crystal Phase Transition Using Coarse-Grained Gay-Berne Anisotropic Potential

Gay-Berne (GB) potential is regarded as an accurate model in the simulation of anisotropic particles, especially for liquid crystal (LC) mesogens. However, its computational complexity leads to an extremely time-consuming process for large systems. Here, we developed a GPU-accelerated molecular dynamics (MD) simulation with coarse-grained GB potential implemented in GALAMOST package to investigate the LC phase transitions for mesogens in small molecules, main-chain or side-chain polymers. For identical mesogens in three different molecules, on cooling from fully isotropic melts, the small molecules form a single-domain smectic-B phase, while the main-chain LC polymers prefer a single-domain nematic phase as a result of connective restraints in neighboring mesogens. The phase transition of side-chain LC polymers undergoes a two-step process: nucleation of nematic islands and formation of multi-domain nematic texture. The particular behavior originates in the fact that the rotational orientation of the mesogenes is hindered by the polymer backbones. Both the global distribution and the local orientation of mesogens are critical for the phase transition of anisotropic particles. Furthermore, compared with the MD simulation in LAMMPS, our GPU-accelerated code is about 4 times faster than the GPU version of LAMMPS and at least 200 times faster than the CPU version of LAMMPS. This study clearly shows that GPU-accelerated MD simulation with GB potential in GALAMOST can efficiently handle systems with anisotropic particles and interactions, and accurately explore phase differences originated from molecular structures.     

The interplay between viscoelastic and thermodynamic properties determines the birefringence of F-actin gels

F-actin gels of increasing concentrations (25-300 microM) display in vitro a progressive onset of birefringence due to orientational ordering of actin filaments. At F-actin concentrations <100 microM, this birefringence can be erased and restored at will by sonication and gentle flow, respectively. Hence, the orientational ordering does not result from a thermodynamic transition to a nematic phase but instead is due to mechanical stresses stored in the gels. In contrast, at F-actin concentrations > or =100 microM, gels display spontaneous birefringence recovery, at rest, which is the sign of true nematic ordering, in good agreement with statistical physics models of the isotropic/nematic transition. Well-aligned samples of F-actin gels could be produced and their small-angle x-ray scattering patterns are quite anisotropic. These patterns show no sign of filament positional short-range order and could be modeled by averaging the form factor with the Maier-Saupe nematic distribution function. The derived nematic order parameter S of the gels ranged from S = 0.7 at 300 microM to S = 0.4 at 25 microM. Both birefringence and small-angle x-ray scattering data indicate that, even in absence of cross-linking proteins, spontaneous cooperative alignment of actin filaments may arise in motile regions of living cells where F-actin concentrations can reach values of a few 100 microM.     

On the theory of the nematic phase and its possible relation to the mitotic spindle structure

Onsager's method of studying the nematic phase is developed for general molecular interactions. It is shown that the symmetry of the molecule helps determine the type of transition that occurs in passing from the isotropic phase to the anisotropic phase. The possible relation between the nematic phase and spindle structure is briefly discussed.     

Effect of Interfaces on Small, Starved Marine Bacteria

The copiotrophic marine Vibrio sp. strain DW1, shown previously in batch culture to increase in numbers at the onset of starvation and then to form viable small cells with low endogenous respiration, appears to have a survival advantage at interfaces. Vibrio sp. strain DW1 behaved differently at interfaces compared with the aqueous phase under starvation conditions: (i) small cells were observed at an air-water interface without nutrients, (ii) nutrients added to the air-water interface quickly produced larger cells at the surface, (iii) motility persisted many hours longer at the solid-water interface of a dialysis membrane in a microchamber at the onset of starvation, and (iv) regrowth and division at the solid-liquid interface occurred quickly and at nutrient concentrations too low to permit growth in the aqueous phase. It was concluded that, if small starved cells from copiotrophic bacteria can reach an interface, additional survival mechanisms become available to them: (i) interfaces constitute areas of favorable nutrient conditions, and (ii) interfaces lacking a sufficient amount of nutrient, nevertheless, trigger cells to become smaller, thus increasing their surface/volume ratio and the packing density.     

Imaging of anisotropic cellulose suspensions using environmental scanning electron microscopy

The effect of concentration on anisotropic phase behavior of acid-hydrolyzed cellulose suspensions has been examined using conventional polarizing microscopy and the novel technique of environmental scanning electron microscopy (ESEM). Microcrystalline cellulose dispersed in water formed biphasic suspensions in a narrow concentration range, 4-12 wt % for a suspension pH of 4, where the upper and lower phases were isotropic and anisotropic (chiral nematic), respectively. It is known from previous work that within the biphasic regime total suspension concentration affects only the volume fractions of the two phases, not phase concentration or interfacial packing. As the total suspension concentration surpassed the upper critical limit (c), however, a single anisotropic phase of increasing concentration was observed. It was evident from polarizing microscopy that the chiral nematic pitch of the anisotropic phase decreased with increasing concentration, which has been attributed to a reduction in the electrostatic double layer thickness of the individual rods, thus increasing intermolecular interactions. Chiral nematic textures were also visible using ESEM. This technique has the advantage of studying individual rod orientation within the liquid crystalline phase as it permits the high resolution of electron microscopy to be applied to hydrated samples in their natural state. To our knowledge this is the first time such lyotropic systems have been observed using electron microscopy.     

Cholesteryl ester hydrolysis in J774 macrophages occurs in the cytoplasm and lysosomes

The relationship of cholesteryl ester hydrolysis to the physical state of the cholesteryl ester in J774 murine macrophages was explored in cells induced to store cholesteryl esters either in anisotropic (ordered) inclusions or isotropic (liquid) inclusions. In contrast to other cell systems, the rate of cholesteryl ester hydrolysis was faster in cells containing anisotropic inclusions than in cells containing isotropic inclusions. Two contributing factors were identified. Kinetic analyses of the rates of hydrolysis are consistent with a substrate competition by co-deposited triglyceride in cells with isotropic inclusions. In addition, hydrolysis of cholesteryl esters in cells with anisotropic droplets is mediated by both cytoplasmic and lysosomal lipolytic enzymes, as shown by using the lysosomotropic agent, chloroquine, and an inhibitor of neutral cholesteryl ester hydrolase, umbelliferyl diethylphosphate. In cells containing anisotropic inclusions, hydrolysis was partially inhibited by incubation in media containing either chloroquine or umbelliferyl diethylphosphate. Together, chloroquine and umbelliferyl diethylphosphate completely inhibited hydrolysis. However, when cells containing isotropic inclusions were incubated with umbelliferyl diethylphosphate, cholesteryl ester hydrolysis was completely inhibited, but chloroquine had no effect. Transmission electron microscopy demonstrated a primarily lysosomal location for lipid droplets in cells with anisotropic droplets and both non-lysosomal and lysosomal populations of lipid droplets in cells with isotropic droplets.These results support the conclusion that there is a lysosomal component to the hydrolysis of stored cholesteryl esters in foam cells.     

In vitro chiral nematic ordering of chitin crystallites

Microfibrillar fragments of purified crab and shrimp chitin were prepared by hydrolysis in 3 HCl at its boiling point (104°C). After removal of the acid by centrifugal washing and dialysis, an ultrasound treatment converts the residual product to a colloidal suspension stabilized by NH₃⁺ charges. When dewatered to a critical concentration, spontaneous formation of a two-phase equilibrium system occurs. The upper phase (lower concentration) is isotropic and the lower phase is anisotropic. The latter displays chiral nematic order and dries to a solid film which mimics the helicoid organization characteristic of the chitin microfibrils in the cuticle of arthropods.     

Bacterial Motility Reveals Unknown Molecular Organization

The water solubility of lyotropic liquid crystals (LCs) makes them very attractive to study the behavior of biological microorganisms in an environment where local symmetry is broken (as often encountered in nature). Several recent studies have shown a dramatic change in the behavior of flagellated bacteria when swimming in solutions of the lyotropic LC disodium cromoglycate (DSCG). In this study, the movements of Escherichia coli bacteria in DSCG-water solutions of different concentrations are observed to improve our understanding of this phenomenon. In addition, the viscosity of DSCG aqueous solutions is measured as a function of concentration at room temperature. We also experimentally identify a previously undescribed isotropic pretransition zone where bacteria start sticking to each other and to surfaces. Simple estimations show that the unbalanced osmotic pressure induced depletion force might be responsible for this sticking phenomenon. An estimate of the bacteria propulsive force and the DSCG aggregates length (versus concentration) are calculated from the measured viscosity of the medium. All these quantities are found to undergo a strong increase in the pretransition zone, starting at a threshold concentration of 6 ± 1 wt % DSCG that is well below the known isotropic-LC transition (∼10 wt %). This study also shines light on the motility of flagellated bacteria in realistic environments, and it opens new avenues for interesting applications such as the use of motile microorganisms to probe the physical properties of their host or smart bandages that could guide bacteria out of wounds.     

Nematic Ordering Pattern Formation in the Process of Self-Organization of Microtubules in a Gravitational Field

Papaseit et al. (Proc. Natl. Acad. Sci. U.S.A. 97, 8364, 2000) showed the decisive role of gravity in the formation of patterns by assemblies of microtubules in vitro. By virtue of a functional scaling, the free energy for MT systems in a gravitational field was constructed. The influence of the gravitational field on MT’s self-organization process, that can lead to the isotropic to nematic phase transition, is the focus of this paper. A coupling of a concentration gradient with orientational order characteristic of nematic ordering pattern formation is the new feature emerging in the presence of gravity. The concentration range corresponding to a phase coexistence region increases with increasing g or MT concentration. Gravity facilitates the isotropic to nematic phase transition leading to a significantly broader transition region. The phase transition represents the interplay between the growth in the isotropic phase and the precipitation into the nematic phase. We also present and discuss the numerical results obtained for local MT concentration change with the height of the vessel, order parameter and phase transition properties. An erratum to this article can be found at     

Adaptive generation of multimaterial grids from imaging data for biomedical Lagrangian fluid–structure simulations

Spatial discretization of complex imaging- derived fluid–solid geometries, such as the cardiac environment, is a critical but often overlooked challenge in biomechanical computations. This is particularly true in problems with Lagrangian interfaces, where the fluid and solid phases share a common interface geometrically. For simplicity and better accuracy, it is also highly desirable for the two phases to have a matching surface mesh at the interface between them. We outline a method for solving this problem, and illustrate the approach with a 3D fluid–solid mesh of the mouse heart. An MRI dataset of a perfusion-fixed mouse heart with 50μm isotropic resolution was semi-automatically segmented using a customized multimaterial connected-threshold approach that divided the volume into non-overlapping regions of blood, tissue, and background. Subsequently a multimaterial marching cubes algorithm was applied to the segmented data to produce two detailed, compatible isosurfaces, one for blood and one for tissue. Both isosurfaces were simultaneously smoothed with a multimaterial smoothing algorithm that exactly conserves the volume for each phase. Using these two isosurfaces, we developed and applied novel automated meshing algorithms to generate anisotropic hybrid meshes on arbitrary biological geometries with the number of layers and the desired element anisotropy for each phase as the only input parameters. Since our meshes adapt to the local feature sizes and include boundary layer prisms, they are more efficient and accurate than non-adaptive, isotropic meshes, and the fluid–structure interaction computations will tend to have relative error equilibrated over the whole mesh.     

Enhanced Invasion of Metastatic Cancer Cells via Extracellular Matrix Interface

Cancer cell invasion is a major component of metastasis and is responsible for extensive cell diffusion into and major destruction of tissues. Cells exhibit complex invasion modes, including a variety of collective behaviors. This phenomenon results in the structural heterogeneity of the extracellular matrix (ECM) in tissues. Here, we systematically investigated the environmental heterogeneity facilitating tumor cell invasion via a combination of in vitro cell migration experiments and computer simulations. Specifically, we constructed an ECM microenvironment in a microfabricated biochip and successfully created a three-dimensional (3D) funnel-like matrigel interface inside. Scanning electron microscopy demonstrated that the interface was at the interior defects of the nano-scale molecular anisotropic orientation and the localized structural density variations in the matrigel. Our results, particularly the correlation of the collective migration pattern with the geometric features of the funnel-like interface, indicate that this heterogeneous in vitro ECM structure strongly guides and promotes aggressive cell invasion in the rigid matrigel space. A cellular automaton model was proposed based on our experimental observations, and the associated quantitative analysis indicated that cell invasion was initiated and controlled by several mechanisms, including microenvironment heterogeneity, long-range cell-cell homotype and gradient-driven directional cellular migration. Our work shows the feasibility of constructing a complex and heterogeneous in vitro 3D ECM microenvironment that mimics the in vivo environment. Moreover, our results indicate that ECM heterogeneity is essential in controlling collective cell invasive behaviors and therefore determining metastasis efficiency.     

Enantioselective separations using chiral supported liquid crystalline membranes

Porous and nonporous supported liquid crystalline membranes were produced by impregnating porous cellulose nitrate supports with cholesteric liquid crystal (LC) materials consisting of 4-cyano-4'-pentylbiphenyl (5CB) mixed with a cholesterol-based dopant (cholesteryl oleyl carbonate [COC], cholesteryl nonanoate [CN], or cholesteryl chloride [CC]). The membranes exhibit selectivity for R-phenylglycine and R-1-phenylethanol because of increased interactions between the S enantiomers and the left-handed cholesteric phase. The selectivity of both phenylglycine and 1-phenylethanol in 5CB/CN membranes decreases with effective pore diameter while the permeabilities increase, as expected. Phenylglycine, which is insoluble in the LC phase, exhibits no transport in the nonporous (completely filled) membranes; however, 1-phenylethanol, which is soluble in the LC phase, exhibits transport but negligible enantioselectivity. The enantioselectivity for 1-phenylethanol was higher (1.20 in 5CB/COC and 5CB/CN membranes) and the permeability was lower in the cholesteric phase than in the isotropic phase. Enantioselectivity was also higher in the 5CB/COC cholesteric phase than in the nematic phase of undoped 5CB (1.03). Enantioselectivity in the cholesteric phase of 5CB doped with CC (1.1), a dopant lacking hydrogen bonding groups, was lower than in the 5CB/COC phases. Finally, enantioselectivity increases with the dopant concentration up to a plateau value at approximately 17 mol%.     

Investigations of Nematic-Isotropic Transition by Means of Constant Pressure Molecular Dynamics Simulations

Abstract

Constant pressure molecular dynamics simulations, which secure the system to be under hydrostatic pressure, are used to simulate the behavior of liquid crystals consisting of anisotropic molecules with both translational and orientational freedom. In order to investigate to what extent can the properties known to real liquid crystalline phases be explained by the anisotropy of the shape of the molecules alone, the molecular dynamic (MD) simulation uses purely repulsive short-range pair potentials representing soft spherocylinders. A clear change in the microscopic as well as the macroscopic physical properties are observed near the phase transition from nematic liquid crystal to isotropic liquid.     

Oscillatory Mechanisms in Pairs of Neurons Connected with Fast Inhibitory Synapses

We study dynamical mechanisms underlying oscillatory behavior in reciprocal inhibitory pairs of neurons, using a two-dimensionalcell model. We introduce one-and-two dimensional phase portraits to illustratethe behaviors, thus reducing the study of dynamical mechanisms to planar geometrical properties. We examined whether other mechanisms besides the escape and release mechanisms (Wang and Rinzel, 1992) might be needed for some cases of reciprocal inhibition, and show that, within the confines of a simple two-dimensional cell model, escape and releaseare sufficient for all cases. We divided the behaviors of a singlecell into six different types and examined the joint behaviors arising from every combination of pairs of cells with behaviors drawn from thesesix types. For the case of two quiescent cells or two cells eachhaving plateau potentials, bifurcation diagrams demonstrate therelations between synaptic threshold and synaptic strength necessaryfor oscillations by escape, oscillations by release, ornetwork-generated plateau potentials. Thus we clarify therelationship between plateau potentials and oscillations in a cell.Using the two dimensional cell model we examine 1:N beating betweencells and find that our simple model displays many of the essentialdynamical properties displayed by more sophisticated models, some ofwhich relate to thalamocortical spindling.     

Modeling cell apoptosis for simulating three-dimensional multicellular morphogenesis based on a reversible network reconnection framework

Morphogenesis in multicellular organisms is accompanied by apoptotic cell behaviors: cell shrinkage and cell disappearance. The mechanical effects of these behaviors are spatiotemporally regulated within multicellular dynamics to achieve proper tissue sizes and shapes in three-dimensional (3D) space. To analyze 3D multicellular dynamics, 3D vertex models have been suggested, in which a reversible network reconnection (RNR) model has successfully expressed 3D cell rearrangements during large deformations. To analyze the effects of apoptotic cell behaviors on 3D multicellular morphogenesis, we modeled cell apoptosis based on the RNR model framework. Cell shrinkage was modeled by the potential energy as a function of individual cell times during the apoptotic phase. Cell disappearance was modeled by merging neighboring polyhedrons at their boundary surface according to the topological rules of the RNR model. To establish that the apoptotic cell behaviors could be expressed as modeled, we simulated morphogenesis driven by cell apoptosis in two types of tissue topology: 3D monolayer cell sheet and 3D compacted cell aggregate. In both types of tissue topology, the numerical simulations successfully illustrated that cell aggregates gradually shrank because of successive cell apoptosis. During tissue shrinkage, the number of cells in aggregates decreased while maintaining individual cell size and shape. Moreover, in case of localizing apoptotic cells within a part of the 3D monolayer cell aggregate, the cell apoptosis caused the global tissue bending by pulling on surrounding cells. In case of localizing apoptotic cells on the surface of the 3D compacted cell aggregate, the cell apoptosis caused successive, directional cell rearrangements from the inside to the surface. Thus, the proposed model successfully provided a basis for expressing apoptotic cell behaviors during 3D multicellular morphogenesis based on an RNR model framework.     

Tensorial elastic network model for protein dynamics: Integration of the anisotropic network model with bond‐bending and twist elasticities

We present a tensorial elastic network model (TNM) to describe the equilibrium fluctuations of proteins near their native fold structure. The model combines the anisotropic network model (ANM), bond bending elasticity, and backbone twist elasticity, and can predict both the isotropic fluctuations, similar to the Gaussian network model (GNM), and anisotropic fluctuations, similar to the ANM. TNM performs equally well for B‐factor predictions as GNM and predicts the anisotropy of B‐factors better than ANM. The model also outperforms the ANM in its predictability of the complete anisotropic displacement parameters. Proteins 2012; © 2012 Wiley Periodicals, Inc.     

Liquid crystal-based detection of thrombin coupled to interactions between a polyelectrolyte and a phospholipid monolayer

Herein, we describe a real-time, label-free biosensing strategy for thrombin detection that uses the orientational properties of nematic liquid crystals (LCs) and the interactions between a polyelectrolyte and a phospholipid monolayer. The imaging principle is based on the disruption of the orientation of 4-cyano-4′-pentylbiphenyl by reorganized 1,2-dioleoyl-sn-glycero-3-phospho-rac-(1-glycerol) sodium salt (DOPG) at the aqueous/LC interface. Positively charged, multiple arginine peptides (poly-l-arginine hydrochloride) interacted with negatively charged DOPG at the aqueous/LC interface, which caused reorganization of the phospholipid layer and induced an orientational transition of LCs from a homeotropic to a planar state. As a result, a dark to bright shift in the optical response was observed. Thrombin cleaves poly-l-arginine hydrochloride into peptides. Thus, when thrombin was added, the optical signals generated by the LCs reverted from bright to dark because of the weakened ability of the fragments to induce electrostatic interactions. The limit of detection of the LC-based sensor was 0.25 ng/mL (6.7 pM) thrombin, and the sensor was fully reusable. The detection limit of our LC-based interface sensor is 600 times lower than that of a previously reported enzyme-linked aptamer assay for the detection of thrombin. Thus, we have established a new, simple thrombin biosensor with high sensitivity and low interference.     

Second-order susceptibilities of anisotropic liquid-crystal thin films

Using optical second-harmonic generation, we have measured the anisotropic orientation behavior of liquid crystal (LC) molecules in thin films. A high degree of in-plane anisotropy was obtained by depositing the LC films onto rubbed polyimide-covered substrates. The analysis with respect to the relevant second-order susceptibility tensor components is based upon the experimental results obtained in transmission geometry using different angles of incidence for the separation of isotropic and anisotropic contributions.     

F-actin, a model polymer for semiflexible chains in dilute, semidilute, and liquid crystalline solutions.

Single actin filaments were analyzed in solutions ranging from dilute (0.2 microgram/ml), where filaments interact only with solvent, to concentrations (4.0 mg/ml) at which F-actin forms a nematic phase. A persistence length of approximately 1.8 microns and an average length of approximately 22 microns (Kaufmann et al., 1992) identify actin as a model for studying the dynamics of semiflexible polymers. In dilute solutions the filaments exhibit thermal bending undulations in addition to diffusive motion. At higher semidilute concentrations (1.4 mg/ml) three-dimensional reconstructions of confocal images of fluorescently labeled filaments in a matrix of unlabeled F-actin reveal steric interactions between filaments, which account for the viscoelastic behavior of these solutions. The restricted undulations of these labeled chains reveal the virtual tube formed around a filament by the surrounding actin. The average tube diameter <a> scales with monomer concentration c as <a> varies; is directly proportional to c-(0.5 +/- 0.15). The diffusion of filaments in semidilute solutions (c = (0.1-2.0) mg/ml) is dominated by diffusion along the filament contour (reptation), and constraint release by remodeling of the surrounding filaments is rare. The self-diffusion coefficient D parallel along the tube decreases linearly with the chain length for semidilute solutions. For concentrations > 2.5 mg/ml a transition occurs from an isotropic entangled phase to a coexistence between isotropic and nematic domains. Analysis of the molecular motions of filaments suggests that the filaments in the aligned domains are in thermal equilibrium and that the diffusion coefficient parallel to the director D parallel is nearly independent of filament length. We also report the novel direct observation of u-shaped defects, called hairpins, in the nematic domains.