Parallel Simulated Annealing Using an Adaptive Resampling Interval

Classification-Based Inference of Dynamical Models of Gene Regulatory Networks

Cell-fate decisions during development are controlled by densely interconnected gene regulatory networks (GRNs) consisting of many genes. Inferring and predictively modeling these GRNs is crucial for understanding development and other physiological processes. Gene circuits, coupled differential equations that represent gene product synthesis with a switch-like function, provide a biologically realistic framework for modeling the time evolution of gene expression. However, their use has been limited to smaller networks due to the computational expense of inferring model parameters from gene expression data using global non-linear optimization. Here we show that the switch-like nature of gene regulation can be exploited to break the gene circuit inference problem into two simpler optimization problems that are amenable to computationally efficient supervised learning techniques. We present FIGR (Fast Inference of Gene Regulation), a novel classification-based inference approach to determining gene circuit parameters. We demonstrate FIGR’s effectiveness on synthetic data generated from random gene circuits of up to 50 genes as well as experimental data from the gap gene system of Drosophila melanogaster, a benchmark for inferring dynamical GRN models. FIGR is faster than global non-linear optimization by a factor of 600 and its computational complexity scales much better with GRN size. On a practical level, FIGR can accurately infer the biologically realistic gap gene network in under a minute on desktop-class hardware instead of requiring hours of parallel computing. We anticipate that FIGR would enable the inference of much larger biologically realistic GRNs than was possible before.

Development of a globally optimised model of the cerebral arteries

The cerebral arteries are difficult to reproduce from first principles, featuring interwoven territories, and intricate layers of grey and white matter with differing metabolic demand. The aim of this study was to identify the ideal configuration of arteries required to sustain an entire brain hemisphere based on minimisation of the energy required to supply the tissue. The 3D distribution of grey and white matter within a healthy human brain was first segmented from magnetic resonance images. A novel simulated annealing algorithm was then applied to determine the optimal configuration of arteries required to supply brain tissue. The model was validated through comparison of this ideal, entirely optimised, brain vasculature with the structure and properties of real arteries. This analysis established that the human cerebral vasculature is highly optimised; closely resembling the most energy efficient arrangement of vessels. In addition to local adherence to fluid dynamical optimisation principles, the optimised vasculature reproduced expected brain perfusion territories, featuring well-defined boundaries between anterior, middle and posterior regions. This validated brain vascular model and algorithm can be used for patient-specific modelling of stroke and cerebral haemodynamics, identification of sub-optimal conditions associated with vascular disease, and optimising vascular structures for tissue engineering applications and artificial organ design.