Stop-and-go movements of plant Golgi stacks are mediated by the acto-myosin system

The Golgi apparatus in plant cells consists of a large number of independent Golgi stack/trans-Golgi network/Golgi matrix units that appear to be randomly distributed throughout the cytoplasm. To study the dynamic behavior of these Golgi units in living plant cells, we have cloned a cDNA from soybean (Glycine max), GmMan1, encoding the resident Golgi protein alpha-1,2 mannosidase I. The predicted protein of approximately 65 kD shows similarity of general structure and sequence (45% identity) to class I animal and fungal alpha-1,2 mannosidases. Expression of a GmMan1::green fluorescent protein fusion construct in tobacco (Nicotiana tabacum) Bright Yellow 2 suspension-cultured cells revealed the presence of several hundred to thousands of fluorescent spots. Immuno-electron microscopy demonstrates that these spots correspond to individual Golgi stacks and that the fusion protein is largely confined to the cis-side of the stacks. In living cells, the stacks carry out stop-and-go movements, oscillating rapidly between directed movement and random "wiggling." Directed movement (maximal velocity 4.2 microm/s) is related to cytoplasmic streaming, occurs along straight trajectories, and is dependent upon intact actin microfilaments and myosin motors, since treatment with cytochalasin D or butanedione monoxime blocks the streaming motion. In contrast, microtubule-disrupting drugs appear to have a small but reproducible stimulatory effect on streaming behavior. We present a model that postulates that the stop-and-go motion of Golgi-trans-Golgi network units is regulated by "stop signals" produced by endoplasmic reticulum export sites and locally expanding cell wall domains to optimize endoplasmic reticulum to Golgi and Golgi to cell wall trafficking.     

Pattern-dependent response modulations in motion-sensitive visual interneurons--a model study

Even if a stimulus pattern moves at a constant velocity across the receptive field of motion-sensitive neurons, such as lobula plate tangential cells (LPTCs) of flies, the response amplitude modulates over time. The amplitude of these response modulations is related to local pattern properties of the moving retinal image. On the one hand, pattern-dependent response modulations have previously been interpreted as 'pattern-noise', because they deteriorate the neuron's ability to provide unambiguous velocity information. On the other hand, these modulations might also provide the system with valuable information about the textural properties of the environment. We analyzed the influence of the size and shape of receptive fields by simulations of four versions of LPTC models consisting of arrays of elementary motion detectors of the correlation type (EMDs). These models have previously been suggested to account for many aspects of LPTC response properties. Pattern-dependent response modulations decrease with an increasing number of EMDs included in the receptive field of the LPTC models, since spatial changes within the visual field are smoothed out by the summation of spatially displaced EMD responses. This effect depends on the shape of the receptive field, being the more pronounced--for a given total size--the more elongated the receptive field is along the direction of motion. Large elongated receptive fields improve the quality of velocity signals. However, if motion signals need to be localized the velocity coding is only poor but the signal provides--potentially useful--local pattern information. These modelling results suggest that motion vision by correlation type movement detectors is subject to uncertainty: you cannot obtain both an unambiguous and a localized velocity signal from the output of a single cell. Hence, the size and shape of receptive fields of motion sensitive neurons should be matched to their potential computational task.     

Imaging of calcium wave propagation in guinea-pig ventricular cell pairs by confocal laser scanning microscopy.

We describe here the use of a confocal laser scanning microscope for imaging fast dynamic changes of the intracellular calcium ion concentration ([Ca2+]i) in isolated ventricular cell pairs. The scanning apparatus of our system, paired galvanometer mirrors, can perform narrow band scanning of an area of interest at a high temporal resolution of less than 70 msec per image. The actual [Ca2+]i is obtained directly through the fluorescence intensity of injected fluo-3, which responds to changes of [Ca2+]i in optically sectioned unit volumes of the cell. Images of the calcium wave obtained during propagation between paired cells revealed that the wavefront is constant in shape and propagates at constant velocity without any delay at the cell-to-cell junction. The confocal laser scanning microscope with depth-discriminating ability is a valuable tool for taking pictures of the sequence of biological events in living cells.     

Combining multiple laser scans of spotted microarrays by means of a two-way ANOVA model

MOTIVATION: Assessment of gene expression on spotted microarrays is based on measurement of fluorescence intensity emitted by hybridized spots. Unfortunately, quantifying fluorescence intensity from hybridized spots does not always correctly reflect gene expression level. Low gene expression levels produce low fluorescence intensities which tend to be confounded with the local background while high gene expression levels produce high fluorescence intensities which rapidly reach the saturation level. Most algorithms that combine data acquired at different voltages of the photomultiplier tube (PMT) assume that a change in scanner setting transforms the intensity measurements by a multiplicative constant. METHODS AND RESULTS: In this paper we introduce a new model of spot foreground intensity which integrates a PMT voltage independent scanner optical bias. This new model is used to implement a "Combining Multiple Scan using a Two-way ANOVA" (CMS2A) method, which is based on a maximum likelihood estimation of the scanner optical bias. After having computed scanner bias, coefficients of the two-way ANOVA model are used for correcting the saturated spots intensities obtained at high PMT voltage by using their counterpart values at lower PMT voltages. The method was compared to state-of-the-art multiple scan algorithms, using data generated from the MAQC study. CMS2A produced fold-changes that were highly correlated with qPCR fold-changes. As the scanner optical bias is accurately estimated within CMS2A, this method allows also avoiding fold-change compression biases whatever the value of this optical bias.     

Heterogeneous transportation of alpha1B-adrenoceptor in living cells

The heterogeneous motion of alpha(1B)-adrenoceptor (alpha(1B)-AR) was visualized in living cells with BODIPY-labeled antagonist of AR by single molecule fluorescence microscopy at high spatial resolution. The moving trajectory was reconstructed by precise localization (better than 20 nm) with a least-square fit of a two-dimensional Gaussian point spread function to each single spot. Trajectory analysis revealed two apparent groups of movements: directed motion and hindered motion. The directed motion had speeds higher than 0.1 mum/s. The histogram of diffusion coefficients of the hindered motion showed distinction between the cell membrane and the cytoplasm: the diffusion coefficient was lower near the cell membrane than in the internal cytoplasm, suggesting that alpha(1B)-AR was located or trapped in different networks, which was consistent with the natural distribution of cytoskeleton in living cells. These results suggested that the heterogeneity in the motion of alpha(1B)-AR in living cell might be associated with different localizations of cell skeleton proteins in the cell, which could provide molecular insight of AR regulation in living cells.     

CompuCell3D Simulations Reproduce Mesenchymal Cell Migration on Flat Substrates

Mesenchymal cell crawling is a critical process in normal development, in tissue function, and in many diseases. Quantitatively predictive numerical simulations of cell crawling thus have multiple scientific, medical, and technological applications. However, we still lack a low-computational-cost approach to simulate mesenchymal three-dimensional (3D) cell crawling. Here, we develop a computationally tractable 3D model (implemented as a simulation in the CompuCell3D simulation environment) of mesenchymal cells crawling on a two-dimensional substrate. The Fürth equation, the usual characterization of mean-squared displacement (MSD) curves for migrating cells, describes a motion in which, for increasing time intervals, cell movement transitions from a ballistic to a diffusive regime. Recent experiments have shown that for very short time intervals, cells exhibit an additional fast diffusive regime. Our simulations’ MSD curves reproduce the three experimentally observed temporal regimes, with fast diffusion for short time intervals, slow diffusion for long time intervals, and intermediate time -interval-ballistic motion. The resulting parameterization of the trajectories for both experiments and simulations allows the definition of time- and length scales that translate between computational and laboratory units. Rescaling by these scales allows direct quantitative comparisons among MSD curves and between velocity autocorrelation functions from experiments and simulations. Although our simulations replicate experimentally observed spontaneous symmetry breaking, short-timescale diffusive motion, and spontaneous cell-motion reorientation, their computational cost is low, allowing their use in multiscale virtual-tissue simulations. Comparisons between experimental and simulated cell motion support the hypothesis that short-time actomyosin dynamics affects longer-time cell motility. The success of the base cell-migration simulation model suggests its future application in more complex situations, including chemotaxis, migration through complex 3D matrices, and collective cell motion.     

Influence of individual cell motility on the 2D front roughness dynamics of tumour cell colonies

The dynamics of in situ 2D HeLa cell quasi-linear and quasi-radial colony fronts in a standard culture medium is investigated. For quasi-radial colonies, as the cell population increased, a kinetic transition from an exponential to a constant front average velocity regime was observed. Special attention was paid to individual cell motility evolution under constant average colony front velocity looking for its impact on the dynamics of the 2D colony front roughness. From the directionalities and velocity components of cell trajectories in colonies with different cell populations, the influence of both local cell density and cell crowding effects on individual cell motility was determined. The average dynamic behaviour of individual cells in the colony and its dependence on both local spatio-temporal heterogeneities and growth geometry suggested that cell motion undergoes under a concerted cell migration mechanism, in which both a limiting random walk-like and a limiting ballistic-like contribution were involved. These results were interesting to infer how biased cell trajectories influenced both the 2D colony spreading dynamics and the front roughness characteristics by local biased contributions to individual cell motion. These data are consistent with previous experimental and theoretical cell colony spreading data and provide additional evidence of the validity of the Kardar-Parisi-Zhang equation, within a certain range of time and colony front size, for describing the dynamics of 2D colony front roughness.     

Magnetic particle motions within living cells. Measurement of cytoplasmic viscosity and motile activity.

Submicrometer magnetic particles, ingested by cells and monitored via the magnetic fields they generate, provide an alternative to optical microscopy for probing movement and viscosity of living cytoplasm, and can be used for cells both in vitro and in vivo. We present methods for preparing lung macrophages tagged with magnetic particles for magnetometric study. Interpretation of the data involves fitting experimental remanent-field decay curves to nonlinear mechanistic models of intracellular particle motion. The model parameters are sensitive to mobility and apparent cytoplasmic viscosity experienced by particle-containing organelles. We present results of parameter estimation for intracellular particle behavior both within control cells and after (a) variable magnetization duration, (b) incubation with cytochalasin D, and (c) particle twisting by external fields. Magnetometric analysis showed cytoplasmic elasticity, dose-dependent motion inhibition by cytochalasin D, and a shear-thinning apparent viscosity.     

恒定磁场出生前暴露对小鼠胚胎发育的影响

目的:探讨不同强度恒定磁场对小鼠胚胎发育的影响。方法:将52只受孕母鼠分别暴露于不同强度的恒定磁场下(0.04T、0.08T、0.12T),17只怀孕雌鼠作为对照,孕18日时处死母鼠,解剖后记录仔鼠总数,活胎数、吸收胎数,并称重胎重,测量胎仔身长、尾长以及计算心脏、肝、脑、肾脏与体重的比值。结果:在0.08T及0.12T磁场暴露下,平均活胎数较对照组明显减少。平均吸收胎数明显增加,与对照组相比差异有显著性(P<0.05)。胎鼠尾长、心、肝、脑、肾与体重比值各组间无显著性差异,但体重与身长在0.12T磁场环境中要较另外三组显著增加(P<0.05)。0.12T磁场暴露下畸形胎仔增加,其中外表畸形较另外三组差异有显著性。结论:磁场对小鼠胚胎发育有一定的影响,并且其影响与磁场强度相关。     

Intracellular fluorescent probe concentrations by confocal microscopy.

A general method is described that takes advantage of the optical sectioning properties of a confocal microscope to enable measurement of both absolute and relative concentrations of fluorescent molecules inside cells. For compartments within cells that are substantially larger than the point spread function, the fluorescence intensity is simply proportional to the concentration of the fluorophore. For small compartments, the fluorescence intensity is diluted by contributions from regions outside the compartment. Corrections for this dilution can be estimated via calibrations that are based on the intensity distribution found in a computationally synthesized model for a cell or organelle that has been blurred by convolution with the microscope point spread function. The method is illustrated with four test cases: estimation of intracellular concentration of a fluorescent calcium indicator; estimation of the relative distribution between the neurite and soma of a neuronal cell of the InsP3 receptor on the endoplasmic reticulum; estimation of the distribution of the bradykinin receptor along the surface of a neuronal cell; and relative distribution of a potentiometric dye between the mitochondria and cytosol as a means of assaying mitochondrial membrane potential.     

Oscillating focus of SopA associated with filamentous structure guides partitioning of F plasmid

The F plasmid is actively partitioned to daughter cells by the sopABC gene. To elucidate the partitioning mechanisms, we simultaneously analysed movements of the plasmid and the SopA ATPase in single living cells. SopA, which is a putative motor protein assembled densely near nucleoid borders and formed a single discrete focus associated with less dense filamentous distribution along the long axis of the cell. The dense SopA focus oscillates between cell poles. The direction of the plasmid motion switches as the SopA focus switches its position. The velocity of the plasmid motion stays constant while it oscillates moving towards the SopA focus. The low density filamentous distribution of SopA persisted throughout the SopA oscillation. The focus associated with filamentous distribution of SopA was also observed in a cell without nucleoid. The SopA filament may guide the movement of the plasmid as a railway track and lead it to cell quarters.     

Wide-Field Multi-Parameter FLIM: long-term minimal invasive observation of proteins in living cells

Time-domain Fluorescence Lifetime Imaging Microscopy (FLIM) is a remarkable tool to monitor the dynamics of fluorophore-tagged protein domains inside living cells. We propose a Wide-Field Multi-Parameter FLIM method (WFMP-FLIM) aimed to monitor continuously living cells under minimum light intensity at a given illumination energy dose. A powerful data analysis technique applied to the WFMP-FLIM data sets allows to optimize the estimation accuracy of physical parameters at very low fluorescence signal levels approaching the lower bound theoretical limit. We demonstrate the efficiency of WFMP-FLIM by presenting two independent and relevant long-term experiments in cell biology: 1) FRET analysis of simultaneously recorded donor and acceptor fluorescence in living HeLa cells and 2) tracking of mitochondrial transport combined with fluorescence lifetime analysis in neuronal processes.     

Molecular dynamics of individual alpha-helices of bacteriorhodopsin in dimyristol phosphatidylocholine. I. Structure and dynamics.

Understanding the role of the lipid bilayer in membrane protein structure and dynamics is needed for tertiary structure determination methods. However, the molecular details are not well understood. Molecular dynamics computer calculations can provide insight into these molecular details of protein:lipid interactions. This paper reports on 10 simulations of individual alpha-helices in explicit lipid bilayers. The 10 helices were selected from the bacteriorhodopsin structure as representative alpha-helical membrane folding components. The bilayer is constructed of dimyristoyl phosphatidylcholine molecules. The only major difference between simulations is the primary sequence of the alpha-helix. The results show dramatic differences in motional behavior between alpha-helices. For example, helix A has much smaller root-mean-squared deviations than does helix D. This can be understood in terms of the presence of aromatic residues at the interface for helix A that are not present in helix D. Additional motions are possible for the helices that contain proline side chains relative to other amino acids. The results thus provide insight into the types of motion and the average structures possible for helices within the bilayer setting and demonstrate the strength of molecular simulations in providing molecular details that are not directly visualized in experiments.     

A theory on the determination of 3D motion and 3D structure from features

The present paper proposes a mathematical theory and a method of recognition of both the 3D structure and the motion of a moving object from its monocular image. Initially, characteristic features are extracted from the 2D perspective image of the object. Because motion of the object induces a change in its 2D perspective image, it also induces a change in the features which depends on the 3D structure and the velocity of the object. This suggests the possibility of detecting the 3D structure and the motion directly from the features and their changing rate, without the need for calculating optical flows. An analysis is made of the relation between the 3D rigid motion of a surface element and the change in local linear features. From this relation, a method is proposed for calculating the velocity of and the normal to the surface element without considering any correspondence of points. An optical flow can also be calculated by this method. Two simple computer simulations are provided.     

Electro-orientation of ellipsoidal erythrocytes. Theory and experiment.

The frequency-dependent orientation of human and llama erythrocytes suspended in isotonic solutions and subjected to linearly polarized electric fields is examined. Human erythrocytes may be represented as oblate spheroids (3.9:3.9:1.1 microns) with two distinguishable orientations, while the llama cells are approximated as ellipsoids with three distinct axes (4.0:2.0:1.1 microns). Under appropriate experimental conditions, both orientations of the human cells and all three orientations of the llama cells are observed. A theoretical cell model which accounts for the membrane as a thin confocal layer of ideal capacitance is used to predict the orientational spectra. The predicted spectra compare favorably in frequency range and orientational sequence with experimental data. Estimates for cell internal conductivity and permittivity are obtained by adjusting the values of these important parameters to achieve the closet fit of the theoretical curves to the data. By the use of this method, the internal conductivity of llama erythrocytes is estimated to be 0.26 S/m (+/- 20%), while the effective internal dielectric constant and conductivity of Euglena gracilis are estimated to be 120 (+/- 10%) and 0.43 S/m (+/- 20%), respectively.     

Contact response of cells can mediate morphogenetic pattern formation

Theories of morphogenetic pattern formation have included Turing's chemical prepatterns, mechanochemical interactions, cell sorting, and other mechanisms involving guided motion or signalling of cells. Many of these theories presuppose long-range cellular communication or other controls such as chemical concentration fields. However, the possibility that direct interactions between cells can lead to order and structure has not been seriously investigated in mathematical models. In this paper we consider this possibility, with emphasis on cells that reorient and align with each other when they come into contact. We show that such contact responses can account for the formation of multicellular patterns called parallel arrays. These patterns typically occur in tissue cultures of fibroblasts, and consist of clusters of cells sharing a common axis of orientation. Using predictions of a mathematical model and computer simulations of cell motion and interactions we show that contact responses alone, in the absence of other global controls, can promote the formation of these patterns. We suggest other situations in which patterns may result from direct cellular communication. Previous theories of morphogenesis are briefly reviewed and compared with this proposed mechanism.     

The impact of fabrication parameters and substrate stiffness in direct writing of living constructs

Biomolecules and living cells can be printed in high‐resolution patterns to fabricate living constructs for tissue engineering. To evaluate the impact of processing cells with rapid prototyping (RP) methods, we modeled the printing phase of two RP systems that use biomaterial inks containing living cells: a high‐resolution inkjet system (BioJet) and a lower‐resolution nozzle‐based contact printing system (PAM²). In the first fabrication method, we reasoned that cell damage occurs principally during drop collision on the printing surface, in the second we hypothesize that shear stresses act on cells during extrusion (within the printing nozzle). The two cases were modeled changing the printing conditions: biomaterial substrate stiffness and volumetric flow rate, respectively, in BioJet and PAM². Results show that during inkjet printing impact energies of about 10^?8 J are transmitted to cells, whereas extrusion energies of the order of 10^?11 J are exerted in direct printing. Viability tests of printed cells can be related to those numerical simulations, suggesting a threshold energy of 10^?9 J to avoid permanent cell damage. To obtain well‐defined living constructs, a combination of these methods is proposed for the fabrication of scaffolds with controlled 3D architecture and spatial distribution of biomolecules and cells. © 2012 American Institute of Chemical Engineers Biotechnol. Prog., 2012     

A model for encoding of magnetic field intensity by magnetite-based magnetoreceptor cells

A conceptual model is proposed for the encoding of magnetic field intensity from the motion of a chain of single-domain magnetite crystals which is located within a receptor cell, connected at one end to the cell membrane, and linked by cytoskeletal filaments to an array of mechanically gated ion channels centred on the end of the chain. In this arrangement, the physical links between the chain and ion channels will restrict the motion of the magnetite chain in response to the external magnetic field to a narrow cone with its axis through the point where the chain is attached to the membrane. The motion of the chain in the presence of an external magnetic field and thermal agitation will open a varying number of channels, causing the membrane potential to oscillate about some mean value that depends on the component of magnetic intensity oriented perpendicular to the cell membrane. The model permits estimation of magnetic intensity by integration of the motion of the magnetite chain over an area of the cell membrane, explains a number of results from physiological recordings in birds and fish, and makes testable predictions for future experimental studies. The model also provides a mechanism at the cellular level for a constant value of the Weber fraction (the ratio of the threshold sensitivity to a stimulus and the magnitude of that stimulus) for the magnetic sense but requires a separate gain control mechanism for modulation of sensitivity over a range of background fields. If magnetic field detection and encoding works as proposed in the model, the magnetoreceptor system may also be able to reconstruct the magnetic field vector using information about the vertical and horizontal axes from the eyes, gravity detectors, or both.     

A framework for three-dimensional simulation of morphogenesis

We present COMPUCELL3D, a software framework for three-dimensional simulation of morphogenesis in different organisms. COMPUCELL3D employs biologically relevant models for cell clustering, growth, and interaction with chemical fields. COMPUCELL3D uses design patterns for speed, efficient memory management, extensibility, and flexibility to allow an almost unlimited variety of simulations. We have verified COMPUCELL3D by building a model of growth and skeletal pattern formation in the avian (chicken) limb bud. Binaries and source code are available, along with documentation and input files for sample simulations, at http:// compucell.sourceforge.net.