Random subwindows and extremely randomized trees for image classification in cell biology

Background: With the improvements in biosensors and high-throughput image acquisition technologies, life science laboratories are able to perform an increasing number of experiments that involve the generation of a large amount of images at different imaging modalities/scales. It stresses the need for computer vision methods that automate image classification tasks. Results: We illustrate the potential of our image classification method in cell biology by evaluating it on four datasets of images related to protein distributions or subcellular localizations, and red-blood cell shapes. Accuracy results are quite good without any specific pre-processing neither domain knowledge incorporation. The method is implemented in Java and available upon request for evaluation and research purpose. Conclusion: Our method is directly applicable to any image classification problems. We foresee the use of this automatic approach as a baseline method and first try on various biological image classification problems

Dynamics of defect formation

A dynamic symmetry-breaking transition with noise and inertia is analyzed. Exact solution of the linearized equation that describes the critical region allows precise calculation (exponent and prefactor) of the number of defects produced as a function of the rate of increase of the critical parameter. The procedure is valid in both the overdamped and underdamped limits. In one space dimension, we perform quantitative comparison with numerical simulations of the nonlinear nonautonomous stochastic partial differential equation and report on signatures of underdamped dynamics.Comment: 4 pages, LaTeX, 4 figures. Submitted to Physical Revie

Cortical Factor Feedback Model for Cellular Locomotion and Cytofission

Eukaryotic cells can move spontaneously without being guided by external cues. For such spontaneous movements, a variety of different modes have been observed, including the amoeboid-like locomotion with protrusion of multiple pseudopods, the keratocyte-like locomotion with a widely spread lamellipodium, cell division with two daughter cells crawling in opposite directions, and fragmentations of a cell to multiple pieces. Mutagenesis studies have revealed that cells exhibit these modes depending on which genes are deficient, suggesting that seemingly different modes are the manifestation of a common mechanism to regulate cell motion. In this paper, we propose a hypothesis that the positive feedback mechanism working through the inhomogeneous distribution of regulatory proteins underlies this variety of cell locomotion and cytofission. In this hypothesis, a set of regulatory proteins, which we call cortical factors, suppress actin polymerization. These suppressing factors are diluted at the extending front and accumulated at the retracting rear of cell, which establishes a cellular polarity and enhances the cell motility, leading to the further accumulation of cortical factors at the rear. Stochastic simulation of cell movement shows that the positive feedback mechanism of cortical factors stabilizes or destabilizes modes of movement and determines the cell migration pattern. The model predicts that the pattern is selected by changing the rate of formation of the actin-filament network or the threshold to initiate the network formation

A multiscale hybrid model for pro-angiogenic calcium signals in a vascular endothelial cell

Cytosolic calcium machinery is one of the principal signaling mechanisms by which endothelial cells (ECs) respond to external stimuli during several biological processes, including vascular progression in both physiological and pathological conditions. Low concentrations of angiogenic factors (such as VEGF) activate in fact complex pathways involving, among others, second messengers arachidonic acid (AA) and nitric oxide (NO), which in turn control the activity of plasma membrane calcium channels. The subsequent increase in the intracellular level of the ion regulates fundamental biophysical properties of ECs (such as elasticity, intrinsic motility, and chemical strength), enhancing their migratory capacity. Previously, a number of continuous models have represented cytosolic calcium dynamics, while EC migration in angiogenesis has been separately approached with discrete, lattice-based techniques. These two components are here integrated and interfaced to provide a multiscale and hybrid Cellular Potts Model (CPM), where the phenomenology of a motile EC is realistically mediated by its calcium-dependent subcellular events. The model, based on a realistic 3-D cell morphology with a nuclear and a cytosolic region, is set with known biochemical and electrophysiological data. In particular, the resulting simulations are able to reproduce and describe the polarization process, typical of stimulated vascular cells, in various experimental conditions.Moreover, by analyzing the mutual interactions between multilevel biochemical and biomechanical aspects, our study investigates ways to inhibit cell migration: such strategies have in fact the potential to result in pharmacological interventions useful to disrupt malignant vascular progressio

Modelling Cell Polarization Driven by Synthetic Spatially Graded Rac Activation

The small GTPase Rac is known to be an important regulator of cell polarization, cytoskeletal reorganization, and motility of mammalian cells. In recent microfluidic experiments, HeLa cells endowed with appropriate constructs were subjected to gradients of the small molecule rapamycin leading to synthetic membrane recruitment of a Rac activator and direct graded activation of membrane-associated Rac. Rac activation could thus be triggered independent of upstream signaling mechanisms otherwise responsible for transducing activating gradient signals. The response of the cells to such stimulation depended on exceeding a threshold of activated Rac. Here we develop a minimal reaction-diffusion model for the GTPase network alone and for GTPase-phosphoinositide crosstalk that is consistent with experimental observations for the polarization of the cells. The modeling suggests that mutual inhibition is a more likely mode of cell polarization than positive feedback of Rac onto its own activation. We use a new analytical tool, Local Perturbation Analysis, to approximate the partial differential equations by ordinary differential equations for local and global variables. This method helps to analyze the parameter space and behaviour of the proposed models. The models and experiments suggest that (1) spatially uniform stimulation serves to sensitize a cell to applied gradients. (2) Feedback between phosphoinositides and Rho GTPases sensitizes a cell. (3) Cell lengthening/flattening accompanying polarization can increase the sensitivity of a cell and stabilize an otherwise unstable polarization

Catalyst composition and impurity-dependent kinetics of nanowire heteroepitaxy.

The mechanisms and kinetics of axial Ge-Si nanowire heteroepitaxial growth based on the tailoring of the Au catalyst composition via Ga alloying are studied by environmental transmission electron microscopy combined with systematic ex situ CVD calibrations. The morphology of the Ge-Si heterojunction, in particular, the extent of a local, asymmetric increase in nanowire diameter, is found to depend on the Ga composition of the catalyst, on the TMGa precursor exposure temperature, and on the presence of dopants. To rationalize the findings, a general nucleation-based model for nanowire heteroepitaxy is established which is anticipated to be relevant to a wide range of material systems and device-enabling heterostructures.S.H. acknowledges funding from ERC grant InsituNANO (No. 279342). A.D.G. acknowledges funding from the Marshall Aid Commemoration Commission and the National Science Foundation. C.D. acknowledges funding from the Royal Society. A portion of the research was also performed using EMSL, a national scientific user facility sponsored by the Department of Energy’s (DOE) Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory (PNNL). PNNL is operated by Battelle for the U.S. DOE under Contract DE-AC05-76RL01830. We gratefully acknowledge the use of facilities within the LeRoy Eyring Center for Solid State Science at Arizona State University. This work was performed in part at CINT, a U.S. DOE, Office of Science User Facility. The research was funded in part by the Laboratory Directed Research and Development Program at LANL, an affirmative action equal opportunity employer operated by Los Alamos National Security, LLC, for the National Nuclear Security Administration of the U.S. DOE under Contract DE-AC52-06NA25396.This document is the Accepted Manuscript version of a Published Work that appeared in final form in ACS Nano, copyright © American Chemical Society after peer review and technical editing by the publisher. To access the final edited and published work see http://pubs.acs.org/doi/abs/10.1021/nn402208p. Gamalski AD, Perea DE, Yoo J, Li N, Olszta MJ, Colby R, Schreiber DK, Ducati C, Picraux ST, Hofmann S, ACS Nano 2013, 7 (9), 7689–7697, doi:10.1021/nn402208

The Roots of Bioinformatics in Theoretical Biology

From the late 1980s onward, the term “bioinformatics” mostly has been used to refer to computational methods for comparative analysis of genome data. However, the term was originally more widely defined as the study of informatic processes in biotic systems. In this essay, I will trace this early history (from a personal point of view) and I will argue that the original meaning of the term is re-emerging

Factors limiting the transmission of HIV mutations conferring drug resistance: fitness costs and genetic bottlenecks

Transmission of HIV strains with drug-resistance mutations (DRMs) causes public health problems in resource-rich countries. We use a stochastic model, with data from viral competition experiments, to analyze the effect of fitness costs (FCs) and genetic bottlenecks on limiting transmission of 10 clinically significant DRMs. Transmission of DRMs with low FCs (∼0.2%) is similar to wild-type; transmission chains last ∼8 generations causing clusters of ∼60 infected individuals. Genetic bottlenecks substantially limit transmission of DRMs with moderately high FCs (∼0.6%); chains last ∼1–3 generations with transmission clusters of 2–7. Transmission of DRMs with extremely high FCs (>6%) only occurs from ∼5% of index cases. DRMs can revert to wild-type and remain as minority strains, within treatment-naïve individuals, undetectable by current resistance assays. We calculate, based on assay sensitivity, the length of time each DRM is detectable within individuals. Taken together, our results imply a hidden epidemic of transmitted resistance may exist

Theoretical Model for Cellular Shapes Driven by Protrusive and Adhesive Forces

The forces that arise from the actin cytoskeleton play a crucial role in determining the cell shape. These include protrusive forces due to actin polymerization and adhesion to the external matrix. We present here a theoretical model for the cellular shapes resulting from the feedback between the membrane shape and the forces acting on the membrane, mediated by curvature-sensitive membrane complexes of a convex shape. In previous theoretical studies we have investigated the regimes of linear instability where spontaneous formation of cellular protrusions is initiated. Here we calculate the evolution of a two dimensional cell contour beyond the linear regime and determine the final steady-state shapes arising within the model. We find that shapes driven by adhesion or by actin polymerization (lamellipodia) have very different morphologies, as observed in cells. Furthermore, we find that as the strength of the protrusive forces diminish, the system approaches a stabilization of a periodic pattern of protrusions. This result can provide an explanation for a number of puzzling experimental observations regarding cellular shape dependence on the properties of the extra-cellular matrix