Models of collective cell motion for cell populations with different aspect ratio : diffusion, proliferation and travelling waves

Continuum, partial differential equation models are often used to describe the collective motion of cell populations, with various types of motility represented by the choice of diffusion coefficient, and cell proliferation captured by the source terms. Previously, the choice of diffusion coefficient has been largely arbitrary, with the decision to choose a particular linear or nonlinear form generally based on calibration arguments rather than making any physical connection with the underlying individual-level properties of the cell motility mechanism. In this work we provide a new link between individual-level models, which account for important cell properties such as varying cell shape and volume exclusion, and population-level partial differential equation models. We work in an exclusion process framework, considering aligned, elongated cells that may occupy more than one lattice site, in order to represent populations of agents with different sizes. Three different idealizations of the individual-level mechanism are proposed, and these are connected to three different partial differential equations, each with a different diffusion coefficient; one linear, one nonlinear and degenerate and one nonlinear and nondegenerate. We test the ability of these three models to predict the population level response of a cell spreading problem for both proliferative and nonproliferative cases. We also explore the potential of our models to predict long time travelling wave invasion rates and extend our results to two dimensional spreading and invasion. Our results show that each model can accurately predict density data for nonproliferative systems, but that only one does so for proliferative systems. Hence great care must be taken to predict density data for with varying cell shape.

Numerical simulation of red blood cells’ motion : a review

The micro-circulation of blood plays an important role in human body by providing oxygen and nutrients to the cells and removing carbon dioxide and wastes from the cells. This process is greatly affected by the rheological properties of the Red Blood Cells (RBCs). Changes in the rheological properties of the RBCs are caused by certain human diseases such as malaria and sickle cell diseases. Therefore it is important to understand the motion and deformation mechanism of RBCs in order to diagnose and treat this kind of diseases. Although, many methods have been developed to explore the behavior of the RBCs in micro-channels, they could not explain the deformation mechanism of the RBCs properly. Recently developed Particle Methods are employed to explain the RBCs’ behavior in micro-channels more comprehensively. The main objective of this study is to critically analyze the present methods, used to model the RBC behavior in micro-channels, in order to develop a computationally efficient particle based model to describe the complete behavior of the RBCs in micro-channels accurately and comprehensively

Simulating droplet motion on virtual leaf surfaces

A curvilinear thin film model is used to simulate the motion of droplets on a virtual leaf surface, with a view to better understand the retention of agricultural sprays on plants. The governing model, adapted from Roy et al. (2002 J. Fluid Mech. 454, 235–261) with the addition of a disjoining pressure term, describes the gravity- and curvature driven flow of a small droplet on a complex substrate: a cotton leaf reconstructed from digitized scan data. Coalescence is the key mechanism behind spray coating of foliage, and our simulations demonstrate that various experimentally observed coalescence behaviours can be reproduced qualitatively. By varying the contact angle over the domain, we also demonstrate that the presence of a chemical defect can act as an obstacle to the droplet’s path, causing break-up. In simulations on the virtual leaf, it is found that the movement of a typical spray size droplet is driven almost exclusively by substrate curvature gradients. It is not until droplet mass is sufficiently increased via coalescence that gravity becomes the dominating force.

Non-linear resonance in strings under narrow band random excitation, part I: Planar response and stability

Non-linear planar response of a string to planar narrow band random excitation is investigated in this paper. A response equation for the mean square deflection σ2 is obtained under a single mode approximation by using the equivalent linearization technique. It is shown that the response is triple valued, as in the case of harmonic excitation, if the centre frequency of excitation Ω lies in a certain specified range. The triple valued response occurs only if the excitation bandwidth β is smaller than a critical value βcrit which is a monotonically increasing function of the intensity of excitation. An approximate method of investigating the almost sure asymptotic stability of the solution is presented and regions of instability in the Ω-σ2 plane have been charted. It is shown that planar response can become unstable either due to an unbounded growth of the in-plane component of motion or due to a spontaneous appearance of an out-of-plane component.

Uniformity of film thickness on rotating planetary planar substrates

For the purposes of obtaining a number of components with nearly identical thickness distributions over the substrate area and of minimizing the inhomogeneities of the film, it is logical to presume that a substrate rotating on its own axis and revolving around another axis will give more uniformity in film thickness than a substrate only revolving around one axis. In relation to the practical applications, an investigation has been undertaken to study the refinement that can be achieved by using a planar planetary substrate holder. It is shown theoretically that the use of the planetary substrate holder under ideal conditions of source and geometry does not offer any further improvement in uniformity of thickness over the conventional rotary work-holder. It is also shown that the geometrical parameters alone have little influence over the uniformity achieved on a planetary substrate, because of the complex cyclidal motion of any point on it. However, for any given geometry, a non-integral speed ratio of the planetary substrate and the work-holder shows considerably less variation in thickness over the substrate area.

Mechanism of RNA Translocation by a Viral Packaging Motor

Molecular motors are proteins that convert chemical energy into mechanical work. The viral packaging ATPase P4 is a hexameric molecular motor that translocates RNA into preformed viral capsids. P4 belongs to the ubiquitous class of hexameric helicases. Although its structure is known, the mechanism of RNA translocation remains elusive. Here we present a detailed kinetic study of nucleotide binding, hydrolysis, and product release by P4. We propose a stochastic-sequential cooperative model to describe the coordination of ATP hydrolysis within the hexamer. In this model the apparent cooperativity is a result of hydrolysis stimulation by ATP and RNA binding to neighboring subunits rather than cooperative nucleotide binding. Simultaneous interaction of neighboring subunits with RNA makes the otherwise random hydrolysis sequential and processive. Further, we use hydrogen/deuterium exchange detected by high resolution mass spectrometry to visualize P4 conformational dynamics during the catalytic cycle. Concerted changes of exchange kinetics reveal a cooperative unit that dynamically links ATP binding sites and the central RNA binding channel. The cooperative unit is compatible with the structure-based model in which translocation is effected by conformational changes of a limited protein region. Deuterium labeling also discloses the transition state associated with RNA loading which proceeds via opening of the hexameric ring. Hydrogen/deuterium exchange is further used to delineate the interactions of the P4 hexamer with the viral procapsid. P4 associates with the procapsid via its C-terminal face. The interactions stabilize subunit interfaces within the hexamer. The conformation of the virus-bound hexamer is more stable than the hexamer in solution, which is prone to spontaneous ring openings. We propose that the stabilization within the viral capsid increases the packaging processivity and confers selectivity during RNA loading. Finally, we use single molecule techniques to characterize P4 translocation along RNA. While the P4 hexamer encloses RNA topologically within the central channel, it diffuses randomly along the RNA. In the presence of ATP, unidirectional net movement is discernible in addition to the stochastic motion. The diffusion is hindered by activation energy barriers that depend on the nucleotide binding state. The results suggest that P4 employs an electrostatic clutch instead of cycling through stable, discrete, RNA binding states during translocation. Conformational changes coupled to ATP hydrolysis modify the electrostatic potential inside the central channel, which in turn biases RNA motion in one direction. Implications of the P4 model for other hexameric molecular motors are discussed.

A real-time algorithm for simulation of flexible objects and hyper-redundant manipulators

Flexible objects such as a rope or snake move in a way such that their axial length remains almost constant. To simulate the motion of such an object, one strategy is to discretize the object into large number of small rigid links connected by joints. However, the resulting discretised system is highly redundant and the joint rotations for a desired Cartesian motion of any point on the object cannot be solved uniquely. In this paper, we revisit an algorithm, based on the classical tractrix curve, to resolve the redundancy in such hyper-redundant systems. For a desired motion of the `head' of a link, the `tail' is moved along a tractrix, and recursively all links of the discretised objects are moved along different tractrix curves. The algorithm is illustrated by simulations of a moving snake, tying of knots with a rope and a solution of the inverse kinematics of a planar hyper-redundant manipulator. The simulations show that the tractrix based algorithm leads to a more `natural' motion since the motion is distributed uniformly along the entire object with the displacements diminishing from the `head' to the `tail'.

Use of 2-pyrimidineamidoxime to generate polynuclear homo-/heterometallic assemblies: synthesis, crystal structures and magnetic study with theoretical investigations on the exchange mechanism

Three new transition metal complexes using 2-pyrimidineamidoxime (pmadH(2)) as multidentate chelating and/or bridging ligand have been synthesized and characterized. The ligand pmadH(2) has two potential bridging functional groups mu-O and mu-(N-O)] and consequently shows several coordination modes. While a polymeric 1D Cu-II complex Cu(pmadH(2))(2)(NO3)](NO3) (1) was obtained upon treatment of Cu(NO3)(2)center dot 3H(2)O with pmadH(2) at room temperature in the absence of base, a high temperature reaction in the presence of base yielded a tetranuclear Cu-II-complex Cu-4(pmad)(2)(pmadH)(2)(NO3)](NO3)(H2O) (2). One of the Cu-II centers is in a square pyramidal environment while the other three are in a square planar geometry. Reaction of the same ligand with an equimolar mixture of both Cu(NO3)(2)center dot 3H(2)O and NiCl2 center dot 6H(2)O yielded a tetranuclear heterometallic (Cu2Ni2II)-Ni-II complex Cu2Ni2(pmad)(2)(pmadH)(2)Cl-2]center dot H2O (3) containing both square planar (Ni-II) and square pyramidal (Cu-II) metal centers. Complexes 1-3 represent the first examples of polynuclear metal complexes of 2-pyrimidineamidoxime. The analysis of variable temperature magnetic susceptibility data of 2 reveals that both ferromagnetic and antiferromagnetic interactions exist in this complex (J(1) = +10.7 cm(-1) and J(2) = -2.7 cm(-1) with g = 2.1) leading to a resultant ferromagnetic behavior. Complex 3 shows expected antiferromagnetic interaction between two Cu-II centers through -N-O- bridging pathway with J(1) = -3.4 cm(-1) and g = 2.08. DFT calculations have been used to corroborate the magnetic results.

Self-consistent microscopic treatment of the effects of self-motion of the probe on ionic and dipolar solvation dynamics

A simple but self-consistent microscopic theory for the time dependent solvation energy of both ions and dipoles is presented which includes, for the first time, the details of the self-motion of the probe on its own solvation dynamics. The theory leads to several interesting predictions. The most important of them is that, for dipolar solvation, both the rotational and the translational motions of the dipolar solute probe can significantly accelerate the rate of solvation. In addition, the rotational self-motion of the solute can also give rise to an additional mechanism of nonexponentiality in solvation time correlation functions in otherwise slow liquids. A comparison between the present theoretical predictions and the recent experimental studies of Maroncelli et al. on solvation dynamics of aniline in l-propanol seems to indicate that the said experiments have missed the initial solvent response up to about 45 ps. After mapping the experimental results on the redefined time scale, the theoretical results can explain the experimental results for solvation of aniline in 1-propanol very well. For ionic solvation, the translational motion is significant for light solutes only. For example, for Li+ in water, translational motion speeds up the solvation by about 20%. The present theory demonstrates that in dipolar solvation the partial quenching of the self-motion due to the presence of specific solute-solvent interactions (such as H-bonding) may lead to a much slower solvation than that when the self-motion is present. This point has been discussed. In addition, we present the theoretical results for solvation of aniline in propylene carbonate, Here, the solvation is predicted to be complete within 15-20 ps.

Synthesis of workspaces of planar manipulators with arbitrary topology using shape representation and simulated annealing

This paper presents a general methodology for the synthesis of the external boundary of the workspaces of a planar manipulator with arbitrary topology. Both the desired workspace and the manipulator workspaces are identified by their boundaries and are treated as simple closed polygons. The paper introduces the concept of best match configuration and shows that the corresponding transformation can be obtained by using the concept of shape normalization available in image processing literature. Introduction of the concept of shape in workspace synthesis allows highly accurate synthesis with fewer numbers of design variables. This paper uses a new global property based vector representation for the shape of the workspaces which is computationally efficient because six out of the seven elements of this vector are obtained as a by-product of the shape normalization procedure. The synthesis of workspaces is formulated as an optimization problem where the distance between the shape vector of the desired workspace and that of the workspace of the manipulator at hand are minimized by changing the dimensional parameters of the manipulator. In view of the irregular nature of the error manifold, the statistical optimization procedure of simulated annealing has been used. A number of worked-out examples illustrate the generality and efficiency of the present method. (C) 1998 Elsevier Science Ltd. All rights reserved.

An Improved Compact Compliant Mechanism for an External Pipe-Crawler

An external pipe-crawling device presented in this paper aids the inspection of pipes in hazardous environments and areas inaccessible to humans. The principal component of our design, which uses inchworm type motion, is a compliant ring mechanism actuated using shape memory alloy (SMA) wire. It was fabricated and tested and was reported in our earlier work. But this device had a drawback of low crawling speed (about 1 mm/min) owing to the delay in heating and cooling of the SMA strips in the linear actuation. Additionally, that design also had the difficulties of mounting on pipes with closed ends, large radial span, and the need for housing for electrical insulation and guiding of the SMA wire. In this paper we present a compact design that overcomes the difficulties of the earlier design. In particular, we present a compact compliant mechanism with two halves so as to enable mounting and un-mounting on any closed or open pipe. Another feature is the presence of insulation and guiding of the SMA wire without housing. This design results in a reduction of the radial span of the ring from 22 mm to 12 mm, and the stiffness of the mechanism and the SMA wire are matched. An SMA helical spring is to used in the place of an SMA strip to increase the crawling speed of the device. A microcontroller-based circuitry is also fitted to cyclically.activate the SMA wires and springs.

Whole arm manipulation using closed and hybrid planar chains: A kinematic study

This paper studies planar whole arm manipulation of a circular object using closed loop and hybrid manipulators. The manipulation is simple with fewer degrees of actuation than the task space. This is an useful operation if the initial and final positions of the object are on the same surface. Closed loop manipulator is a 4/5 bar mechanism. In hybrid manipulators a open loop manipulator with 3/4 links is attached to the floating link of 4/5 bar mechanism. The mobility analysis is carried out to find the connectivity of the object with reference to frame. The manipulation (forward kinematics) starts from a given configuration of the object and the manipulator. In hybrid manipulators determination of initial configuration involves inverse kinematics of open loop manipulator. The input joint velocities are used to demonstrate the manipulation. Conditions are specified for prehensile manipulation.

Barrier crossing by a flexible long chain molecule - The kink mechanism

The problem and related earlier work All the above problems involve the passage of a long chain molecule, through a region in space, where the free energy per segment is higher, thus effectively presenting a barrier for the motion of the molecule. This is what we refer to as the Kramers proble...

Simulation of length-preserving motions of flexible one dimensional oObjects using optimization

During the motion of one dimensional flexible objects such as ropes, chains, etc., the assumption of constant length is realistic. Moreover,their motion appears to be naturally minimizing some abstract distance measure, wherein the disturbance at one end gradually dies down along the curve defining the object. This paper presents purely kinematic strategies for deriving length-preserving transformations of flexible objects that minimize appropriate ‘motion’. The strategies involve sequential and overall optimization of the motion derived using variational calculus. Numerical simulations are performed for the motion of a planar curve and results show stable converging behavior for single-step infinitesimal and finite perturbations 1 as well as multi-step perturbations. Additionally, our generalized approach provides different intuitive motions for various problem-specific measures of motion, one of which is shown to converge to the conventional tractrix-based solution. Simulation results for arbitrary shapes and excitations are also included.