Abstract: Bacterial communities such as biofilms are widely recognised as being important for survival and persistence of bacteria in harsh environments. Mechanistic models of biofilm growth indicate that the way in which the surface is seeded can effect the morphology of simulated biofilms. Experimental studies indicate that genes which are important for chemotaxis also influence biofilm formation, perhaps by influencing aggregation on a surface. Understanding aggregation and microcolony formation could therefore help clarify factors influencing biofilm formation as well as understanding how groups may influence the fitness of bacteria. In this paper I develop an Individual Based Model to examine how different behaviours involved in microcolony formation on a surface determines patterns of group sizes, and link patterns to bacterial fitness. I also provide a method for comparing data with model hypotheses to identify bacterial behaviours in experimental systems.

Abstract: Multivariate polynomial dynamical systems over finite fields have been studied in several contexts, including engineering and mathematical biology. An important problem is to construct models of such systems from a partial specification of dynamic properties, e.g., from a collection of state transition measurements. Here, we consider static models, which are directed graphs that represent the causal relationships between system variables, so-called wiring diagrams. This paper contains an algorithm which computes all possible minimal wiring diagrams for a given set of state transition measurements. The paper also contains several statistical measures for model selection. The algorithm uses primary decomposition of monomial ideals as the principal tool. An application to the reverse-engineering of a gene regulatory network is included. The algorithm and the statistical measures are implemented in Macaulay2 and are available from the authors.

Abstract: We introduce a family of positive definite kernels specifically optimized for the manipulation of 3D structures of molecules with kernel methods. The kernels are based on the comparison of the three-points pharmacophores present in the 3D structures of molecul es, a set of molecular features known to be particularly relevant for virtual screening applications. We present a computationally demanding exact implementation of these kernels, as well as fast approximations related to the classical fingerprint-based approa ches. Experimental results suggest that this new approach outperforms state-of-the-art algorithms based on the 2D structure of mol ecules for the detection of inhibitors of several drug targets.

Abstract: This paper focuses on the study of certain classes of Boolean functions that have appeared in several different contexts. Nested canalyzing functions have been studied recently in the context of Boolean network models of gene regulatory networks. In the same context, polynomial functions over finite fields have been used to develop network inference methods for gene regulatory networks. Finally, unate cascade functions have been studied in the design of logic circuits and binary decision diagrams. This paper shows that the class of nested canalyzing functions is equal to that of unate cascade functions. Furthermore, it provides a description of nested canalyzing functions as a certain type of Boolean polynomial function. Using the polynomial framework one can show that the class of nested canalyzing functions, or, equivalently, the class of unate cascade functions, forms an algebraic variety which makes their analysis amenable to the use of techniques from algebraic geometry and computational algebra. As a corollary of the functional equivalence derived here, a formula in the literature for the number of unate cascade functions provides such a formula for the number of nested canalyzing functions.

Abstract: Finding coarse-grained, low-dimensional descriptions is an important task in the analysis of complex, stochastic models of gene regulatory networks. This task involves (a) identifying observables that best describe the state of these complex systems and (b) characterizing the dynamics of the observables. In a previous paper [13], we assumed that good observables were known a priori, and presented an equation-free approach to approximate coarse-grained quantities (i.e, effective drift and diffusion coefficients) that characterize the long-time behavior of the observables. Here we use diffusion maps [9] to extract appropriate observables ("reduction coordinates") in an automated fashion; these involve the leading eigenvectors of a weighted Laplacian on a graph constructed from network simulation data. We present lifting and restriction procedures for translating between physical variables and these data-based observables. These procedures allow us to perform equation-free coarse-grained, computations characterizing the long-term dynamics through the design and processing of short bursts of stochastic simulation initialized at appropriate values of the data-based observables.

Abstract: Canalization of genetic regulatory networks has been argued to be favored by evolutionary processes due to the stability that it can confer to phenotype expression. We explore whether a significant amount of canalization and partial canalization can arise in purely random networks in the absence of evolutionary pressures. We use a mapping of the Boolean functions in the Kauffman N-K model for genetic regulatory networks onto a k-dimensional Ising hypercube to show that the functions can be divided into different classes strictly due to geometrical constraints. The classes can be counted and their properties determined using results from group theory and isomer chemistry. We demonstrate that partially canalized functions completely dominate all possible Boolean functions, particularly for higher k. This indicates that partial canalization is extremely common, even in randomly chosen networks, and has implications for how much information can be obtained in experiments on native state genetic regulatory networks.

Abstract: An analysis of data from 212 permanent sample plots provided no evidence of any decline in rainforest productivity after three cycles of selection logging in the tropical rainforests of north Queensland. Relative productivity was determined as the difference between observed diameter increments and increments predicted from a diameter increment function which incorporated tree size, stand density and site quality. Analyses of variance and regression analyses revealed no significant decline in productivity after repeated harvesting. There is evidence to support the assertion that if any permanent productivity decline exists, it does not exceed six per cent per harvest.

Abstract: We present a maskless microscope projection lithography system (MPLS), in which photomasks have been replaced by a Digital Micromirror Device type spatial light modulator (DMD, Texas Instruments). Employing video projector technology high resolution patterns, designed as bitmap images on the computer, are displayed using a micromirror array consisting of about 786000 tiny individually addressable tilting mirrors. The DMD, which is located in the image plane of an infinity corrected microscope, is projected onto a substrate placed in the focal plane of the microscope objective. With a 5x(0.25 NA) Fluar microscope objective, a fivefold reduction of the image to a total size of 9 mm2 and a minimum feature size of 3.5 microns is achieved. Our system can be used in the visible range as well as in the near UV (with a light intensity of up to 76 mW/cm2 around the 365 nm Hg-line). We developed an inexpensive and simple method to enable exact focusing and controlling of the image quality of the projected patterns. Our MPLS has originally been designed for the light-directed in situ synthesis of DNA microarrays. One requirement is a high UV intensity to keep the fabrication process reasonably short. Another demand is a sufficient contrast ratio over small distances (of about 5 microns). This is necessary to achieve a high density of features (i.e. separated sites on the substrate at which different DNA sequences are synthesized in parallel fashion) while at the same time the number of stray light induced DNA sequence errors is kept reasonably small. We demonstrate the performance of the apparatus in light-directed DNA chip synthesis and discuss its advantages and limitations.

Abstract: We have developed an alignment-free method that calculates phylogenetic distances using a maximum likelihood approach for a model of sequence change on patterns that are discovered in unaligned sequences. To evaluate the phylogenetic accuracy of our method, and to conduct a comprehensive comparison of existing alignment-free methods (freely available as Python package decaf+py at this http URL), we have created a dataset of reference trees covering a wide range of phylogenetic distances. Amino acid sequences were evolved along the trees and input to the tested methods; from their calculated distances we infered trees whose topologies we compared to the reference trees. We find our pattern-based method statistically superior to all other tested alignment-free methods on this dataset. We also demonstrate the general advantage of alignment-free methods over an approach based on automated alignments when sequences violate the assumption of collinearity. Similarly, we compare methods on empirical data from an existing alignment benchmark set that we used to derive reference distances and trees. Our pattern-based approach yields distances that show a linear relationship to reference distances over a substantially longer range than other alignment-free methods. The pattern-based approach outperforms alignment-free methods and its phylogenetic accuracy is statistically indistinguishable from alignment-based distances.

Abstract: A class of models of biological population and communities with a singular equilibrium at the origin is analyzed; it is shown that these models can possess a dynamical regime of deterministic extinction, which is crucially important from the biological standpoint. This regime corresponds to the presence of a family of homoclinics to the origin, so-called elliptic sector. The complete analysis of possible topological structures in a neighborhood of the origin, as well as asymptotics to orbits tending to this point, is given. An algorithmic approach to analyze system behavior with parameter changes is presented. The developed methods and algorithm are applied to existing mathematical models of biological systems. In particular, we analyze a model of anticancer treatment with oncolytic viruses, a parasite-host interaction model, and a model of Chagas' disease.

Abstract: Motivated by applications in systems biology, we seek a probabilistic framework based on Markov processes to represent intracellular processes. We review the formal relationships between different stochastic models referred to in the systems biology literature. As part of this review, we present a novel derivation of the differential Chapman-Kolmogorov equation for a general multidimensional Markov process made up of both continuous and jump processes. We start with the definition of a time-derivative for a probability density but place no restrictions on the probability distribution, in particular, we do not assume it to be confined to a region that has a surface (on which the probability is zero). In our derivation, the master equation gives the jump part of the Markov process while the Fokker-Planck equation gives the continuous part. We thereby sketch a {}``family tree'' for stochastic models in systems biology, providing explicit derivations of their formal relationship and clarifying assumptions involved.

Abstract: In this paper, we propose a mean-field model which attempts to bridge the gap between random Boolean networks and more realistic stochastic modeling of genetic regulatory networks. The main idea of the model is to replace all regulatory interactions to any one gene with an average or effective interaction, which takes into account the repression and activation mechanisms. We find that depending on the set of regulatory parameters, the model exhibits rich nonlinear dynamics. The model also provides quantitative support to the earlier qualitative results obtained for random Boolean networks.

Abstract: We present a dynamical model of lipoprotein metabolism derived by combining a cascading process in the blood stream and cellular level regulatory dynamics. We analyse the existence and stability of equilibria and show that this low-dimensional, nonlinear model exhibits bistability between a low and a high cholesterol state. A sensitivity analysis indicates that the intracellular concentration of cholesterol is robust to parametric variations while the plasma cholesterol can vary widely. We show how the dynamical response to time-dependent inputs can be used to diagnose the state of the system. We also establish the connection between parameters in the system and medical and genetic conditions.

Abstract: We model a plant leaf by using two-dimensional domain walls with internal structures Such domain walls can be found as soliton solutions in field theory describing magnetic materials The radiation scattered by such domain walls behaves quite similar to the spectral reflectance of plant leaves The model nicely simulates the spectral reflectance of a plant leaf as a function of the wavelength

Abstract: Kolmogorov complexity has inspired several alignment-free distance measures, based on the comparison of lengths of compressions, which have been applied successfully in many areas. One of these measures, the socalled Universal Similarity Metric (USM), has been used by Krasnogor and Pelta to compare simple protein contact maps, showing that it yielded good clustering on four small datasets. We report an extensive test of this metric using a much larger and representative protein dataset: the domain dataset used by Sierk and Pearson to evaluate seven protein structure comparison methods and two protein sequence comparison methods. One result is that Krasnogor-Pelta method has less domain discriminant power than any one of the methods considered by Sierk and Pearson when using these simple contact maps. In another test, we found that the USM based distance has low agreement with the CATH tree structure for the same benchmark of Sierk and Pearson. In any case, its agreement is lower than the one of a standard sequential alignment method, SSEARCH. Finally, we manually found lots of small subsets of the database that are better clustered using SSEARCH than USM, to confirm that Krasnogor-Pelta’s conclusions were based on datasets that were too small.

Abstract: Affected relatives are essential for pedigree linkage analysis, however, they cause a violation of the independent sample assumption in case-control association studies. To avoid the correlation between samples, a common practice is to take only one affected sample per pedigree in association analysis. Although several methods exist in handling correlated samples, they are still not widely used in part because these are not easily implemented, or because they are not widely known. We advocate the effective sample size method as a simple and accessible approach for case-control association analysis with correlated samples. This method modifies the chi-square test statistic, p-value, and 95% confidence interval of the odds-ratio by replacing the apparent number of allele or genotype counts with the effective ones in the standard formula, without the need for specialized computer programs. We present a simple formula for calculating effective sample size for many types of relative pairs and relative sets. For allele frequency estimation, the effective sample size method captures the variance inflation exactly. For genotype frequency, simulations showed that effective sample size provides a satisfactory approximation. A gene which is previously identified as a type 1 diabetes susceptibility locus, the interferon-induced helicase gene (IFIH1), is shown to be significantly associated with rheumatoid arthritis when the effective sample size method is applied. This significant association is not established if only one affected sib per pedigree were used in the association analysis. Relationship between the effective sample size method and other methods -- the generalized estimation equation, variance of eigenvalues for correlation matrices, and genomic controls -- are discussed.

Abstract: The organization and mining of malaria genomic and post-genomic data is highly motivated by the necessity to predict and characterize new biological targets and new drugs. Biological targets are sought in a biological space designed from the genomic data from Plasmodium falciparum, but using also the millions of genomic data from other species. Drug candidates are sought in a chemical space containing the millions of small molecules stored in public and private chemolibraries. Data management should therefore be as reliable and versatile as possible. In this context, we examined five aspects of the organization and mining of malaria genomic and post-genomic data: 1) the comparison of protein sequences including compositionally atypical malaria sequences, 2) the high throughput reconstruction of molecular phylogenies, 3) the representation of biological processes particularly metabolic pathways, 4) the versatile methods to integrate genomic data, biological representations and functional profiling obtained from X-omic experiments after drug treatments and 5) the determination and prediction of protein structures and their molecular docking with drug candidate structures. Progresses toward a grid-enabled chemogenomic knowledge space are discussed.

Abstract: In this paper the existence and unicity of a stable periodic orbit is proven, for a class of piecewise affine differential equations in dimension 3 or more, provided their interaction structure is a negative feedback loop. It is also shown that the same systems converge toward a unique stable equilibrium point in dimension 2. This extends a theorem of Snoussi, which showed the existence of these orbits only. The considered class of equations is usually studied as a model of gene regulatory networks. It is not assumed that all decay rates are identical, which is biologically irrelevant, but has been done in the vast majority of previous studies. Our work relies on classical results about fixed points of monotone, concave operators acting on positive variables. Moreover, the used techniques are very likely to apply in more general contexts, opening directions for future work.

Abstract: Spatially distributed problems are often approximately modelled in terms of partial differential equations (PDEs) for appropriate coarse-grained quantities (e.g. concentrations). The derivation of accurate such PDEs starting from finer scale, atomistic models, and using suitable averaging, is often a challenging task; approximate PDEs are typically obtained through mathematical closure procedures (e.g. mean-field approximations). In this paper, we show how such approximate macroscopic PDEs can be exploited in constructing preconditioners to accelerate stochastic simulations for spatially distributed particle-based process models. We illustrate how such preconditioning can improve the convergence of equation-free coarse-grained methods based on coarse timesteppers. Our model problem is a stochastic reaction-diffusion model capable of exhibiting Turing instabilities.

Abstract: Scaling analysis of heart rate time series has emerged as an useful tool for assessment of autonomic cardiac control. We investigate the heart rate time series of ten athletes (five males and five females), by applying detrended fluctuation analysis (DFA). High resolution ECGs are recorded under standardized resting conditions over 30 minutes and subsequently heart rate time series are extracted and artefacts filtered. We find three distinct regions of scale-invariance, which correspond to the well-known VLF, LF, and HF bands in the power spectra of heart rate variability. The scaling exponents alpha are alphaHF: 1.15 [0.96-1.22], alphaLF: 0.68 [0.57-0.84], alphaVLF: 0.83[0.82-0.99]; p<10^-5). In conclusion, DFA scaling exponents of heart rate time series should be fitted to the VLF, LF, and HF ranges, respectively.

Abstract: Scientists investigate the dynamics of complex systems with quantitative models, employing them to synthesize knowledge, to explain observations, and to forecast future system behavior. Complete specification of systems is impossible, so models must be simplified abstractions. Thus, the art of modeling involves deciding which system elements to include and determining how they should be represented. We view modeling as search through a space of candidate models that is guided by model objectives, theoretical knowledge, and empirical data. In this contribution, we introduce a method for representing process-based models that facilitates the discovery of models that explain observed behavior. This representation casts dynamic systems as interacting sets of processes that act on entities. Using this approach, a modeler first encodes relevant ecological knowledge into a library of generic entities and processes, then instantiates these theoretical components, and finally assembles candidate models from these elements. We illustrate this methodology with a model of the Ross Sea ecosystem.

Abstract: Fish farms represent a growing source of disturbance to shallow benthic ecosystems like seagrass meadows. Despite some existing insights on the mechanisms underlying decline, efficient tools to quantitatively predict the response of benthic communities to fish farm effluents have not yet been developed. We explored relationships of fish farm organic and nutrient input rates to the sediments with population dynamics of the key seagrass species (Posidonia oceanica) in deep meadows growing around four Mediterranean Sea bream and Sea bass fish farms. We performed 2 annual shoot censuses on permanent plots at increasing distance from cages. Before each census we measured sedimentation rates adjacent to the plots using benthic sediment traps. High shoot mortality rates were recorded near the cages, up to 20 times greater than at control sites. Recruitment rates remained similar to undisturbed meadows and could not compensate mortality, leading to rapid seagrass decline within the first 100 meters from cages. Seagrass mortality increased with total (R2= 0.47, p< 0.0002), organic matter (R2= 0.36, p= 0.001), nitrogen (R2= 0.34, p= 0.002) and phosphorus (R2= 0.58, p< 3 x 10-5) sedimentation rates. P. oceanica decline accelerated above a phosphorus loading threshold of 50 mg m-2 day-1. Benthic sedimentation rates seem a powerful predictor of seagrass mortality from fish farming, integrating local hydrodynamics, waste effluents variability and several environmental mechanisms, fuelled by organic inputs and leading to seagrass loss. Coupling direct measurements of benthic sedimentation rates with dynamics of key species is proposed as an efficient way to predict and minimize fish farm impacts to benthic communities.

Abstract: Statistical properties of interbeat intervals cascade are evaluated by considering the joint probability distribution $P(\Delta x_2,\tau_2;\Delta x_1,\tau_1)$ for two interbeat increments $\Delta x_1$ and $\Delta x_2$ of different time scales $\tau_1$ and $\tau_2$. We present evidence that the conditional probability distribution $P(\Delta x_2,\tau_2|\Delta x_1,\tau_1)$ may obey a Chapman-Kolmogorov equation. The corresponding Kramers-Moyal (KM) coefficients are evaluated. It is shown that while the first and second KM coefficients, i.e., the drift and diffusion coefficients, take on well-defined and significant values, the higher-order coefficients in the KM expansion are very small. As a result, the joint probability distributions of the increments in the interbeat intervals obey a Fokker-Planck equation. The method provides a novel technique for distinguishing the two classes of subjects in terms of the drift and diffusion coefficients, which behave differently for two classes of the subjects, namely, healthy subjects and those with congestive heart failure.

Abstract: This paper presents efficient algorithms for solving the problem of aligning a protein structure template to a query amino-acid sequence, known as protein threading problem. We consider the problem as a special case of graph matching problem. We give formal graph and integer programming models of the problem. After studying the properties of these models, we propose two kinds of Lagrangian relaxation for solving them. We present experimental results on real life instances showing the efficiency of our approaches.

Abstract: Through this paper, we call for a distributed, internet-based collaboration to address one of the worst plagues of our present world, malaria. The spirit is a non-proprietary peer-production of information-embedding goods. And we propose to use the grid technology to enable such a world wide "open source" like collaboration. The first step towards this vision has been achieved during the summer on the EGEE grid infrastructure where 46 million ligands were docked for a total amount of 80 CPU years in 6 weeks in the quest for new drugs.

Abstract: We study a classical two-state stochastic system in a sea of substrates and products (absorbing states), which can be interpreted as a single Michaelis-Menten catalyzing enzyme or as a channel on a cell surface. We introduce a novel general method and use it to derive the expression for the full counting statistics of transitions among the absorbing states. For the evolution of the system under a periodic perturbation of the kinetic rates, the latter contains a term with a purely geometrical (the Berry phase) interpretation. This term gives rise to a pump current between the absorbing states, which is due entirely to the stochastic nature of the system. We calculate the first two cumulants of this current, and we argue that it is observable experimentally.

Abstract: A new type of microfluidic system for biological cell manipulation, a CMOS/microfluidic hybrid, is demonstrated. The hybrid system starts with a custom-designed CMOS (complementary metal-oxide semiconductor) chip fabricated in a semiconductor foundry using standard integration circuit technology. A microfluidic channel is post-fabricated on top of the CMOS chip to provide biocompatible environment. The motion of individual biological cells that are tagged with magnetic beads is directly controlled by the CMOS chip that generates localized magnetic filed patterns using an on-chip array of micro-electromagnets. The speed and the programmability of the CMOS chip further allow for the dynamic reconfiguration of the magnetic fields, substantially increasing the manipulation capability of the hybrid system. The concept of a hybrid system is verified by simultaneously manipulating individual biological cells with microscopic resolution. A new operation protocol that exploits the fast speed of electronics to trap and move a large number of cells with less power consumption is also demonstrated. Combining the advantages of microelectronics, the CMOS/microfluidic hybrid approach presents a new model for a multifunctional lab-on-a chip for biological and medical applications.

Abstract: Erikkson showed that singular value decomposition(SVD) of flattenings determined a partition of a phylogenetic tree to be a split. In this paper, based on his work, we develop new statistically consistent algorithms fit for grid computing to construct a phylogenetic tree by computing SVD of flattenings with the small fixed number of rows.

Abstract: We present a novel method for finding low dimensional views of high dimensional data: Targeted Projection Pursuit. The method proceeds by finding projections of the data that best approximate a target view. Two versions of the method are introduced; one version based on Procrustes analysis and one based on a single layer perceptron. These versions are capable of finding orthogonal or non-orthogonal projections respectively. The method is quantitatively and qualitatively compared with other dimension reduction techniques. It is shown to find two-dimensional views that display the classification of cancers from gene expression data with a visual separation equal to, or better than, existing dimension reduction techniques.