  The day we understand the time evolution of subcellular elements at a level
of detail comparable to physical systems governed by Newton's laws of motion
seems far away. Even so, quantitative approaches to cellular dynamics add to
our understanding of cell biology, providing data-guided frameworks that allow
us to develop better predictions about and methods for control over specific
biological processes and system-wide cell behavior. In this paper we describe
an approach to optimizing the use of transcription factors in the context of
cellular reprogramming. We construct an approximate model for the natural
evolution of a synchronized population of fibroblasts, based on data obtained
by sampling the expression of some 22,083 genes at several times along the cell
cycle. (These data are based on a colony of cells that have been cell cycle
synchronized) In order to arrive at a model of moderate complexity, we cluster
gene expression based on the division of the genome into topologically
associating domains (TADs) and then model the dynamics of the expression levels
of the TADs. Based on this dynamical model and known bioinformatics, we develop
a methodology for identifying the transcription factors that are the most
likely to be effective toward a specific cellular reprogramming task. The
approach used is based on a device commonly used in optimal control. From this
data-guided methodology, we identify a number of validated transcription
factors used in reprogramming and/or natural differentiation. Our findings
highlight the immense potential of dynamical models models, mathematics, and
data guided methodologies for improving methods for control over biological
processes.
