In silico estimates of the free energy rates in growing tumor spheroids

The physics of solid tumor growth can be considered at three distinct size scales: the tumor scale, the cell-extracellular matrix (ECM) scale and the sub-cellular scale. In this paper we consider the tumor scale in the interest of eventually developing a system-level understanding of the progression of cancer. At this scale, cell populations and chemical species are best treated as concentration fields that vary with time and space. The cells have chemo-mechanical interactions with each other and with the ECM, consume glucose and oxygen that are transported through the tumor, and create chemical byproducts. We present a continuum mathematical model for the biochemical dynamics and mechanics that govern tumor growth. The biochemical dynamics and mechanics also engender free energy changes that serve as universal measures for comparison of these processes. Within our mathematical framework we therefore consider the free energy inequality, which arises from the first and second laws of thermodynamics. With the model we compute preliminary estimates of the free energy rates of a growing tumor in its pre-vascular stage by using currently available data from single cells and multicellular tumor spheroids.Comment: 27 pages with 5 figures and 2 tables. Figures and tables appear at the end of the pape

Onset of collective and cohesive motion

We study the onset of collective motion, with and without cohesion, of groups of noisy self-propelled particles interacting locally. We find that this phase transition, in two space dimensions, is always discontinuous, including for the minimal model of Vicsek et al. [Phys. Rev. Lett. {\bf 75},1226 (1995)] for which a non-trivial critical point was previously advocated. We also show that cohesion is always lost near onset, as a result of the interplay of density, velocity, and shape fluctuations.Comment: accepted for publication in Phys. Rev. Let

Hydrodynamic fluctuations and instabilities in ordered suspensions of self-propelled particles

We construct the hydrodynamic equations for {\em suspensions} of self-propelled particles (SPPs) with spontaneous orientational order, and make a number of striking, testable predictions:(i) SPP suspensions with the symmetry of a true {\em nematic} are {\em always} absolutely unstable at long wavelengths.(ii) SPP suspensions with {\em polar}, i.e., head-tail {\em asymmetric}, order support novel propagating modes at long wavelengths, coupling orientation, flow, and concentration. (iii) In a wavenumber regime accessible only in low Reynolds number systems such as bacteria, polar-ordered suspensions are invariably convectively unstable.(iv) The variance in the number N of particles, divided by the mean , diverges as $^{2/3}$ in polar-ordered SPP suspensions.Comment: submitted to Phys Rev Let

Actin Fusion Proteins Alter the Dynamics of Mechanically Induced Cytoskeleton Rearrangement

Mechanical forces can regulate various functions in living cells. The cytoskeleton is a crucial element for the transduction of forces in cell-internal signals and subsequent biological responses. Accordingly, many studies in cellular biomechanics have been focused on the role of the contractile acto-myosin system in such processes. A widely used method to observe the dynamic actin network in living cells is the transgenic expression of fluorescent proteins fused to actin. However, adverse effects of GFP-actin fusion proteins on cell spreading, migration and cell adhesion strength have been reported. These shortcomings were shown to be partly overcome by fusions of actin binding peptides to fluorescent proteins. Nevertheless, it is not understood whether direct labeling by actin fusion proteins or indirect labeling via these chimaeras alters biomechanical responses of cells and the cytoskeleton to forces. We investigated the dynamic reorganization of actin stress fibers in cells under cyclic mechanical loading by transiently expressing either egfp-Lifeact or eyfp-actin in rat embryonic fibroblasts and observing them by means of live cell microscopy. Our results demonstrate that mechanically-induced actin stress fiber reorganization exhibits very different kinetics in EYFP-actin cells and EGFP-Lifeact cells, the latter showing a remarkable agreement with the reorganization kinetics of non-transfected cells under the same experimental conditions

Active Brownian Particles. From Individual to Collective Stochastic Dynamics

We review theoretical models of individual motility as well as collective dynamics and pattern formation of active particles. We focus on simple models of active dynamics with a particular emphasis on nonlinear and stochastic dynamics of such self-propelled entities in the framework of statistical mechanics. Examples of such active units in complex physico-chemical and biological systems are chemically powered nano-rods, localized patterns in reaction-diffusion system, motile cells or macroscopic animals. Based on the description of individual motion of point-like active particles by stochastic differential equations, we discuss different velocity-dependent friction functions, the impact of various types of fluctuations and calculate characteristic observables such as stationary velocity distributions or diffusion coefficients. Finally, we consider not only the free and confined individual active dynamics but also different types of interaction between active particles. The resulting collective dynamical behavior of large assemblies and aggregates of active units is discussed and an overview over some recent results on spatiotemporal pattern formation in such systems is given.Comment: 161 pages, Review, Eur Phys J Special-Topics, accepte

Gene Expression Variability within and between Human Populations and Implications toward Disease Susceptibility

Variations in gene expression level might lead to phenotypic diversity across individuals or populations. Although many human genes are found to have differential mRNA levels between populations, the extent of gene expression that could vary within and between populations largely remains elusive. To investigate the dynamic range of gene expression, we analyzed the expression variability of ∼18, 000 human genes across individuals within HapMap populations. Although ∼20% of human genes show differentiated mRNA levels between populations, our results show that expression variability of most human genes in one population is not significantly deviant from another population, except for a small fraction that do show substantially higher expression variability in a particular population. By associating expression variability with sequence polymorphism, intriguingly, we found SNPs in the untranslated regions (5′ and 3′UTRs) of these variable genes show consistently elevated population heterozygosity. We performed differential expression analysis on a genome-wide scale, and found substantially reduced expression variability for a large number of genes, prohibiting them from being differentially expressed between populations. Functional analysis revealed that genes with the greatest within-population expression variability are significantly enriched for chemokine signaling in HIV-1 infection, and for HIV-interacting proteins that control viral entry, replication, and propagation. This observation combined with the finding that known human HIV host factors show substantially elevated expression variability, collectively suggest that gene expression variability might explain differential HIV susceptibility across individuals

The role of noise and positive feedback in the onset of autosomal dominant diseases

<p>Abstract</p> <p>Background</p> <p>Autosomal dominant (AD) diseases result when a single mutant or non-functioning gene is present on an autosomal chromosome. These diseases often do not emerge at birth. There are presently two prevailing theories explaining the expression of AD diseases. One explanation originates from the Knudson two-hit theory of hereditary cancers, where loss of heterozygosity or occurrence of somatic mutations impairs the function of the wild-type copy. While these somatic second hits may be sufficient for stable disease states, it is often difficult to determine if their occurrence necessarily marks the initiation of disease progression. A more direct consequence of a heterozygous genetic background is haploinsufficiency, referring to a lack of sufficient gene function due to reduced wild-type gene copy number; however, haploinsufficiency can involve a variety of additional mechanisms, such as noise in gene expression or protein levels, injury and second hit mutations in other genes. In this study, we explore the possible contribution to the onset of autosomal dominant diseases from intrinsic factors, such as those determined by the structure of the molecular networks governing normal cellular physiology.</p> <p>Results</p> <p>First, simple models of single gene insufficiency using the positive feedback loops that may be derived from a three-component network were studied by computer simulation using Bionet software. The network structure is shown to affect the dynamics considerably; some networks are relatively stable even when large stochastic variations in are present, while others exhibit switch-like dynamics. In the latter cases, once the network switches over to the disease state it remains in that state permanently. Model pathways for two autosomal dominant diseases, AD polycystic kidney disease and mature onset diabetes of youth (MODY) were simulated and the results are compared to known disease characteristics.</p> <p>Conclusions</p> <p>By identifying the intrinsic mechanisms involved in the onset of AD diseases, it may be possible to better assess risk factors as well as lead to potential new drug targets. To illustrate the applicability of this study of pathway dynamics, we simulated the primary pathways involved in two autosomal dominant diseases, Polycystic Kidney Disease (PKD) and mature onset diabetes of youth (MODY). Simulations demonstrate that some of the primary disease characteristics are consistent with the positive feedback - stochastic variation theory presented here. This has implications for new drug targets to control these diseases by blocking the positive feedback loop in the relevant pathways.</p

Effects of Transcriptional Pausing on Gene Expression Dynamics

Stochasticity in gene expression affects many cellular processes and is a source of phenotypic diversity between genetically identical individuals. Events in elongation, particularly RNA polymerase pausing, are a source of this noise. Since the rate and duration of pausing are sequence-dependent, this regulatory mechanism of transcriptional dynamics is evolvable. The dependency of pause propensity on regulatory molecules makes pausing a response mechanism to external stress. Using a delayed stochastic model of bacterial transcription at the single nucleotide level that includes the promoter open complex formation, pausing, arrest, misincorporation and editing, pyrophosphorolysis, and premature termination, we investigate how RNA polymerase pausing affects a gene's transcriptional dynamics and gene networks. We show that pauses' duration and rate of occurrence affect the bursting in RNA production, transcriptional and translational noise, and the transient to reach mean RNA and protein levels. In a genetic repressilator, increasing the pausing rate and the duration of pausing events increases the period length but does not affect the robustness of the periodicity. We conclude that RNA polymerase pausing might be an important evolvable feature of genetic networks

Genomic Restructuring in the Tasmanian Devil Facial Tumour: Chromosome Painting and Gene Mapping Provide Clues to Evolution of a Transmissible Tumour

Devil facial tumour disease (DFTD) is a fatal, transmissible malignancy that threatens the world's largest marsupial carnivore, the Tasmanian devil, with extinction. First recognised in 1996, DFTD has had a catastrophic effect on wild devil numbers, and intense research efforts to understand and contain the disease have since demonstrated that the tumour is a clonal cell line transmitted by allograft. We used chromosome painting and gene mapping to deconstruct the DFTD karyotype and determine the chromosome and gene rearrangements involved in carcinogenesis. Chromosome painting on three different DFTD tumour strains determined the origins of marker chromosomes and provided a general overview of the rearrangement in DFTD karyotypes. Mapping of 105 BAC clones by fluorescence in situ hybridisation provided a finer level of resolution of genome rearrangements in DFTD strains. Our findings demonstrate that only limited regions of the genome, mainly chromosomes 1 and X, are rearranged in DFTD. Regions rearranged in DFTD are also highly rearranged between different marsupials. Differences between strains are limited, reflecting the unusually stable nature of DFTD. Finally, our detailed maps of both the devil and tumour karyotypes provide a physical framework for future genomic investigations into DFTD

A Methodological Framework for the Reconstruction of Contiguous Regions of Ancestral Genomes and Its Application to Mammalian Genomes

The reconstruction of ancestral genome architectures and gene orders from homologies between extant species is a long-standing problem, considered by both cytogeneticists and bioinformaticians. A comparison of the two approaches was recently investigated and discussed in a series of papers, sometimes with diverging points of view regarding the performance of these two approaches. We describe a general methodological framework for reconstructing ancestral genome segments from conserved syntenies in extant genomes. We show that this problem, from a computational point of view, is naturally related to physical mapping of chromosomes and benefits from using combinatorial tools developed in this scope. We develop this framework into a new reconstruction method considering conserved gene clusters with similar gene content, mimicking principles used in most cytogenetic studies, although on a different kind of data. We implement and apply it to datasets of mammalian genomes. We perform intensive theoretical and experimental comparisons with other bioinformatics methods for ancestral genome segments reconstruction. We show that the method that we propose is stable and reliable: it gives convergent results using several kinds of data at different levels of resolution, and all predicted ancestral regions are well supported. The results come eventually very close to cytogenetics studies. It suggests that the comparison of methods for ancestral genome reconstruction should include the algorithmic aspects of the methods as well as the disciplinary differences in data aquisition