The tree of life : the powerful mathematical idea at the heart of evolution

The way in which all living things — from microscopic bacteria to human beings to giant sequoia — are related to one another follows a deep and unexpected mathematical pattern, which we now know as the 'tree of life'. This pattern was discovered by naturalists over centuries, but it was Charles Darwin and Alfred Russel Wallace who realised that it held the key to understanding the origin and diversity of life on Eart

Tuning Genetic Clocks Employing DNA Binding Sites

Periodic oscillations play a key role in cell physiology from the cell cycle to circadian clocks. The interplay of positive and negative feedback loops among genes and proteins is ubiquitous in these networks. Often, delays in a negative feedback loop and/or degradation rates are a crucial mechanism to obtain sustained oscillations. How does nature control delays and kinetic rates in feedback networks? Known mechanisms include proper selection of the number of steps composing a feedback loop and alteration of protease activity, respectively. Here, we show that a remarkably simple means to control both delays and effective kinetic rates is the employment of DNA binding sites. We illustrate this design principle on a widely studied activator-repressor clock motif, which is ubiquitous in natural systems. By suitably employing DNA target sites for the activator and/or the repressor, one can switch the clock “on” and “off” and precisely tune its period to a desired value. Our study reveals a design principle to engineer dynamic behavior in biomolecular networks, which may be largely exploited by natural systems and employed for the rational design of synthetic circuits.United States. Air Force Office of Scientific Research (Grant FA9550-09-1-0211)National Science Foundation (U.S.). (Communication and Information Foundations) (Grant 1058127

Modeling of the Bacterial Mechanism of Methicillin-Resistance by a Systems Biology Approach

BACKGROUND: A microorganism is a complex biological system able to preserve its functional features against external perturbations and the ability of the living systems to oppose to these external perturbations is defined "robustness". The antibiotic resistance, developed by different bacteria strains, is a clear example of robustness and of ability of the bacterial system to acquire a particular functional behaviour in response to environmental changes. In this work we have modeled the whole mechanism essential to the methicillin-resistance through a systems biology approach. The methicillin is a beta-lactamic antibiotic that act by inhibiting the penicillin-binding proteins (PBPs). These PBPs are involved in the synthesis of peptidoglycans, essential mesh-like polymers that surround cellular enzymes and are crucial for the bacterium survival. METHODOLOGY: The network of genes, mRNA, proteins and metabolites was created using CellDesigner program and the data of molecular interactions are stored in Systems Biology Markup Language (SBML). To simulate the dynamic behaviour of this biochemical network, the kinetic equations were associated with each reaction. CONCLUSIONS: Our model simulates the mechanism of the inactivation of the PBP by methicillin, as well as the expression of PBP2a isoform, the regulation of the SCCmec elements (SCC: staphylococcal cassette chromosome) and the synthesis of peptidoglycan by PBP2a. The obtained results by our integrated approach show that the model describes correctly the whole phenomenon of the methicillin resistance and is able to respond to the external perturbations in the same way of the real cell. Therefore, this model can be useful to develop new therapeutic approaches for the methicillin control and to understand the general mechanism regarding the cellular resistance to some antibiotics

Population distribution analyses reveal a hierarchy of molecular players underlying parallel endocytic pathways.

Single-cell-resolved measurements reveal heterogeneous distributions of clathrin-dependent (CD) and -independent (CLIC/GEEC: CG) endocytic activity in Drosophila cell populations. dsRNA-mediated knockdown of core versus peripheral endocytic machinery induces strong changes in the mean, or subtle changes in the shapes of these distributions, respectively. By quantifying these subtle shape changes for 27 single-cell features which report on endocytic activity and cell morphology, we organize 1072 Drosophila genes into a tree-like hierarchy. We find that tree nodes contain gene sets enriched in functional classes and protein complexes, providing a portrait of core and peripheral control of CD and CG endocytosis. For 470 genes we obtain additional features from separate assays and classify them into early- or late-acting genes of the endocytic pathways. Detailed analyses of specific genes at intermediate levels of the tree suggest that Vacuolar ATPase and lysosomal genes involved in vacuolar biogenesis play an evolutionarily conserved role in CG endocytosis

Structural Analysis of Biodiversity

Large, recently-available genomic databases cover a wide range of life forms, suggesting opportunity for insights into genetic structure of biodiversity. In this study we refine our recently-described technique using indicator vectors to analyze and visualize nucleotide sequences. The indicator vector approach generates correlation matrices, dubbed Klee diagrams, which represent a novel way of assembling and viewing large genomic datasets. To explore its potential utility, here we apply the improved algorithm to a collection of almost 17000 DNA barcode sequences covering 12 widely-separated animal taxa, demonstrating that indicator vectors for classification gave correct assignment in all 11000 test cases. Indicator vector analysis revealed discontinuities corresponding to species- and higher-level taxonomic divisions, suggesting an efficient approach to classification of organisms from poorly-studied groups. As compared to standard distance metrics, indicator vectors preserve diagnostic character probabilities, enable automated classification of test sequences, and generate high-information density single-page displays. These results support application of indicator vectors for comparative analysis of large nucleotide data sets and raise prospect of gaining insight into broad-scale patterns in the genetic structure of biodiversity

Optimal Resting-Growth Strategies of Microbial Populations in Fluctuating Environments

Bacteria spend most of their lifetime in non-growing states which allow them to survive extended periods of stress and starvation. When environments improve, they must quickly resume growth to maximize their share of limited nutrients. Cells with higher stress resistance often survive longer stress durations at the cost of needing more time to resume growth, a strong disadvantage in competitive environments. Here we analyze the basis of optimal strategies that microorganisms can use to cope with this tradeoff. We explicitly show that the prototypical inverse relation between stress resistance and growth rate can explain much of the different types of behavior observed in stressed microbial populations. Using analytical mathematical methods, we determine the environmental parameters that decide whether cells should remain vegetative upon stress exposure, downregulate their metabolism to an intermediate optimum level, or become dormant. We find that cell-cell variability, or intercellular noise, is consistently beneficial in the presence of extreme environmental fluctuations, and that it provides an efficient population-level mechanism for adaption in a deteriorating environment. Our results reveal key novel aspects of responsive phenotype switching and its role as an adaptive strategy in changing environments

Automatic Design of Synthetic Gene Circuits through Mixed Integer Non-linear Programming

Automatic design of synthetic gene circuits poses a significant challenge to synthetic biology, primarily due to the complexity of biological systems, and the lack of rigorous optimization methods that can cope with the combinatorial explosion as the number of biological parts increases. Current optimization methods for synthetic gene design rely on heuristic algorithms that are usually not deterministic, deliver sub-optimal solutions, and provide no guaranties on convergence or error bounds. Here, we introduce an optimization framework for the problem of part selection in synthetic gene circuits that is based on mixed integer non-linear programming (MINLP), which is a deterministic method that finds the globally optimal solution and guarantees convergence in finite time. Given a synthetic gene circuit, a library of characterized parts, and user-defined constraints, our method can find the optimal selection of parts that satisfy the constraints and best approximates the objective function given by the user. We evaluated the proposed method in the design of three synthetic circuits (a toggle switch, a transcriptional cascade, and a band detector), with both experimentally constructed and synthetic promoter libraries. Scalability and robustness analysis shows that the proposed framework scales well with the library size and the solution space. The work described here is a step towards a unifying, realistic framework for the automated design of biological circuits

In-Silico Patterning of Vascular Mesenchymal Cells in Three Dimensions

Cells organize in complex three-dimensional patterns by interacting with proteins along with the surrounding extracellular matrix. This organization provides the mechanical and chemical cues that ultimately influence a cell's differentiation and function. Here, we computationally investigate the pattern formation process of vascular mesenchymal cells arising from their interaction with Bone Morphogenic Protein-2 (BMP-2) and its inhibitor, Matrix Gla Protein (MGP). Using a first-principles approach, we derive a reaction-diffusion model based on the biochemical interactions of BMP-2, MGP and cells. Simulations of the model exhibit a wide variety of three-dimensional patterns not observed in a two-dimensional analysis. We demonstrate the emergence of three types of patterns: spheres, tubes, and sheets, and show that the patterns can be tuned by modifying parameters in the model such as the degradation rates of proteins and chemotactic coefficient of cells. Our model may be useful for improved engineering of three-dimensional tissue structures as well as for understanding three dimensional microenvironments in developmental processes.National Institutes of Health (U.S.) (GM69811)United States. Dept. of Energy (DOE CSGF fellowship