Nathan Weinstein

Finding Missing Interactions of the Arabidopsis thaliana Root Stem Cell Niche Gene Regulatory Network.

Over the last few decades, the Arabidopsis thaliana root stem cell niche (RSCN) has become a model system for the study of plant development and stem cell niche dynamics. Currently, many of the molecular mechanisms involved in RSCN maintenance and development have been described. A few years ago, we published a gene regulatory network (GRN) model integrating this information. This model suggested that there were missing components or interactions. Upon updating the model, the observed stable gene configurations of the RSCN could not be recovered, indicating that there are additional missing components or interactions in the model. In fact, due to the lack of experimental data, GRNs inferred from published data are usually incomplete. However, predicting the location and nature of the missing data is a not trivial task. Here, we propose a set of procedures for detecting and predicting missing interactions in Boolean networks. We used these procedures to predict putative missing interactions in the A. thaliana RSCN network model. Using our approach, we identified three necessary interactions to recover the reported gene activation configurations that have been experimentally uncovered for the different cell types within the RSCN: (1) a regulation of PHABULOSA to restrict its expression domain to the vascular cells, (2) a self-regulation of WOX5, possibly by an indirect mechanism through the auxin signaling pathway, and (3) a positive regulation of JACKDAW by MAGPIE. The procedures proposed here greatly reduce the number of possible Boolean functions that are biologically meaningful and experimentally testable and that do not contradict previous data. We believe that these procedures can be used on any Boolean network. However, because the procedures were designed for the specific case of the RSCN, formal demonstrations of the procedures should be shown in future efforts.

Building Qualitative Models of Plant Regulatory Networks with SQUAD.

There is a wealth of molecular information that has enhanced the reconstruction of regulatory networks and pathways in plants, specially in Arabidopsis thaliana due to its central role as a model system. Most of this information, however, has proven to be insufficient to construct quantitative models due to the absence of sufficiently accurate measurements of kinetic constants. Hence, there have been efforts to develop qualitative methodologies that permit the use of the vast information regarding the regulatory interactions among molecules.

A network model for the specification of vulval precursor cells and cell fusion control in Caenorhabditis elegans

The vulva of Caenorhabditis elegans has been long used as an experimental model of cell differentiation and organogenesis. While it is known that the signaling cascades of Wnt, Ras/MAPK, and NOTCH interact to form a molecular network, there is no consensus regarding its precise topology and dynamical properties. We inferred the molecular network, and developed a multivalued synchronous discrete dynamic model to study its behavior. The model reproduces the patterns of activation reported for the following types of cell: vulval precursor, first fate, second fate, second fate with reversed polarity, third fate, and fusion fate. We simulated the fusion of cells, the determination of the first, second, and third fates, as well as the transition from the second to the first fate. We also used the model to simulate all possible single loss- and gain-of-function mutants, as well as some relevant double and triple mutants. Importantly, we associated most of these simulated mutants to multivulva, vulvaless, egg-laying defective, or defective polarity phenotypes. The model shows that it is necessary for RAL-1 to activate NOTCH signaling, since the repression of LIN-45 by RAL-1 would not suffice for a proper second fate determination in an environment lacking DSL ligands. We also found that the model requires the complex formed by LAG-1, LIN-12, and SEL-8 to inhibit the transcription of eff-1 in second fate cells. Our model is the largest reconstruction to date of the molecular network controlling the specification of vulval precursor cells and cell fusion control in C. elegans. According to our model, the process of fate determination in the vulval precursor cells is reversible, at least until either the cells fuse with the ventral hypoderm or divide, and therefore the cell fates must be maintained by the presence of extracellular signals.

The combination of the functionalities of feedback circuits is determinant for the attractors’ number and size in pathway-like Boolean networks

Molecular regulation was initially assumed to follow both a unidirectional and a hierarchical organization forming pathways. Regulatory processes, however, form highly interlinked networks with non-hierarchical and non-unidirectional structures that contain statistically overrepresented circuits or motifs. Here, we analyze the behavior of pathways containing non-unidirectional (i.e. bidirectional) and non-hierarchical interactions that create motifs. In comparison with unidirectional and hierarchical pathways, our pathways have a high diversity of behaviors, characterized by the size and number of attractors. Motifs have been studied individually showing that feedback circuit motifs regulate the number and size of attractors. It is less clear what happens in molecular networks that usually contain multiple feedbacks. Here, we find that the way feedback circuits couple to each other (i.e., the combination of the functionalities of feedback circuits) regulate both the number and size of the attractors. We show that the different expected results of epistasis analysis (a method to infer regulatory interactions) are produced by many non-hierarchical and non-unidirectional structures. Thus, these structures cannot be correctly inferred by epistasis analysis. Finally, we show that the combinations of functionalities, combined with other network properties, allow for a better characterization of regulatory structures.

A model of the regulatory network involved in the control of the cell cycle and cell differentiation in the Caenorhabditis elegans vulva

BackgroundThere are recent experimental reports on the cross-regulation between molecules involved in the control of the cell cycle and the differentiation of the vulval precursor cells (VPCs) of Caenorhabditis elegans. Such discoveries provide novel clues on how the molecular mechanisms involved in the cell cycle and cell differentiation processes are coordinated during vulval development. Dynamic computational models are helpful to understand the integrated regulatory mechanisms affecting these cellular processes.ResultsHere we propose a simplified model of the regulatory network that includes sufficient molecules involved in the control of both the cell cycle and cell differentiation in the C. elegans vulva to recover their dynamic behavior. We first infer both the topology and the update rules of the cell cycle module from an expected time series. Next, we use a symbolic algorithmic approach to find which interactions must be included in the regulatory network. Finally, we use a continuous-time version of the update rules for the cell cycle module to validate the cyclic behavior of the network, as well as to rule out the presence of potential artifacts due to the synchronous updating of the discrete model. We analyze the dynamical behavior of the model for the wild type and several mutants, finding that most of the results are consistent with published experimental results.ConclusionsOur model shows that the regulation of Notch signaling by the cell cycle preserves the potential of the VPCs and the three vulval fates to differentiate and de-differentiate, allowing them to remain completely responsive to the concentration of LIN-3 and lateral signal in the extracellular microenvironment.

The combination of the functionalities of feedback circuits is determinant for the number and size of attractors of molecular networks

Molecular regulation was initially assumed to follow both a unidirectional and a hierarchical organization forming pathways. Regulatory processes, however, form highly interlinked networks with non-hierarchical and non-unidirectional structures that contain statistically overrepresented circuits (motifs). Here, we analyze the behavior of pathways containing non-hierarchical and non-unidirectional interactions that create motifs. In comparison with unidirectional and hierarchical pathways, our pathways have a high diversity of behaviors, characterized by the size and number of attractors. Motifs have been studied individually showing that feedback circuit motifs regulate the number and size of attractors. It is less clear what happens in molecular networks that usually contain multiple feedbacks. Here, we find that the way feedback circuits couple to each other (i.e., the combination of the functionalities of feedback circuits) regulate both the precise number and size of the attractors. We show that the different sets of expected results of epistasis analysis (a method to infer regulatory interactions) are produced by many non-hierarchical and non-unidirectional structures. Thus, these structures cannot be correctly inferred by epistasis analysis. Finally, we show that the structures producing the epistasis results have remarkably similar sets of combinations of functionalities, that combined with other network properties could greatly improve epistasis analysis.

Organogenesis of the C. elegans Vulva and Control of Cell Fusion

The vulva of Caenorhabditis elegans is widely used as a paradigm for the study of organogenesis and is composed of seven toroids, formed by the migration of cells and the formation of homotypic contacts. Five of the toroids contain two or four nuclei and cell membrane fusion is one of the main driving forces during the morphogenesis of the vulva. The network of genes involved in the control of cell fusion during the formation of the vulva must determine which cells fuse and when. Especially during the formation of the vulval toroids, when those cells that fuse to form each ring, must not fuse with the neighbor cells, which form other separate rings. This is achieved through very fine control on the expression and function of several key genes.

A Network Model to Explore the Effect of the Micro-environment on Endothelial Cell Behavior during Angiogenesis

Angiogenesis is an important adaptation mechanism of the blood vessels to the changing requirements of the body during development, aging, and wound healing. Angiogenesis allows existing blood vessels to form new connections or to reabsorb existing ones. Blood vessels are composed of a layer of endothelial cells (ECs) covered by one or more layers of mural cells (smooth muscle cells or pericytes). We constructed a computational Boolean model of the molecular regulatory network involved in the control of angiogenesis. Our model includes the ANG/TIE, HIF, AMPK/mTOR, VEGF, IGF, FGF, PLCγ/Calcium, PI3K/AKT, NO, NOTCH, and WNT signaling pathways, as well as the mechanosensory components of the cytoskeleton. The dynamical behavior of our model recovers the patterns of molecular activation observed in Phalanx, Tip, and Stalk ECs. Furthermore, our model is able to describe the modulation of EC behavior due to extracellular micro-environments, as well as the effect due to loss- and gain-of-function mutations. These properties make our model a suitable platform for the understanding of the molecular mechanisms underlying some pathologies. For example, it is possible to follow the changes in the activation patterns caused by mutations that promote Tip EC behavior and inhibit Phalanx EC behavior, that lead to the conditions associated with retinal vascular disorders and tumor vascularization. Moreover, the model describes how mutations that promote Phalanx EC behavior are associated with the development of arteriovenous and venous malformations. These results suggest that the network model that we propose has the potential to be used in the study of how the modulation of the EC extracellular micro-environment may improve the outcome of vascular disease treatments.

Computational Simulation of the Hemodynamic Behavior of a Blood Vessel Network

During development, blood vessel networks adapt to gradual changes in the oxygen required by surrounding tissue, shear stress, and mechanical stretch. The possible adaptations include remodeling the vessel network and thickening the walls of blood vessels. However, the treatment of several vascular diseases including cerebral arteriovenous malformations, arteriosclerosis, aneurysms, and vascular retinal disorders, may lead to abrupt changes that could produce hemorrhage or other problems. Modeling the hemodynamic behavior of a blood vessel network may help assess or even diminish the risks associated with each treatment. In this work, we briefly describe the radiological studies available to study the anatomy and hemodynamics of a patient. We then describe the segmentation, smoothing, healing, skeletonyzation, and meshing processes that are needed to obtain an initial model for the numerical simulations. Additionally, we state some important concepts about blood rheology and blood vessel elasticity. Further, we include a system of equations to describe the interaction between flowing blood and the elastic blood vessels.

Modeling the Blood Vessels of the Brain

The results described in this work are part of a larger project. The long term goal of this project is to help physicians predict the hemodynamic changes, and associated risks, caused by different treatment options for brain arteriovenous malformations. First, we need to build a model of the vascular architecture of each specific patient. Our approach to build these models is described in this work. Later we will use the model of the vascular architecture to simulate the velocity and pressure gradients of the blood flowing within the vessels, and the stresses on the blood vessel walls, before and after treatment. We are developing a computer program to describe each blood vessel as a parametric curve, where each point within this curve includes a normal vector that points in the opposite direction of the pressure gradient. The shape of the cross section of the vessel in each point is described as an ellipse. Our program is able to describe the geometry of a blood vessel using as an input a cloud of dots. The program allows us to model any blood vessel, and other tubular structures.