Global parameter estimation for thermodynamic models of transcriptional regulation.

Deciphering the mechanisms involved in gene regulation holds the key to understanding the control of central biological processes, including human disease, population variation, and the evolution of morphological innovations. New experimental techniques including whole genome sequencing and transcriptome analysis have enabled comprehensive modeling approaches to study gene regulation. In many cases, it is useful to be able to assign biological significance to the inferred model parameters, but such interpretation should take into account features that affect these parameters, including model construction and sensitivity, the type of fitness calculation, and the effectiveness of parameter estimation. This last point is often neglected, as estimation methods are often selected for historical reasons or for computational ease. Here, we compare the performance of two parameter estimation techniques broadly representative of local and global approaches, namely, a quasi-Newton/Nelder-Mead simplex (QN/NMS) method and a covariance matrix adaptation-evolutionary strategy (CMA-ES) method. The estimation methods were applied to a set of thermodynamic models of gene transcription applied to regulatory elements active in the Drosophila embryo. Measuring overall fit, the global CMA-ES method performed significantly better than the local QN/NMS method on high quality data sets, but this difference was negligible on lower quality data sets with increased noise or on data sets simplified by stringent thresholding. Our results suggest that the choice of parameter estimation technique for evaluation of gene expression models depends both on quality of data, the nature of the models [again, remains to be established] and the aims of the modeling effort.

Bio-inspired compensatory strategies for damage to flapping robotic propulsors.

Natural swimmers and flyers can fully recover from catastrophic propulsor damage by altering stroke mechanics: some fish can lose even 76% of their propulsive surface without loss of thrust. We consider applying these principles to enable robotic flapping propulsors to autonomously repair functionality. However, direct transference of these alterations from an organism to a robotic flapping propulsor may be suboptimal owing to irrelevant evolutionary pressures. Instead, we use machine learning techniques to compare these alterations with those optimal for a robotic system. We implement an online artificial evolution with hardware-in-the-loop, performing experimental evaluations with a flexible plate. To recoup thrust, the learned strategy increased amplitude, frequency and angle of attack (AOA) amplitude, and phase-shifted AOA by approximately 110°. Only amplitude increase is reported by most fish literature. When recovering side force, we find that force direction is correlated with AOA. No clear amplitude or frequency trend is found, whereas frequency increases in most insect literature. These results suggest that how mechanical flapping propulsors most efficiently adjust to damage may not align with natural swimmers and flyers.

Advancing X-ray scattering metrology using inverse genetic algorithms.

We compare the speed and effectiveness of two genetic optimization algorithms to the results of statistical sampling via a Markov chain Monte Carlo algorithm to find which is the most robust method for determining real space structure in periodic gratings measured using critical dimension small angle X-ray scattering. Both a covariance matrix adaptation evolutionary strategy and differential evolution algorithm are implemented and compared using various objective functions. The algorithms and objective functions are used to minimize differences between diffraction simulations and measured diffraction data. These simulations are parameterized with an electron density model known to roughly correspond to the real space structure of our nanogratings. The study shows that for X-ray scattering data, the covariance matrix adaptation coupled with a mean-absolute error log objective function is the most efficient combination of algorithm and goodness of fit criterion for finding structures with little foreknowledge about the underlying fine scale structure features of the nanograting.